Local topological analysis at the distal radius by HR-pQCT: Application to in vivo bone microarchitecture and fracture assessment in the OFELY study

Fixed and Relative Positioning of Scans for High Resolution Peripheral Quantitative Computed Tomography

INTRODUCTION

High resolution peripheral quantitative computed tomography (HR-pQCT) imaging protocol requires defining where to position the ∼1 cm thick scan along the bone length. Discrepancies between the use of two positioning methods, the relative and fixed offset, may be problematic in the comparison between studies and participants. This study investigated how bone landmarks scale linearly with length and how this scaling affects both positioning methods aimed at providing a consistent anatomical location for scan acquisition.

METHODS

Using CT images of the radius (N = 25) and tibia (N = 42), 10 anatomical landmarks were selected along the bone length. The location of these landmarks was converted to a percent length along the bone, and the variation in their location was evaluated across the dataset. The absolute location of the HR-pQCT scan position using both offset methods was identified for all bones and converted to a percent length position relative to the HR-pQCT reference line for comparison. A secondary analysis of the location of the scan region specifically within the metaphysis was explored at the tibia.

RESULTS

The location of landmarks deviated from a linear relationship across the dataset, with a range of 3.6 % at the radius sites, and 4.5 % at the tibia sites. The consequent variation of the position of the scan at the radius was 0.6 % and 0.3 %, and at the tibia 2.4 % and 0.5 %, for the fixed and relative offset, respectively. The position of the metaphyseal junction with the epiphysis relative to the scan position was poorly correlated to bone length, with R2 = 0.06 and 0.37, for the fixed and relative offset respectively.

CONCLUSION

The variation of the scan position by either method is negated by the intrinsic variation of the bone anatomy with respect both to total bone length as well as the metaphyseal region. Therefore, there is no clear benefit of either offset method. However, the lack of difference due to the inherent variation in the underlying anatomy implies that it is reasonable to compare studies even if they are using different positioning methods.

The clinical application of high-resolution peripheral computed tomography (HR-pQCT) in adults: state of the art and future directions

High-resolution peripheral computed tomography (HR-pQCT) was developed to image bone microarchitecture in vivo at peripheral skeletal sites. Since the introduction of HR-pQCT in 2005, clinical research to gain insight into pathophysiology of skeletal fragility and to improve prediction of fractures has grown. Meanwhile, the second-generation HR-pQCT device has been introduced, allowing novel applications such as hand joint imaging, assessment of subchondral bone and cartilage thickness in the knee, and distal radius fracture healing. This article provides an overview of the current clinical applications and guidance on interpretation of results, as well as future directions. Specifically, we provide an overview of (1) the differences and reference data for HR-pQCT variables by age, sex, and race/ethnicity; (2) fracture risk prediction using HR-pQCT; (3) the ability to monitor response of anti-osteoporosis therapy with HR-pQCT; (4) the use of HR-pQCT in patients with metabolic bone disorders and diseases leading to secondary osteoporosis; and (5) novel applications of HR-pQCT imaging. Finally, we summarize the status of the application of HR-pQCT in clinical practice and discuss future directions. From the clinical perspective, there are both challenges and opportunities for more widespread use of HR-pQCT. Assessment of bone microarchitecture by HR-pQCT improves fracture prediction in mostly normal or osteopenic elderly subjects beyond DXA of the hip, but the added value is marginal. The prospects of HR-pQCT in clinical practice need further study with respect to medication effects, metabolic bone disorders, rare bone diseases, and other applications such as hand joint imaging and fracture healing. The mostly unexplored potential may be the differentiation of patients with only moderately low BMD but severe microstructural deterioration, which would have important implications for the decision on therapeutical interventions.

Characterising variability and regional correlations of microstructure and mechanical competence of human tibial trabecular bone: An in-vivo HR-pQCT study

OBJECTIVE

Quantifying spatial distribution of trabecular bone mechanical competence and microstructure is important for early diagnosis of skeletal disorders and potential risk of fracture. The objective of this study was to determine a spatial distribution of trabecular mechanical and morphological properties in human distal tibia and examine the contribution of regional variability of trabecular microarchitecture to mechanical competence.

METHODS

A total of 340 representative volume elements at five anatomic regions of trabecular bone - anterior, posterior, lateral, medial and centre - from ten white European-origin postmenopausal women were studied. Region-specific trabecular parameters such as trabecular volume fraction, trabecular thickness, trabecular number, trabecular surface area, trabecular separation, plate-like structure fraction and finite element analysis of trabecular stiffness were determined based on in-vivo high resolution peripheral quantitative computed tomographic (HR-pQCT) images of distal tibiae from ten postmenopausal women. Mean values were compared using analysis of variance. The correlations between morphological parameters and stiffness were calculated.

