Diabetes Risk Factors and Bone Microarchitecture as Assessed by High-Resolution Peripheral Quantitative Computed Tomography in Adults With Long-Standing Type 1 Diabetes

OBJECTIVE

To determine whether type 1 diabetes and its complications are associated with bone geometry and microarchitecture.

RESEARCH DESIGN AND METHODS

This cross-sectional study was embedded in a long-term observational study. High-resolution peripheral quantitative computed tomography (HR-pQCT) scans of the distal radius and distal and diaphyseal tibia were performed in a subset of 183 participants with type 1 diabetes from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study and 94 control participants without diabetes. HbA1c, skin advanced glycation end products (AGEs), and diabetes-related complications were assessed in EDIC participants with >30 years of follow-up.

RESULTS

Compared with control participants (aged 60 ± 8 years, 65% female), EDIC participants (aged 60 ± 7 years, diabetes duration 38 ± 5 years, 51% female) had lower total bone mineral density (BMD) at the distal radius (-7.9% [95% CI -15.2%, -0.6%]; P = 0.030) and distal tibia (-11.3% [95% CI -18.5%, -4.2%]; P = 0.001); larger total area at all sites (distal radius 4.7% [95% CI 0.5%, 8.8%; P = 0.030]; distal tibia 5.9% [95% CI 2.1%, 9.8%; P = 0.003]; diaphyseal tibia 3.4% [95% CI 0.8%, 6.1%; P = 0.011]); and poorer radius trabecular and cortical microarchitecture. Estimated failure load was similar between the two groups. Among EDIC participants, higher HbA1c, AGE levels, and macroalbuminuria were associated with lower total BMD. Macroalbuminuria was associated with larger total area and lower cortical thickness at the distal radius. Higher HbA1c and AGE levels and lower glomerular filtration rate, peripheral neuropathy, and retinopathy were associated with deficits in trabecular microarchitecture.

CONCLUSIONS

Type 1 diabetes is associated with lower BMD, larger bone area, and poorer trabecular microarchitecture. Among participants with type 1 diabetes, suboptimal glycemic control, AGE accumulation, and microvascular complications are associated with deficits in bone microarchitecture and lower BMD.

Erosion Identification in Metacarpophalangeal Joints in Rheumatoid Arthritis using High-Resolution Peripheral Quantitative Computed Tomography

Bone erosions are a pathological feature of several forms of inflammatory arthritis including rheumatoid arthritis (RA). The increased presence and size of erosions are associated with poor outcomes, joint function, and disease progression. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides unparalleled in vivo visualization of bone erosions. However, at this resolution, discontinuities in the cortical shell (cortical breaks) that are associated with normal physiological processes and pathology are also visible. The Study grouP for xtrEme Computed Tomography in Rheumatoid Arthritis previously used a consensus process to develop a definition of pathological erosion in HR-pQCT: a cortical break detected in at least two consecutive slices, in at least two perpendicular planes, non-linear in shape, with underlying trabecular bone loss. However, despite the availability of a consensus definition, erosion identification is a demanding task with challenges in inter-rater variability. The purpose of this work is to introduce a training tool to provide users with guidance on identifying pathological cortical breaks on HR-pQCT images for erosion analysis. The protocol presented here uses a custom-built module (Bone Analysis Module (BAM) - Training), implemented as an extension to an open-source image processing software (3D Slicer). Using this module, users can practice identifying erosions and compare their results to erosions annotated by expert rheumatologists.

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.

Interactions Between HR-pQCT Bone Density and D3 Cr Muscle Mass (or HR-pQCT Bone Structure and HR-pQCT Muscle Density) in Predicting Fractures: The Osteoporotic Fractures in Men Study

We examined if an interaction exists between bone and muscle in predicting fractures in older men. Prospective data from the Osteoporotic Fractures in Men study was used to build Cox proportional hazards models. Predictors included HR-pQCT total volumetric BMD (Tt.BMD), trabecular BMD (Tb.BMD), cortical BMD (Ct.BMD) and cortical area (Ct.Ar) at distal radius/tibia, HR-pQCT muscle volume and density (diaphyseal tibia), D3 -creatine dilution (D3 Cr) muscle mass, and grip strength and leg force, analyzed as continuous variables and as quartiles. Incident fractures were self-reported every 4 months via questionnaires and centrally adjudicated by physician review of radiology reports. Potential confounders (demographics, comorbidities, lifestyle factors, etc.) were considered. A total of 1353 men (mean age 84.2 ± 4.0 years, 92.7% white) were followed for 6.03 ± 2.11 years. In the unadjusted (continuous) model, there were no interactions (p > 0.05) between any muscle variable (D3 Cr muscle mass, muscle volume, muscle density, grip strength or leg force) and Tt.BMD at distal radius/tibia for fractures (all: n = 182-302; nonvertebral: n = 149-254; vertebral: n = 27-45). No consistent interactions were observed when interchanging Tt.BMD for Tb.BMD/Ct.BMD or for Ct.Ar (bone structure) at the distal radius/tibia in the unadjusted (continuous) models. Compared with men in quartiles (Q) 2-4 of D3 Cr muscle mass and Q2-4 of distal tibia Tt.BMD, men in Q1 of both had increased risk for all fractures (hazard ratio (HR) = 2.00; 95% confidence interval [CI] 1.24-3.23, p = 0.005) and nonvertebral fractures (HR = 2.10; 95% CI 1.25-3.52, p < 0.001) in the multivariable-adjusted model. Confidence intervals overlapped (p > 0.05) when visually inspecting other quartile groups in the multivariable-adjusted model. In this prospective cohort study of older men, there was no consistent interactions between bone and muscle variables on fracture risk. Larger sample sizes and longer follow-up may be needed to clarify if there is an interaction between bone and muscle on fracture risk in men. © 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).

A Laplace-Hamming Binarization Approach for Second-Generation HR-pQCT Rescues Fine Feature Segmentation

Although second-generation high-resolution peripheral quantitative computed tomography (XCTII) provides the highest-resolution in vivo bone microstructure assessment, the manufacturer's standard image processing protocol omits fine features in both trabecular and cortical compartments. To optimize fine structure segmentation, we implemented a binarization approach based on a Laplace-Hamming (LH) segmentation and documented the reproducibility and accuracy of XCTII structure segmentation using both the standard Gaussian-based binarization and the proposed LH segmentation approach. To evaluate reproducibility, 20 volunteers (9 women, 11 men; aged 23-75 years) were recruited, and three repeat scans of the radii and tibias were acquired using the manufacturer's standard in vivo protocol. To evaluate accuracy, cadaveric structure phantoms (14 radii, 6 tibias) were scanned on XCTII using the same standard in vivo protocol and on μCT at 24.5 μm resolution. XCTII images were analyzed twice-first, with the manufacturer's standard patient evaluation protocol and, second, with the proposed LH segmentation approach. The LH approach rescued fine features evident in the grayscale images but omitted or overrepresented (thickened) by the standard approach. The LH approach significantly reduced error in trabecular volume fraction (BV/TV) and thickness (Tb.Th) compared with the standard approach; however, higher error was introduced for trabecular separation (Tb.Sp). The LH approach improved the correlation between XCTII and μCT for cortical porosity (Ct.Po) and significantly reduced error in cortical pore diameter (Ct.Po.Dm) compared with the standard approach. The LH approach resulted in improved precision compared with the standard approach for BV/TV, Tb.Th, Ct.Po, and Ct.Po.Dm at the radius and for Ct.Po at the tibia. Our results suggest that the proposed LH approach produces substantially improved binary masks, reduces proportional bias, and provides greater accuracy and reproducibility in important outcome metrics, all due to more accurate segmentation of the fine features in both trabecular and cortical compartments. © 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).

Analysis of Hand Joint Space Morphology in Women and Men with Hereditary Hemochromatosis

Hereditary hemochromatosis (HH) causes unbalanced iron deposition in many organs including the joints leading to severe cartilage loss and bone damage in the metacarpophalangeal joints (MCPJ). High-resolution peripheral quantitative computed tomography (HR-pQCT) and its joint space width (JSW) quantification algorithm quantifies in vivo 3D joint morphology. We therefore aimed to (i) determine feasibility and performance of the JSW algorithm in HH, (ii) quantify joint space morphology, and (iii) investigate the relationship between morphological and clinical parameters in HH. Here, we performed an exploratory study on 24 HH patients and sex- and age-matched controls using HR-pQCT imaging of MCPJ. Mineralized bone structure was automatically segmented from the grayscale image data and periosteal surface bone masks and joint space masks were generated. Mean, minimal, and maximal joint space width (JSW; JSW.MIN; JSW.MAX), JSW heterogeneity (JSW.SD), JSW asymmetry (JSW.AS), and joint space volume (JSV) were computed. Demographics and, for HH patients, disease-specific parameters were recorded. Segmentation of JS was very good with 79.7% of MCPJs successfully segmented at first attempt and 20.3% requiring semi-manual correction. HH men showed larger JSV at all MCPs (+ 25.4% < JSV <  + 41.8%, p < 0.05), larger JSW.MAX at MCP 3-4 (+ 14%, 0.006 < p < 0.062), and wider JSW (+ 13%, p = 0.043) at MCP 4 relative to HH women. Compared to controls, both HH men and HH women showed larger JSW.AS and smaller JSW.MIN at all MCP levels, reaching significance for HH men at MCP 2 and 3 (JSW.AS: + 323% < JSW.AS <  + 359%, 0.020 < p < 0.043; JSW.MIN: - 216% < JSW.MIN < - 225%, p < 0.043), and for women at MCP 3 (JSW.AS: + 180%, p = 0.025; JSW.MIN: - 41.8%, p = 0.022). Time since HH diagnosis was correlated positively with MCP 4 JSW.AS and JSW.SD (0.463 < ρ < 0.499, p < 0.040), and the number of phlebotomies since diagnosis was correlated with JSW.SD at all MCPs (0.432 < ρ < 0.535, p < 0.050). HR-pQCT-based JSW quantification in MCPJ of HH patients is feasible, performs well even in narrow JS, and allows to define the microstructural joint burden of HH.

