Harmonizing finite element modelling for non-invasive strength estimation by high-resolution peripheral quantitative computed tomography

Fracture risk based on high-resolution peripheral quantitative computed tomography measures does not vary with age in older adults-the bone microarchitecture international consortium prospective cohort study

Fracture risk increases with lower areal bone mineral density (aBMD); however, aBMD-related estimate of risk may decrease with age. This may depend on technical limitations of 2-dimensional (2D) dual energy X-ray absorptiometry (DXA) which are reduced with 3D high-resolution peripheral quantitative computed tomography (HR-pQCT). Our aim was to examine whether the predictive utility of HR-pQCT measures with fracture varies with age. We analyzed associations of HR-pQCT measures at the distal radius and distal tibia with two outcomes: incident fractures and major osteoporotic fractures. We censored follow-up time at first fracture, death, last contact or 8 years after baseline. We estimated hazard ratios (HR) and 95%CI for the association between bone traits and fracture incidence across age quintiles. Among 6835 men and women (ages 40-96) with at least one valid baseline HR-pQCT scan who were followed prospectively for a median of 48.3 months, 681 sustained fractures. After adjustment for confounders, bone parameters at both the radius and tibia were associated with higher fracture risk. The estimated HRs for fracture did not vary significantly across age quintiles for any HR-pQCT parameter measured at either the radius or tibia. In this large cohort, the homogeneity of the associations between the HR-pQCT measures and fracture risk across age groups persisted for all fractures and for major osteoporotic fractures. The patterns were similar regardless of the HR-pQCT measure, the type of fracture, or the statistical models. The stability of the associations between HR-pQCT measures and fracture over a broad age range shows that bone deficits or low volumetric density remain major determinants of fracture risk regardless of age group. The lower risk for fractures across measures of aBMD in older adults in other studies may be related to factors which interfere with DXA but not with HR-pQCT measures.

Bone remodeling and responsiveness to mechanical stimuli in individuals with type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) has been linked to increased osteocyte apoptosis, local accumulation of mineralized lacunar spaces, and microdamage suggesting an impairment of the mechanoregulation network in affected individuals. Diabetic neuropathy might exacerbate this dysfunction through direct effects on bone turnover, and indirect effects on balance, muscle strength, and gait. However, the in vivo effects of impaired bone mechanoregulation on bone remodeling in humans remain underexplored. This longitudinal cohort study assessed consenting participants with T1DM and varying degree of distal symmetric sensorimotor polyneuropathy (T1DM, n = 20, median age 46.5 yr, eight female) and controls (CTRL; n = 9, median age 59.0 yr, four female) at baseline and 4-yr follow-up. Nerve conduction in participants with T1DM was tested using DPNCheck and bone remodeling was quantified with longitudinal high-resolution peripheral quantitative-computed tomography (HR-pQCT, 82 μm) at the standard distal sites. Local trabecular bone formation (Tb.F) and resorption (Tb.R) sites were captured by implementing 3D rigid image registration of HR-pQCT images, and the mechanical environment across the bone microarchitecture at these sites was simulated using micro-finite element analysis. We calculated odds ratios to determine the likelihood of bone formation (ORF) and resorption (ORR) with increasing/decreasing strain in percent as markers for mechanoregulation. At the distal radius, Tb.F was 47% lower and Tb.R was 59% lower in T1DM participants compared with CTRL (P < .05). Tb.F correlated positively with nerve conduction amplitude (R = 0.69, P < .05) in participants with T1DM and negatively with glycated hemoglobin (HbA1c) (R = -0.45, P < .05). Additionally, ORF was 34% lower and ORR was 18% lower in T1DM compared with CTRL (P < .05). Our findings represent in vivo evidence suggesting that bone remodeling in individuals with T1DM is in a state of low responsiveness to mechanical stimuli, resulting in impaired bone formation and resorption rates; these correlate to the degree of neuropathy and level of diabetes control.

Recommendations for High-resolution Peripheral Quantitative Computed Tomography Assessment of Bone Density, Microarchitecture, and Strength in Pediatric Populations

PURPOSE OF REVIEW

The purpose of this review is to summarize current approaches and provide recommendations for imaging bone in pediatric populations using high-resolution peripheral quantitative computed tomography (HR-pQCT).

RECENT FINDINGS

Imaging the growing skeleton is challenging and HR-pQCT protocols are not standardized across centers. Adopting a single-imaging protocol for all studies is unrealistic; thus, we present three established protocols for HR-pQCT imaging in children and adolescents and share advantages and disadvantages of each. Limiting protocol variation will enhance the uniformity of results and increase our ability to compare study results between different research groups. We outline special cases along with tips and tricks for acquiring and processing scans to minimize motion artifacts and account for growing bone. The recommendations in this review are intended to help researchers perform HR-pQCT imaging in pediatric populations and extend our collective knowledge of bone structure, architecture, and strength during the growing years.

Comparison of Bone Quality Among Winter Endurance Athletes with and Without Risk Factors for Relative Energy Deficiency in Sport (REDs): A Cross-Sectional Study

Relative Energy Deficiency in Sport (REDs) is a syndrome describing the relationship between prolonged and/or severe low energy availability and negative health and performance outcomes. The high energy expenditures incurred during training and competition put endurance athletes at risk of REDs. The objective of this study was to investigate differences in bone quality in winter endurance athletes classified as either low-risk versus at-risk for REDs. Forty-four participants were recruited (M = 18; F = 26). Bone quality was assessed at the distal radius and tibia using high resolution peripheral quantitative computed tomography (HR-pQCT), and at the hip and spine using dual X-ray absorptiometry (DXA). Finite element analysis was used to estimate bone strength. Participants were grouped using modified criteria from the REDs Clinical Assessment Tool Version 1. Fourteen participants (M = 3; F = 11), were classified as at-risk of REDs (≥ 3 risk factors). Measured with HR-pQCT, cortical bone area (radius) and bone strength (radius and tibia) were 6.8%, 13.1% and 10.3% lower (p = 0.025, p = 0.033, p = 0.027) respectively, in at-risk compared with low-risk participants. Using DXA, femoral neck areal bone density was 9.4% lower in at-risk compared with low-risk participants (p = 0.005). At-risk male participants had 21.9% lower femoral neck areal bone density (via DXA) than low-risk males (p = 0.020) with no significant differences in females. Overall, 33.3% of athletes were at-risk for REDs and had lower bone quality than those at low-risk.

A Fracture Risk Assessment Tool for High Resolution Peripheral Quantitative Computed Tomography