RESULTS

Significant regional variation in trabecular microarchitecture of the human distal tibia was observed (p < 0.05), with up to 106% differences between lowest (central and anterior) and highest (medial and posterior) regions. Higher proportion of plate-like trabecular morphology (63% and 53%) was found in medial and posterior regions in the distal tibia. Stiffness estimated from finite element models also differed significantly (p < 0.05), with stiffness being 4.5 times higher in the highest (medial) than lowest (central) regions. The bone volume fraction was the strongest correlate of stiffness in all regions.

CONCLUSION

A novel finding of this study is the fact that significant regional variation of stiffness derived from two-phased FEA model with individual trabecula representation correlated highly to regional morphology obtained from in-vivo HR-pQCT images at the distal tibia. The correlations between regional morphological parameters and mechanical competence of trabecular bone were consistent at all regions studied, with regional BV/TV showing the highest correlation. The method developed for regional analysis of trabecular mechanical competence may offer a better insight into the relationship between mechanical behaviour and microstructure of bone. The findings provide evidence needed to further justify a larger-cohort feasibility study for early detection of bone degenerative diseases: examining regional variations in mechanical competence and trabecular specifications may allow better understanding of fracture risks in addition to others contributing factors.

Bone Microarchitecture Assessed by HR-pQCT as Predictor of Fracture Risk in Postmenopausal Women: The OFELY Study

Several cross-sectional studies have shown that impairment of bone microarchitecture contributes to skeletal fragility. The aim of this study was to prospectively investigate the prediction of fracture (Fx) by bone microarchitecture assessed by high-resolution peripheral computed tomography (HR- pQCT) in postmenopausal women. We measured microarchitecture at the distal radius and tibia with HR-pQCT in the OFELY study, in addition to areal BMD with dual-energy X-ray absorptiometry (DXA) in 589 women, mean ± SD age 68 ± 9 years. During a median [IQ] 9.4 [1.0] years of follow-up, 135 women sustained an incident fragility Fx, including 81 women with a major osteoporotic Fx (MOP Fx). After adjustment for age, women who sustained Fx had significantly lower total and trabecular volumetric densities (vBMD) at both sites, cortical parameters (area and thickness at the radius, vBMD at the tibia), trabecular number (Tb.N), connectivity density (Conn.D), stiffness, and estimated failure load at both sites, compared with control women. After adjustment for age, current smoking, falls, prior Fx, use of osteoporosis-related drugs, and total hip BMD, each quartile decrease of several baseline values of bone microarchitecture at the radius was associated with significant change of the risk of Fx (HR of 1.39 for Tb.BMD [p = 0.001], 1.32 for Tb.N [p = 0.01], 0.76 for Tb.Sp.SD [p = 0.01], 1.49 [p = 0.01] for Conn.D, and 1.27 for stiffness [p = 0.02]). At the tibia, the association remained significant for stiffness and failure load in the multivariate model for all fragility Fx and for Tt.BMD, stiffness, and failure load for MOP Fx. We conclude that impairment of bone microarchitecture-essentially in the trabecular compartment of the radius-predict the occurrence of incident fracture in postmenopausal women. This assessment may play an important role in identifying women at high risk of fracture who could not be adequately detected by BMD measurement alone, to benefit from a therapeutic intervention. © 2017 American Society for Bone and Mineral Research.

A Comparison of Peripheral Imaging Technologies for Bone and Muscle Quantification: a Mixed Methods Clinical Review

PURPOSE OF REVIEW

Bone and muscle peripheral imaging technologies are reviewed for their association with fractures and frailty. A narrative systematized review was conducted for bone and muscle parameters from each imaging technique. In addition, meta-analyses were performed across all bone quality parameters.

RECENT FINDINGS

The current body of evidence for bone quality's association with fractures is strong for (high-resolution) peripheral quantitative computed tomography (pQCT), with trabecular separation (Tb.Sp) and integral volumetric bone mineral density (vBMD) reporting consistently large associations with various fracture types across studies. Muscle has recently been linked to fractures and frailty, but the quality of evidence remains weaker from studies of small sample sizes. It is increasingly apparent that musculoskeletal tissues have a complex relationship with interrelated clinical endpoints such as fractures and frailty. Future studies must concurrently address these relationships in order to decipher the relative importance of one causal pathway from another.

