Quantification of lower leg arterial calcifications by high-resolution peripheral quantitative computed tomography

Vascular calcifications and bone health seem to be etiologically linked via common risk factors such as aging and subclinical chronic inflammation. Epidemiologic studies have shown significant associations between low bone mineral density (BMD), fragility fractures and calcifications of the coronary arteries and the abdominal aorta. In the last decade, high-resolution peripheral quantitative computed tomography (HR-pQCT) has emerged as in-vivo research tool for the assessment of peripheral bone geometry, density, and microarchitecture. Although vascular calcifications are frequently observed as incidental findings in HR-pQCT scans, they have not yet been incorporated into quantitative HR-pQCT analyses. We developed a semi-automated algorithm to quantify lower leg arterial calcifications (LLACs), captured by HR-pQCT. The objective of our study was to determine validity and reliability of the LLAC measure. HR-pQCT scans were downscaled to a voxel size of 250μm. After subtraction of bone volumes from the scans, LLACs were detected and contoured by a semi-automated, dual-threshold seed-point segmentation. LLAC mass (in mg hydroxyapatite; HA) was calculated as the product of voxel-based calcification volume (mm(3)) and mean calcification density (mgHA/cm(3))/1000. To determine validity, we compared LLACs to coronary artery calcifications (CACs), as quantified by multi-detector computed tomography (MDCT) and Agatston scoring in forty-six patients on chronic hemodialysis. Moreover, we investigated associations of LLACs with age, time on dialysis, type-2 diabetes mellitus, history of stroke, and myocardial infarction. In a second step, we determined intra- and inter-reader reliability of the LLAC measure. In the validity study, LLACs were present (>0mgHA) in 76% of patients, 78% of patients had CACs (>0mgHA). Median LLAC was 6.65 (0.08-24.40)mgHA and median CAC as expressed by Agatston score was 266.3 (15.88-1877.28). We found a significant positive correlation between LLAC and CAC (rho=0.6; p<0.01). Dialysis patients with type-2 diabetes mellitus (DM; 35%) and history of stroke (13%) had higher median LLAC than patients without those conditions (DM 20.0 fold greater, p=0.006; Stroke 5.1 fold greater, p=0.047). LLAC was positively correlated with time on dialysis (rho=0.337, p=0.029), there was a trend towards a positive association of LLAC and age (rho=0.289, p=0.053). The reliability study yielded excellent intra- and inter-reader agreement of the LLAC measure (intra-reader ICC=0.999, 95% CI=0.998-1.000; inter-reader ICC=0.998, 95% CI=0.994-0.999). Our study indicates that the LLAC measure has good validity and excellent reliability. The use of HR-pQCT for the simultaneous evaluation of arterial calcifications, peripheral bone geometry, bone density, and bone microarchitecture should facilitate future research on osteo-vascular interactions and potential associations with cardiovascular events.

Postmenopausal women treated with combination parathyroid hormone (1-84) and ibandronate demonstrate different microstructural changes at the radius vs. tibia: the PTH and Ibandronate Combination Study (PICS)

SUMMARY

In postmenopausal women receiving combination parathyroid hormone (PTH) (1-84) therapy and ibandronate, we evaluated bone microarchitecture and biomechanics using high-resolution peripheral quantitative computed tomography (HR-pQCT). Cortical and trabecular changes were different at the nonweight-bearing radius vs. the weight-bearing tibia, with more favorable overall changes at the tibia.

INTRODUCTION

PTH therapy and bisphosphonates decrease fracture risk in postmenopausal osteoporosis, but their effects on bone microstructure and strength have not been fully characterized, particularly during combination therapy. PTH increases trabecular bone mineral density (BMD) substantially but may decrease cortical BMD, possibly by stimulating intracortical remodeling. We evaluated bone microarchitecture and biomechanics with HR-pQCT at the radius (a nonweight-bearing site) and tibia (weight bearing) in women receiving combination PTH(1-84) and ibandronate.

METHODS

Postmenopausal women with low bone mass (n = 43) were treated with 6 months of PTH(1-84) (100 μg/day), either as one 6- or two 3-month courses, in combination with ibandronate (150 mg/month) over 2 years. HR-pQCT was performed before and after therapy.

RESULTS

Because changes in HR-pQCT parameters did not differ between treatment arms, groups were pooled into one cohort for analysis. Trabecular BMD increased at both radius and tibia (p < 0.01 for each). Cortical thickness and BMD decreased at the radius (p < 0.01), consistent with changes in dual-energy X-ray absorptiometry, while these parameters did not change at the tibia (p ≤ 0.02 for difference between radius and tibia). In contrast, cortical porosity increased at the tibia (p < 0.01) but not radius. Stiffness and failure load decreased at the radius (p < 0.0001) but did not change at the tibia.

CONCLUSIONS

Cortical and trabecular changes in response to the PTH/ibandronate treatment combinations utilized in this study were different at the nonweight-bearing radius vs. the weight-bearing tibia, with more favorable overall changes at the tibia. Our findings support the possibility that weight bearing may optimize the effects of osteoporosis therapy.

Heterogeneity of bone microstructure in the femoral head in patients with osteoporosis: an ex vivo HR-pQCT study

INTRODUCTION

Trabecular bone in the femoral head has a complicated and heterogeneous structure with few studies having analyzed heterogeneity in this structure quantitatively. We analyze trabecular bone microstructure in the femoral head with osteoporosis (OP) using high resolution peripheral quantitative CT (HR-pQCT) to investigate its regional characteristics.

METHODS

Fifteen femoral heads extracted from female OP patients with femoral neck fracture (85 ± 7, 67-94 years) were scanned by HR-pQCT at 41 μm voxel size. The femoral head was segmented into 15 regions (3 longitudinal regions: superior, center, and inferior, and 5 axial subregions: center, medial, lateral, anterior, posterior). Of these 15 regions, five were excluded due to overlap with the fracture site, leaving a total of 10 regions of cancellous bone microstructures to be quantitatively assessed using the following parameters: bone volume fraction, trabecular thickness, number, separation, connectivity density, structure model index, and degree and orientation of anisotropy. These parameters were compared among each region.

RESULTS

Trabecular bone at the center, superior, and supero-posterior regions of the femoral head had higher bone volume, trabecular number, thickness, narrower bone marrow spaces, higher connectivity and anisotropy, and more plate-like structure. This plate-like structure ran supero-inferiorly and antero-posteriorly at the superior and center regions. Bone volume at the anterior, posterior, and medial regions was almost half of the central and superior regions.

CONCLUSION

Significant heterogeneity of the trabecular bone microstructure in the OP femoral head was showed quantitatively in this study. These data offer new insight into bone microstructural anatomy and may prove to provide useful information on clinical medicine such as hip surgeries.

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

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

Computational identification and quantification of trabecular microarchitecture classes by 3-D texture analysis-based clustering

High resolution peripheral quantitative computed tomography (HR-pQCT) permits the non-invasive assessment of cortical and trabecular bone density, geometry, and microarchitecture. Although researchers have developed various post-processing algorithms to quantify HR-pQCT image properties, few of these techniques capture image features beyond global structure-based metrics. While 3D-texture analysis is a key approach in computer vision, it has been utilized only infrequently in HR-pQCT research. Motivated by high isotropic spatial resolution and the information density provided by HR-pQCT scans, we have developed and evaluated a post-processing algorithm that quantifies microarchitecture characteristics via texture features in HR-pQCT scans. During a training phase in which clustering was applied to texture features extracted from each voxel of trabecular bone, three distinct clusters, or trabecular microarchitecture classes (TMACs) were identified. These TMACs represent trabecular bone regions with common texture characteristics. The TMACs were then used to automatically segment the voxels of new data into three regions corresponding to the trained cluster features. Regional trabecular bone texture was described by the histogram of relative trabecular bone volume covered by each cluster. We evaluated the intra-scanner and inter-scanner reproducibility by assessing the precision errors (PE), intra class correlation coefficients (ICC) and Dice coefficients (DC) of the method on 14 ultradistal radius samples scanned on two HR-pQCT systems. DC showed good reproducibility in intra-scanner set-up with a mean of 0.870±0.027 (no unit). Even in the inter-scanner set-up the ICC showed high reproducibility, ranging from 0.814 to 0.964. In a preliminary clinical test application, the TMAC histograms appear to be a good indicator, when differentiating between postmenopausal women with (n=18) and without (n=18) prevalent fragility fractures. In conclusion, we could demonstrate that 3D-texture analysis and feature clustering seems to be a promising new HR-pQCT post-processing tool with good reproducibility, even between two different scanners.

Multicenter precision of cortical and trabecular bone quality measures assessed by high-resolution peripheral quantitative computed tomography

High-resolution peripheral quantitative computed tomography (HR-pQCT) has recently been introduced as a clinical research tool for in vivo assessment of bone quality. The utility of this technology to address important skeletal health questions requires translation to standardized multicenter data pools. Our goal was to evaluate the feasibility of pooling data in multicenter HR-pQCT imaging trials. Reproducibility imaging experiments were performed using structure and composition-realistic phantoms constructed from cadaveric radii. Single-center precision was determined by repeat scanning over short-term (<72 hours), intermediate-term (3-5 months), and long-term intervals (28 months). Multicenter precision was determined by imaging the phantoms at nine different HR-pQCT centers. Least significant change (LSC) and root mean squared coefficient of variation (RMSCV) for each interval and across centers was calculated for bone density, geometry, microstructure, and biomechanical parameters. Single-center short-term RMSCVs were <1% for all parameters except cortical thickness (Ct.Th) (1.1%), spatial variability in cortical thickness (Ct.Th.SD) (2.6%), standard deviation of trabecular separation (Tb.Sp.SD) (1.8%), and porosity measures (6% to 8%). Intermediate-term RMSCVs were generally not statistically different from short-term values. Long-term variability was significantly greater for all density measures (0.7% to 2.0%; p < 0.05 versus short-term) and several structure measures: cortical thickness (Ct.Th) (3.4%; p < 0.01 versus short-term), cortical porosity (Ct.Po) (15.4%; p < 0.01 versus short-term), and trabecular thickness (Tb.Th) (2.2%; p < 0.01 versus short-term). Multicenter RMSCVs were also significantly higher than short-term values: 2% to 4% for density and micro-finite element analysis (µFE) measures (p < 0.0001), 2.6% to 5.3% for morphometric measures (p < 0.001), whereas Ct.Po was 16.2% (p < 0.001). In the absence of subject motion, multicenter precision errors for HR-pQCT parameters were generally less than 5%. Phantom-based multicenter precision was comparable to previously reported in in vivo single-center precision errors, although this was approximately two to five times worse than ex vivo short-term precision. The data generated from this study will contribute to the future design and validation of standardized procedures that are broadly translatable to multicenter study designs.