Most fracture risk assessment tools use clinical risk factors combined with bone mineral density (BMD) to improve assessment of osteoporosis; however, stratifying fracture risk remains challenging. This study developed a fracture risk assessment tool that uses information about volumetric bone density and three-dimensional structure, obtained using high-resolution peripheral quantitative compute tomography (HR-pQCT), to provide an alternative approach for patient-specific assessment of fracture risk. Using an international prospective cohort of older adults (n = 6802) we developed a tool to predict osteoporotic fracture risk, called μFRAC. The model was constructed using random survival forests, and input predictors included HR-pQCT parameters summarizing BMD and microarchitecture alongside clinical risk factors (sex, age, height, weight, and prior adulthood fracture) and femoral neck areal BMD (FN aBMD). The performance of μFRAC was compared to the Fracture Risk Assessment Tool (FRAX) and a reference model built using FN aBMD and clinical covariates. μFRAC was predictive of osteoporotic fracture (c-index = 0.673, p < 0.001), modestly outperforming FRAX and FN aBMD models (c-index = 0.617 and 0.636, respectively). Removal of FN aBMD and all clinical risk factors, except age, from μFRAC did not significantly impact its performance when estimating 5-year and 10-year fracture risk. The performance of μFRAC improved when only major osteoporotic fractures were considered (c-index = 0.733, p < 0.001). We developed a personalized fracture risk assessment tool based on HR-pQCT that may provide an alternative approach to current clinical methods by leveraging direct measures of bone density and structure. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Precision of bone mechanoregulation assessment in humans using longitudinal high-resolution peripheral quantitative computed tomography in vivo

Local mechanical stimuli in the bone microenvironment are essential for the homeostasis and adaptation of the skeleton, with evidence suggesting that disruption of the mechanically-driven bone remodelling process may lead to bone loss. Longitudinal clinical studies have shown the combined use of high-resolution peripheral quantitative computed tomography (HR-pQCT) and micro-finite element analysis can be used to measure load-driven bone remodelling in vivo; however, quantitative markers of bone mechanoregulation and the precision of these analyses methods have not been validated in human subjects. Therefore, this study utilised participants from two cohorts. A same-day cohort (n = 33) was used to develop a filtering strategy to minimise false detections of bone remodelling sites caused by noise and motion artefacts present in HR-pQCT scans. A longitudinal cohort (n = 19) was used to develop bone imaging markers of trabecular bone mechanoregulation and characterise the precision for detecting longitudinal changes in subjects. Specifically, we described local load-driven formation and resorption sites independently using patient-specific odds ratios (OR) and 99 % confidence intervals. Conditional probability curves were computed to link the mechanical environment to the remodelling events detected on the bone surface. To quantify overall mechanoregulation, we calculated a correct classification rate measuring the fraction of remodelling events correctly identified by the mechanical signal. Precision was calculated as root-mean-squared averages of the coefficient of variation (RMS-SD) of repeated measurements using scan-rescan pairs at baseline combined with a one-year follow-up scan. We found no significant mean difference (p < 0.01) between scan-rescan conditional probabilities. RMS-SD was 10.5 % for resorption odds, 6.3 % for formation odds, and 1.3 % for correct classification rates. Bone was most likely to be formed in high-strain and resorbed in low-strain regions for all participants, indicating a consistent, regulated response to mechanical stimuli. For each percent increase in strain, the likelihood of bone resorption decreased by 2.0 ± 0.2 %, and the likelihood of bone formation increased by 1.9 ± 0.2 %, totalling 38.3 ± 1.1 % of strain-driven remodelling events across the entire trabecular compartment. This work provides novel robust bone mechanoregulation markers and their precision for designing future clinical studies.

Hip Fractures in Older Adults Are Associated With the Low Density Bone Phenotype and Heterogeneous Deterioration of Bone Microarchitecture

Femoral neck areal bone mineral density (FN aBMD) is a key determinant of fracture risk in older adults; however, the majority of individuals who have a hip fracture are not considered osteoporotic according to their FN aBMD. This study uses novel tools to investigate the characteristics of bone microarchitecture that underpin bone fragility. Recent hip fracture patients (n = 108, 77% female) were compared with sex- and age-matched controls (n = 216) using high-resolution peripheral quantitative computed tomography (HR-pQCT) imaging of the distal radius and tibia. Standard morphological analysis of bone microarchitecture, micro-finite element analysis, and recently developed techniques to identify void spaces in bone microarchitecture were performed to evaluate differences between hip fracture patients and controls. In addition, a new approach for phenotyping bone microarchitecture was implemented to evaluate whether hip fractures in males and females occur more often in certain bone phenotypes. Overall, hip fracture patients had notable deterioration of bone microarchitecture and reduced bone mineral density compared with controls, especially at weight-bearing sites (tibia and femoral neck). Hip fracture patients were more likely to have void spaces present at either site and had void spaces that were two to four times larger on average when compared with non-fractured controls (p < 0.01). Finally, bone phenotyping revealed that hip fractures were significantly associated with the low density phenotype (p < 0.01), with the majority of patients classified in this phenotype (69%). However, female and male hip fracture populations were distributed differently across the bone phenotype continuum. These findings highlight how HR-pQCT can provide insight into the underlying mechanisms of bone fragility by using information about bone phenotypes and identification of microarchitectural defects (void spaces). The added information suggests that HR-pQCT can have a beneficial role in assessing the severity of structural deterioration in bone that is associated with osteoporotic hip fractures. © 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).

High performance multi-platform computing for large-scale image-based finite element modeling of bone

BACKGROUND

Image-based finite element (FE) modeling of bone is a non-invasive method to estimate bone stiffness and strength. High-resolution imaging data as input allows for inclusion of bone microarchitecture but results in large amounts of data unsuitable for traditional FE solvers. Bone-specific mesh-free solvers have been developed over the past 20 years to improve on memory efficiency in simulated bone loading applications. The objective of this study was to provide linear performance benchmarking for a bone-specific, mesh-free solver (FAIM) using µCT and HR-pQCT image data on Mac, Linux, and Windows operating systems using both single- and multi-thread CPU and GPU processing.

METHODS

The focus is on the linear gradient-descent solver using standardized uniaxial loading of bone models from µCT, and first- and second-generation HR-pQCT scans of the radius and tibia. Convergence, speedup, memory, and batch performance tests were completed using CPUs and GPUs on three laboratory-based systems with Windows, Linux, and Mac operating systems.

RESULTS

Although varying by system and model size, time-per-iteration was as low as 0.03 s when an HR-pQCT-based radius model (6.45 million DOF) was solved with 3 GPUs. Strong scaling was achieved with GPU and CPU parallel processing, with strong parallel efficiencies when models were solved using 3 GPUs or ≤ 10 CPU threads. Errors in force, strain energy density, and Von Mises stress were as low as 0.1% when a convergence tolerance of 10^-3 or smaller was used.

CONCLUSION

The results of this study indicate that to maximize computational efficiency and minimize model solution times using FAIM software under the standardized tested conditions using µCT, XCT1 and XCT2 HR-pQCT image data, convergence tolerance set to 10^-4, and 10 threads or 2 GPUs are sufficient for efficient solution times. Less strict convergence tolerances will improve solution times but will introduce more error in the outcome measures.

HR-pQCT parameters of the distal radius correlate with the bending strength of the radial diaphysis