Biochemical Markers of Bone Turnover in Older Adults With Type 1 Diabetes

CONTEXT

Type 1 diabetes (T1D) is characterized by high fracture risk, yet little is known regarding diabetes-related mechanisms or risk factors.

OBJECTIVE

Determine whether glycemic control, advanced glycation end products (AGEs), and microvascular complications are associated with bone turnover markers among older T1D adults.

DESIGN

Cross-sectional.

SETTING

Epidemiology of Diabetes Interventions and Complications study (6 of 27 clinical centers).

PARTICIPANTS

232 T1D participants followed for >30 years.

EXPOSURES

Glycemic control ascertained as concurrent and cumulative hemoglobin A1c (HbA1c); kidney function, by estimated glomerular filtration rates (eGFR); and AGEs, by skin intrinsic fluorescence.

MAIN OUTCOME MEASURES

Serum procollagen 1 intact N-terminal propeptide (PINP), bone-specific alkaline phosphatase (bone ALP), serum C-telopeptide (sCTX), tartrate-resistant acid phosphatase 5b (TRACP5b), and sclerostin.

RESULTS

Mean age was 59.6 ± 6.8 years, and 48% were female. In models with HbA1c, eGFR, and AGEs, adjusted for age and sex, higher concurrent HbA1c was associated with lower PINP [β -3.4 pg/mL (95% CI -6.1, -0.7), P = 0.015 for each 1% higher HbA1c]. Lower eGFR was associated with higher PINP [6.9 pg/mL (95% CI 3.8, 10.0), P < 0.0001 for each -20 mL/min/1.73 m2 eGFR], bone ALP [1.0 U/L (95% CI 0.2, 1.9), P = 0.011], sCTX [53.6 pg/mL (95% CI 32.6, 74.6), P < 0.0001], and TRACP5b [0.3 U/L (95% CI 0.1, 0.4), P = 0.002]. However, AGEs were not associated with any bone turnover markers in adjusted models. HbA1c, eGFR, and AGEs were not associated with sclerostin levels.

CONCLUSIONS

Among older adults with T1D, poor glycemic control is a risk factor for reduced bone formation, while reduced kidney function is a risk factor for increased bone resorption and formation.

Augmenting Osteoporosis Imaging with Machine Learning

PURPOSE OF REVIEW

In this paper, we discuss how recent advancements in image processing and machine learning (ML) are shaping a new and exciting era for the osteoporosis imaging field. With this paper, we want to give the reader a basic exposure to the ML concepts that are necessary to build effective solutions for image processing and interpretation, while presenting an overview of the state of the art in the application of machine learning techniques for the assessment of bone structure, osteoporosis diagnosis, fracture detection, and risk prediction.

RECENT FINDINGS

ML effort in the osteoporosis imaging field is largely characterized by "low-cost" bone quality estimation and osteoporosis diagnosis, fracture detection, and risk prediction, but also automatized and standardized large-scale data analysis and data-driven imaging biomarker discovery. Our effort is not intended to be a systematic review, but an opportunity to review key studies in the recent osteoporosis imaging research landscape with the ultimate goal of discussing specific design choices, giving the reader pointers to possible solutions of regression, segmentation, and classification tasks as well as discussing common mistakes.

CT muscle density, D3Cr muscle mass and body fat associations with physical performance, mobility outcomes and mortality risk in older men

BACKGROUND

Muscle mass declines with age, while body adiposity increases. Sarcopenic obesity has been proposed to be particularly deleterious. However, previous methods for estimating muscle mass have been inadequate, and the relative contributions of total body fat vs. muscle fat to adverse outcomes have been unclear.

METHODS

In a large cohort of older men (N= 1017), we measured muscle mass (D3 creatine dilution), muscle density (high resolution peripheral computed tomography in the diaphyseal tibia) as a proxy of muscle fat, and total body fat (dual energy x-ray absorptiometry). We examined their associations with physical performance (walking speed, grip strength, chair stand time), the risk of mobility outcomes (mobility limitations, mobility disability), and the risk of death over ~5 years.

RESULTS

In combined models, lower muscle mass and muscle density were independently associated with worse physical performance and the risk of adverse outcomes, while total body fat was minimally related to physical performance and not related to mobility outcomes or mortality. For example, the relative risks for mortality per 1 standardized unit increase in muscle density, muscle mass, and total body fat were 0.84 (95% CI: 0.74, 0.70), 0.70 (0.57, 0.86), and 0.90 (0.64, 1.25), respectively.

CONCLUSIONS

Muscle mass and muscle density were associated with physical performance and adverse outcomes, and had independent, additive effects. There was little additional contribution of total body fat.

Differences in bone mineral density and morphometry measurements by fixed versus relative offset methods in high-resolution peripheral quantitative computed tomography

INTRODUCTION

High-resolution peripheral quantitative computed tomography (HR-pQCT), which enables in vivo analysis of bone morphometry, is widely used in osteoporosis research. The scan position is usually determined by the fixed offset method; however, there are concerns that the scan position can become relatively proximal if limb length is short. The present study compared bone mineral density and morphometry measured using the fixed and relative offset methods, in which the scan position is determined based on the lengths of the forearm and lower leg, and investigated factors responsible for measurement differences between the two methods.

METHODS

A total of 150 healthy Japanese subjects, comprising 75 men and 75 women, with a mean age of 45.1 years, were enrolled in this study. The distal radius and tibia were scanned using the fixed and relative offset methods; the fixed offset method involved scanning the radius and tibia at 9 mm and 22 mm, respectively, proximal to their distal articular surfaces. By contrast, the relative offset method entailed scanning the radius at 4% of the forearm length and the tibia at 7.3% of the lower leg length, proximal to their respective distal articular surfaces. The percent overlap between the scan positions of the two methods was measured using the scout views. Measurement values obtained with the two methods were compared. The correlation between the differences in the values among the two methods and forearm length, lower leg length, and body height was examined.

RESULTS

The subjects had a mean height of 164.3 ± 14.3 cm, mean forearm length of 252.9 ± 17.3 mm, and mean lower leg length of 346.7 ± 22.3 mm. The mean percent overlap was 85.0 ± 9.1% (59.2-99.6%) for the radius and 79.8 ± 12.5% (48.3-99.8%) for the tibia. Fixed offset scanning yielded higher total volumetric bone mineral density (Tt.vBMD) and cortical vBMD (Ct.vBMD) and greater cortical thickness (Ct.Th) (all p < 0.001). The differences between the two methods in terms of Tt.vBMD, Ct.vBMD and Ct.Th were significantly greater with shorter forearm length, lower leg length, and body height (radius: 0.51 < |r| < 0.63, tibia: 0.61 < |r| < 0.95).

CONCLUSION

Measurements of bone mineral density and morphometry obtained using the fixed offset method differed from those obtained using the relative offset method, which takes body size into account. Shorter body height, forearm length, and lower leg length were found to correlate with greater measurement differences. In populations with smaller stature, use of the fixed offset method results in relatively proximal images; thus, caution should be exercised when comparing groups of different height.

Longitudinal Evolution of Bone Microarchitecture and Bone Strength in Type 2 Diabetic Postmenopausal Women With and Without History of Fragility Fractures-A 5-Year Follow-Up Study Using High Resolution Peripheral Quantitative Computed Tomography

Introduction

Diabetic bone disease is characterized by an increased fracture risk which may be partly attributed to deficits in cortical bone quality such as higher cortical porosity. However, the temporal evolution of bone microarchitecture, strength, and particularly of cortical porosity in diabetic bone disease is still unknown. Here, we aimed to prospectively characterize the 5-year changes in bone microarchitecture, strength, and cortical porosity in type 2 diabetic (T2D) postmenopausal women with (DMFx) and without history of fragility fractures (DM) and to compare those to nondiabetic fracture free controls (Co) using high resolution peripheral quantitative computed tomography (HR-pQCT).

Methods

Thirty-two women underwent baseline HR-pQCT scanning of the ultradistal tibia and radius and a FU-scan 5 years later. Bone microarchitectural parameters, including cortical porosity, and bone strength estimates via µFEA were calculated for each timepoint and annualized. Linear regression models (adjusted for race and change in BMI) were used to compare the annualized percent changes in microarchitectural parameters between groups.

Results

At baseline at the tibia, DMFx subjects exhibited the highest porosity of the three groups (66.3% greater Ct.Po, 71.9% higher Ct.Po.Volume than DM subjects, p < 0.022). Longitudinally, porosity increased significantly over time in all three groups and at similar annual rates, while DMFx exhibited the greatest annual decreases in bone strength indices (compared to DM 4.7× and 6.7× greater decreases in failure load [F] and stiffness [K], p < 0.025; compared to Co 14.1× and 22.2× greater decreases in F and K, p < 0.020).

Conclusion

Our data suggest that despite different baseline levels in cortical porosity, T2D women with and without fractures experienced long-term porosity increases at a rate similar to non-diabetics. However, the annual loss in bone strength was greatest in T2D women with a history of a fragility fractures. This suggests a potentially non-linear course of cortical porosity development in T2D bone disease: major porosity may develop early in the course of disease, followed by a smaller steady annual increase in porosity which in turn can still have a detrimental effect on bone strength-depending on the amount of early cortical pre-damage.

Global and Spatial Compartmental Interrelationships of Bone Density, Microstructure, Geometry and Biomechanics in the Distal Radius in a Colles' Fracture Study Using HR-pQCT

BACKGROUND

Bone parameters derived from HR-pQCT have been investigated on a parameter-by-parameter basis for different clinical conditions. However, little is known regarding the interrelationships of bone parameters and the spatial distribution of these interrelationships. In this work: 1) we investigate compartmental interrelationships of bone parameters; 2) assess the spatial distribution of interrelationships of bone parameters; and 3) compare interrelationships of bone parameters between postmenopausal women with and without a recent Colles' fracture.