Quantitative and semiquantitative bone erosion assessment on high-resolution peripheral quantitative computed tomography in rheumatoid arthritis

OBJECTIVE

To develop novel quantitative and semiquantitative bone erosion measures at metacarpophalangeal (MCP) and wrist joints in patients with rheumatoid arthritis (RA) using high-resolution peripheral quantitative computed tomography (HR-pQCT), and to correlate these measurements with disease duration and bone marrow edema (BME) patterns derived from magnetic resonance imaging (MRI).

METHODS

Sixteen patients with RA and 7 healthy subjects underwent hand and wrist HR-pQCT and 3-Tesla MRI. Bone erosions of the MCP2, MCP3, and distal radius were evaluated by measuring maximal erosion dimension on axial slices, which is a simple and fast measurement, and then were graded (grades 0-3) based on the maximal dimension. Correlation coefficients were calculated between (1) sum maximal dimensions, highest grades, and sum grades of bone erosions; (2) erosion measures and the clinical evaluation; (3) erosion measures and BME volume in distal radius.

RESULTS

The inter- and intrareader agreements of maximal erosion dimensions were excellent (intraclass correlation coefficients 0.89, 0.99, and root mean square error 9.4%, 4.7%, respectively). Highest grades and sum grades were significantly correlated to sum maximal dimensions of all erosions. Number of erosions, sum maximal erosion dimensions, highest grades, and sum grades correlated significantly with disease duration. Number of erosions, sum maximal dimensions, and erosion grading of the distal radius correlated significantly with BME volume.

CONCLUSION

HR-pQCT provides a sensitive method with high reader agreement in assessment of structural bone damage in RA. The good correlation of erosion measures with disease duration as well as BME volume suggests that they could become feasible measures of erosions in RA.

Increased cortical porosity in type 2 diabetic postmenopausal women with fragility fractures

The primary goal of this study was to assess peripheral bone microarchitecture and strength in postmenopausal women with type 2 diabetes with fragility fractures (DMFx) and to compare them with postmenopausal women with type 2 diabetics without fractures (DM). Secondary goals were to assess differences in nondiabetic postmenopausal women with fragility fractures (Fx) and nondiabetic postmenopausal women without fragility fractures (Co), and in DM and Co women. Eighty women (mean age 61.3 ± 5.7 years) were recruited into these four groups (DMFx, DM, Fx, and Co; n = 20 per group). Participants underwent dual-energy X-ray absorptiometry (DXA) and high-resolution peripheral quantitative computed tomography (HR-pQCT) of the ultradistal and distal radius and tibia. In the HR-pQCT images volumetric bone mineral density and cortical and trabecular structure measures, including cortical porosity, were calculated. Bone strength was estimated using micro-finite element analysis (µFEA). Differential strength estimates were obtained with and without open cortical pores. At the ultradistal and distal tibia, DMFx had greater intracortical pore volume (+52.6%, p = 0.009; +95.4%, p = 0.020), relative porosity (+58.1%, p = 0.005; +87.9%, p = 0.011) and endocortical bone surface (+10.9%, p = 0.031; +11.5%, p = 0.019) than DM. At the distal radius DMFx had 4.7-fold greater relative porosity (p < 0.0001) than DM. At the ultradistal radius, intracortical pore volume was significantly higher in DMFx than DM (+67.8%, p = 0.018). DMFx also displayed larger trabecular heterogeneity (ultradistal radius: +36.8%, p = 0.035), and lower total and cortical BMD (ultradistal tibia: -12.6%, p = 0.031; -6.8%, p = 0.011) than DM. DMFx exhibited significantly higher pore-related deficits in stiffness, failure load, and cortical load fraction at the ultradistal and distal tibia, and the distal radius than DM. Comparing nondiabetic Fx and Co, we only found a nonsignificant trend with increase in pore volume (+38.9%, p = 0.060) at the ultradistal radius. The results of our study suggest that severe deficits in cortical bone quality are responsible for fragility fractures in postmenopausal diabetic women.

Age- and gender-related differences in cortical geometry and microstructure: Improved sensitivity by regional analysis

OBJECTIVE

While the importance of cortical structure quantification is increasingly underscored by recent literature, conventional analysis techniques obscure potentially important regional variations in cortical structure. The objective of this study was to characterize the spatial variability in cortical geometry and microstructure at the distal radius and tibia using high resolution peripheral quantitative computed tomography (HR-pQCT). We show that spatially-resolved analysis is able to identify cortical sub-regions with increased sensitivity to the effects of gender and aging.

METHODS

HR-pQCT scans of 146 volunteers (92 female/54 male) spanning a wide range of ages (20-78years) were analyzed. For each subject, radius and tibia scans were obtained using a clinical HR-pQCT system. Measures describing geometry (cortical bone thickness (Ct.Th)), microstructure (porosity (Ct.Po), pore diameter (Ct.Po.Dm), and pore size heterogeneity (Ct.Po.Dm SD)), and cortical bone density were calculated from the image data. Biomechanical parameters describing load and stress distribution were calculated using linear finite element analysis. Cortical quadrants were defined based on anatomic axes to quantify regional parameter variation. Subjects were categorized by gender, and age, and menopausal status for analysis.

RESULTS

Significant regional variation was found in all geometric and microstructural parameters in both the radius and tibia. In general, the radius showed more pronounced and significant variations in all parameters as compared with the tibia. At both sites, Ct.Po displayed the greatest regional variations. Correlation coefficients for Ct.Po and Ct.Th with respect to load and stress distribution provided evidence of an association between regional cortical structure and biomechanics in the tibia. Comparing women to men, differences in Ct.Po were most pronounced in the anterior quadrant of the radius (36% lower in women (p<0.01)) and the posterior quadrant of the tibia (27% lower in women (p<0.01)). Comparing elderly to young women, differences in Ct.Po were most pronounced in the lateral quadrant of the radius (328% higher in elderly women (p<0.001)) and the anterior quadrant of the tibia (433% higher in elderly women (p<0.001)). Comparing elderly to young men, the most pronounced age differences were found in the anterior radius (205% higher in elderly men, (p<0.001)) and the anterior tibia (190% higher in elderly men (p<0.01)). All subregional Ct.Po differences provided greater sensitivity to gender and age effects than those based on the global means.

CONCLUSION

These results show significant regional variation in all geometric and microarchitectural parameters studied in both the radius and tibia. Quantification of region-specific parameters provided increased sensitivity in the analysis of age- and gender-related differences, in many cases providing statistically significant differentiation of groups where conventional global analysis failed to detect differences. These results suggest that regional analysis may be important in studies of disease and therapeutic effects, particularly where microstructural parameters based on global analyses have thus far failed to identify a response in bone quality.

The effect of voxel size on high-resolution peripheral computed tomography measurements of trabecular and cortical bone microstructure

PURPOSE

Accurate quantification of bone microstructure plays a significant role in understanding bone mechanics and response to disease or treatment. High-resolution peripheral quantitative computed tomography (HR-pQCT) allows for the quantification of trabecular and cortical structure in vivo, with the capability of generating images at multiple voxel sizes (41, 82, and 123 μm). The aim of this study was to characterize the effect of voxel size on structural measures of trabecular and cortical bone and to determine accuracy in reference to micro-CT ([micro sign]CT), the gold standard for bone microstructure quantification.

METHODS

Seventeen radii from human cadaver specimens were imaged at each HR-pQCT voxel size and subsequently imaged using [micro sign]CT. Bone density and microstructural assessment was performed in both the trabecular and cortical compartments, including cortical porosity quantification. Two distinct analysis techniques were applied to the 41 μm HR-pQCT data: the standard clinical indirect analysis and a direct analysis requiring no density or structural model assumptions. Analysis parameters were adjusted to enable segmentation and structure extraction at each voxel size.

RESULTS

For trabecular microstructural measures, the 41 μm HR-pQCT data displayed the strongest correlations and smallest errors compared to [micro sign]CT data. The direct analysis technique applied to the 41 μm data yielded an additional improvement in accuracy, especially for measures of trabecular thickness. The 123 μm data performed poorly, with all microstructural measures either having moderate or nonsignificant correlations with [micro sign]CT data. Trabecular densitometric measures showed strong correlations to [micro sign]CT data across all voxel sizes. Cortical thickness was strongly correlated with [micro sign]CT values across all HR-pQCT voxel sizes. The accuracy of cortical porosity parameters was highly dependent on voxel size; again, the 41 μm data was most strongly correlated. Measures of cortical density and pore diameter at all HR-pQCT voxel sizes had either weak or nonsignificant correlations.

CONCLUSIONS

This study demonstrates the effect of voxel size on the accuracy of HR-pQCT measurements of trabecular and cortical microstructure and presents parameters for HR-pQCT analysis at nonstandard resolutions. For all parameters measured, correlations were strongest at 41 μm. Weak correlations for porosity measures indicate that a better understanding of pore structure and resolution dependence is needed.