High resolution, peripheral quantitative computed tomography (HR-pQCT) scanners can now characterize an individual's trabecular architecture, cortical structure, and volumetric bone mineral density at a nominal resolution of 61 μm. While predictions of failure load of the distal radius and tibial diaphysis in compression by finite element analysis (FEA) of HR-pQCT scans have been validated against mechanical tests of cadaveric bones in compression, namely for images with nominal resolutions of 82 μm and 165 μm, the HR-pQCT parameters that best predict bending strength of cortical bone remain unknown. Therefore, we scanned cadaveric forearms from 31 elderly donors (Female: 72.8 ± 8.8 years and Male: 72.1 ± 6.3 years), and then loaded the radial diaphysis to failure in three-point bending after denuding each bone (38 in total). The cortical parameters had stronger correlations with ultimate moment than the trabecular parameters such that cortical area and estimated failure load of the distal radius had the highest Spearman correlation coefficients (r = 0.89 and r = 0.81, respectively, p < 0.0001). Despite being a known determinant of bone strength, cortical porosity of the distal radius did not correlate with ultimate moment (p = 0.8537). In multivariate linear regressions with section modulus (SM) of the radial diaphysis as one of two predictors of bending strength, cortical area and cortical thickness were each significant contributors to the prediction of ultimate moment. Their contribution was one-half and one-third, respectively, of the contribution from SM. None of the HR-pQCT parameters were strongly correlated with post-yield displacement, an indicator of bone brittleness. In support of HR-pQCT imaging of the distal radius to identify individuals with osteoporosis, the present study found that parameters of the cortex and failure load predictions by linear FEA are strongly related to the bending strength of cortical bone.

Incomplete recovery of bone strength and trabecular microarchitecture at the distal tibia 1 year after return from long duration spaceflight

Determining the extent of bone recovery after prolonged spaceflight is important for understanding risks to astronaut long-term skeletal health. We examined bone strength, density, and microarchitecture in seventeen astronauts (14 males; mean 47 years) using high-resolution peripheral quantitative computed tomography (HR-pQCT; 61 μm). We imaged the tibia and radius before spaceflight, at return to Earth, and after 6- and 12-months recovery and assessed biomarkers of bone turnover and exercise. Twelve months after flight, group median tibia bone strength (F.Load), total, cortical, and trabecular bone mineral density (BMD), trabecular bone volume fraction and thickness remained − 0.9% to − 2.1% reduced compared with pre-flight (p ≤ 0.001). Astronauts on longer missions (> 6-months) had poorer bone recovery. For example, F.Load recovered by 12-months post-flight in astronauts on shorter (< 6-months; − 0.4% median deficit) but not longer (− 3.9%) missions. Similar disparities were noted for total, trabecular, and cortical BMD. Altogether, nine of 17 astronauts did not fully recover tibia total BMD after 12-months. Astronauts with incomplete recovery had higher biomarkers of bone turnover compared with astronauts whose bone recovered. Study findings suggest incomplete recovery of bone strength, density, and trabecular microarchitecture at the weight-bearing tibia, commensurate with a decade or more of terrestrial age-related bone loss.

3D Image Registration Marginally Improves the Precision of HR-pQCT Measurements Compared to Cross-Sectional-Area Registration in Adults With Osteogenesis Imperfecta

Repositioning error in longitudinal high-resolution peripheral-quantitative computed tomography (HR-pQCT) imaging can lead to different bone volumes being assessed over time. To identify the same bone volumes at each time point, image registration is used. While cross-sectional area image registration corrects axial misalignment, 3D registration additionally corrects rotations. Other registration methods involving matched angle analysis (MA) or boundary transformations (3D-TB) can be used to limit interpolation error in 3D-registering micro-finite-element data. We investigated the effect of different image registration methods on short-term in vivo precision in adults with osteogenesis imperfecta, a collagen-related genetic disorder resulting in low bone mass, impaired quality, and increased fragility. The radii and tibiae of 29 participants were imaged twice on the same day with full repositioning. We compared the precision error of different image registration methods for density, microstructural, and micro-finite-element outcomes with data stratified based on anatomical site, motion status, and scanner generation. Regardless of the stratification, we found that image registration improved precision for total and trabecular bone mineral densities, trabecular and cortical bone mineral contents, area measurements, trabecular bone volume fraction, separation, and heterogeneity, as well as cortical thickness and perimeter. 3D registration marginally outperformed cross-sectional area registration for some outcomes, such as trabecular bone volume fraction and separation. Similarly, precision of micro-finite-element outcomes was improved after image registration, with 3D-TB and MA methods providing greatest improvements. Our regression model confirmed the beneficial effect of image registration on HR-pQCT precision errors, whereas motion had a detrimental effect on precision even after image registration. Collectively, our results indicate that 3D registration is recommended for longitudinal HR-pQCT imaging in adults with osteogenesis imperfecta. Since our precision errors are similar to those of healthy adults, these results can likely be extended to other populations, although future studies are needed to confirm this. © 2022 American Society for Bone and Mineral Research (ASBMR).

Response to High-Dose Vitamin D Supplementation Is Specific to Imaging Modality and Skeletal Site

High-dose vitamin D supplementation (4000 or 10,000 IU/d) in vitamin D-sufficient individuals results in a dose-dependent decrease in radius and tibia total bone mineral density (Tt.BMD) compared with 400 IU/d. This exploratory analysis examined whether the response to high-dose vitamin D supplementation depends on imaging modality and skeletal site. Participants were aged 55 to 70 years, not osteoporotic, with serum 25(OH)D 30 to 125 nM. Participants' radius and tibia were scanned on high-resolution peripheral quantitative computed tomography (HR-pQCT) to measure Tt.BMD, trabecular bone volume fraction (Tb.BV/TV), trabecular separation (Tb.Sp), cortical thickness (Ct.Th), and finite element analysis (FEA) estimated failure load. Three-dimensional image registration was used. Dual-energy X-ray absorptiometry (DXA) scans of the hip, spine, and radius measured areal BMD (aBMD) and trabecular bone score (TBS). Constrained linear mixed-effects models determined treatment group-by-time and treatment group-by-time-by-sex interactions. The treatment group-by-time interaction previously observed for radial Tt.BMD was observed at both ultradistal (UD, p < 0.001) and 33% (p < 0.001) aBMD sites. However, the treatment group-by-time-by-sex interaction observed for radial Tt.BMD was not observed with aBMD at either the UD or 33% site, and the 4000 and 400 groups did not differ. Registered radial FEA results mirrored Tt.BMD. An increase in Tb.Sp and decrease in Ct.Th underpinned dose-dependent changes in radial BMD and strength. We observed no effects in DXA-based aBMD at the hip or spine or TBS. At the tibia, we observed a time-by-treatment group effect for Tb.BV/TV. Given that DXA measures at the radius did not detect sex differences or differences between the 4000 and 400 groups, HR-pQCT at the radius may be more sensitive for examining bone changes after vitamin D supplementation. Although DXA did not reveal treatment effects at the hip or spine, whether that is a true skeletal site difference or a lack of modality sensitivity remains unclear. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Reference data and calculators for second-generation HR-pQCT measures of the radius and tibia at anatomically standardized regions in White adults

High-resolution peripheral quantitative computed tomography (HR-pQCT) is a powerful tool to assess bone health. To determine how an individual's or population of interest's HR-pQCT outcomes compare to expected, reference data are required. This study provides reference data for HR-pQCT measures acquired in a population of White adults.

PURPOSE