METHODS

Images from the unaffected radius in fracture cases (n=84), and from the non-dominant radius of controls (n=98) were obtained using HR-pQCT. Trabecular voxel-based maps of local bone volume fraction (L.Tb.BV/TV), homogenized volumetric bone mineral density (H.Tb.BMD), homogenized μFEA-derived strain energy density (H.Tb.SED), and homogenized inter-trabecular distances (H.Tb.1/N) were generated; as well as surface-based maps of apparent cortical bone thickness (Surf.app.Ct.Th), porosity-weighted cortical bone thickness (Surf.Ct.SIT), mean cortical BMD (Surf.Ct.BMD), and mean cortical SED (Surf.Ct.SED). Anatomical correspondences across the parametric maps in the study were established via spatial normalization to a common template. Mean values of the parametric maps before spatial normalization were used to assess compartmental Spearman's rank partial correlations of bone parameters (e.g., between H.Tb.BMD and L.Tb.BV/TV or between Surf.Ct.BMD and Surf.app.Ct.Th). Spearman's rank partial correlations were also assessed for each voxel and vertex of the spatially normalized parametric maps, thus generating maps of Spearman's rank partial correlation coefficients. Correlations were performed independently within each group, and compared between groups using the Fisher's Z transformation.

RESULTS

All within-group global trabecular and cortical Spearman's rank partial correlations were significant; and the correlations of H.Tb.BMD-L.Tb.BV/TV, H.Tb.BMD-H.Tb.1/N, L.Tb.BV/TV-H.Tb.1/N, Surf.Ct.BMD-Surf.Ct.SED and Surf.Ct.SIT-Surf.Ct.SED were significantly different between controls and fracture cases. The spatial analyses revealed significant heterogeneous voxel- and surface-based correlation coefficient maps across the distal radius for both groups; and the correlation maps of H.Tb.BMD-L.Tb.BV/TV, H.Tb.BMD-H.Tb.1/N, L.Tb.BV/TV-H.Tb.1/N, H.Tb.1/N-H.Tb.SED and Surf.app.Ct.Th - Surf.Ct.SIT yielded small clusters of significant correlation differences between groups.

DISCUSSION

The heterogeneous spatial distribution of interrelationships of bone parameters assessing density, microstructure, geometry and biomechanics, along with their global and local differences between controls and fracture cases, may help us further understand different bone mechanisms of bone fracture.

Reliability and Change in Erosion Measurements by High-resolution Peripheral Quantitative Computed Tomography in a Longitudinal Dataset of Rheumatoid Arthritis Patients

OBJECTIVE

The aim of this multireader exercise was to assess the reliability and change over time of erosion measurements in patients with rheumatoid arthritis (RA) using high-resolution peripheral quantitative computed tomography (HR-pQCT).

METHODS

HR-pQCT scans of 23 patients with RA were assessed at baseline and 12 months. Four experienced readers examined the dorsal, palmar, radial, and ulnar surfaces of the metacarpal head (MH) and phalangeal base (PB) of the second and third digits, blinded to time order. In total, 368 surfaces (23 patients´ 16 surfaces) were evaluated per timepoint to characterize cortical breaks as pathological (erosion) or physiological, and to quantify erosion width and depth. Reliability was evaluated by intraclass correlation coefficients (ICC), percentage agreement, and Light k; change over time was defined by means ± SD of erosion numbers and dimensions.

RESULTS

ICC for the mean measurements of width and depth of the pathological breaks ranged between 0.819-0.883, and 0.771-0.907, respectively. Most physiological cortical breaks were found at the palmar PB, whereas most pathological cortical breaks were located at the radial MH. There was a significant increase in both the numbers and the dimensions of erosions between baseline and follow-up (P = 0.0001 for erosion numbers, width, and depth in axial plane; P = 0.001 for depth in perpendicular plane).

CONCLUSION

This exercise confirmed good reliability of HR-pQCT erosion measurements and their ability to detect change over time.

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.

High-Resolution Peripheral Quantitative Computed Tomography for Bone Evaluation in Inflammatory Rheumatic Disease

High resolution peripheral quantitative computed tomography (HR-pQCT) is a 3-dimensional imaging modality with superior sensitivity for bone changes and abnormalities. Recent advances have led to increased use of HR-pQCT in inflammatory arthritis to report quantitative volumetric measures of bone density, microstructure, local anabolic (e.g., osteophytes, enthesiophytes) and catabolic (e.g., erosions) bone changes and joint space width. These features may be useful for monitoring disease progression, response to therapy, and are responsive to differentiating between those with inflammatory arthritis conditions and healthy controls. We reviewed 69 publications utilizing HR-pQCT imaging of the metacarpophalangeal (MCP) and/or wrist joints to investigate arthritis conditions. Erosions are a marker of early inflammatory arthritis progression, and recent work has focused on improvement and application of techniques to sensitively identify erosions, as well as quantifying erosion volume changes longitudinally using manual, semi-automated and automated methods. As a research tool, HR-pQCT may be used to detect treatment effects through changes in erosion volume in as little as 3 months. Studies with 1-year follow-up have demonstrated progression or repair of erosions depending on the treatment strategy applied. HR-pQCT presents several advantages. Combined with advances in image processing and image registration, individual changes can be monitored with high sensitivity and reliability. Thus, a major strength of HR-pQCT is its applicability in instances where subtle changes are anticipated, such as early erosive progression in the presence of subclinical inflammation. HR-pQCT imaging results could ultimately impact decision making to uptake aggressive treatment strategies and prevent progression of joint damage. There are several potential areas where HR-pQCT evaluation of inflammatory arthritis still requires development. As a highly sensitive imaging technique, one of the major challenges has been motion artifacts; motion compensation algorithms should be implemented for HR-pQCT. New research developments will improve the current disadvantages including, wider availability of scanners, the field of view, as well as the versatility for measuring tissues other than only bone. The challenge remains to disseminate these analysis approaches for broader clinical use and in research.

Objective measures of moderate to vigorous physical activity are associated with higher distal limb bone strength among elderly men

Our aim was to determine the association between objectively measured physical activity (PA) and bone strength of the distal limbs among older men. We studied 994 men from the MrOS cohort study (mean age 83.9) who had repeat (Year 7 and 14) 5-day activity assessment with at least 90% wear time (SenseWearPro3 Armband) and Year 14 measures using high resolution peripheral quantitative tomography (HR-pQCT) (Scanco). Total energy expenditure (TEE), total steps per day, peak cadence (mean of top 30 steps/min over 24 h) and time spent in a given level of activity: sedentary (reference, <1.5 metabolic equivalents of task [METs]), light (1.5 to <3 METs), or moderate to vigorous physical activity(MVPA: ≥3 METs) were calculated as mean over the two time points. Estimated failure load was determined from HR-pQCT data using finite element analysis. We used standardized variables and adjusted for potential confounders using linear regression. The means ±SDs for daily activity were: 2338 ± 356 kcal/d [TEE]; 5739 ± 2696 steps/day [step count], 60 ± 20 cpm [peak cadence], 67 ± 28 min/d [light activity], and 85 ± 52 min/d [MVPA]. Higher TEE, step count, and peak cadence were each associated with higher failure load of the distal radius (effect sizes respectively: 0.13 [95% CI: 0.05, 0.20], 0.11 [95% CI: 0.04, 0.18], and 0.08 [95% CI: 0.01, 0.15]) and higher failure load of the distal tibia (effect sizes respectively 0.21 [95% CI: 0.13, 0.28], 0.19 [95% CI: 0.13, 0.26], 0.19 [95% CI, 0.13, 0.25]). Time spent in MVPA vs. time sedentary was related to bone strength at both sites after adjustment, whereas time spent in light activity vs. time sedentary was not. TEE was associated with compartmental area and BMD parameters at distal tibia, but only area parameters at the distal radius. In summary, MVPA over a 7-year period of time may have a modest association with bone strength and geometry among older men.

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.

Consensus approach for 3D joint space width of metacarpophalangeal joints of rheumatoid arthritis patients using high-resolution peripheral quantitative computed tomography

Background

Joint space assessment for rheumatoid arthritis (RA) by ordinal conventional radiographic scales is susceptible to floor and ceiling effects. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides superior resolution, and may detect earlier changes. The goal of this work was to compare existing 3D methods to calculate joint space width (JSW) metrics in human metacarpophalangeal (MCP) joints with HR-pQCT and reach consensus for future studies. Using the consensus method, we established reproducibility with repositioning as well as feasibility for use in second-generation HR-pQCT scanners.

Methods

Three published JSW methods were compared using datasets from individuals with RA from three research centers. A SPECTRA consensus method was developed to take advantage of strengths of the individual methods. Using the SPECTRA method, reproducibility after repositioning was tested and agreement between scanner generations was also established.

Results

When comparing existing JSW methods, excellent agreement was shown for JSW minimum and mean (ICC 0.987-0.996) but not maximum and volume (ICC 0.000-0.897). Differences were identified as variations in volume definitions and algorithmic differences that generated high sensitivity to boundary conditions. The SPECTRA consensus method reduced this sensitivity, demonstrating good scan-rescan reliability (ICC >0.911) except for minimum JSW (ICC 0.656). There was strong agreement between results from first- and second-generation HR-pQCT (ICC >0.833).

Conclusions

The SPECTRA consensus method combines unique strengths of three independently-developed algorithms and leverages underlying software updates to provide a mature analysis to measure 3D JSW. This method is robust with respect to repositioning and scanner generations, suggesting its suitability for detecting change.

Hip Fracture Discrimination Based on Statistical Multi-parametric Modeling (SMPM)