Vascular patterning and permeability in prostate cancer models with differing osteogenic properties

Bone metastasis is a major cause of morbidity and mortality in prostate cancer. However, the lack of clinically relevant models hinders our understanding of the disease as well as development of effective therapies and imaging approaches. We used noninvasive MRI, histology and micro CT to further characterize the newly established prostate cancer bone metastases-derived model MDA-PCa-118b, and to compare it to the well-established PC-3MM2 model with regard to bone structure and vascular patterning. The PC-3MM2 model is highly osteolytic whereas the MDA-PCa-118b model shows a robust osteoblastic reaction, as often seen in clinical cases. Macromolecular contrast enhanced MRI revealed differences in vascular permeability patterns, which appeared peripheral for PC-3MM2 and nodular for MDA-PCa-118b, matching the microscopic cellular composition of each model: PC-3MM2 exclusively recruits endothelial cells to form thin tumor-core blood vessels and enlarged, leaky peripheral vessels, whereas MDA-PCa-118b also recruits bone-forming cells and pericytes such that small tumor nests are encircled with leaky vessels and embedded in bone-like tissue dotted with pericyte-covered vessels. Despite these structural differences, vascular permeability was reduced in both tumor models by either imatinib or SU10944 treatment. This study highlights the importance of clinically relevant osteogenic models of human prostate cancer and the value of such models not only in enhancing our understanding of tumorigenesis, metastasis and response to therapy, but also for development of appropriate methods for noninvasive imaging of these processes.

Noninvasive imaging of bone microarchitecture

The noninvasive quantification of peripheral compartment-specific bone microarchitecture is feasible with high-resolution peripheral quantitative computed tomography (HR-pQCT) and high-resolution magnetic resonance imaging (HR-MRI). In addition to classic morphometric indices, both techniques provide a suitable basis for virtual biomechanical testing using finite element (FE) analyses. Methodical limitations, morphometric parameter definition, and motion artifacts have to be considered to achieve optimal data interpretation from imaging studies. With increasing availability of in vivo high-resolution bone imaging techniques, special emphasis should be put on quality control including multicenter, cross-site validations. Importantly, conclusions from interventional studies investigating the effects of antiosteoporotic drugs on bone microarchitecture should be drawn with care, ideally involving imaging scientists, translational researchers, and clinicians.

Morphology of the human vertebral endplate

It is presumed that poor intervertebral disc cell nutrition is a contributing factor in degeneration, and is exacerbated by vertebral endplate sclerosis. Yet, quantitative relationships between endplate morphology and degeneration are unavailable. We investigated how endplate bone microstructure relates to indices of disc degeneration, such as morphologic grade, proteoglycan content, and cell density. Intervertebral core samples [n = 96, 14 subjects, L1-L5 level, ages 35-85 (64 ± 16 years), degeneration grade 1 (n = 4), grade 2 (n = 32), grade 3 (n = 44), grade 4 (n = 10), grade 5 (n = 6)] that included subchondral bone, cartilage endplate, and adjacent nucleus were harvested from human cadaveric lumbar spines. The morphology of the vertebral endplate was analyzed using µCT and the adjacent nucleus tissue was collected for biochemical and cellular analyses. Relationships between vertebral endplate morphology and adjacent disc degeneration were analyzed. Contrary to the prevailing notion, vertebral endplate porosity increased between 50% and 130% and trabecular thickness decreased by between 20% and 50% with advancing disc degeneration (p < 0.05). We also observed that nucleus cell density increased (R(2)  = 0.33, p < 0.05) and proteoglycan content decreased (R(2)  = 0.47, p < 0.05) as the endplate became more porous. Our data suggest that endplate sclerosis is not a fundamental factor contributing to disc degeneration. Rather, the opposite was observed in our samples, as the endplate became progressively more porous with age and degeneration. Since ischemic disc cell behavior is commonly associated with degenerative change, this may be related to other factors such as the quality of vertebral capillaries, as opposed to decreased permeability of intervening tissues.

Visual grading of motion induced image degradation in high resolution peripheral computed tomography: impact of image quality on measures of bone density and micro-architecture

Motion artifacts are a common finding during high-resolution peripheral quantitative computed tomography (HR-pQCT) image acquisitions. To date it is not clear (i) when to repeat an acquisition, (ii) when to exclude a motion-degraded dataset post hoc, and (iii) how motion induced artifacts impact measures of trabecular and cortical parameters. In this study we present inter- and intra-observer reproducibility of a qualitative image quality grading score and report the prevalence of repeat acquisitions in our population. Finally the errors in bone density and micro-architectural parameters estimated from repeat acquisitions with and without motion degradation are presented. The relationship between these errors and the image quality grade is evaluated for each parameter. Repeat acquisitions performed due to operator-observed motion in the reconstructed image occurred for 22.7% of the exams (29.7% radius, 15.7% tibia). Of this subset, 88 exams with repeat acquisitions had at least one acquisition graded 1 (best quality). In this subset, the percent differences in bone density and micro-architecture measures tended to increase as the relative image quality decreased. Micro-architectural parameters were more sensitive to motion compared to geometric and densitometric parameters. These results provide estimates of the error in bone quality measures due to motion artifacts and provide an initial framework for developing standardized quality control criteria for cross-sectional and longitudinal HR-pQCT studies.

Imaging longitudinal changes in articular cartilage and bone following doxycycline treatment in a rabbit anterior cruciate ligament transection model of osteoarthritis

OBJECTIVE

The development of osteoarthritis following traumatic anterior cruciate ligament (ACL) injury is well established. However, few reliable indicators of early osteoarthritic changes have been established, which has limited the development of effective therapies. T(1ρ) and T(2) mapping techniques have the ability to provide highly accurate and quantitative measurements of articular cartilage degeneration in vivo. Relating these cartilaginous changes to high-resolution bone-densitometric evaluations of the late-stage osteoarthritic bone is crucial in elucidating the mechanisms of development of traumatic osteoarthritis (OA) and potential therapies for early- or late-stage intervention.

METHODS

Twelve rabbits were monitored with in vivo magnetic resonance imaging (MRI) scans following ACL transection surgery with a contralateral leg sham operation. Six of the rabbits were treated with oral doxycycline for the duration of the experiment. At 12 weeks, the excised knees from three animals from each group (n=6 overall) were subjected to micro-computed tomography (CT) analysis.

RESULTS

Consistent with previous studies, initial elevations in T(1ρ) and T(2) values in ACL-transected animals were observed with relative normalization towards values see in sham-operated legs over the 12-week study. This biphasic pattern could hold diagnostic potential to differentiate osteoarthritic cartilage by tracking the relative proportions of T(1ρ) and T(2) values as they rise with inflammation then fall as collagen and proteoglycan loss leads to further dehydration. The addition of doxycycline resulted in inconclusive, yet potentially interesting, cartilaginous changes in several compartments of the rabbit legs. Micro-CT studies demonstrated decreased bone densitometrics in ACL-transected knees. Correlation studies suggest that the cartilaginous changes may be associated with some aspects of bony change and the development of OA.

CONCLUSION

We conclude that there are definite relationships between cartilaginous changes as seen on MRI and late-stage microstructural bony changes after traumatic ACL injury in rabbits. In addition, doxycycline may show promise in mitigating early-stage cartilage damage that may serve to lessen late-stage osteoarthritic changes. This study demonstrates the ability to track OA progression and therapeutic efficacy with imaging modalities in vivo.

Trabecular reorganization in consecutive iliac crest biopsies when switching from bisphosphonate to strontium ranelate treatment

BACKGROUND

Several agents are available to treat osteoporosis while addressing patient-specific medical needs. Individuals' residual risk to severe fracture may require changes in treatment strategy. Data at osseous cellular and microstructural levels due to a therapy switch between agents with different modes of action are rare. Our study on a series of five consecutively taken bone biopsies from an osteoporotic individual over a six-year period analyzes changes in cellular characteristics, bone microstructure and mineralization caused by a therapy switch from an antiresorptive (bisphosphonate) to a dual action bone agent (strontium ranelate).

METHODOLOGY/PRINCIPAL FINDINGS

Biopsies were progressively taken from the iliac crest of a female patient. Four biopsies were taken during bisphosphonate therapy and one biopsy was taken after one year of strontium ranelate (SR) treatment. Furthermore, serum bone markers and dual x-ray absorptiometry measurements were acquired. Undecalcified histology was used to assess osteoid parameters and bone turnover. Structural indices and degree of mineralization were determined using microcomputed tomography, quantitative backscattered electron imaging, and combined energy dispersive x-ray/µ-x-ray-fluorescence microanalysis.

CONCLUSIONS/SIGNIFICANCE

Microstructural data revealed a notable increase in bone volume fraction after one year of SR treatment compared to the bisphosphonate treatment period. Indices of connectivity density, structure model index and trabecular bone pattern factor were predominantly enhanced indicating that the architectural transformation from trabecular rods to plates was responsible for the bone volume increase and less due to changes in trabecular thickness and number. Administration of SR following bisphosphonates led to a maintained mineralization profile with an uptake of strontium on the bone surface level. Reactivated osteoclasts designed tunneling, hook-like intratrabecular resorption sites. The appearance of tunneling resorption lacunae and the formation of both mini-modeling units and osteon-like structures within increased plate-like cancellous bone mass provides additional information on the mechanisms of strontium ranelate following bisphosphonate treatment, which may deserve special attention when monitoring a treatment switch.

Does vertebral bone marrow fat content correlate with abdominal adipose tissue, lumbar spine bone mineral density, and blood biomarkers in women with type 2 diabetes mellitus?

PURPOSE

To compare vertebral bone marrow fat content quantified with proton MR spectroscopy ((1)H-MRS) with the volume of abdominal adipose tissue, lumbar spine volumetric bone mineral density (vBMD), and blood biomarkers in postmenopausal women with and without type 2 diabetes mellitus (T2DM).