To provide age- and sex-specific reference data for high-resolution peripheral quantitative computed tomography (HR-pQCT) measures of the distal and diaphyseal radius and tibia acquired using a second-generation scanner and percent-of-length offsets proximal from the end of the bone.

METHODS

Data were acquired in White adults (aged 18-80 years) living in the Midwest region of the USA. HR-pQCT scans were performed at the 4% distal radius, 30% diaphyseal radius, 7.3% distal tibia, and 30% diaphyseal tibia. Centile curves were fit to the data using the LMS approach.

RESULTS

Scans of 867 females and 317 males were included. The fitted centile curves reveal HR-pQCT differences between ages, sexes, and sites. They also indicate differences when compared to data obtained by others using fixed length offsets. Excel-based calculators based on the current data were developed and are provided to enable computation of subject-specific percentiles, z-scores, and t-scores and to plot an individual's outcomes on the fitted curves. In addition, regression equations are provided to convert estimated failure load acquired with the conventional criteria utilized with first-generation scanners and those specifically developed for second-generation scanners.

CONCLUSION

The current study provides unique data and resources. The combination of the reference data and calculators provide clinicians and investigators an ability to assess HR-pQCT outcomes in an individual or population of interest, when using the described scanning and analysis procedure. Ultimately, the expectation is these data will be expanded over time so the wealth of information HR-pQCT provides becomes increasingly interpretable and utilized.

The bone strain index predicts fragility fractures. The OFELY study

Recently, the bone strain index (BSI), a new index of bone strength based on a finite element model (FEA) from dual X-ray absorptiometry (DXA), has been developed. BSI represents the average equivalent strain inside the bone, assuming that a higher strain level (high BSI) indicates a condition of higher risk. Our study aimed to analyze the relationship between BSI and age, BMI and areal BMD in pre- and postmenopausal women and to prospectively investigate fracture prediction (Fx) by BSI in postmenopausal women. Methods. At the 14th annual follow-up of the OFELY study, BSI was measured at spine (Spine BSI) and femoral scans (Neck and Total Hip BSI), in addition to areal BMD with DXA (Hologic QDR 4500) in 846 women, mean (SD) age 60 yr (15). The FRAX® (fracture risk assessment tool) for major osteoporotic fractures (MOF) was calculated with FN areal BMD (aBMD) at baseline; incident fragility fractures were annually registered until January 2016. Results. In premenopausal women (n = 261), Neck and Total Hip BSI were slightly negatively correlated with age (Spearman r = -0.13 and -0.15 respectively, p = 0.03), whereas all BSIs were positively correlated with BMI (r = +0.20 to 0.37, p < 0.01) and negatively with BMD (r = -0.69 to -0.37, p < 0.0001). In postmenopausal women (n = 585), Neck and Total Hip BSI were positively correlated with age (Spearman r = +0.26 and +0.31 respectively, p < 0.0001), whereas Spine BSI was positively correlated with BMI (r = +0.22, p < 0.0001) and all BSIs were negatively correlated with BMD (r = -0.81 to -0.60, p < 0.0001). During a median [IQ] 9.3 [1.0] years of follow-up, 133 postmenopausal women reported an incident fragility Fx, including 80 women with a major osteoporotic Fx (MOF) and 26 women with clinical vertebral Fx (VFx). Each SD increase of BSI value was associated with a significant increase of the risk of all fragility Fx with an age-adjusted HR of 1.23 for Neck BSI (p = 0.02); 1.27 for Total Hip BSI (p = 0.004) and 1.35 for Spine BSI (p < 0.0001). After adjustment for FRAX®, the association remained statistically significant for Total Hip BSI (HR 1.24, p = 0.02 for all fragility Fx; 1.31, p = 0.01 for MOF) and Spine BSI (HR 1.33, p < 0.0001 for all fragility Fx; 1.33, p = 0.005 for MOF; 1.67, p = 0.002 for clinical VFx). In conclusion, spine and femur BSI, an FEA DXA derived index, predict incident fragility fracture in postmenopausal women, regardless of FRAX®.

Bone Microarchitecture Phenotypes Identified in Older Adults Are Associated With Different Levels of Osteoporotic Fracture Risk

Prevalence of osteoporosis is more than 50% in older adults, yet current clinical methods for diagnosis that rely on areal bone mineral density (aBMD) fail to detect most individuals who have a fragility fracture. Bone fragility can manifest in different forms, and a "one-size-fits-all" approach to diagnosis and management of osteoporosis may not be suitable. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides additive information by capturing information about volumetric density and microarchitecture, but interpretation is challenging because of the complex interactions between the numerous properties measured. In this study, we propose that there are common combinations of bone properties, referred to as phenotypes, that are predisposed to different levels of fracture risk. Using HR-pQCT data from a multinational cohort (n = 5873, 71% female) between 40 and 96 years of age, we employed fuzzy c-means clustering, an unsupervised machine-learning method, to identify phenotypes of bone microarchitecture. Three clusters were identified, and using partial correlation analysis of HR-pQCT parameters, we characterized the clusters as low density, low volume, and healthy bone phenotypes. Most males were associated with the healthy bone phenotype, whereas females were more often associated with the low volume or low density bone phenotypes. Each phenotype had a significantly different cumulative hazard of major osteoporotic fracture (MOF) and of any incident osteoporotic fracture (p < 0.05). After adjustment for covariates (cohort, sex, and age), the low density followed by the low volume phenotype had the highest association with MOF (hazard ratio = 2.96 and 2.35, respectively), and significant associations were maintained when additionally adjusted for femoral neck aBMD (hazard ratio = 1.69 and 1.90, respectively). Further, within each phenotype, different imaging biomarkers of fracture were identified. These findings suggest that osteoporotic fracture risk is associated with bone phenotypes that capture key features of bone deterioration that are not distinguishable by aBMD. © 2021 American Society for Bone and Mineral Research (ASBMR).

Influence of loading conditions in finite element analysis assessed by HR-pQCT on ex vivo fracture prediction

Many fractures occur in individuals with normal areal Bone Mineral Density (aBMD) measured by Dual X-ray Absorptiometry (DXA). High Resolution peripheral Quantitative Computed Tomography (HR-pQCT) allows for non-invasive evaluation of bone stiffness and strength through micro finite element (μFE) analysis at the tibia and radius. These μFE outcomes are strongly associated with fragility fractures but do not provide clear enhancement compared with DXA measurements. The objective of this study was to establish whether a change in loading conditions in standard μFE analysis assessed by HR-pQCT enhance the discrimination of low-trauma fractured radii (n = 11) from non-fractured radii (n = 16) obtained experimentally throughout a mechanical test reproducing a forward fall. Micro finite element models were created using HR-pQCT images, and linear analyses were performed using four different types of loading conditions (axial, non-axial with two orientations and torsion). No significant differences were found between the failure load assessed with the axial and non-axial models. The different loading conditions tested presented the same area under the receiver operating characteristic (ROC) curves of 0.79 when classifying radius fractures with an accuracy of 81.5%. In comparison, the area under the curve (AUC) is 0.77 from DXA-derived ultra-distal aBMD of the forearm with an accuracy of 85.2%. These results suggest that the restricted HR-pQCT scanned region seems not sensitive to loading conditions for the prediction of radius fracture risk based on ex vivo experiments (n = 27).