Studies using quantitative computed tomography (QCT) and data-driven image analysis techniques have shown that trabecular and cortical volumetric bone mineral density (vBMD) can improve the hip fracture prediction of dual-energy X-ray absorptiometry areal BMD (aBMD). Here, we hypothesize that (1) QCT imaging features of shape, density and structure derived from data-driven image analysis techniques can improve the hip fracture discrimination of classification models based on mean femoral neck aBMD (Neck.aBMD), and (2) that data-driven cortical bone thickness (Ct.Th) features can improve the hip fracture discrimination of vBMD models. We tested our hypotheses using statistical multi-parametric modeling (SMPM) in a QCT study of acute hip fracture of 50 controls and 93 fragility fracture cases. SMPM was used to extract features of shape, vBMD, Ct.Th, cortical vBMD, and vBMD in a layer adjacent to the endosteal surface to develop hip fracture classification models with machine learning logistic LASSO. The performance of these classification models was evaluated in two aspects: (1) their hip fracture classification capability without Neck.aBMD, and (2) their capability to improve the hip fracture classification of the Neck.aBMD model. Assessments were done with 10-fold cross-validation, areas under the receiver operating characteristic curve (AUCs), differences of AUCs, and the integrated discrimination improvement (IDI) index. All LASSO models including SMPM-vBMD features, and the majority of models including SMPM-Ct.Th features performed significantly better than the Neck.aBMD model; and all SMPM features significantly improved the hip fracture discrimination of the Neck.aBMD model (Hypothesis 1). An interesting finding was that SMPM-features of vBMD also captured Ct.Th patterns, potentially explaining the superior classification performance of models based on SMPM-vBMD features (Hypothesis 2). Age, height and weight had a small impact on model performances, and the model of shape, vBMD and Ct.Th consistently yielded better performances than the Neck.aBMD models. Results of this study clearly support the relevance of bone density and quality on the assessment of hip fracture, and demonstrate their potential on patient and healthcare cost benefits.

Structural Changes over a Short Period Are Associated with Functional Assessments in Rheumatoid Arthritis

Objective.

To investigate the correlation between changes in radiological quantitative assessment with changes in clinical and functional assessment from baseline to 3 months in patients with rheumatoid arthritis (RA).

Methods.

Twenty-eight patients with RA [methotrexate (MTX) and anti-tumor necrosis factor–α (TNF-α) group with high disease activity (n = 18); and MTX group with low disease activity (n = 10)] underwent assessments at baseline and 3 months: clinical [28-joint count Disease Activity Score (DAS28)], functional [Health Assessment Questionnaire (HAQ) and Michigan Hand Outcome Questionnaire (MHQ)], and imaging-based [3 Tesla magnetic resonance imaging (MRI) and high-resolution peripheral quantitative computed tomography (HR-pQCT)]. MR images were evaluated semiquantitatively [RA MRI scoring (RAMRIS)] and quantitatively for the volume of synovitis and bone marrow edema (BME) lesions. Erosion volumes were measured using HR-pQCT.

Results.

After 3 months, the anti-TNF-α group demonstrated an improvement in disease activity through DAS28, HAQ, and MHQ. MRI showed significant decreases in synovitis and BME volume for the anti-TNF-α group, and significant increases in the MTX group. HR-pQCT showed significant decreases in bone erosion volume for the anti-TNF-α group, and significant increases in the MTX group. No significance was observed using RAMRIS. Changes in synovitis, BME, and erosion volumes, but not RAMRIS, were significantly correlated with changes in DAS28, HAQ, and MHQ.

Conclusion.

Quantitative measures were more sensitive than semiquantitative grading when evaluating structural and inflammatory changes with treatment, and were associated with patient clinical and functional outcomes. Multimodality imaging with 3T MRI and HR-pQCT may provide promising biomarkers that help determine disease progression and therapy response.

Accelerated Bone Loss in Older Men: Effects on Bone Microarchitecture and Strength

Accelerated bone loss (ABL) shown on routine dual-energy X-ray absorptiometry (DXA) may be accompanied by microarchitectural changes, increased cortical porosity, and lower bone strength. To test this hypothesis, we performed a cross-sectional study and used high-resolution peripheral quantitative computed tomography (HR-pQCT) scans (Scanco Medical AG, Brüttisellen, Switzerland) to measure estimated bone strength and microarchitecture in the distal radius and distal and diaphyseal tibia. We studied 1628 men who attended the year 14 exam of the Osteoporotic Fractures in Men (MrOS) study. We retrospectively characterized areal bone mineral density (aBMD) change from the year 7 to year 14 exam in three categories: "accelerated" loss, ≥10% loss at either the total hip or femoral neck (n = 299, 18.4%); "expected" loss, <10% (n = 1061, 65.2%), and "maintained" BMD, ≥0% (n = 268, 16.5%). The ABL cut-off was a safety alert established for MrOS. We used regression models to calculate adjusted mean HR-pQCT parameters in men with ABL, expected loss, or maintained BMD. Men who experienced ABL were older and had a lower body mass index and aBMD and experienced greater weight loss compared with other men. Total volumetric BMD and trabecular and cortical volumetric BMD were lower in men with ABL compared with the expected or maintained group. Men with ABL had significantly lower trabecular bone volume fraction (BV/TV), fewer trabeculae, and greater trabecular separation at both the distal radius and tibia than men with expected loss or who maintained aBMD, all p trend <0.001. Men with ABL had lower cortical thickness and lower estimated bone strength, but there was no difference in cortical porosity except at the tibia diaphyseal site. In summary, men with ABL have lower estimated bone strength, poorer trabecular microarchitecture, and thinner cortices than men without ABL but have similar cortical porosity. These impairments may lead to an increased risk of fracture. © 2018 American Society for Bone and Mineral Research.

Weight loss in men in late life and bone strength and microarchitecture: a prospective study

Weight loss in men in late life was associated with lower bone strength. In contrast, weight gain was not associated with a commensurate increase in bone strength. Future studies should measure concurrent changes in weight and parameters of bone strength and microarchitecture and evaluate potential causal pathways underlying these associations.

INTRODUCTION

Our aim was to determine associations of weight loss with bone strength and microarchitecture.

METHODS

We used data from 1723 community-dwelling men (mean age 84.5 years) who attended the MrOS study Year (Y) 14 exam and had high-resolution peripheral quantitative computed tomography (HR-pQCT) scans at ≥ 1 skeletal sites (distal tibia, distal radius, or diaphyseal tibia). Weight change from Y7 to Y14 exams (mean 7.3 years between exams) was classified as moderate weight loss (loss ≥ 10%), mild weight loss (loss 5 to < 10%), stable weight (< 5% change), or weight gain (gain ≥ 5%). Mean HR-pQCT parameters (95%CI) were calculated by weight change category using linear regression models adjusted for age, race, site, health status, body mass index, limb length, and physical activity. The primary outcome measure was estimated failure load.

RESULTS

There was a nonlinear association of weight change with failure load at each skeletal site with different associations for weight loss vs. weight gain (p < 0.03). Failure load and total bone mineral density (BMD) at distal sites were lower with greater weight loss with 7.0-7.6% lower failure loads and 4.3-5.8% lower BMDs among men with moderate weight loss compared to those with stable weight (p < 0.01, both comparisons). Cortical, but not trabecular, BMDs at distal sites were lower with greater weight loss. Greater weight loss was associated with lower cortical thickness at all three skeletal sites.

CONCLUSION

Weight loss in men in late life is associated with lower peripheral bone strength and total BMD with global measures reflecting cortical but not trabecular parameters.

Volumetric Bone Mineral Density and Failure Load of Distal Limbs Predict Incident Clinical Fracture Independent HR-pQCT BMD and Failure Load Predicts Incident Clinical Fracture of FRAX and Clinical Risk Factors Among Older Men

Our objective was to determine the associations of peripheral bone strength and microarchitecture with incident clinical and major osteoporotic fracture among older men after adjusting for major clinical risk factors. We used a prospective cohort study design with data from 1794 men (mean age 84.4 years) in the Osteoporotic Fractures in Men (MrOS) study. Eligible men attended the year 14 visit, had high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of the distal radius and distal or diaphyseal tibia, DXA measured BMD, and were followed for mean 1.7 years for incident fracture. Failure load was estimated using finite element analysis. We used Cox proportional hazards models with standardized HR-pQCT parameters as exposure variables. Primary outcome was clinical fracture (n = 108). Covariates included either Fracture Risk Assessment Tool (FRAX) major osteoporotic fracture probability calculated with BMD (FRAX-BMD), or individual clinical risk factors (CRF) including age, total hip BMD, race, falls, and prevalent fracture after age 50 years. Lower failure load was associated with higher risk of incident clinical fracture and incident major osteoporotic fracture. For clinical fracture with FRAX-BMD adjustment, the associations ranged from hazard ratio (HR) 1.58 (95% CI, 1.25 to 2.01) to 2.06 (95% CI, 1.60 to 2.66) per SD lower failure load at the diaphyseal tibia and distal radius. These associations were attenuated after adjustment for individual CRFs, but remained significant at the distal sites. Associations of volumetric BMD with these outcomes were similar to those for failure load. At the distal radius, lower trabecular BMD, number, and thickness, and lower cortical BMD, thickness, and area were all associated with higher risk of clinical fracture, but cortical porosity was not. Among community-dwelling older men, HR-pQCT measures including failure load, volumetric BMD, and microstructure parameters at peripheral sites (particularly distal radius) are robust independent predictors of clinical and major osteoporotic fracture. © 2018 American Society for Bone and Mineral Research.

High dairy protein intake is associated with greater bone strength parameters at the distal radius and tibia in older men: a cross-sectional study

Dairy protein but not plant protein was associated with bone strength of the radius and tibia in older men. These results are consistent with previous results in women and support similar findings related to fracture outcomes. Bone strength differences were largely due to thickness and area of the bone cortex.

INTRODUCTION

Our objective was to determine the association of protein intake by source (dairy, non-dairy animal, plant) with bone strength and bone microarchitecture among older men.

METHODS