MATERIALS AND METHODS

Thirteen postmenopausal women with T2DM and 13 age- and body mass index-matched healthy controls were included in this study. All subjects underwent (1)H-MRS of L1-L3 to quantify vertebral bone marrow fat content (FC) and unsaturated lipid fraction (ULF). Quantitative computed tomography (QCT) was performed to assess vBMD of L1-L3. The volumes of abdominal subcutaneous/visceral/total adipose tissue were determined from the QCT images and adjusted for abdominal body volume (SAT(adj)/VAT(adj)/TAT(adj)). Fasting blood tests included plasma glucose and HbA1c.

RESULTS

Mean FC showed an inverse correlation with vBMD (r = -0.452; P < 0.05) in the whole study population. While mean FC was similar in the diabetic women and healthy controls (69.3 ± 7.5% versus 67.5 ± 6.1%; P > 0.05), mean ULF was significantly lower in the diabetic group (6.7 ± 1.0% versus 7.9 ± 1.6%; P < 0.05). SAT(adj) and TAT(adj) correlated significantly with mean FC in the whole study population (r = 0.538 and r = 0.466; P < 0.05). In contrast to the control group, significant correlations of mean FC with VAT(adj) and HbA1c were observed in the diabetic group (r = 0.642 and r = 0.825; P < 0.05).

CONCLUSION

This study demonstrated that vertebral bone marrow fat content correlates significantly with SAT(adj), TAT(adj), and lumbar spine vBMD in postmenopausal women with and without T2DM, but with VAT(adj) and HbA1c only in women with T2DM.

Teriparatide in bisphosphonate-resistant osteoporosis: microarchitectural changes and clinical results after 6 and 18 months

A number of osteoporotic patients under bisphosphonate treatment present persistent fragility fractures and bone loss despite good compliance. The objective of this 18-month prospective study was to investigate the effect of teriparatide [rhPTH(1-34)] in 25 female osteoporotics who were inadequate responders to oral bisphosphonates and to correlate microarchitectural changes in three consecutive iliac crest biopsies measured by micro-computed tomography (μCT) with bone mineral density (BMD) and bone serum markers. Scanned biopsies at baseline (M0), 6 months (M6), and 18 months (M18) demonstrated early significant (P < 0.01) increases in bone volume per tissue volume (+34%) and trabecular number (+14%) at M6 with only moderate changes in most μCT structural parameters between M6 and M18. μCT-measured bone tissue density was significantly decreased at M18, expressing an overall lower degree of tissue mineralization characteristic for new bone formation despite unchanged trabecular thickness due to increased intratrabecular tunneling at M18. μCT results were consistent with serum bone turnover markers, reaching maximal levels of bone alkaline phosphatase and serum β-crosslaps at M6, with subsequent decline until M18. BMD assessed by DXA demonstrated persistent increases at the lumbar spine until M12, whereas no significant change was observed at the hip. Type (alendronate/risedronate) and duration (3.5 ± 4 years) of prior bisphosphonate treatment did not influence outcome on μCT, BMD, or bone marker results. The overall results indicate a positive ceiling effect of teriparatide on bone microarchitecture and bone markers after 6 and 12 months for lumbar spine BMD, with no additional gain until M18 in bisphosphonate nonresponders.

High-resolution computed tomography for clinical imaging of bone microarchitecture

BACKGROUND

The role of bone structure, one component of bone quality, has emerged as a contributor to bone strength. The application of high-resolution imaging in evaluating bone structure has evolved from an in vitro technology for small specimens to an emerging clinical research tool for in vivo studies in humans. However, many technical and practical challenges remain to translate these techniques into established clinical outcomes.

QUESTIONS/PURPOSES

We reviewed use of high-resolution CT for evaluating trabecular microarchitecture and cortical ultrastructure of bone specimens ex vivo, extension of these techniques to in vivo human imaging studies, and recent studies involving application of high-resolution CT to characterize bone structure in the context of skeletal disease.

METHODS

We performed the literature review using PubMed and Google Scholar. Keywords included CT, MDCT, micro-CT, high-resolution peripheral CT, bone microarchitecture, and bone quality.

RESULTS

Specimens can be imaged by micro-CT at a resolution starting at 1 μm, but in vivo human imaging is restricted to a voxel size of 82 μm (with actual spatial resolution of ~ 130 μm) due to technical limitations and radiation dose considerations. Presently, this mode is limited to peripheral skeletal regions, such as the wrist and tibia. In contrast, multidetector CT can assess the central skeleton but incurs a higher radiation burden on the subject and provides lower resolution (200-500 μm).

CONCLUSIONS

CT currently provides quantitative measures of bone structure and may be used for estimating bone strength mathematically. The techniques may provide clinically relevant information by enhancing our understanding of fracture risk and establishing the efficacy of antifracture for osteoporosis and other bone metabolic disorders.

Ultrashort echo time MRI of cortical bone at 7 tesla field strength: a feasibility study

PURPOSE

To implement and examine the feasibility of a three-dimensional (3D) ultrashort TE (UTE) sequence on a 7 Tesla (T) clinical MR scanner in comparison with 3T MRI at high isotropic resolution.

MATERIALS AND METHODS

Using an in-house built saddle coil at both field strengths we have imaged mid-diaphysial sections of five fresh cadaveric specimens of the distal tibia. An additional in vivo scan was performed at 7 Tesla using a quadrature knee coil.

RESULTS

Using the same type of saddle coil at both field strengths, a significant increase in SNR at 7T compared with 3T (factor 1.7) was found. Significantly shorter T2* values were found at the higher field strength (T2* = 552.2 ± 126 μs at 7T versus T2* = 1163 ± 391 μs at 3T).

CONCLUSION

UHF MRI at 7T has great potential for imaging tissues with short T2.

Quantitative characterization of subject motion in HR-pQCT images of the distal radius and tibia

Image quality degradation due to subject motion is a common artifact affecting in vivo high-resolution peripheral quantitative computed tomography (HR-pQCT) of bones. These artifacts confound the accuracy and reproducibility of bone density, geometry, and cortical and trabecular structure measurements. Observer-based systems for grading image quality and criteria for deciding when to repeat an acquisition and post hoc data quality control remain highly subjective and non-standardized. This study proposes an objective, quantitative technique for measuring subject motion in HR-pQCT acquisitions from raw projection data, using image similarity measures applied to parallelized projections at 0° and 180°. A total of 88 HR-pQCT exams with repeated acquisitions of the distal radius (N = 54) or distal tibia (N = 34) of 49 women (age = 59 ± 14 year) and 3 men (46 ± 2 year) were retrospectively evaluated. All images were graded from 1 (no visible motion artifacts) to 5 (severe motion artifacts) according to the manufacturer-suggested image quality grading system. In addition, to serve as the reference case without motion artifacts, two cadaveric wrist and two ankle specimens were imaged twice with repositioning. The motion-induced error was calculated as the percent difference in each bone parameter for the paired scans with and without visually apparent motion artifacts. Quantitative motion estimates (QMEs) for each motion-degraded scan were calculated using two different image similarity measures: sum of squared differences (SSD) and normalized cross-correlation (NCC). The mean values of QME(SSD) and QME(NCC) increased with the image quality grade for both radius and tibia. Quality grades were differentiated between grades 2 and 3 using QME(SSD), but not with QME(NCC), in addition to between grades 4 and 5. Both QMEs correlated significantly to the motion-induced errors in the measurements and their empirical relationship was derived. Subject motion had greater impact on the precision of trabecular structure indices than on the densitometric indices. The results of this study may provide a basis for establishing a threshold for motion artifacts in accordance to the study design as well as a standardized quality control protocol across operators and imaging centers.

Longitudinal evaluation of the effects of alendronate on MRI bone microarchitecture in postmenopausal osteopenic women

We evaluated longitudinal effects of alendronate on MRI-based trabecular bone structure parameters derived from dual thresholding and fuzzy clustering (BE-FCM) trabecular bone segmentation. Treatment effects were observed in the distal tibia after 24 months. The BE-FCM method increased correlations to HR-pQCT-based parameters.

INTRODUCTION

High-resolution magnetic resonance imaging (MRI) allows for non-invasive bone microarchitecture analysis. The goal of this study was to examine the potential of MRI-based trabecular bone structure parameters to monitor effects of alendronate in humans in vivo, and to compare the results to HR-pQCT and DXA measurements.

MATERIALS AND METHODS

Postmenopausal osteopenic women were divided into alendronate treatment and control groups, and imaged at baseline, 12 months, and 24 months (n = 52 at baseline) using 3T MRI, HR-pQCT, and DXA. Image acquisition sites included distal tibia (MRI and HR-pQCT), distal radius (MRI, DXA, and HR-pQCT), and the proximal femur (MRI and DXA). Two different regions of interest were evaluated. One contained the trabecular bone region within the entire MRI acquisition, and the second contained a subregion matched to the region contained in the HR-pQCT acquisition. The trabecular bone was segmented using two different methods; dual thresholding and BE-FCM. Trabecular bone structure parameters included bone volume fraction (BV/TV), number (Tb.N), spacing (Tb.Sp), and thickness (Tb.Th), along with seven geodesic topological analysis (GTA) parameters. Longitudinal changes and correlations to HR-pQCT and DXA measurements were evaluated.

RESULTS

Apparent Tb.N and four GTA parameters showed treatment effects (p < 0.05) in the distal tibia after 24 months in the entire MRI region using BE-FCM, as well as Tb.N using dual thresholding. No treatment effects after 24 months were observed in the HR-pQCT or in MRI analysis for the HR-pQCT-matched regions. Apparent BV/TV and Tb.N from BE-FCM had significantly higher correlations to HR-pQCT values compared to those derived from thresholding.

CONCLUSIONS