Bone Phenotyping Approaches in Human, Mice and Zebrafish - Expert Overview of the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal traits TranslatiOnal NEtwork")

A synoptic overview of scientific methods applied in bone and associated research fields across species has yet to be published. Experts from the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal Traits translational Network") Working Group 2 present an overview of the routine techniques as well as clinical and research approaches employed to characterize bone phenotypes in humans and selected animal models (mice and zebrafish) of health and disease. The goal is consolidation of knowledge and a map for future research. This expert paper provides a comprehensive overview of state-of-the-art technologies to investigate bone properties in humans and animals - including their strengths and weaknesses. New research methodologies are outlined and future strategies are discussed to combine phenotypic with rapidly developing -omics data in order to advance musculoskeletal research and move towards "personalised medicine".

Bone density, microarchitecture and strength in elite figure skaters is discipline dependent

OBJECTIVES

In elite figure skaters, to determine if there was a difference in volumetric bone mineral density and bone strength between 1) figure skaters and population-based normative data, 2) single or pair skaters and ice dancers, and 3) the landing and takeoff legs.

DESIGN

Cross-sectional.

METHODS

Figure skaters had their non-dominant distal radius and bilateral tibia scanned using high-resolution peripheral quantitative computed tomography. Volumetric bone mineral density was determined at the total, cortical and trabecular compartments, and finite element analysis estimated bone strength. Normative data was used to compare the total bone mineral density of figure skaters to a population-based cohort. Independent t-tests compared differences between skating discipline, and paired t-tests compared skeletal parameters for the landing and takeoff leg.

RESULTS

Twenty elite skaters (mean age 22 ± 6.2; female = 11, male = 9) completed scans. Compared with the general population, the mean percentile rank for skaters' total volumetric bone mineral density was below normal at the radius (27th percentile) and normal at the tibia (54th percentile). Single or pair skaters had more robust bone in the landing compared with their takeoff leg. Specifically, the landing leg had higher total bone mineral density (2.8%) and trabecular bone mineral density (6.5%), and superior bone strength (8.5%) than the takeoff leg (p < 0.05).

CONCLUSIONS

Volumetric bone mineral density and strength differences in figure skaters were discipline dependent. Side-to-side differences were observed in single and pair skaters where the landing leg is denser, larger and stronger than the takeoff leg.

Using 3D image registration to maximize the reproducibility of longitudinal bone strength assessment by HR-pQCT and finite element analysis

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D image registration to account for misalignment between images. This reduced variability in longitudinal FE estimates and improved our ability to measure in vivo changes in HR-pQCT studies of bone strength.

INTRODUCTION

We developed and validated a finite element (FE) approach for longitudinal high-resolution peripheral quantitative computed tomography (HR-pQCT) studies using 3D rigid-body registration (3DR) to maximize reproducibility by accounting for misalignment between images.

METHODS

In our proposed approach, we used the full common bone volume defined by 3DR to estimate standard FE parameters. Using standard HR-pQCT imaging protocols, we validated the 3DR approach with ex vivo samples of the distal radius (n = 10, four repeat scans) by assessing whether 3DR can reduce measurement variability from repositioning error. We used in vivo data (n = 40, five longitudinal scans) to assess the sensitivity of 3DR to detect changes in bone strength at the distal radius by the standard deviation of the rate of change (σ), where the ideal value of σ is minimized to define true change. FE estimates by 3DR were compared to estimates by no registration (NR) and slice-matching (SM).

RESULTS

Group-wise comparisons of ex vivo variation (CV_RMS, %) found that FE measurement precision was improved by SM (CV_RMS < 0.80%) and 3DR (CV_RMS < 0.62%) compared to NR (CV_RMS~2%), and 3DR was advantageous as repositioning error increased. Longitudinal in vivo reproducibility was minimized by 3DR for failure load estimates (σ = 0.008 kN/month).

CONCLUSION

Although 3D registration cannot negate motion artifacts, it plays an important role in detecting and reducing variability in FE estimates for longitudinal HR-pQCT data and is well suited for estimating effects of interventions in in vivo longitudinal studies of bone strength.

Restoration of Stiffness During Fracture Healing at the Distal Radius, Using HR-pQCT and Finite Element Methods

Finite element analysis (FE) coupled with high-resolution peripheral quantitative computed tomography (HR-pQCT) allows for noninvasive in vivo assessment of fracture stiffness at peripheral locations including the distal radius. Previous studies have reported the ability of FE analysis to capture significant longitudinal changes in fracture stiffness. We hypothesized that continuum-based FE methods are necessary to capture significant changes in FE-estimated stiffness in men and women, with closed reductions and casting, over the course of their fracture healing process. The primary aim of the study was to evaluate the performance of 3 micro-FE (μFE) methods, 2 density-based (continuum) methods, and a homogeneous method. A total of 30 participants with stable distal radius fractures completed follow-ups at 2, 4, 6, 8, 12, and 26 weeks postfracture. Participants had their fractured wrist scanned using HR-pQCT at each follow-up; the contralateral wrist was also scanned at the initial assessment to represent baseline conditions. Images were used to generate continuum and homogeneous µFE models. Uniaxial compression and torsional tests were completed, with apparent stiffness determined as the primary outcome measure. Stiffness of the fractured wrist was compared to stiffness of the uninjured contralateral wrist to quantify the change in stiffness. Days since fracture significantly predicted change in stiffness for continuum and homogeneous µFE methods (p < 0.05). Continuum µFE methods appeared to account for partially mineralized tissues, resulting in a graduated recovery of stiffness (1% per week). Homogeneous µFE methods were more sensitive to stages of healing progression, resulting in a faster recovery of stiffness (3% per week). Our findings demonstrate the capability of µFE to capture the restoration of stiffness at the fractured side to prefracture stiffness in men and women, up to 6 months postfracture.

Multisite longitudinal calibration of HR-pQCT scanners and precision in osteogenesis imperfecta

BACKGROUND

For high-resolution peripheral quantitative computed tomography (HR-pQCT) to be used in longitudinal multi-center studies to assess disease and treatment effects, data must be aggregated across multiple timepoints and scanners. This requires an understanding of the factors contributing to scanner precision, and multi-scanner cross-calibration procedures, especially for clinical populations with severe phenotypes, like osteogenesis imperfecta (OI).

METHODS

To address this, we first evaluated single- and multi-center short- and long-term precision errors of standard HR-pQCT parameters. Two imaging phantoms were circulated among 13 sites (7 XtremeCT and 6 XtremeCT2) and scanned in triplicate at 3 timepoints/site. Additionally, duplicate in vivo radial and tibial scans were acquired in 29 individuals with OI. Secondly, we investigated subject- and scanner-related factors that contribute to precision errors using regression analysis. Thirdly, we proposed a reference site selection criterion for multisite cross-calibration and demonstrated the external validity of phantom-based calibrations.

RESULTS