We used data from 1016 men (mean 84.3 years) who attended the Year 14 exam of the Osteoporotic Fractures in Men (MrOS) study, completed a food frequency questionnaire (500-5000 kcal/day), were not taking androgen or androgen agonists, and had high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of the distal radius and distal or diaphyseal tibia. Protein was expressed as percentage of total energy intake (TEI); mean ± SD for TEI = 1548 ± 607 kcal/day and for total protein = 16.2 ± 2.9%TEI. We used linear regression with standardized HR-pQCT parameters as dependent variables and adjusted for age, limb length, center, education, race/ethnicity, marital status, smoking, alcohol intake, physical activity level, corticosteroids use, supplement use (calcium and vitamin D), and osteoporosis medications.

RESULTS

Higher dairy protein intake was associated with higher estimated failure load at the distal radius and distal tibia [radius effect size = 0.17 (95% CI 0.07, 0.27), tibia effect size = 0.13 (95% CI 0.03, 0.23)], while higher non-dairy animal protein was associated with higher failure load at only the distal radius. Plant protein intake was not associated with failure load at any site.

CONCLUSION

The association between protein intake and bone strength varied by source of protein. These results support a link between dairy protein intake and skeletal health, but an intervention study is needed to evaluate causality.

Assessment of 3-month changes in bone microstructure under anti-TNFα therapy in patients with rheumatoid arthritis using high-resolution peripheral quantitative computed tomography (HR-pQCT)

Background

Although one study showed minimal progression of erosions in patients with rheumatoid arthritis (RA) one year after TNFα inhibition therapy, no studies have investigated very early bone changes after initiation of anti-TNFα treatment. We investigated the effects of 3-month anti-TNFα treatment on bone erosion progression and bone microarchitecture in RA patients using high-resolution peripheral quantitative computed tomography (HR-pQCT).

Methods

Patients with RA (n = 27) (17 in the anti-TNFα and 10 in the MTX-only group) underwent assessment of disease activity score in 28 joints (DAS-28), radiographs, 3-T magnetic resonance imaging (MRI) and HR-pQCT of metacarpophalangeal and wrist joints at baseline and 3 months. HR-pQCT-derived erosion volume, joint volume/width and bone microarchitecture were computed and joint destruction was assessed using Sharp and RAMRIS scorings on radiographs and MRI, respectively.

Results

Overall, 73 erosions were identified by HR-pQCT at baseline. Over 3 months, the anti-TNFα group had decreased mean erosion volume; increased erosion volume was observed in one clinical non-responder. The MTX-only group in contrast, trended toward increasing erosion volume despite low disease activity. In the anti-TNFα group, joint-space width and volume of MCP joints decreased significantly and was positively correlated with erosion volume changes (R² = 0.311, p = 0.013; R² = 0.527, p = 0.003, respectively). In addition, erosion volume changes were significantly negatively correlated with changes in trabecular bone mineral density (R² = 0.353, p = 0.020) in this group. We observed significant correlation between percentage change in erosion volume and change in DAS-28 erythrocyte sedimentation rate and C-reactive protein CRP scores (R² = 0.558, p < 0.001; R² = 0.745, p < 0.001, respectively) in all patients.

Conclusions

Using HR-pQCT, our data suggest that anti-TNFα treatment prevents erosion progression and deterioration of bone microarchitecture within the first 3 months of treatment, one patient not responding to treatment, had significant progression of bone erosions within this short time period. Patients with low disease activity scores (<3.2) can have continuous HR-pQCT-detectable progression of erosive disease with MTX treatment only. HR-pQCT can be a sensitive, powerful tool to quantify bone changes and monitor RA treatment short term (such as 3 months).

The SPECTRA Collaboration OMERACT Special Interest Group: Current Research and Future Directions

OBJECTIVE

High-resolution peripheral quantitative computed tomography (HR-pQCT) has the potential to improve radiographic progression determination in clinical trials and longitudinal observational studies. The goal of this work was to describe the current state of research presented at Outcome Measures in Rheumatology (OMERACT) 2016 and ensuing future directions outlined during discussion among attendees.

METHODS

At OMERACT 2016, SPECTRA (Study grouP for xtrEme-Computed Tomography in Rheumatoid Arthritis) introduced efforts to (1) validate the HR-pQCT according to OMERACT guidelines, focusing on rheumatoid arthritis (RA), and (2) find alternatives for automated joint space width (JSW) analysis. The Special Interest Group (SIG) was presented to patient research partners, physicians/researchers, and SIG leaders followed by a 40-min discussion on future directions.

RESULTS

A consensus definition for RA erosion using HR-pQCT was demonstrated through a systematic literature review and a Delphi exercise. Histopathology and perfusion studies were presented that analyzed the true characteristics of cortical breaks in HR-pQCT images, and to provide criterion validity. Results indicate that readers were able to discriminate between erosion and small vascular channels. Moderate reliability (ICC 0.206-0.871) of direct erosion size measures was shown, which improved (> 0.9) only when experienced readers were considered. Quantification of erosion size was presented for scoring, direct measurement, and volumetric approaches, as well as a reliability exercise for direct measurement. Three methods for JSW measurement were compared, all indicating excellent reproducibility with differences at the extremes (i.e., near-zero and joint edge thickness).

CONCLUSION

Initial reports on HR-pQCT are promising; however, to consider its use in clinical trials and longitudinal observational studies, it is imperative to assess the responsiveness of erosion measurement quantification.

The comparability of HR-pQCT bone measurements is improved by scanning anatomically standardized regions

We investigated the sensitivity of distal bone density, structure, and strength measurements by high-resolution peripheral quantitative computed tomography (HR-pQCT) to variability in limb length. Our results demonstrate that HR-pQCT should be performed at a standard %-of-total-limb-length to avoid substantial measurement bias in population study comparisons and the evaluation of individual skeletal status in a clinical context.

INTRODUCTION

High-resolution peripheral quantitative computed tomography (HR-pQCT) measures of bone do not account for anatomic variability in bone length: a 1-cm volume is acquired at a fixed offset from an anatomic landmark. Our goal was to evaluate HR-pQCT measurement variability introduced by imaging fixed vs. proportional volumes and to propose a standard protocol for relative anatomic positioning.

METHODS

Double-length (2-cm) scans were acquired in 30 adults. We compared measurements from 1-cm sub-volumes located at the default fixed offset, and the average %-of-length offset. The average position corresponded to 4.0% ± 1.1 mm for radius, and 7.2% ± 2.2 mm for tibia. We calculated the RMS difference in bone parameters and T-scores to determine the measurement variability related to differences in limb length. We used anthropometric ratios to estimate the mean limb length for published HR-pQCT reference data, and then calculated mean %-of-length offsets.

RESULTS

Variability between fixed vs. relative scan positions was highest in the radius, and for cortical bone in general (RMS difference Ct.Th = 19.5%), while individuals had T-score differentials as high as +3.0 SD (radius Ct.BMD). We estimated that average scan position for published HR-pQCT reference data corresponded to 4.0% at the radius, and 7.3% at tibia.

CONCLUSION

Variability in limb length introduces significant bias to HR-pQCT measures, confounding cross-sectional analyses and limiting the clinical application for individual assessment of skeletal status. We propose to standardize scan positioning using 4.0 and 7.3% of total bone length for the distal radius and tibia, respectively.

Operator variability in scan positioning is a major component of HR-pQCT precision error and is reduced by standardized training

In this study, we determined that operator positioning precision contributes significant measurement error in high-resolution peripheral quantitative computed tomography (HR-pQCT). Moreover, we developed software to quantify intra- and inter-operator variability and demonstrated that standard positioning training (now available as a web-based application) can significantly reduce inter-operator variability.

INTRODUCTION

HR-pQCT is increasingly used to assess bone quality, fracture risk, and anti-fracture interventions. The contribution of the operator has not been adequately accounted in measurement precision. Operators acquire a 2D projection ("scout view image") and define the region to be scanned by positioning a "reference line" on a standard anatomical landmark. In this study, we (i) evaluated the contribution of positioning variability to in vivo measurement precision, (ii) measured intra- and inter-operator positioning variability, and (iii) tested if custom training software led to superior reproducibility in new operators compared to experienced operators.

METHODS

To evaluate the operator in vivo measurement precision, we compared precision errors calculated in 64 co-registered and non-co-registered scan-rescan images. To quantify operator variability, we developed software that simulates the positioning process of the scanner's software. Eight experienced operators positioned reference lines on scout view images designed to test intra- and inter-operator reproducibility. Finally, we developed modules for training and evaluation of reference line positioning. We enrolled six new operators to participate in a common training, followed by the same reproducibility experiments performed by the experienced group.

RESULTS

In vivo precision errors were up to threefold greater (Tt.BMD and Ct.Th) when variability in scan positioning was included. The inter-operator precision errors were significantly greater than the short-term intra-operator precision (p < 0.001). New trained operators achieved comparable intra-operator reproducibility to experienced operators and lower inter-operator reproducibility (p < 0.001). Precision errors were significantly greater for the radius than for the tibia.

CONCLUSION