This study demonstrates the influence of computational methods and region of interest definitions on measurements of trabecular bone structure, and the feasibility of MRI-based quantification of longitudinal changes in bone microarchitecture due to bisphosphonate therapy. The results suggest that there may be a need to reevaluate the current standard HR-pQCT region definition for increased treatment sensitivity.

A longitudinal HR-pQCT study of alendronate treatment in postmenopausal women with low bone density: Relations among density, cortical and trabecular microarchitecture, biomechanics, and bone turnover

The goal of this study was to characterize longitudinal changes in bone microarchitecture and function in women treated with an established antifracture therapeutic. In this double-blind, placebo-controlled pilot study, 53 early postmenopausal women with low bone density (age = 56 ± 4 years; femoral neck T-score = -1.5 ± 0.6) were monitored by high-resolution peripheral quantitative computed tomography (HR-pQCT) for 24 months following randomization to alendronate (ALN) or placebo (PBO) treatment groups. Subjects underwent annual HR-pQCT imaging of the distal radius and tibia, dual-energy X-ray absorptiometry (DXA), and determination of biochemical markers of bone turnover (BSAP and uNTx). In addition to bone density and microarchitecture assessment, regional analysis, cortical porosity quantification, and micro-finite-element analysis were performed. After 24 months of treatment, at the distal tibia but not the radius, HR-pQCT measures showed significant improvements over baseline in the ALN group, particularly densitometric measures in the cortical and trabecular compartments and endocortical geometry (cortical thickness and area, medullary area) (p < .05). Cortical volumetric bone mineral density (vBMD) in the tibia alone showed a significant difference between treatment groups after 24 months (p < .05); however, regionally, significant differences in Tb.vBMD, Tb.N, and Ct.Th were found for the lateral quadrant of the radius (p < .05). Spearman correlation analysis revealed that the biomechanical response to ALN in the radius and tibia was specifically associated with changes in trabecular microarchitecture (|ρ| = 0.51 to 0.80, p < .05), whereas PBO progression of bone loss was associated with a broad range of changes in density, geometry, and microarchitecture (|ρ| = 0.56 to 0.89, p < .05). Baseline cortical geometry and porosity measures best predicted ALN-induced change in biomechanics at both sites (ρ > 0.48, p < .05). These findings suggest a more pronounced response to ALN in the tibia than in the radius, driven by trabecular and endocortical changes.

Variations in morphological and biomechanical indices at the distal radius in subjects with identical BMD

Determination of osteoporotic status is based primarily on areal bone mineral density (aBMD) obtained through dual X-ray absorptiometry (DXA). However, many fractures occur in patients with T-scores above the WHO threshold of osteoporosis, in part because DXA measures are insensitive to biomechanically important alterations in bone quality. The goal of this study was to determine--within groups of subjects with identical radius aBMD values--the extant variation in densitometric, geometric, microstructural, and biomechanical parameters. High resolution peripheral quantitative computed tomography (HR-pQCT) and DXA radius data from males and females spanning large ranges in age, osteoporotic status, and anthropometrics were compiled. 262 distal radius datasets were processed for this study. HR-pQCT scans were analyzed according to the manufacturer's standard clinical protocol to quantify densitometric, geometric, and microstructural indices. Micro-finite element analysis was performed to calculate biomechanical indices. Factor of risk of wrist fracture was calculated. Simulated aBMD calculated from HR-pQCT data was used to group scans for evaluation of variation in quantified indices. Indices reflecting the greatest variation within aBMD level were BMD in the central portion of the trabecular compartment (max CV 142), trabecular heterogeneity (max CV 90), and intra-cortical porosity (max CV 151). Of the biomechanical indices, cortical load fraction had the greatest variation (max CV 38). Substantial variations in indices reflecting density, structure, and biomechanical competence exist among subjects with identical aBMD levels. Overlap of these indices among osteoporotic status groups reflects the reported incidence of osteoporotic fracture in subjects classified as osteopenic or normal.

High-resolution peripheral quantitative computed tomographic imaging of cortical and trabecular bone microarchitecture in patients with type 2 diabetes mellitus

CONTEXT

Cross-sectional epidemiological studies have found that patients with type 2 diabetes mellitus (T2DM) have a higher incidence of certain fragility fractures despite normal or elevated bone mineral density (BMD).

OBJECTIVE

In this study, high-resolution peripheral quantitative computed tomography was applied to characterize cortical and trabecular microarchitecture and biomechanics in the peripheral skeleton of female patients with T2DM.

DESIGN AND SETTING

A cross-sectional study was conducted in patients with T2DM recruited from a diabetic outpatient clinic.

PARTICIPANTS

Elderly female patients (age, 62.9 ± 7.7 yr) with a history of T2DM (n = 19) and age- and height-matched controls (n = 19) were recruited.

OUTCOME MEASURES

Subjects were imaged using high-resolution peripheral quantitative computed tomography at the distal radius and tibia. Quantitative measures of volumetric (BMD), cross-sectional geometry, trabecular and cortical microarchitecture were calculated. Additionally, compressive mechanical properties were determined by micro-finite element analysis.

RESULTS

Compared to the controls, the T2DM cohort had 10% higher trabecular volumetric BMD (P < 0.05) adjacent to the cortex and higher trabecular thickness in the tibia (13.8%; P < 0.05). Cortical porosity differences alone were consistent with impaired bone strength and were significant in the radius (>+50%; P < 0.05), whereas pore volume approached significance in the tibia (+118%; P = 0.1).

CONCLUSION

The results of this pilot investigation provide a potential explanation for the inability of standard BMD measures to explain the elevated fracture incidence in patients with T2DM. The findings suggest that T2DM may be associated with impaired resistance to bending loads due to inefficient redistribution of bone mass, characterized by loss of intracortical bone offset by an elevation in trabecular bone density.

Reproducibility of direct quantitative measures of cortical bone microarchitecture of the distal radius and tibia by HR-pQCT

Quantitative cortical microarchitectural end points are important for understanding structure-function relations in the context of fracture risk and therapeutic efficacy. This technique study details new image-processing methods to automatically segment and directly quantify cortical density, geometry, and microarchitecture from HR-pQCT images of the distal radius and tibia. An automated segmentation technique was developed to identify the periosteal and endosteal margins of the distal radius and tibia and detect intracortical pore space morphologically consistent with Haversian canals. The reproducibility of direct quantitative cortical bone indices based on this method was assessed in a pooled data set of 56 subjects with two repeat acquisitions for each site. The in vivo precision error was characterized using root mean square coefficient of variation (RMSCV%) from which the least significant change (LSC) was calculated. Bland-Altman plots were used to characterize bias in the precision estimates. The reproducibility of cortical density and cross-sectional area measures was high (RMSCV <1% and <1.5%, respectively) with good agreement between young and elder medians. The LSC for cortical porosity (Ct.Po) was somewhat smaller in the radius (0.58%) compared with the distal tibia (0.84%) and significantly different between young and elder medians in the distal tibia (LSC: 0.75% vs. 0.92%, p<0.001). The LSC for pore diameter and distribution (Po.Dm and Po.Dm.SD) ranged between 15 and 23 microm. Bland-Altman analysis revealed moderate bias for integral measures of area and volume but not for density or microarchitecture. This study indicates that HR-pQCT measures of cortical bone density and architecture can be measured in vivo with high reproducibility and limited bias across a biologically relevant range of values. The results of this study provide informative data for the design of future clinical studies of bone quality.

High-resolution imaging techniques for the assessment of osteoporosis

The importance of assessing the bone's microarchitectural make-up in addition to its mineral density in the context of osteoporosis has been emphasized in several publications. The high spatial resolution required to resolve the bone's microstructure in a clinically feasible scan time is challenging. At present, the best suited modalities meeting these requirements in vivo are high-resolution peripheral quantitative imaging (HR-pQCT) and magnetic resonance imaging (MRI). Whereas HR-pQCT is limited to peripheral skeleton regions like the wrist and ankle, MRI can also image other sites like the proximal femur but usually with lower spatial resolution. In addition, multidetector computed tomography has been used for high-resolution imaging of trabecular bone structure; however, the radiation dose is a limiting factor. This article provides an overview of the different modalities, technical requirements, and recent developments in this emerging field. Details regarding imaging protocols as well as image postprocessing methods for bone structure quantification are discussed.

Regional variations of gender-specific and age-related differences in trabecular bone structure of the distal radius and tibia

Regional variation in trabecular structure across axial sections is often obscured by the conventional global analysis, which takes an average value for the entire trabecular compartment. The objective of this study is to characterize spatial variability in trabecular structure within a cross-section at the distal radius and tibia, and gender and age effects using in vivo high-resolution peripheral quantitative computed tomography (HR-pQCT). HR-pQCT images of the distal radius and tibia were acquired from 146 healthy individuals aged 20-78 years. Trabecular bone volume fraction (BV/TV), number (Tb.N), thickness (Tb.Th), separation (Tb.Sp), and heterogeneity (Tb.1/N.SD) were obtained in a total of 11 regions-the entire trabecular compartment (the global means), inner, outer, and eight defined subregions. Regional variations were examined with respect to the global means, and compared between women and men, and between young (20-29 years old) and elderly (65-79 years old) adults. Substantial regional variations in trabecular bone structure at the distal radius and tibia were revealed (e.g. BV/TV varied -40% to +57% and -59% to +100% of the global means, respectively, for elderly women). The inner-lateral (IL) subregion had low BV/TV, Tb.N, and Tb.Th, and low Tb.Sp and Tb.1/N.SD at both sites; the opposite was true in the outer-anterior (OA) subregion at the distal radius and the outer-medial (OM) and -posterior (OP) subregions at the distal tibia. Gender differences were most pronounced in the inner-anterior (IA) subregion compared to the other regions or the global mean differences at both sites. Trabecular structure associated with age and differed between young and elderly adults predominantly in the inner-posterior (IP) subregion at the distal radius and in the IL and IA subregions at the distal tibia; on the other hand, it remained unchanged in the OA subregion at the distal radius and in the OM subregion at the distal tibia for both women and men. This study demonstrated that not only the conventional global analysis can obscure regional differences, but also assuming bone status from that of smaller subregion may introduce a confounding sampling error. Therefore, a combined approach of investigating the entire region, each subregion, and the cortical compartment may offer more complete information.