Our results show excellent short-term single-site precision in both phantoms (CV % < 0.5%) and in density, microarchitecture and finite element parameters of OI participants (CV % = 0.75 to 1.2%). In vivo reproducibility significantly improved with (i) cross sectional area image registration versus no registration and (ii) scans with no motion artifacts. While reproducibility was similar across OI subtypes and anatomical sites, XtremeCT2 scanners achieved ~2.5% better precision than XtremeCT for trabecular parameters. Finally, we demonstrate that multisite longitudinal precision errors resulting from inconsistencies between scanners can be partially corrected through scanner cross-calibration.

CONCLUSIONS

This study is the first to assess long-term reproducibility and cross-calibration in a study using first and second generation HR-pQCT scanners. The results presented in this context provide timely guidelines for future use of this powerful clinical imaging modality in multi-center longitudinal clinical trials.

A new approach for quantifying localized bone loss by measuring void spaces

The application of high resolution peripheral quantitative computed tomography (HR-pQCT) for the study of bone health has provided valuable insight into the role bone microarchitecture has in determining bone strength and fracture risk. However, conventional density and morphological parameters struggle to distinguish whether localized bone loss is present, visible as heterogeneous deterioration in the trabecular network. This is because current HR-pQCT parameters quantify a global average of properties in the cortical or trabecular compartment. This study proposes a new metric we term "void space" that segments volumes of localized deterioration in the trabecular bone network from HR-pQCT scans and quantifies void space as the void space to total volume ratio (VS/TV, %). A simple and fully automated protocol for segmenting and quantifying void space in HR-pQCT scans is presented, along with the assessment of accuracy, precision, and cross-calibration between generations of HR-pQCT systems. Finally, prevalence of void space and the association with standard HR-pQCT parameters is demonstrated using a large population-based cohort (n = 1236). Void space detection was found to be highly reproducible (accuracy >95%, least significant change <1.76% VS/TV) and correlation between scanner generations was strong (R² = 0.87), although the first generation system struggled to identify small voids. Assessment of void space prevalence in the population-based cohort revealed that void spaces were more common in females than males, prevalence increased with age, and void spaces were typically systemic (occurring at both scan sites rather than only one). A comparison of group-wise differences between participants with and without void space demonstrated that individuals with void spaces had significantly worse trabecular properties for both sexes and at both scan sites. Specifically, the median trabecular bone mineral density, bone volume fraction, and trabecular number were below the 25th percentile of the population, while trabecular separation and inhomogeneity were above the 75th percentile of the population in participants with void spaces. Cortical bone characteristics did not differ between participants with and without void spaces. When the void space region was excluded from morphological analysis so that only the remaining "functional bone" was considered, trabecular properties of participants with void spaces were greatly improved, especially for those who were the greatest outliers. Void space is an intuitive morphological parameter that captures localized deterioration in the trabecular bone network, and has the potential to provide valuable insight into the assessment of bone fragility.

Discriminating value of HR-pQCT for fractures in women with similar FRAX scores: A substudy of the FRISBEE cohort

Areal bone mineral density (aBMD) has a low sensitivity to identify women at high fracture risk. The FRAX algorithm, by combining several clinical risk factors, might improve fracture prediction compared to aBMD alone. Several micro-architectural and biomechanical parameters which can be measured by high-resolution peripheral quantitative computed tomography (HR-pQCT) are associated with fracture risk. HR-pQCT in combination or not with finite element analysis (FEA) may be used to improve bone strength prediction. Our aim was to assess whether HR-pQCT measurements (densities, cortical and trabecular microarchitecture, biomechanical proprieties assessed by FEA) had an added value in predicting fractures in a subgroup of women belonging to the Belgian FRISBEE cohort. One hundred nineteen women who sustained a fracture (aged 60 to 85 years) during the initial follow-up of our cohort had a radius and tibia examination by HR-pQCT and were compared with controls matched for their FRAX score at baseline. We found that low distal radius total (OR = 1.41 [1.07-1.86] per SD, p < 0.05) and trabecular densities (OR = 1.45 [1.10-1.90], p < 0.01), trabecular number (OR = 1.32 [1.01-1.72], p < 0.05), intra individual distribution of separation (OR = 0.73 [0.54-0.99], p < 0.05) as several FEA parameters were significantly associated with fractures. At the distal tibia, impaired cortical density (OR = 1.32 [1.03-1.70] per SD, p < 0.05) and thickness (OR = 1.29 [1.01-1.63], p < 0.05) and apparent modulus (OR = 1.30 [1.01-1.66], p < 0.05) were significantly correlated with fractures. A low ultra distal radial aBMD (UDR) measured at the time of HR-pQCT was significantly associated with fractures (OR = 1.67 [1.22-2.28], p < 0.01). Women from both groups were followed further after the realization of the HR-pQCT and 46 new fractures were registered. In this second part of the study, low UDR aBMD (OR = 1.66 [1.18-2.35], p < 0.01), total (OR = 1.48 [1.08-2.03], p < 0.05), cortical (OR = 1.40 [1.04-1.87], p < 0.05) and trabecular (OR = 1.37 [1.01-1.85], p < 0.05) densities or apparent modulus (OR = 1.49 [1.07-2.05], p < 0.05) at the radius were associated with a significant increase of fracture risk. At the tibia, only the cortical density was significantly associated with the fracture risk (OR = 1.34 [1.02-2.76], p < 0.05). These results confirm the interest of HR-pQCT measurements for the evaluation of fracture risk, also in women matched for their baseline FRAX score. They also highlight that UDR aBMD contains pertinent information.

Improvements in radiographic and clinical assessment of distal radius fracture healing by FE-estimated bone stiffness

Bone strength determined from finite element (FE) modelling provides an estimate of fracture healing progression following a distal radius fracture (DRF), but how these measures relate to patient-reported outcomes and functional outcomes remains unknown. We hypothesized that changes in bone stiffness and bone mineral density measured using high-resolution peripheral quantitative computed tomography (HR-pQCT) are associated with clinically available measures of functional and patient-reported outcomes. We also aimed to identify which clinical outcome measures best predict fracture stiffness and could therefore be used to inform cast removal.

Participants (n = 30) with stable distal radius fractures were followed for two week intervals from the time of fracture until two months post-fracture, then at three months and six months post-fracture. At each follow-up, participants underwent clinical, radiographic, and functional assessments, as well as had their fractured wrist scanned using HR-pQCT. Recovery of bone stiffness during fracture healing was determined from micro-FE (μFE) models generated from HR-pQCT image data.

During the DRF healing process, significant longitudinal changes were found in μFE-estimated stiffness, patient-reported outcomes, grip strength, range of motion (ROM), tenderness, number of cortices healed based on radiographs, and fracture line visibility (p < 0.05); however, no significant change was detected in HR-pQCT based total bone mineral density. Patient-reported outcomes, such as the Patient-Rated Wrist Evaluation (PRWE) and the Quick Disability of the Arm, Shoulder and Hand (QuickDASH) questionnaire, correlated strongly with μFE-estimated stiffness (0.61 ≥ r_m ≥ 0.66). Based on μFE-estimated stiffness, PRWE and QuickDASH are the best predictors of stiffness recovery (p < 0.05) and may be used to guide duration of cast immobilization in the clinical setting.

•