Operator reference line positioning contributes significantly to in vivo measurement precision and is significantly greater for multi-operator datasets. Inter-operator variability can be significantly reduced using a systematic training platform, now available online ( http://webapps.radiology.ucsf.edu/refline/ ).

Serum miRNA Signatures Are Indicative of Skeletal Fractures in Postmenopausal Women With and Without Type 2 Diabetes and Influence Osteogenic and Adipogenic Differentiation of Adipose Tissue-Derived Mesenchymal Stem Cells In Vitro

Standard DXA measurements, including Fracture Risk Assessment Tool (FRAX) scores, have shown limitations in assessing fracture risk in Type 2 Diabetes (T2D), underscoring the need for novel biomarkers and suggesting that other pathomechanisms may drive diabetic bone fragility. MicroRNAs (miRNAs) are secreted into the circulation from cells of various tissues proportional to local disease severity and were recently found to be crucial to bone homeostasis and T2D. Here, we studied, if and which circulating miRNAs or combinations of miRNAs can discriminate best fracture status in a well-characterized study of diabetic bone disease and postmenopausal osteoporosis (n = 80 postmenopausal women). We then tested the most discriminative and most frequent miRNAs in vitro. Using miRNA-qPCR-arrays, we showed that 48 miRNAs can differentiate fracture status in T2D women and that several combinations of four miRNAs can discriminate diabetes-related fractures with high specificity and sensitivity (area under the receiver-operating characteristic curve values [AUCs], 0.92 to 0.96; 95% CI, 0.88 to 0.98). For the osteoporotic study arm, 23 miRNAs were fracture-indicative and potential combinations of four miRNAs showed AUCs from 0.97 to 1.00 (95% CI, 0.93 to 1.00). Because a role in bone homeostasis for those miRNAs that were most discriminative and most present among all miRNA combinations had not been described, we performed in vitro functional studies in human adipose tissue-derived mesenchymal stem cells to investigate the effect of miR-550a-5p, miR-188-3p, and miR-382-3p on osteogenesis, adipogenesis, and cell proliferation. We found that miR-382-3p significantly enhanced osteogenic differentiation (p < 0.001), whereas miR-550a-5p inhibited this process (p < 0.001). Both miRNAs, miR-382-3p and miR-550a-5p, impaired adipogenic differentiation, whereas miR-188-3p did not exert an effect on adipogenesis. None of the miRNAs affected significantly cell proliferation. Our data suggest for the first time that miRNAs are linked to fragility fractures in T2D postmenopausal women and should be further investigated for their diagnostic potential and their detailed function in diabetic bone. © 2016 American Society for Bone and Mineral Research.

Statistical Parametric Mapping of HR-pQCT Images: A Tool for Population-Based Local Comparisons of Micro-Scale Bone Features

HR-pQCT enables in vivo multi-parametric assessments of bone microstructure in the distal radius and distal tibia. Conventional HR-pQCT image analysis approaches summarize bone parameters into global scalars, discarding relevant spatial information. In this work, we demonstrate the feasibility and reliability of statistical parametric mapping (SPM) techniques for HR-pQCT studies, which enable population-based local comparisons of bone properties. We present voxel-based morphometry (VBM) to assess trabecular and cortical bone voxel-based features, and a surface-based framework to assess cortical bone features both in cross-sectional and longitudinal studies. In addition, we present tensor-based morphometry (TBM) to assess trabecular and cortical bone structural changes. The SPM techniques were evaluated based on scan-rescan HR-pQCT acquisitions with repositioning of the distal radius and distal tibia of 30 subjects. For VBM and surface-based SPM purposes, all scans were spatially normalized to common radial and tibial templates, while for TBM purposes, rescans (follow-up) were spatially normalized to their corresponding scans (baseline). VBM was evaluated based on maps of local bone volume fraction (BV/TV), homogenized volumetric bone mineral density (vBMD), and homogenized strain energy density (SED) derived from micro-finite element analysis; while the cortical bone framework was evaluated based on surface maps of cortical bone thickness, vBMD, and SED. Voxel-wise and vertex-wise comparisons of bone features were done between the groups of baseline and follow-up scans. TBM was evaluated based on mean square errors of determinants of Jacobians at baseline bone voxels. In both anatomical sites, voxel- and vertex-wise uni- and multi-parametric comparisons yielded non-significant differences, and TBM showed no artefactual bone loss or apposition. The presented SPM techniques demonstrated robust specificity thus warranting their application in future clinical HR-pQCT studies.

Bone microstructure in men assessed by HR-pQCT: Associations with risk factors and differences between men with normal, low, and osteoporosis-range areal BMD

Purpose

The primary objective of this study was to analyze the relationships between bone microstructure and strength, and male osteoporosis risk factors including age, body mass index, serum 25-hydroxyvitamin D level, and testosterone level. A secondary objective was to compare microstructural and strength parameters between men with normal, low, and osteoporosis-range areal bone mineral density (aBMD).

Methods

Seventy-eight healthy male volunteers (mean age 62.4 ± 7.8 years, range 50–84 years) were recruited. The participants underwent dual-energy X-ray absorptiometry (DXA) and high-resolution peripheral quantitative computed tomography (HR-pQCT) of the ultra-distal radius and tibia. From the HR-pQCT images, volumetric bone mineral density (BMD) and cortical and trabecular bone microstructure were evaluated, and bone strength and cortical load fraction (Ct.LF) were estimated using micro-finite element analysis (μFEA).

Results

Age was more strongly correlated with bone microstructure than other risk factors. Age had significant positive correlations with cortical porosity at both ultra-distal radius and tibia (r = 0.36, p = 0.001, and r = 0.47, p < 0.001, respectively). At the tibia, age was negatively correlated with cortical BMD, whereas it was positively correlated with trabecular BMD. In μFEA, age was negatively correlated with Ct.LF, although not with bone strength. Compared with men with normal aBMD, men with low or osteoporosis-range aBMD had significantly poor trabecular bone microstructure and lower bone strength at the both sites, while there was no significant difference in cortical bone.

Conclusions

Cortical bone microstructure was negatively affected by aging, and there was a suggestion that the influence of aging may be particularly important at the weight-bearing sites.

•

Bone microstructure of the ultra-distal radius and tibia was analyzed by HR-pQCT.

•

The subjects consisted of 78 healthy male volunteers ranging from 50 to 84 years.

•

Cortical bone was more impaired with age compared with trabecular bone.

•

Cortical porosity was increased with age, particularly at the ultra-distal tibia.

Kartogenin treatment prevented joint degeneration in a rodent model of osteoarthritis: A pilot study

Osteoarthritis (OA) is a major degenerative joint disease characterized by progressive loss of articular cartilage, synovitis, subchondral bone changes, and osteophyte formation. Currently there is no treatment for OA except temporary pain relief and end-stage joint replacement surgery. We performed a pilot study to determine the effect of kartogenin (KGN, a small molecule) on both cartilage and subchondral bone in a rat model of OA using multimodal imaging techniques. OA was induced in rats (OA and KGN treatment group) by anterior cruciate ligament transection (ACLT) surgery in the right knee joint. Sham surgery was performed on the right knee joint of control group rats. KGN group rats received weekly intra-articular injection of 125 μM KGN 1 week after surgery until week 12. All rats underwent in vivo magnetic resonance imaging (MRI) at 3, 6, and 12 weeks after surgery. Quantitative MR relaxation measures (T_1ρ and T₂ ) were determined to evaluate changes in articular cartilage. Cartilage and bone turnover markers (COMP and CTX-I) were determined at baseline, 3, 6, and 12 weeks. Animals were sacrificed at week 12 and the knee joints were removed for micro-computed tomography (micro-CT) and histology. KGN treatment significantly lowered the T_1ρ and T₂ relaxation times indicating decreased cartilage degradation. KGN treatment significantly decreased COMP and CTX-I levels indicating decreased cartilage and bone turnover rate. KGN treatment also prevented subchondral bone changes in the ACLT rat model of OA. Thus, kartogenin is a potential drug to prevent joint deterioration in post-traumatic OA. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1780-1789, 2016.

Determining Metacarpophalangeal Flexion Angle Tolerance for Reliable Volumetric Joint Space Measurements by High-resolution Peripheral Quantitative Computed Tomography

Objective.

The position-dependence of a method to measure the joint space of metacarpophalangeal (MCP) joints using high-resolution peripheral quantitative computed tomography (HR-pQCT) was studied.

Methods.

Cadaveric MCP were imaged at 7 flexion angles between 0 and 30 degrees. The variability in reproducibility for mean, minimum, and maximum joint space widths and volume measurements was calculated for increasing degrees of flexion.

Results.

Root mean square coefficient of variance values were < 5% under 20 degrees of flexion for mean, maximum, and volumetric joint spaces. Values for minimum joint space width were optimized under 10 degrees of flexion.

Conclusion.

MCP joint space measurements should be acquired at < 10 degrees of flexion in longitudinal studies.

Cortical bone laminar analysis reveals increased midcortical and periosteal porosity in type 2 diabetic postmenopausal women with history of fragility fractures compared to fracture-free diabetics

We investigated the characteristics and spatial distribution of cortical bone pores in postmenopausal women with type 2 diabetes (T2D). High porosity in the midcortical and periosteal layers in T2D subjects with fragility fractures suggests that these cortical zones might be particularly susceptible to T2D-induced toxicity and may reflect cortical microangiopathy.

INTRODUCTION

Elevated cortical porosity is regarded as one of the main contributors to the high skeletal fragility in T2D. However, to date, it remains unclear if diabetic cortical porosity results from vascular cortical changes or from an expansion in bone marrow space. Here, we used a novel cortical laminar analysis technique to investigate the characteristics and spatial radial distribution of cortical pores in a T2D group with prior history of fragility fractures (DMFx, assigned high-risk group) and a fracture-free T2D group (DM, assigned low-risk group) and to compare their results to non-diabetic controls with (Fx) and without fragility fractures (Co).

METHODS

Eighty postmenopausal women (n = 20/group) underwent high-resolution peripheral quantitative computed tomography (HR-pQCT) of the distal tibia and radius. Cortical bone was divided into three layers of equal width including an endosteal, midcortical, and periosteal layer. Within each layer, total pore area (TPA), total pore number (TPN), and average pore area (APA) were calculated. Statistical analysis employed Mann-Whitney tests and ANOVA with post hoc tests.

RESULTS

Compared to the DM group, DMFx subjects exhibited +90 to +365 % elevated global porosity (p = 0.001). Cortical laminar analysis revealed that this increased porosity was for both skeletal sites confined to the midcortical layer, followed by the periosteal layer (midcortical +1327 % TPA, p ≤ 0.001, periosteal +634 % TPA, p = 0.002), and was associated in both layers and skeletal sites with high TPN (+430 % TPN, p < 0.001) and high APA (+71.5 % APA, p < 0.001).

CONCLUSION

High porosity in the midcortical and periosteal layers in the high-risk T2D group suggests that these cortical zones might be particularly susceptible to T2D-induced toxicity and may reflect cortical microangiopathy.

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.

Novel anthropomorphic hip phantom corrects systemic interscanner differences in proximal femoral vBMD

Quantitative computed tomography (QCT) is increasingly used in osteoporosis studies to assess volumetric bone mineral density (vBMD), bone quality and strength. However, QCT is confronted by technical issues in the clinical research setting, such as potentially confounding effects of body size on vBMD measurements and lack of standard approaches to scanner cross-calibration, which affects measurements of vBMD in multicenter settings. In this study, we addressed systematic inter-scanner differences and subject-dependent body size errors using a novel anthropomorphic hip phantom, containing a calibration hip to estimate correction equations, and a contralateral test hip to assess the quality of the correction. We scanned this phantom on four different scanners and we applied phantom-derived corrections to in vivo images of 16 postmenopausal women scanned on two scanners. From the phantom study, we found that vBMD decreased with increasing phantom size in three of four scanners and that inter-scanner variations increased with increasing phantom size. In the in vivo study, we observed that inter-scanner corrections reduced systematic inter-scanner mean vBMD differences but that the inter-scanner precision error was still larger than expected from known intra-scanner precision measurements. In conclusion, inter-scanner corrections and body size influence should be considered when measuring vBMD from QCT images.

Automatic multi-parametric quantification of the proximal femur with quantitative computed tomography

BACKGROUND

Quantitative computed tomography (QCT) imaging is the basis for multiple assessments of bone quality in the proximal femur, including volumetric bone mineral density (vBMD), tissue volume, estimation of bone strength using finite element modeling (FEM), cortical bone thickness, and computational-anatomy-based morphometry assessments.

METHODS

Here, we present an automatic framework to perform a multi-parametric QCT quantification of the proximal femur. In this framework, the proximal femur is cropped from the bilateral hip scans, segmented using a multi-atlas based segmentation approach, and then assigned volumes of interest through the registration of a proximal femoral template. The proximal femur is then subjected to compartmental vBMD, compartmental tissue volume, FEM bone strength, compartmental surface-based cortical bone thickness, compartmental surface-based vBMD, local surface-based cortical bone thickness, and local surface-based cortical vBMD computations. Consequently, the template registrations together with vBMD and surface-based cortical bone parametric maps enable computational anatomy studies. The accuracy of the segmentation was validated against manual segmentations of 80 scans from two clinical facilities, while the multi-parametric reproducibility was evaluated using repeat scans with repositioning from 22 subjects obtained on CT imaging systems from two manufacturers.

RESULTS

Accuracy results yielded a mean dice similarity coefficient of 0.976±0.006, and a modified Haussdorf distance of 0.219±0.071 mm. Reproducibility of QCT-derived parameters yielded root mean square coefficients of variation (CVRMS) between 0.89-1.66% for compartmental vBMD; 0.20-1.82% for compartmental tissue volume; 3.51-3.59% for FEM bone strength; 1.89-2.69% for compartmental surface-based cortical bone thickness; and 1.08-2.19% for compartmental surface-based cortical vBMD. For local surface-based assessments, mean CVRMS were between 3.45-3.91% and 2.74-3.15% for cortical bone thickness and vBMD, respectively.

CONCLUSIONS

The automatic framework presented here enables accurate and reproducible QCT multi-parametric analyses of the proximal femur. Our subjects were elderly, with scans obtained across multiple clinical sites and manufacturers, thus documenting its value for clinical trials and other multi-site studies.

Spatial distribution of intracortical porosity varies across age and sex

Cortical bone porosity is a major determinant of strength, stiffness, and fracture toughness of cortical tissue. The goal of this work was to investigate changes in spatial distribution and microstructure of cortical porosity associated with aging in men and women. The specific aims were to: 1) develop an automated technique for spatial analysis of cortical microstructure based on HR-pQCT data, and; 2) apply this technique to explore sex- and age-specific spatial distribution and microstructure of porosity within the cortex. We evaluated HR-pQCT images of the distal tibia from a cross-sectional cohort of 145 individuals, characterizing detectable pores as being in the endosteal, midcortical, or periosteal layers of the cortex. Metrics describing porosity, pore number, and pore size were quantified within each layer and compared across sexes, age groups, and cortical layers. The elderly cohort (65-78 years, n=22) displayed higher values than the young cohort (20-29 years, n=29) for all parameters both globally and within each layer. While all three layers displayed significant age-related porosity increases, the greatest difference in porosity between the young and elderly cohort was in the midcortical layer (+344%, p<0.001). Similarly, the midcortical layer reflected the greatest differences between young and elderly cohorts in both pore number (+243%, p<0.001) and size (+28%, p<0.001). Females displayed greater age-related changes in porosity and pore number than males. Females and males displayed comparable small to non-significant changes with age in pore size. In summary, considerable variability exists in the spatial distribution of detectable cortical porosity at the distal tibia, and this variability is dependent on age and sex. Intracortical pore distribution analysis may ultimately provide insight into both mechanisms of pore network expansion and biomechanical consequences of pore distribution.