Higher doses of bisphosphonates further improve bone mass, architecture, and strength but not the tissue material properties in aged rats

We report the results of a series of experiments designed to determine the effects of ibandronate (Ibn) and risedronate (Ris) on a number of bone quality parameters in aged osteopenic rats to explain how bone material and bone mass may be affected by the dose of bisphosphonates (BP) and contribute to their anti-fracture efficacy. Eighteen-month old female rats underwent either ovariectomy or sham surgery. The ovariectomized (OVX) groups were left untreated for 2 months to develop osteopenia. Treatments started at 20 months of age as follows: sham and OVX control (treated with saline), OVX + risedronate 30 and 90 (30 or 90 microg/kg/dose), and OVX + ibandronate 30 and 90 (30 or 90 microg/kg/dose). The treatments were given monthly for 4 months by subcutaneous injection. At sacrifice at 24 months of age the 4th lumbar vertebra was used for microCT scans (bone mass, architecture, and degree of mineralization of bone, DMB) and histomorphometry, and the 6th lumbar vertebra, tibia, and femur were collected for biomechanical testing to determine bone structural and material strength, cortical fracture toughness, and tissue elastic modulus. The compression testing of the vertebral bodies (LVB6) was simulated using finite-element analysis (FEA) to also estimate the bone structural stiffness. Both Ibn and Ris dose-dependently increased bone mass and improved vertebral bone microarchitecture and mechanical properties compared to OVX control. Estimates of vertebral maximum stress from FEA were correlated with vertebral maximum load (r=0.5, p<0.001) and maximum stress (r=0.4, p<0.005) measured experimentally. Tibial bone bending modulus and cortical strength increased compared to OVX with both BP but no dose-dependent effect was observed. DMB and elastic modulus of trabecular bone were improved with Ibn 30 compared to OVX but were not affected in other BP-treated groups. DMB of tibial cortical bone showed no change with BP treatments. The fracture toughness examined in midshaft femurs did not change with BP even with the higher doses. In summary, the anti-fracture efficacy of BP is largely due to their preservation of bone mass and while the higher doses further improve the bone structural properties do not improve the localized bone material characteristics such as tissue strength, elastic modulus, and cortical toughness.

Contribution of the intra-specimen variations in tissue mineralization to PTH- and raloxifene-induced changes in stiffness of rat vertebrae

The intra-specimen spatial variation in mineralization of bone tissue can be changed by drug treatments that alter bone remodeling. However, the contribution of such changes to the overall biomechanical effect of a treatment on bone strength is not known. To provide insight into this issue, we used a rat model to determine the effects of ovariectomy, parathyroid hormone, and raloxifene (vs. sham) on the contribution of spatial variations in mineralization to treatment-induced changes in vertebral stiffness. Mineral density was measured from 6-microm voxel-sized quantitative micro-CT scans. Whole-vertebral and trabecular stiffness values were estimated using finite element analysis of these micro-CT scans, first including all intra-specimen variations in mineral density in the model and then excluding such variations by using a specimen-specific average density throughout each specimen. As expected, we found appreciable effects of treatment on overall bone stiffness, the effect being greater for the trabecular compartment (up to 52% reduction vs. sham, p<0.0001) than the whole vertebra (p=0.055). Intra-specimen mean mineralization was not changed with treatment but the intra-specimen variation in mineralization was, although the effect was small (4%). Intra-specimen spatial variations in mineralization accounted for 10-12% and 5-6% of overall stiffness of the trabecular compartment and whole vertebral body, respectively. However, after accounting for all treatment effects on bone geometry and trabecular microstructure, any treatment effects due to changes in mineralization were negligible (<2%), although statistically detectable (p<0.02). We conclude that, despite a role in the general biomechanical behavior of bone, the spatial variations in tissue mineralization, as measured by quantitative micro-CT, did not appreciably contribute to ovariectomy-, PTH-, or raloxifene-induced changes in stiffness of the whole bone or the trabecular compartment in these rat vertebrae.

Automated simulation of areal bone mineral density assessment in the distal radius from high-resolution peripheral quantitative computed tomography

SUMMARY

An automated image processing method is presented for simulating areal bone mineral density measures using high-resolution peripheral quantitative computed tomography (HR-pQCT) in the ultra-distal radius. The accuracy of the method is validated against clinical dual X-ray absorptiometry (DXA). This technique represents a useful reference to gauge the utility of novel 3D quantification methods applied to HR-pQCT in multi-center clinical studies and potentially negates the need for separate forearm DXA measurements.

INTRODUCTION

Osteoporotic status is primarily assessed by measuring areal bone mineral density (aBMD) using 2D dual X-ray absorptiometry (DXA). However, this technique does not sufficiently explain bone strength and fracture risk. High-resolution peripheral quantitative computed tomography (HR-pQCT) has been introduced as a method to quantify 3D bone microstructure and biomechanics. In this study, an automated method is proposed to simulate aBMD measures from HR-pQCT distal radius images.

METHODS

A total of 117 subject scans were retrospectively analyzed from two clinical bone quality studies. The distal radius was imaged by HR-pQCT and DXA on one of two devices (Hologic or Lunar). Areal BMD was calculated by simulation from HR-pQCT images (aBMD(sim)) and by standard DXA analysis (aBMD(dxa)).

RESULTS

The reproducibility of the simulation technique was 1.1% (root mean-squared coefficient of variation). HR-pQCT-based aBMD(sim) correlated strongly to aBMD(dxa) (Hologic: R (2) = 0.82, Lunar: R (2) = 0.87), though aBMD(sim) underestimated aBMD(dxa) for both DXA devices (p < 0.0001). Finally, aBMD(sim) predicted aBMD at the proximal femur and lumbar spine with equal power compared to aBMD(dxa).

CONCLUSION

The results demonstrate that aBMD can be simulated from HR-pQCT images of the distal radius. This approach has the potential to serve as a surrogate forearm aBMD measure for clinical HR-pQCT studies when axial bone mineral density values are not required.

Accuracy of volumetric bone mineral density measurement in high-resolution peripheral quantitative computed tomography

Accurate bone mineral density (BMD) quantification is critical in clinical assessment of fracture risk and in the research of age-, disease-, and treatment-related musculoskeletal changes. The development of high-resolution peripheral quantitative computed tomography (HR-pQCT) imaging has made possible in vivo assessment of compartmental volumetric BMD (vBMD) and bone micro-architecture in the distal radius and tibia. HR-pQCT imaging relies on a polychromatic X-ray source and therefore is subject to beam hardening as well as scatter artifacts. In light of these limitations, we hypothesize that the accuracy of HR-pQCT vBMD measurement in the trabecular compartment (vBMD(trab)) is not independent of bone density and geometry, but rather influenced by variations in trabecular bone volume fraction and cortical thickness. The goal of this study, therefore, was to evaluate the accuracy of HR-pQCT vBMD(trab) measurement in the radius and tibia, and to determine the dependence of this measurement on geometric and densitometric parameters. Our approach was to use a series of idealized hydroxyapatite (HA) phantoms with varying densities and geometries to quantify the accuracy of HR-pQCT analysis. Two sets of custom-made HA phantoms designed to mimic the distal tibia and distal radius were manufactured. Geometric and densitometric specifications were based on a dataset of healthy volunteers and osteopenic patients. Multiple beam hardening correction (BHC) algorithms were implemented and evaluated in their ability to reduce measurement error. Substantial errors in measured vBMD(trab) were found. Overestimation of vBMD(trab) increased proportional to cortical shell thickness and decreased proportional to insert density. The most pronounced vBMD(trab) overestimation therefore occurred in the phantoms with the lowest insert densities and highest shell thickness, where error was as high as 20 mg HA/cm3 (33%) in the radius phantom and 25 mg HA/cm(3) (41%) in the tibia phantom. Error in vBMD(trab) propagates to the calculation of micro-architectural measures; 41% error in vBMD(trab) will produce 41% error in volume fraction (BV/TV) and trabecular thickness (Tb.Th), and 5% error in trabecular separation (Tb.Sp). BHC algorithms supplied by the manufacturer failed to eliminate these errors. Our results confirm that geometric and densitometric variations influence the accuracy of HR-pQCT vBMD(trab) measurements, and must be considered when interpreting data across populations or time-points.

Assessment of trabecular bone structure using MDCT: comparison of 64- and 320-slice CT using HR-pQCT as the reference standard

Objectives

The aim of our study was to perform trabecular bone structure analysis with images from 64- and 320-slice multidetector computed tomography (MDCT) and to compare these with high-resolution peripheral computed tomography (HR-pQCT).

Materials and methods

Twenty human cadaver distal forearm specimens were imaged on a 64- and 320-slice MDCT system at 120 kVp, 200 mA and 135 kVp, 400 mA (in-plane pixel size 234 µm; slice thickness 500 µm). HR-pQCT imaging was performed at an isotropic voxel size of 41 µm. Bone volume fraction (BV/TV), trabecular number (Tb.N), thickness (Tb.Th) and separation (Tb.Sp) were computed.

Results

MDCT-derived BV/TV and Tb.Sp were highly correlated (r = 0.92–0.96, p < 0.0001) with the corresponding HR-pQCT parameters. Tb.Th was the only structure measure that did not yield any significant correlation.

Conclusion

The 64- and 320-slice MDCT systems both perform equally well in depicting trabecular bone architecture. However, because of constrained resolutions accurate derivation of trabecular bone measures is limited to only a subset of microarchitectural parameters.