Recovery of fracture stiffness may inform time required for cast immobilization.

•

Patient reported outcomes predict rate of fracture stiffness recovery.

•

Radiographic outcomes correlate weakly with fracture stiffness.

•

Patient reported outcomes may inform duration of cast immobilization.

Validation of HR-pQCT against micro-CT for morphometric and biomechanical analyses: A review

High-resolution peripheral quantitative computed-tomography (HR-pQCT) has the potential to become a powerful clinical assessment and diagnostic tool. Given the recent improvements in image resolution, from 82 to 61 μm, this technology may be used to accurately quantify in vivo bone microarchitecture, a key biomarker of degenerative bone diseases. However, computational methods to assess bone microarchitecture were developed for micro computed tomography (micro-CT), a higher-resolution technology only available for ex vivo studies, and validation of these computational analysis techniques against the gold-standard micro-CT has been inconsistent and incomplete. Herein, we review methods for segmentation of bone compartments and microstructure, quantification of bone morphology, and estimation of mechanical strength using finite-element analysis, highlighting the need throughout for improved standardization across the field. Studies have relied on homogenous datasets for validation, which does not allow for robust comparisons between methods. Consequently, the adaptation and validation of novel segmentation approaches has been slow to non-existent, with most studies still using the manufacturer's segmentation for morphometric analysis despite the existence of better performing alternative approaches. The promising accuracy of HR-pQCT for capturing morphometric indices is overshadowed by considerable variability in outcomes between studies. For finite element analysis (FEA) methods, the use of disparate material models and FEA tools has led to a fragmented ability to assess mechanical bone strength with HR-pQCT. Further, the scarcity of studies comparing 62 μm HR-pQCT to the gold standard micro-CT leaves the validation of this imaging modality incomplete. This review revealed that without standardization, the capabilities of HR-pQCT cannot be adequately assessed. The need for a public, extendable, heterogeneous dataset of HR-pQCT and corresponding gold-standard micro-CT images, which would allow HR-pQCT users to benchmark existing and novel methods and select optimal methods depending on the scientific question and data at hand, is now evident. With more recent advancements in HR-pQCT, the community must learn from its past and provide properly validated technologies to ensure that HR-pQCT can truly provide value in patient diagnosis and care.

Sex- and Site-Specific Reference Data for Bone Microarchitecture in Adults Measured Using Second-Generation HR-pQCT

There are currently no population-based reference data sets available for volumetric bone mineral density and microarchitecture parameters measured using the second-generation high-resolution peripheral quantitative computed tomography (HR-pQCT), yet the technology is rapidly becoming a standard for studies of bone microarchitecture. Although cross-calibrated data sets from the first-generation HR-pQCT have been reported, they are not suitable for second-generation bone microarchitecture properties because of fundamental differences between scanner generations. This study provides site- and sex-specific centile curves across the adult life span for second-generation HR-pQCT properties. A total of 1236 adult participants (768 female and 468 male) from the Calgary area between the ages of 18 and 90 years were scanned at the distal tibia and radius using the second-generation HR-pQCT. Bone densities, microarchitectural properties, and failure load estimated using finite element analysis were determined using standard in vivo protocol. Site- and sex-specific centile curves were generated using the generalized additive models for location, scale, and shape (GAMLSS) method. These data provide reference curves appropriate for predominantly white male and female adults, which can be used as a tool to assess patient- or cohort-specific bone health. © 2020 American Society for Bone and Mineral Research.

HR-pQCT in vivo imaging of periarticular bone changes in chronic inflammatory diseases: Data from acquisition to impact on treatment indications

Imaging is essential for the assessment of bone and inflammatory joint diseases. There are several imaging techniques available that differ regarding resolution, radiation exposure, time expending, precision, cost, availability or ability to predict disease progression. High-resolution peripheral quantitative computed tomography (HR-pQCT) that was introduced in 2004 allows the in vivo evaluation of peripheral bone microarchitecture and demonstrated high precision in assessing bone changes in inflammatory musculoskeletal diseases. This review summarizes the use of HR-pQCT for the evaluation of the hand skeleton in inflammatory joint diseases. We conducted a review of the literature regarding the protocols that involve hand joints assessment and evaluation of bone changes as erosions and osteophytes in chronic inflammatory diseases. Apart from measuring bone density and structure of the radius and the tibia, HR-pQCT has contributed to assessment of bone erosions and osteophytes, considered the hallmark of diseases as rheumatoid arthritis and psoriatic arthritis, respectively. In this way, there are some conventions recently established by rheumatic study groups that we just summarized here in order to standardize HR-pQCT measurements.

Guidelines for the assessment of bone density and microarchitecture in vivo using high-resolution peripheral quantitative computed tomography

INTRODUCTION

The application of high-resolution peripheral quantitative computed tomography (HR-pQCT) to assess bone microarchitecture has grown rapidly since its introduction in 2005. As the use of HR-pQCT for clinical research continues to grow, there is an urgent need to form a consensus on imaging and analysis methodologies so that studies can be appropriately compared. In addition, with the recent introduction of the second-generation HrpQCT, which differs from the first-generation HR-pQCT in scan region, resolution, and morphological measurement techniques, there is a need for guidelines on appropriate reporting of results and considerations as the field adopts newer systems.

METHODS

A joint working group between the International Osteoporosis Foundation, American Society of Bone and Mineral Research, and European Calcified Tissue Society convened in person and by teleconference over several years to produce the guidelines and recommendations presented in this document.

RESULTS

An overview and discussion is provided for (1) standardized protocol for imaging distal radius and tibia sites using HR-pQCT, with the importance of quality control and operator training discussed; (2) standardized terminology and recommendations on reporting results; (3) factors influencing accuracy and precision error, with considerations for longitudinal and multi-center study designs; and finally (4) comparison between scanner generations and other high-resolution CT systems.

CONCLUSION

This article addresses the need for standardization of HR-pQCT imaging techniques and terminology, provides guidance on interpretation and reporting of results, and discusses unresolved issues in the field.

HR-pQCT Measures of Bone Microarchitecture Predict Fracture: Systematic Review and Meta-Analysis

High-resolution peripheral quantitative computed tomography (HR-pQCT) is a noninvasive imaging modality for assessing volumetric bone mineral density (vBMD) and microarchitecture of cancellous and cortical bone. The objective was to (1) assess fracture-associated differences in HR-pQCT bone parameters; and (2) to determine if HR-pQCT is sufficiently precise to reliably detect these differences in individuals. We systematically identified 40 studies that used HR-pQCT (39/40 used XtremeCT scanners) to assess 1291 to 3253 and 3389 to 10,687 individuals with and without fractures, respectively, ranging in age from 10.9 to 84.7 years with no comorbid conditions. Parameters describing radial and tibial bone density, microarchitecture, and strength were extracted and percentage differences between fracture and control subjects were estimated using a random effects meta-analysis. An additional meta-analysis of short-term in vivo reproducibility of bone parameters assessed by XtremeCT was conducted to determine whether fracture-associated differences exceeded the least significant change (LSC) required to discern measured differences from precision error. Radial and tibial HR-pQCT parameters, including failure load, were significantly altered in fracture subjects, with differences ranging from -2.6% (95% confidence interval [CI] -3.4 to -1.9) in radial cortical vBMD to -12.6% (95% CI -15.0 to -10.3) in radial trabecular vBMD. Fracture-associated differences reported by prospective studies were consistent with those from retrospective studies, indicating that HR-pQCT can predict incident fracture. Assessment of study quality, heterogeneity, and publication biases verified the validity of these findings. Finally, we demonstrated that fracture-associated deficits in total and trabecular vBMD and certain tibial cortical parameters can be reliably discerned from HR-pQCT-related precision error and can be used to detect fracture-associated differences in individual patients. Although differences in other HR-pQCT measures, including failure load, were significantly associated with fracture, improved reproducibility is needed to ensure reliable individual cross-sectional screening and longitudinal monitoring. In conclusion, our study supports the use of HR-pQCT in clinical fracture prediction. © 2019 American Society for Bone and Mineral Research.