Quantitative characterization of metacarpal and radial bone in rheumatoid arthritis using high resolution- peripheral quantitative computed tomography

AIM

The objectives of this study were: (i) to develop a standardized method of quantifying bone mineral density (BMD) and microarchitecture in the hand and wrist bones of patients with rheumatoid arthritis (RA) using high resolution- peripheral quantitative computed tomography (HR-pQCT); (ii) to compare quantitative bone parameters between RA and post-menopausal osteopenic (PM-OP) subjects; and (iii) to correlate quantitative bone parameters at the distal radius with those at the metacarpal heads in RA subjects.

METHODS

HR-pQCT imaging of the dominant hand and wrist was performed in 12 female RA patients. BMD and trabecular parameters for the 2-12% head region of the second and third metacarpals were calculated and compared between RA patients and healthy controls. Bone parameters were also calculated for 110 slices of the distal radius in RA patients and compared to data from controls and PM-OP women from a previous study.

RESULTS

Compared to controls, RA patients had significantly decreased BMD, trabecular volume and number, and increased trabecular heterogeneity in the third metacarpal and distal radius. Significantly lower trabecular number and significantly higher ratio of outer annular trabecular BMD to inner trabecular BMD were observed in patients with RA, compared to patients with osteopenia (P < 0.05). Trabecular BMD in the third metacarpal and in the distal radius were significantly correlated (ρ = 0.918, P < 0.0001) in RA patients.

CONCLUSION

This study established a standardized method for quantifying bone density and trabecular properties in the hand and wrist bones of RA patients using HR-pQCT. Deterioration of bone structure in RA patients was found comparable to that in osteopenic women, and trabecular bone loss near affected joints was found to be correlated with bone loss away from joints.

Volumetric femoral BMD, bone geometry, and serum sclerostin levels differ between type 2 diabetic postmenopausal women with and without fragility fractures

While type 2 diabetes (T2D) is associated with higher skeletal fragility, specific risk stratification remains incompletely understood. We found volumetric bone mineral density, geometry, and serum sclerostin differences between low-fracture risk and high-fracture risk T2D women. These features might help identify T2D individuals at high fracture risk in the future.

INTRODUCTION

Diabetic bone disease, an increasingly recognized complication of type 2 diabetes mellitus (T2D), is associated with high skeletal fragility. Exactly which T2D individuals are at higher risk for fracture, however, remains incompletely understood. Here, we analyzed volumetric bone mineral density (vBMD), geometry, and serum sclerostin levels in two specific T2D subsets with different fracture risk profiles. We examined a T2D group with prior history of fragility fractures (DMFx, assigned high-risk group) and a fracture-free T2D group (DM, assigned low-risk group) and compared their results to nondiabetic controls with (Fx) and without fragility fractures (Co).

METHODS

Eighty postmenopausal women (n = 20 per group) underwent quantitative computed tomography (QCT) to compute vBMD and bone geometry of the proximal femur. Additionally, serum sclerostin, vitamin D, parathyroid hormone (PTH), HbA1c, and glomerular filtration rate (GFR) levels were measured. Statistical analyses employed linear regression models.

RESULTS

DMFx subjects exhibited up to 33 % lower femoral neck vBMD than DM subjects across all femoral sites (-19 % ≤ ΔvBMD ≤ -33 %, 0.008 ≤ p ≤0.021). Additionally, DMFx subjects showed significantly thinner cortices (-6 %, p = 0.046) and a trend toward larger bone volume (+10 %, p = 0.055) relative to DM women and higher serum sclerostin levels when compared to DM (+31.4 %, p = 0.013), Fx (+25.2 %, p = 0.033), and control (+22.4 %, p = 0.028) subjects.

CONCLUSION

Our data suggest that volumetric bone parameters by QCT and serum sclerostin levels can identify T2D individuals at high risk of fracture and might therefore show promise as clinical tools for fracture risk assessment in T2D. However, future research is needed to establish diabetes-specific QCT- and sclerostin-reference databases.

Validation of bone marrow fat quantification in the presence of trabecular bone using MRI

BACKGROUND

To validate six-echo, chemical-shift based MRI with T2 * correction for the quantification of bone marrow fat content in the presence of trabecular bone.

METHODS

Ten bone phantoms were made using trabecular bone cores extracted from the distal femur and proximal tibia of 20 human cadaveric knees. Bone marrow was removed from the cores and the marrow spaces were filled with water-fat gelatin to mimic bone marrow of known fat fractions. A chemical-shift based water-fat separation method with T2 * correction was used to generate fat fraction maps. The proton density fat fractions (PDFF) between marrow regions with and without bone were compared with the reference standard of known fat fraction using the squared Pearson correlation coefficient and unpaired t-test.

RESULTS