Assessment of trabecular bone structure of the calcaneus using multi-detector CT: correlation with microCT and biomechanical testing

The prediction of bone strength can be improved when determining bone mineral density (BMD) in combination with measures of trabecular microarchitecture. The goal of this study was to assess parameters of trabecular bone structure and texture of the calcaneus by clinical multi-detector row computed tomography (MDCT) in an experimental in situ setup and to correlate these parameters with microCT (microCT) and biomechanical testing. Thirty calcanei in 15 intact cadavers were scanned using three different protocols on a 64-slice MDCT scanner with an in-plane pixel size of 208 microm and 500 microm slice thickness. Bone cores were harvested from each specimen and microCT images with a voxel size of 16 microm were obtained. After image coregistration, trabecular bone structure and texture were evaluated in identical regions on the MDCT images. After data acquisition, uniaxial compression testing was performed. Significant correlations between MDCT- and microCT-derived measures of bone volume fraction (BV/TV), trabecular thickness (Tb.Th) and trabecular separation (Tb.Sp) were found (range, R(2)=0.19-0.65, p<0.01 or 0.05). The MDCT-derived parameters of volumetric BMD, app. BV/TV, app. Tb.Th and app. Tb.Sp were capable of predicting 60%, 63%, 53% and 25% of the variation in bone strength (p<0.01). When combining those measures with one additional texture index (either GLCM, TOGLCM or MF.euler), prediction of mechanical competence was significantly improved to 86%, 85%, 71% and 63% (p<0.01). In conclusion, this study showed the feasibility of trabecular microarchitecture assessment using MDCT in an experimental setup simulating the clinical situation. Multivariate models of BMD or structural parameters combined with texture indices improved prediction of bone strength significantly and might provide more reliable estimates of fracture risk in patients.

Assessment of bone tissue mineralization by conventional x-ray microcomputed tomography: comparison with synchrotron radiation microcomputed tomography and ash measurements

Assessment of bone tissue mineral density (TMD) may provide information critical to the understanding of mineralization processes and bone biomechanics. High-resolution three-dimensional assessment of TMD has recently been demonstrated using synchrotron radiation microcomputed tomography (SRmuCT); however, this imaging modality is relatively inaccessible due to the scarcity of SR facilities. Conventional desktop muCT systems are widely available and have been used extensively to assess bone microarchitecture. However, the polychromatic source and cone-shaped beam geometry complicate assessment of TMD by conventional muCT. The goal of this study was to evaluate muCT-based measurement of degree and distribution of tissue mineralization in a quantitative, spatially resolved manner. Specifically, muCT measures of bone mineral content (BMC) and TMD were compared to those obtained by SRmuCT and gravimetric methods. Cylinders of trabecular bone were machined from human femoral heads (n = 5), vertebrae (n = 5), and proximal tibiae (n = 4). Cylinders were imaged in saline on a polychromatic muCT system at an isotropic voxel size of 8 microm. Volumes were reconstructed using beam hardening correction algorithms based on hydroxyapatite (HA)-resin wedge phantoms of 200 and 1200 mg HA/cm3. SRmuCT imaging was performed at an isotropic voxel size of 7.50 microm at the National Synchrotron Light Source. Attenuation values were converted to HA concentration using a linear regression derived by imaging a calibration phantom. Architecture and mineralization parameters were calculated from the image data. Specimens were processed using gravimetric methods to determine ash mass and density, muCT-based BMC values were not affected by altering the beam hardening correction. Volume-averaged TMD values calculated by the two corrections were significantly different (p = 0.008) in high volume fraction specimens only, with the 1200 mg HA/cm3 correction resulting in a 4.7% higher TMD value. MuCT and SRmuCT provided significantly different measurements of both BMC and TMD (p < 0.05). In high volume fraction specimens, muCT with 1200 mg HA/cm3 correctionteg resulted in BMC and TMD values 16.7% and 15.0% lower, respectively, than SRmuCT values. In low volume fraction specimens, muCT with 1200 mg HA/cm3 correction resulted in BMC and TMD values 12.8% and 12.9% lower, respectively, than SRmuCT values. MuCT and SRmuCT values were well-correlated when volume fraction groups were considered individually (BMC R2 = 0.97-1.00; TMD R2 = 0.78-0.99). Ash mass and density were higher than the SRmuCT equivalents by 8.6% in high volume fraction specimens and 10.9% in low volume fraction specimens (p < 0.05). BMC values calculated by tomography were highly correlated with ash mass (ash versus muCT R2 = 0.96-1.00; ash versus SRmuCT R2 = 0.99-1.00). TMD values calculated by tomography were moderately correlated with ash density (ash versus muCT R2 = 0.64-0.72; ash versus SRmuCT R2 = 0.64). Spatially resolved comparisons highlighted substantial geometric nonuniformity in the muCT data, which were reduced (but not eliminated) using the 1200 mg HA/cm3 beam hardening correction, and did not exist in the SRmuCT data. This study represents the first quantitative comparison of muCT mineralization evaluation against SRnuCT and gravimetry. Our results indicate that muCT mineralization measures are underestimated but well-correlated with SRmuCT and gravimetric data, particularly when volume fraction groups are considered individually.

In vivo ultra-high-field magnetic resonance imaging of trabecular bone microarchitecture at 7 T

PURPOSE

To investigate the feasibility of 7 T magnetic resonance imaging (MRI) to visualize and quantify trabecular bone structure in vivo by comparison with 3T MRI and in vivo three-dimensional (3D) high-resolution peripheral quantitative computed tomography (HR-pQCT).

MATERIALS AND METHODS

The distal tibiae of 10 healthy volunteers were imaged. Therefore, fully balanced steady state free precession (bSSFP) and spin-echo (bSSSE) pulse sequences were implemented and optimized for 7 T. Structural bone parameters, such as apparent bone-volume over total-volume fraction (app.BV/TV), apparent trabecular plate separation (app.TbSp), apparent trabecular plate thickness (app.TbTh), and apparent trabecular plate number (app.TbN), were derived.

RESULTS

All structural trabecular bone parameters correlated well (r > 0.6) between 7T and 3T, and between 7 T and HR-pQCT (r > 0.69), with the exception of app.TbTh, which correlated modestly (r = 0.41) between field strengths and very low with HR-pQCT (r < 0.16). Regarding absolute values, app.TbN varied only 4% between field strengths, and only 0.6% between 7 T and HR-pQCT. App.TbSp correlated best between 7 T and HR-pQCT (r = 0.89). Using bSSSE, significant smaller trabecular thickness and significant higher trabecular number were found compared to bSSFP.

CONCLUSION

We concluded that imaging and quantification of the trabecular bone architecture at 7 T is feasible and preferably done using bSSSE. There exists great potential for ultra-high-field (UHF) MRI applied to trabecular bone measurements.

Quantitative assessment of bone tissue mineralization with polychromatic micro-computed tomography

Micro-computed tomography (microCT) has become an important tool for morphological characterization of cortical and trabecular bone. Quantitative assessment of bone tissue mineral density (TMD) from microCT images may be possible; however, the methods for calibration and accuracy have not been thoroughly evaluated. This study investigated hydroxyapatite (HA) phantom sampling limitations, short-term reproducibility of phantom measurements, and accuracy of TMD measurements by correlation to ash density. Additionally, the performance of a global and a local threshold for determining TMD was tested. The full length of a commercial density phantom was imaged by microCT, and mean calibration parameters were determined for a volume of interest (VOI) at 10 random positions along the longitudinal axis. Ten different VOI lengths were used (0.9-13 mm). The root mean square error (RMSE) was calculated for each scan length. Short-term reproducibility was assessed by five repeat phantom measurements for three source voltage settings. Accuracy was evaluated by imaging rat cortical bone (n = 16) and bovine trabecular bone (n = 15), followed by ash gravimetry. Phantom heterogeneity was associated with <0.5% RMSE. The coefficient of variation for five repeat measurements was generally <0.25% across all energies and phantom densities. Bone mineral content was strongly correlated to ash weight (R (2) = 1.00 for both specimen groups and both threshold methods). Ash density was well correlated for the trabecular bone specimens (R (2) > 0.80). In cortical bone specimens, the correlation was somewhat weaker when a global threshold was applied (R (2) = 0.67) compared to the local threshold method (R (2) = 0.78).

Assessment of trabecular bone structure comparing magnetic resonance imaging at 3 Tesla with high-resolution peripheral quantitative computed tomography ex vivo and in vivo

In vivo high-resolution peripheral quantitative micro-CT (HR-pQCT) is a new modality for imaging peripheral sites like the distal tibia and the distal radius, providing structural bone parameters. Comparing HR-pQCT with MRI, we found that both modalities are capable of offering meaningful information on trabecular structure.

BACKGROUND

Magnetic resonance imaging (MRI) has emerged as the leading in vivo method for measuring trabecular bone micro-architecture and providing structural information. Recently, an in vivo HR-pQCT modality was introduced for imaging peripheral sites like the distal tibia and the distal radius, providing structural bone parameters. The goal of this work was to compare and evaluate the performances and in vivo capabilities of HR-pQCT in comparison with MRI at 3 Tesla.

METHODS

To this end images of 8 human specimens (5 tibiae and 3 radii) and 11 participants (6 tibia and 5 radii) were acquired with both modalities. Additionally, the radius specimens were scanned with micro-CT (muCT), which was used as a standard of reference. Structural parameters calculated from MRI were compared with results from HR-pQCT images and additionally muCT for the radii specimens.

RESULTS

High correlations (r > 0.7) were found for trabecular number and trabecular spacing between the two modalities in vivo and ex vivo. 2D and 3D analysis revealed high correlations (r > 0.8) in structural bone parameters for all measurements. Using micro-CT as standard of reference both results from QCT and MRI correlated well.