Microarchitecture and Peripheral BMD are Impaired in Postmenopausal White Women With Fracture Independently of Total Hip T-Score: An International Multicenter Study

Because single-center studies have reported conflicting associations between microarchitecture and fracture prevalence, we included high-resolution peripheral quantitative computed tomography (HR-pQCT) data from five centers worldwide into a large multicenter analysis of postmenopausal women with and without fracture. Volumetric BMD (vBMD) and microarchitecture were assessed at the distal radius and tibia in 1379 white postmenopausal women (age 67 ± 8 years); 470 (34%) had at least one fracture including 349 with a major fragility fracture. Age, height, weight, and total hip T-score differed across centers and were employed as covariates in analyses. Women with fracture had higher BMI, were older, and had lower total hip T-score, but lumbar spine T-score was similar between groups. At the radius, total and trabecular vBMD and cortical thickness were significantly lower in fractured women in three out of five centers, and trabecular number in two centers. Similar results were found at the tibia. When data from five centers were combined, however, women with fracture had significantly lower total, trabecular, and cortical vBMD (2% to 7%), lower trabecular number (4% to 5%), and thinner cortices (5% to 6%) than women without fracture after adjustment for covariates. Results were similar at the radius and tibia. Similar results were observed with analysis restricted to major fragility fracture, vertebral and hip fractures, and peripheral fracture (at the radius). When focusing on osteopenic women, each SD decrease of total and trabecular vBMD was associated with a significantly increased risk of major fragility fracture (OR = 1.55 to 1.88, p < 0.01) after adjustment for covariates. Moreover, trabecular architecture modestly improved fracture discrimination beyond peripheral total vBMD. In conclusion, we observed differences by center in the magnitude of fracture/nonfracture differences at both the distal radius and tibia. However, when data were pooled across centers and the sample size increased, we observed significant and consistent deficits in vBMD and microarchitecture independent of total hip T-score in all postmenopausal white women with fracture and in the subgroup of osteopenic women, compared to women who never had a fracture. © 2016 American Society for Bone and Mineral Research.

Age-related changes in bone strength from HR-pQCT derived microarchitectural parameters with an emphasis on the role of cortical porosity

The high resolution peripheral computed tomography (HR-pQCT) technique has seen recent developments with regard to the assessment of cortical porosity. In this study, we investigated the role of cortical porosity on bone strength in a large cohort of women. The distal radius and distal tibia were scanned by HR-pQCT. We assessed bone strength by estimating the failure load by microfinite element analysis (μFEA), with isotropic and homogeneous material properties. We built a multivariate model to predict it, using a few microarchitecture variables including cortical porosity. Among 857 Caucasian women analyzed with μFEA, we found that cortical and trabecular properties, along with the failure load, impaired slightly with advancing age in premenopausal women, the correlations with age being modest, with |rage| ranging from 0.14 to 0.38. After the onset of the menopause, those relationships with age were stronger for most parameters at both sites, with |rage| ranging from 0.10 to 0.64, notably for cortical porosity and failure load, which were markedly deteriorated with increasing age. Our multivariate model using microarchitecture parameters revealed that cortical porosity played a significant role in bone strength prediction, with semipartial r(2)=0.22 only at the tibia in postmenopausal women. In conclusion, in our large cohort of women, we observed a small decline of bone strength at the tibia before the onset of menopause. We also found an age-related increase of cortical porosity at both scanned sites in premenopausal women. In postmenopausal women, the relatively high increase of cortical porosity accounted for the decline in bone strength only at the tibia.

A survey of micro-finite element analysis for clinical assessment of bone strength: the first decade

Micro-Finite Element (micro-FE) analysis is now widely used in biomedical research as a tool to derive bone mechanical properties as they relate to its microstructure. With the development of in vivo high-resolution peripheral quantitative CT (HR-pQCT) scanners, it can now be applied to analyze bone in-vivo in the peripheral skeleton. In this survey, the results of several experimental and clinical studies are summarized that addressed the feasibility of this approach to predict bone strength in-vivo. Specific questions that will be addressed are: how accurate are strength predictions based on micro-FE; how reproducible are the results; and, is it a better predictor of bone fracture risk than DXA based measures? Based on results of experimental studies, it is first concluded that micro-FE based on HR-pQCT images can accurately predict the strength of the distal radius during a fall on the outstretched hand using either linear elastic analysis, implementing a 'Pistoia criterion' or similar criterion in combination with an 'effective' Young's modulus or using non-linear analyses. When evaluating results of clinical reproducibility studies, it is concluded that for single-center studies, errors at the radius are less than 4.4% and 3.7% and at the tibia less than 3.6% and 2.3% for stiffness and strength, respectively. In multicenter trials, however, these errors can be increased by some 1.8% and 1.4% for stiffness and strength, respectively. Finally, based on the results of large cohort studies, it is concluded that micro-FE calculated stiffness better separates cases from controls than bone density parameters for subjects with fragility fractures at any site, but not for subjects with only radius fractures. In this latter case, however, combinations of micro-FE derived parameters can significantly improve the separation.