Strong correlations were found between the known fat fraction and measured PDFF in marrow without trabecular bone (R(2) = 0.99; slope = 0.99, intercept = 0.94) as well as in marrow with trabecular bone (R(2) = 0.97; slope = 1.0, intercept = -3.58). Measured PDFF between regions with and without bone were not significantly different (P = 0.5). However, PDFF was systematically underestimated by -3.2% fat fraction in regions containing trabecular bone.

CONCLUSION

Our implementation of a six-echo chemical-shift based MRI pulse sequence with T2 * correction provided an accurate means of determining fat content in bone marrow in the presence of trabecular bone.

Three-dimensional analysis of subchondral cysts in hip osteoarthritis: an ex vivo HR-pQCT study

INTRODUCTION

Subchondral cysts are deeply related to the pathogenesis of osteoarthritis (OA), but the factors contributing to cyst formation are not well known. A three-dimensional analysis of subchondral cysts at the micro-structural level was conducted using a high-resolution peripheral quantitative CT (HR-pQCT), and their relationships with cartilage attrition and subchondral bone microstructure were investigated.

METHODS

Femoral heads extracted from ten female patients with hip OA were scanned using an HR-pQCT at a voxel size of 41μm. The volume fractions, numbers, and sizes of the cysts were measured in the subchondral bone region under the area of cartilage loss. Furthermore, the areas of cartilage loss, as well as the microstructure of the subchondral bones, were also measured, and their correlations with the cysts were analyzed.

RESULTS

The volume fractions of cysts within subchondral bone regions varied from 2% to 33%, the numbers of cysts varied from 6 to 87, and the sizes varied from 1mm(3) to 657mm(3). There was a positive correlation between the number of cysts and bone volume (r>0.8, p<0.01).

CONCLUSION

The degree of cyst formation showed a wide distribution in number and volume, and there was a close relationship between multiple cyst formation and bone sclerosis, which might be caused by reactive bone formation that occurred around each cyst.

Improved trabecular bone structure of 20-month-old male spontaneously hypertensive rats

A few clinical studies have reported that elderly male participants with hypertensive disease frequently have higher bone mineral density (BMD) than the normotensive participants at several skeletal sites. The detailed mechanism is still unknown; therefore, a study of bone structure and density using the hypertensive animal models could be informative. We used micro-computed tomography to quantitatively evaluate the tibial and 3rd lumbar vertebral bones in the 20-month-old male spontaneous hypertensive rat (SHR). The BMD, volume fraction, and the microarchitecture changes of the SHR were compared to those of same-age normotensive controls (Wistar-Kyoto rat, WKY). We found that in the very old (20 month) male rats, the trabecular bone fraction and microstructure were higher than those in the same-age normotensive controls. The observation of the association of hypertension with BMD and bone strength in hypertensive rats warrants further investigations of bone mass and strength in elderly males with hypertension.

Structural analysis of cortical porosity applied to HR-pQCT data

PURPOSE

The investigation of cortical porosity is an important aspect of understanding biological, pathoetiological, and biomechanical processes occurring within the skeleton. With the emergence of HR-pQCT as a noninvasive tool suitable for clinical use, cortical porosity at appendicular sites can be directly visualized in vivo. The aim of this study was to introduce a novel topological analysis of the cortical pore network for HR-pQCT data and determine the influence of resolution on measures of cortical pore network microstructure and topology.

METHODS

Cadaveric radii were scanned using HR-pQCT at two different voxel sizes (41 and 82 μm) and also using μCT at a voxel size of 18 μm. HR-pQCT and μCT image sets were spatially coregistered. Segmentation and quantification of cortical porosity (Ct.Po) and mean pore diameter (Ct.Po.Dm) were achieved using an established extended cortical analysis technique. Topological classification of individual pores was performed using topology-preserving skeletonization and multicolor dilation algorithms. Based on the pore skeleton topological classification, the following parameters were quantified: total number of planar surface-skeleton canals (N.Slabs), tubular curve-skeleton canals (N.Tubes), and junction elements (N.Junctions), mean slab volume (Slab.Vol), mean tube volume (Tube.Vol), mean slab orientation (Slab.θ), mean tube orientation (Tube.θ), N.Slabs/N.Tubes, and integral (total) slab volume/integral tube volume (iSlab.Vol/iTube.Vol). An in vivo reproducibility study was also conducted to assess short-term precision of the topology parameters. Precision error was characterized using root mean square coefficient of variation (RMSCV%).

RESULTS

Correlations to μCT values for Ct.Po were significant for both the 41 and 82 μm HR-pQCT data (41: r(2) = 0.82, p < 0.001, 82: r(2) = 0.75, p < 0.001). For Ct.Po.Dm, only the 41 μm data were significantly predictive of μCT values (r(2) = 0.72, p < 0.01) Data at both HR-pQCT voxel sizes were strongly predictive of the μCT values for N.Slabs (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.84, p < 0.001), N.Tubes (41: r(2) = 0.94, p < 0.001; 82: r(2) = 0.84, p < 0.001), and N.Junctions (41: r(2) = 0.93, p < 0.001; 82: r(2) = 0.78, p < 0.001), though proportional bias was evident in these correlations. Weak correlations were seen for iSlab.Vol/iTube.Vol at both voxel sizes (41: r(2) = 0.52, p < 0.01; 82: r(2) = 0.39, p < 0.05). Slab.Vol was significantly correlated to μCT data at 41 μm (r(2) = 0.60, p < 0.01) but not at 82 μm, while Tube.Vol was significantly correlated at both voxel sizes (41: r(2) = 0.79, p < 0.001; 82: r(2) = 0.68, p < 0.01). In vivo precision error for these parameters ranged from 2.31 to 9.68 RMSCV%.

CONCLUSIONS

Strong correlations between μCT- and HR-pQCT-derived measurements were found, particularly in HR-pQCT images obtained at 41 μm. These data are in agreement with our previous study investigating the effect of voxel size on standard HR-pQCT metrics of trabecular and cortical microstructure, and extend our previous findings to include topological descriptors of the cortical pore network.

Bone structure and perfusion quantification of bone marrow edema pattern in the wrist of patients with rheumatoid arthritis: a multimodality study

OBJECTIVE

To quantify bone structure and perfusion parameters in regions of bone marrow edema pattern (BMEP), non-edematous bone marrow (NBM), and pannus tissue areas in the wrists of patients with rheumatoid arthritis (RA) using 3-Tesla (3T) magnetic resonance imaging (MRI), and high resolution peripheral quantitative computed tomography (HR-pQCT).

METHODS

Sixteen subjects fulfilling American College of Rheumatology classification were imaged using a HR-pQCT system and a 3T MRI scanner with an 8-channel wrist coil. Coronal T2-weighted and dynamic contrast-enhanced (DCE-MRI) images were acquired. BMEP and pannus tissue areas were segmented semiautomatically in T2-weighted images. NBM areas were placed at a similar distance from the joint space as BMEP regions. MR and HR-pQCT images were registered, and bone variables were calculated within the BMEP and NBM regions. Perfusion parameters in BMEP, pannus tissue, and NBM regions were calculated based on the signal-time curve obtained from DCE-MRI.

RESULTS

Eighteen BMEP areas were segmented, 15 of them presented proximal to pannus-filled erosions. Significant increases in bone density and trabecular thickness and number were observed in all BMEP regions compared to NMB (p < 0.05). Significantly elevated perfusion measures were observed in both BMEP and pannus tissue regions compared to NBM (p < 0.05).

CONCLUSION

BMEP regions showed significantly increased bone density and structures as well as perfusion measures, suggesting bone remodeling and active inflammation. Combining MRI and HR-pQCT provides a powerful multimodality approach for understanding BMEP and erosions, and for potentially identifying novel imaging markers for disease progression in RA.

The influence of disuse on bone microstructure and mechanics assessed by HR-pQCT

Numerous clinical cohorts are exposed to reduced skeletal loading and associated bone loss, including surgical patients, stroke and spinal cord injury victims, and women on bed rest during pregnancy. In this context, understanding disuse-related bone loss is critical to developing interventions to prevent fractures and the associated morbidity, mortality, and cost to the health care system. The aim of this pilot study was to use high-resolution peripheral QCT (HR-pQCT) to examine changes in trabecular and cortical microstructure and biomechanics during a period of non weight bearing (WB) and during recovery following return to normal WB. Surgical patients requiring a 6-week non WB period (n=12, 34.8±7.7 yrs) were scanned at the affected and contralateral tibia prior to surgery, after the 6-week non WB period, and 6 and 13 weeks after returning to full WB. At the affected ultradistal tibia, integral vBMD (including both trabecular and cortical compartments) decreased with respect to baseline (-1.2%), trabecular number increased (+5.6%), while trabecular thickness (-5.4%), separation (-4.6%), and heterogeneity (-7.2%) decreased (all p<0.05). Six weeks after return to full WB, trabecular structure measures reverted to baseline levels. In contrast, integral vBMD continued to decrease after 6 (-2.0%, p<0.05) and 13 weeks (-2.5%, p=0.07) of full WB. At the affected distal site, the disuse period resulted in increased porosity (+16.1%, p<0.005), which remained elevated after 6 weeks (+16.8%, p<0.01) and after 13 weeks (+16.2%, p<0.05). A novel topological analysis applied to the distal tibia cortex demonstrated increased number of canals with surface topology ("slabs" +21.7%, p<0.01) and curve topology ("tubes" +15.0%, p<0.05) as well as increased number of canal junctions (+21.4%, p<0.05) following the disuse period. Porosity increased uniformly through increases in both pore size and number. Finite element analysis at the ultradistal tibia showed decreased stiffness and failure load (-2.8% and -2.4%, p<0.01) following non WB. These biomechanical predictions remained depressed following 6 and 13 weeks of full WB. Finite element analysis at the distal site followed similar trends. Our results suggest that detectable microstructural and biomechanical degradation occurs--particularly within the cortical compartment--as a result of non WB and persists following return to normal loading. A better understanding of these microstructural changes and their short- and long-term influence on biomechanics may have clinical relevance in the context of disuse-related fracture prevention.