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.