A Digital Model to Simulate Effects of Bone Architecture Variations on Texture at Spatial Resolutions of CT, HR-pQCT, and μCT Scanners

The quantification of changes in the trabecular bone structure induced by musculoskeletal diseases like osteoarthritis, osteoporosis, rheumatoid arthritis, and others by means of a texture analysis is a valuable tool which is expected to improve the diagnosis and monitoring of a disease. The reaction of texture parameters on different alterations in the architecture of the fine trabecular network and inherent imaging factors such as spatial resolution or image noise has to be understood in detail to ensure an accurate and reliable determination of the current bone state. Therefore, a digital model for the quantitative analysis of cancellous bone structures was developed. Five parameters were used for texture analysis: entropy, global and local inhomogeneity, local anisotropy, and variogram slope. Various generic structural changes of cancellous bone were simulated for different spatial resolutions. Additionally, the dependence of the texture parameters on tissue mineralization and noise was investigated. The present work explains changes in texture parameter outcomes based on structural changes originating from structure modifications and reveals that a texture analysis could provide useful information for a trabecular bone analysis even at resolutions below the dimensions of single trabeculae.

Fast trabecular bone strength predictions of HR-pQCT and individual trabeculae segmentation-based plate and rod finite element model discriminate postmenopausal vertebral fractures

Although high-resolution peripheral quantitative computed tomography (HR-pQCT) has advanced clinical assessment of trabecular bone microstructure, nonlinear microstructural finite element (µFE) prediction of yield strength using a HR-pQCT voxel model is impractical for clinical use due to its prohibitively high computational costs. The goal of this study was to develop an efficient HR-pQCT-based plate and rod (PR) modeling technique to fill the unmet clinical need for fast bone strength estimation. By using an individual trabecula segmentation (ITS) technique to segment the trabecular structure into individual plates and rods, a patient-specific PR model was implemented by modeling each trabecular plate with multiple shell elements and each rod with a beam element. To validate this modeling technique, predictions by HR-pQCT PR model were compared with those of the registered high-resolution micro-computed tomography (HR-µCT) voxel model of 19 trabecular subvolumes from human cadaveric tibia samples. Both the Young's modulus and yield strength of HR-pQCT PR models strongly correlated with those of µCT voxel models (r²  = 0.91 and 0.86). Notably, the HR-pQCT PR models achieved major reductions in element number (>40-fold) and computer central processing unit (CPU) time (>1200-fold). Then, we applied PR model µFE analysis to HR-pQCT images of 60 postmenopausal women with (n = 30) and without (n = 30) a history of vertebral fracture. HR-pQCT PR model revealed significantly lower Young's modulus and yield strength at the radius and tibia in fracture subjects compared to controls. Moreover, these mechanical measurements remained significantly lower in fracture subjects at both sites after adjustment for areal bone mineral density (aBMD) T-score at the ultradistal radius or total hip. In conclusion, we validated a novel HR-pQCT PR model of human trabecular bone against µCT voxel models and demonstrated its ability to discriminate vertebral fracture status in postmenopausal women. This accurate nonlinear µFE prediction of the HR-pQCT PR model, which requires only seconds of desktop computer time, has tremendous promise for clinical assessment of bone strength.

High-resolution peripheral quantitative computed tomography for the assessment of bone strength and structure: a review by the Canadian Bone Strength Working Group

Bone structure is an integral determinant of bone strength. The availability of high resolution peripheral quantitative computed tomography (HR-pQCT) has made it possible to measure three-dimensional bone microarchitecture and volumetric bone mineral density in vivo, with accuracy previously unachievable and with relatively low-dose radiation. Recent studies using this novel imaging tool have increased our understanding of age-related changes and sex differences in bone microarchitecture, as well as the effect of different pharmacological therapies. One advantage of this novel tool is the use of finite element analysis modelling to non-invasively estimate bone strength and predict fractures using reconstructed three-dimensional images. In this paper, we describe the strengths and limitations of HR-pQCT and review the clinical studies using this tool.