CONCLUSION

Both imaging modalities were found to perform equally well regarding trabecular bone measurements.

Resolution dependence of the non-metric trabecular structure indices

Non-metric indices of topological features of trabecular bone structure, such as structure model index (SMI), connectivity density (Conn.D), and degree of anisotropy (DA), provide unique information relevant to bone quality. With recent technological advancement, in vivo assessment of these indices may be possible from images acquired using high-resolution imaging techniques such as high-resolution peripheral quantitative computed tomography (HR-pQCT). However, more detailed investigation of the dependence of non-metric indices on spatial resolution is needed to determine their applicability. The purpose of this study was to determine whether these three non-metric indices are affected by the spatial resolution of CT images. First, the SMI, Conn.D, and DA were calculated for trabecular bone specimens with varying plate-like and rod-like structures from resampled muCT images across a range of spatial resolutions and compared to the reference values. To account for differences in size across different species and anatomical sites, the results are reported in normalized resolution units. Next, the impact of resolution on the non-metric indices for cores of human distal tibia trabecular bone from clinical HR-pQCT images was evaluated to determine the applicability of the non-metric indices to in vivo imaging. We found that the non-metric indices of trabecular bone structure were affected by spatial resolution of CT images. Particularly, the SMI deviated from the high-resolution muCT reference value depending on the structure type, whether plate-like or rod-like. Both Conn.D and DA were underestimated in the images obtained at an in vivo resolution. It is not trivial to determine absolute threshold for validity of these non-metric indices without considering a specific study design (e.g. relative resolution, the size of the treatment effect to detect, and specimen type). The results of this study provide an upper bound for the accuracy of the non-metric indices under limited resolution scenarios.

The effects of geometric and threshold definitions on cortical bone metrics assessed by in vivo high-resolution peripheral quantitative computed tomography

This study evaluates in vivo methods for calculating cortical thickness (Ct.Th) with respect to sensitivity to tissue-level changes in mineralization and the ability to predict whole-bone mechanical properties. Distal radial and tibial images obtained from normal volunteers using high-resolution peripheral quantitative computed tomography (HR-pQCT) were segmented using three thresholds including the manufacturer default and +/-5% in terms of equivalent mineral density. Ct.Th was determined in two ways: using a direct three-dimensional (3D) method and using an annular method, where cortical bone volume is divided by periosteal surface area. D(comp) (mg HA/cm(3)) was calculated based on the mean density-calibrated linear attenuation values of the cortical compartment. Finite element analysis was performed to evaluate the predictive ability of the annular and direct Ct.Th methods. Using the direct 3D method, a +/-5% change in threshold resulted in a 2% mean difference in Ct.Th for both the radius and tibia. An average difference of 5% was found using the annular method. The change in threshold produced changes in D(comp) ranging 0.50-1.56% for both the tibia and radius. Annular Ct.Th correlated more strongly with whole-bone apparent modulus (R(2)=0.64 vs. R(2)=0.41). Both thickness calculation methods and threshold selection have a direct impact on cortical parameters derived from HR-pQCT images. Indirectly, these results suggest that moderate changes in tissue-level mineralization can affect cortical measures. Furthermore, while the direct 3D Ct.Th method is less sensitive to threshold effects, both methods are moderate predictors of mechanical strength, with the annular method being the stronger correlate.

A local adaptive threshold strategy for high resolution peripheral quantitative computed tomography of trabecular bone

High resolution peripheral quantitative computed tomography (HR-pQCT) is a promising method for detailed in vivo 3D characterization of the densitometric, geometric, and microstructural features of human bone. Currently, a hybrid densitometric, direct, and plate model-based calculation is used to quantify trabecular microstructure. In the present study, this legacy methodology is compared to direct methods derived from a new local thresholding scheme independent of densitometric and model assumptions. Human femoral trabecular bone samples were acquired from patients undergoing hip replacement surgery. HR-pQCT (82 microm isotropic voxels) and micro-tomography (16 microm isotropic voxels) images were acquired. HR-pQCT images were segmented and analyzed in three ways: (1) using the hybrid method provided by the manufacturer based on a fixed global threshold, (2) using direct 3D methods based on the fixed global threshold segmentation, and (3) using direct 3D methods based on a novel local threshold scheme. The results were compared against standard direct 3D indices from microCT analysis. Standard trabecular parameters determined by HR-pQCT correlated strongly to microCT. BV/TV and Tb.Th were significantly underestimated by the hybrid method and significantly overestimated by direct methods based on the global threshold segmentation while the local method yielded optimal intermediate results. The direct-local method also performed favorably for Tb.N (R(2) = 0.85 vs. R(2) = 0.70 for direct-global method) and Tb.Sp (R(2) = 0.93 vs. R(2) = 0.85 for the hybrid method and R(2) = 0.87 for the direct-global method). These results indicate that direct methods, with the aid of advanced segmentation techniques, may yield equivalent or improved accuracy for quantification of trabecular bone microstructure without relying on densitometric or model assumptions.

Wavelet-based characterization of vertebral trabecular bone structure from magnetic resonance images at 3 T compared with micro-computed tomographic measurements

Trabecular bone structure and bone density contribute to the strength of bone and are important in the study of osteoporosis. Wavelets are a powerful tool in characterizing and quantifying texture in an image. The purpose of this study was to validate wavelets as a tool in computing trabecular bone thickness directly from gray-level images. To this end, eight cylindrical cores of vertebral trabecular bone were imaged using 3-T magnetic resonance imaging (MRI) and micro-computed tomography (microCT). Thickness measurements of the trabecular bone from the wavelet-based analysis were compared with standard 2D structural parameters analogous to bone histomorphometry (MR images) and direct 3D distance transformation methods (microCT images). Additionally, bone volume fraction was determined using each method. The average difference in trabecular thickness between the wavelet and standard methods was less than the size of 1 pixel size for both MRI and microCT analysis. A correlation (R) of .94 for microCT measurements and that of .52 for MRI were found for the bone volume fraction. Based on these results, we conclude that wavelet-based methods deliver results comparable with those from established MR histomorphometric measurements. Because the wavelet transform is more robust with respect to image noise and operates directly on gray-level images, it could be a powerful tool for computing structural bone parameters from MR images acquired using high resolution and thus limited signal scenarios.

Evaluation of fetal bone structure and mineralization in IGF-I deficient mice using synchrotron radiation microtomography and Fourier transform infrared spectroscopy

The role of insulin like growth factor-I (IGF-I) during pre-natal development has not been evaluated in detail. However, the high degree of growth retardation and peri-natal mortality in IGF-I deficient mouse models indicates that it plays a critical role during this time. Techniques to assess the structure and quality of bone in small animal fetuses could be beneficial in better understanding its role in bone metabolism and skeletal development. Synchrotron microtomography (SR-microCT) and Fourier transform infrared spectroscopy (FTIR) may provide methods to visualize and quantify differences in the structure and mineral density of bone in small animal fetuses. Tibia and spine from IGF-I deficient and wildtype fetal mice (18th day gestation) were imaged using SR-microCT. Three-dimensional structural indices and the degree of mineralization were determined for each sample. Mineralization was also assessed using FTIR and von Kossa staining. Bone volume was systematically lower in IGF-I -/- animals (tibia: -15%, p<0.05) while both sites were found to have a more rod-like architecture (24%, p<0.05; 113%, p<0.01) and lower trabecular separation (-16%, p<0.05; -21%, p<0.05). These structural results were mostly consistent with those seen in adult models of IGF-I deficiency. The degree of mineralization as measured by SR-microCT was higher in the IGF-I tibial metaphysis (11.7%, p<0.0001), while FTIR of the whole bone showed mineralization to be lower in the knockout group (-11%, p<0.05). Interestingly, von Kossa staining revealed no mineral content in the IGF-I -/- spinal ossification center while SR-microCT clearly indicated the presence of highly attenuating components, if somewhat lower in IGF-I -/- animals (-2.2%, p<0.05). This indicates that IGF-I deficiency is linked to subtle differences in the mineral environment and mineralization progression. The advantages unique to SR-microCT allow for 3D visualization and quantification of pre-natal bone microstructure and mineral density in mice which was not previously possible.

Sensitivity of damage predictions to tissue level yield properties and apparent loading conditions

High-resolution finite element models of trabecular bone failure could be used to augment current techniques for measuring damage in trabecular bone. However, the sensitivity of such models to the assumed tissue yield properties and apparent loading conditions is unknown. The goal of this study was to assess the sensitivity of the amount and mode (tension vs. compression) of tissue level yielding in trabecular bone to these factors. Linear elastic, high-resolution finite element models of nine bovine tibial trabecular bone specimens were used to calculate the fraction of the total tissue volume that exceeded each criterion for apparent level loading to the reported elastic limit in both on-axis and transverse compression and tension, and in shear. Four candidate yield criteria were studied, based on values suggested in the literature. Both the amount and the failure mode of yielded tissue were sensitive to the magnitudes of the tissue yield strains, the degree of tension-compression asymmetry of the yield criterion, and the applied apparent loads. The amount of yielded tissue was most sensitive to the orientation of the applied apparent loading, with the most tissue yielding for loading along the principal trabecular orientation and the least for loading perpendicular to it, regardless of the assumed tissue level yield criterion. Small changes in the magnitudes and the degree of asymmetry of the tissue yield criterion resulted in much larger changes in the amount of yielded tissue in the model. The results indicate that damage predictions based on high-resolution finite element models are highly sensitive to the assumed tissue yield properties. As such, good estimates of these values are needed before high-resolution finite element models can be applied to the study of trabecular bone damage. Regardless of the assumed tissue yield properties, the amount and type of damage that occurs in trabecular bone depends on the relative orientations of the applied apparent loads to the trabecular architecture, and this parameter should be controlled for both experimental and computational damage studies.