The Correction of Systematic Error due to Plaster and Fiberglass Casts on HR-pQCT Bone Parameters Measured In Vivo at the Distal Radius

Due to difficulty assessing healing of distal radius fractures using conventional radiography, there is interest in using high resolution peripheral quantitative computed tomography (HR-pQCT) to track healing at the microarchitectural level. Unfortunately, the plaster-of-Paris and fiberglass casts used to immobilize fractures affect HR-pQCT measurements due to beam hardening, and increased noise. The challenge is compounded because casts have variable thickness, and an individual patient will often have their cast changed 2-3 times during the course of treatment. This study quantifies the effect of casts within a clinically relevant range of thicknesses on measured bone parameters at the distal radius, and establishes conversion equations to correct for systematic error in due to cast presence. Eighteen nonfractured participants were scanned by HR-pQCT in three conditions: no cast, plaster-of-Paris cast, and fiberglass cast. Measured parameters were compared between the baseline scan (no cast) and each cast scan to evaluate if systematic error exists due to cast presence. A linear regression model was used to determine an appropriate conversion for parameters that were found to have systematic error. Plaster-of-Paris casts had a greater range of thicknesses (3.2-9.5 mm) than the fiberglass casts (3.0-5.4 mm), and induced a greater magnitude of systematic error overall. Key parameters of interest were bone mineral density (total, cortical, and trabecular) and trabecular bone volume fraction, all of which were found to have systematic error due to presence of either cast type. Linear correlations between baseline and cast scans for these parameters were excellent (R² > 0.98), and appropriate conversions could be determined within a margin of error less than a ±6% for the plaster-of-Paris cast, and ±4% for the fiberglass cast. We have demonstrated the effects of cast presence on bone microarchitecture measurements, and presented a method to correct for systematic error, in support of future use of HR-pQCT to study fracture healing.

Validation of distal radius failure load predictions by homogenized- and micro-finite element analyses based on second-generation high-resolution peripheral quantitative CT images

This study developed a well-standardized and reproducible approach for micro-finite element (mFE) and homogenized-FE (hFE) analyses that can accurately predict the distal radius failure load using either mFE or hFE models when using the approaches and parameters developed in this study.

INTRODUCTION

Micro-FE analyses based on high-resolution peripheral quantitative CT (HR-pQCT) images are frequently used to predict distal radius failure load. With the introduction of a second-generation HR-pQCT device, however, the default modelling approach no longer provides accurate results. The aim of this study was to develop a well-standardized and reproducible approach for mFE and hFE analyses that can provide precise and accurate results for distal radius failure load predictions based on second-generation HR-pQCT images.

METHODS

Second-generation HR-pQCT was used to scan the distal 20-mm section of 22 cadaver radii. The sections were excised and mechanically tested afterwards. For these sections, mFE and hFE models were made that were used to identify required material parameters by comparing predicted and measured results. Using these parameters, the models were cropped to represent the 10-mm region recommended for clinical studies to test their performance for failure load prediction.

RESULTS

After identification of material parameters, the measured failure load of the 20-mm segments was in good agreement with the results of mFE models (R² = 0.969, slope = 1.035) and hFE models (R² = 0.966, slope = 0.890). When the models were restricted to the clinical region, mFE still accurately predicted the measured failure load (R² = 0.955, slope = 1.021), while hFE predictions were precise but tended to overpredict the failure load (R² = 0.952, slope = 0.780).

CONCLUSIONS

It was concluded that it is possible to accurately predict the distal radius failure load using either mFE or hFE models when using the approaches and parameters developed in this study.

Cortical and trabecular bone microarchitecture as an independent predictor of incident fracture risk in older women and men in the Bone Microarchitecture International Consortium (BoMIC): a prospective study

BACKGROUND

Although areal bone mineral density (aBMD) assessed by dual-energy x-ray absorptiometry (DXA) is the clinical standard for determining fracture risk, most older adults who sustain a fracture have T scores greater than -2·5 and thus do not meet the clinical criteria for osteoporosis. Importantly, bone fragility is due to low BMD and deterioration in bone structure. We assessed whether indices of high-resolution peripheral quantitative CT (HR-pQCT) were associated with fracture risk independently of femoral neck aBMD and the Fracture Risk Assessment Tool (FRAX) score.

METHODS

We assessed participants in eight cohorts from the USA (Framingham, Mayo Clinic), France (QUALYOR, STRAMBO, OFELY), Switzerland (GERICO), Canada (CaMos), and Sweden (MrOS). We used Cox proportional hazard ratios (HRs) to estimate the association between HR-pQCT bone indices (per 1 SD of deficit) and incident fracture, adjusting for age, sex, height, weight, and cohort, and then additionally for femoral neck DXA aBMD or FRAX.

FINDINGS

7254 individuals (66% women and 34% men) were assessed. Mean baseline age was 69 years (SD 9, range 40-96). Over a mean follow-up of 4·63 years (SD 2·41) years, 765 (11%) participants had incident fractures, of whom 633 (86%) had femoral neck T scores greater than -2·5. After adjustment for age, sex, cohort, height, and weight, peripheral skeleton failure load had the greatest association with risk of fracture: tibia HR 2·40 (95% CI 1·98-2·91) and radius 2·13 (1·77-2·56) per 1 SD decrease. HRs for other bone indices ranged from 1·12 (95% CI 1·03-1·23) per 1 SD increase in tibia cortical porosity to 1·58 (1·45-1·72) per 1 SD decrease in radius trabecular volumetric bone density. After further adjustment for femoral neck aBMD or FRAX score, the associations were reduced but remained significant for most bone parameters. A model including cortical volumetric bone density, trabecular number, and trabecular thickness at the distal radius and a model including these indices plus cortical area at the tibia were the best predictors of fracture.

INTERPRETATION

HR-pQCT indices and failure load improved prediction of fracture beyond femoral neck aBMD or FRAX scores alone. Our findings from a large international cohort of men and women support previous reports that deficits in trabecular and cortical bone density and structure independently contribute to fracture risk. These measurements and morphological assessment of the peripheral skeleton might improve identification of people at the highest risk of fracture.

FUNDING

National Institutes of Health National Institute of Arthritis Musculoskeletal and Skin Diseases.