Peter Augat

High-Resolution Magnetic Resonance Imaging: Three-Dimensional Trabecular Bone Architecture and Biomechanical Properties

The purpose of this study was to use high-resolution magnetic resonance (MR) imaging combined with image analysis to investigate the three-dimensional (3D) trabecular structure, anisotropy, and connectivity of human vertebral, femoral, and calcaneal specimens. The goal was to determine whether: (a) MR-derived measures depict known skeletal-site-specific differences in architecture and orientation of trabeculae; (b) 3D architectural parameters combined with bone mineral density (BMD) improve the prediction of the elastic modulus using a fabric tensor formulation; (c) MR-derived 3D architectural parameters combined with BMD improve the prediction of strength using a multiple regression model, and whether these results corresponded to the results obtained using higher resolution depictions of trabecular architecture. A total of 94 specimens (12 x 12 x 12 mm cubes) consisting of trabecular bone only were obtained, of which there were 7 from the calcaneus, 15 from distal femur, 47 from the proximal femur, and 25 from the vertebral bodies. MR images were obtained using a 1.5 Tesla MR scanner at a spatial resolution of 117 x 117 x 300 microm. Additionally, BMD was determined using quantitative computed tomography (QCT), and the specimens were nondestructively tested and the elastic modulus (YM) was measured along three orthogonal axes corresponding to the anatomic superior-inferior (axial), medial-lateral (sagittal), and anterior-posterior (coronal) directions. A subset of the specimens (n=67) was then destructively tested in the superior-inferior (axial) direction to measure the ultimate compressive strength. The MR images were segmented into bone and marrow phases and then analyzed in 3D. Ellipsoids were fitted to the mean intercept lengths, using single value decomposition and the primary orientation of the trabeculae and used to calculate the anisotropy of trabecular architecture. Stereological measures were derived using a previously developed model and measures such as mean trabecular width, spacing, and number were derived. Because the spatial resolution of MR images is comparable to trabecular bone dimensions, these measures may be subject to partial volume effects and were thus treated as apparent measures, such as BV/TV, Tb.Sp, Tb.N, and Tb.Th rather than absolute measures, as would be derived from histomorphometry. In addition, in a subset of specimens, the Euler number per unit volume was determined to characterize the connectivity of the trabecular network. There were significant differences in the BMD, trabecular architectural measures, elastic modulus, and strength at the different skeletal sites. The primary orientation axes for most of the specimens was the anatomic superior-inferior (axial) direction. Using the fabric tensor formulation, in addition to BMD, improved the prediction of YM (SI), while including some of the architectural parameters significantly improved the prediction of strength. In comparing MR-derived 3D measures with those obtained from 20 microm optical images (n=18; 9 vertebrae, 9 femur specimens), good correlations were found for the apparent Tb.Sp and Tb.N, moderate correlation was seen for the apparent BV/TV, and poor correlation was found for the apparent Tb.Th. Using these higher resolution images, the fabric tensor formulation for predicting the elastic modulus also showed improved correlation between the measured and calculated modulus in the axial (SI) direction. In summary, high-resolution MR images may be used to assess 3D architecture of trabecular bone, and the inclusion of some of the 3D architectural measures provides an improved assessment of biomechanical properties. Further studies are clearly warranted to establish the role of architecture in predicting overall bone quality, and the role of trabecular architecture measures in clinical practice. (ABSTRACT TRUNCATED)

In Vivo High Resolution MRI of the Calcaneus: Differences in Trabecular Structure in Osteoporosis Patients

The purpose of this study was to use high resolution (HR) magnetic resonance (MR) images of the calcaneus to investigate the trabecular structure of patients with and without osteoporotic hip fractures and to compare these techniques with bone mineral density (BMD) in differentiating fracture and nonfracture patients. Axial and sagittal HR MR images of the calcaneus were obtained in 50 females (23 postmenopausal patients with osteoporotic hip fractures and 27 postmenopausal controls). A three‐dimensional gradient–echo sequence was used with a slice thickness of 500 μm and in plane resolution of 195 × 195 μm. Texture analysis was performed using morphological features, analogous to standard histomorphometry and fractal dimension. Additionally, BMD measurements of the hip (dual‐energy X‐ray absorptiometry) were obtained in all patients. Significant differences between both patient groups were obtained using morphological parameters and fractal dimension as well as hip BMD (p < 0.05). Odds ratios for the texture parameters apparent (app.) bone volume/total volume and app. trabecular separation were higher than for hip BMD. Receiver operator characteristic values of texture measures and hip BMD were comparable. In conclusion, trabecular structure measures derived from HR MR images of the calcaneus can differentiate between postmenopausal women with and without osteoporotic hip fractures.

Volumetric quantitative computed tomography of the proximal femur : Precision and relation to bone strength

We have developed a three-dimensional computed tomography (CT) scanning and image analysis method for measurement of trabecular and integral bone mineral density (BMD) and geometry in automatically determined femoral-neck and trochanteric subregions of the proximal femur. We measured the correlation of the density and geometry variables to femoral strength assessed in vitro under loading simulating a single-limb condition and a fall to the side. While BMD alone accounted for 48%-77% of the variability in strength for the stance loading configuration, femoral neck cross-sectional area (minCSA) and femoral neck axis length (FNAL) also contributed independently to femoral strength, and a combination of BMD and geometry variables explained 87%-93% of the variance in the data. For the fall loading configuration, trochanteric trabecular BMD alone explained 87% of the variability of strength. The reproducibility in vivo of the technique was assessed in a group of seven postmenopausal women, who underwent repeat scans with repositioning. For trabecular BMD, the precision was 1.1% and 0.6% for the femoral neck and trochanteric subregions, respectively, compared to 3.3% and 1.6% for the corresponding integral envelopes. Thus, trabecular BMD measurements were reproducible and highly correlated to biomechanical strength measurements. These results support further exploration of quantitative CT for assessment of osteoporosis at the proximal femur.

A Comparative Study of Trabecular Bone Properties in the Spine and Femur Using High Resolution MRI and CT

The purpose of this study was to use high resolution (HR) magnetic resonance (MR) and computed tomography (CT) images combined with texture analysis to investigate the trabecular structure of human vertebral and femoral specimens and to compare these techniques with bone mineral density (BMD) in the prediction of bone strength. Twenty‐nine bone cubes were harvested from 12 proximal femur cadaver specimens and 29 from 8 spines. HR MR and CT images were obtained, and texture analysis techniques were used to assess trabecular structure. Additionally, BMD, elastic modulus (EM), and maximum compressive strength were determined. R2 for EM versus texture measures computed in the MR images was higher (R2 = 0.27–0.64, p < 0.01) in the spine than in the femur specimens (R2 = 0.12–0.22, p < 0.05). R2 values were similar in the CT images. R2 for EM versus BMD was 0.66 (p < 0.01) in the spine and 0.61 (p < 0.01) in the femur specimens. In the MR images, texture measures combined with BMD in a multivariate‐regression model significantly increased R2, while improvement was less significant in the CT images. Thus, texture analysis may provide additional information needed to analyze bone strength and quality.

Ultrasound velocity of trabecular cubes reflects mainly bone density and elasticity.

Abstract. Studies have indicated that quantitative ultrasound (QUS) variables may be influenced by the mechanical properties of bone which in turn are determined by bones material and structural properties. However, from these studies it is unclear what role density, elasticity, and structure play in determining velocity. Eighteen defatted, 12-mm cubic trabecular bone specimens were cut from cadaveric specimens. Amplitude-dependent speed of sound (SOS) using a single point QUS system was assessed in three orthogonal axes. Magnetic resonance images were obtained, from which measures of apparent trabeuclar structure were derived. The specimens were nondestructively tested in compression along three orthogonal axes defined by the sides of the cubes. The elastic modulus (in the three directions) and the strength (in one direction) were determined. Trabecular BMD was measured by quantitative computed tomography. SOS varied significantly with direction of measurement, with the highest value in the axial direction (axial:1715 m/s, sagittal: 1662 m/second, and coronal: 1676 m/s). SOS of each of the three axes was generally associated with the various mechanical (r = 0.30–0.87), density (r = 0.81–0.93), and bone structural variables (0.3–0.8). However, after adjusting the SOS correlations by density, only the correlation with elasticity remained significant in the coronal direction. BMD alone explained 88–93% of variance in SOS whereas in the multivariate model, BMD plus elasticity and/or anisotropic variables explained 96–98% of the variance in SOS. Variability of SOS is explained mostly by density and to a small extent by elasticity or anisotropy. Since only 2–6% of the variance of the QUS measurement is not explained by density and elasticity, one could conclude that the remaining variance reflects other properties of bone or perhaps simply measurement error. Evidence that these other properties may be structure related is only found in the anisotropy of QUS parameter.

Trabecular Bone Architecture in the Distal Radius Using Magnetic Resonance Imaging in Subjects with Fractures of the Proximal Femur

Abstract:To determine whether magnetic resonance (MR)-derived measures of trabecular bone architecture in the distal radius are predictive for prevalent hip fractures, 20 subjects with hip fractures and 19 age-matched postmenopausal controls were studied. Bone mineral density (BMD) measures at the hip (dual-energy X-ray absorptiometry, DXA) and the distal radius (peripheral quantitative computed tomography, pQCT) were also obtained. We compared the MR-based structural measures derived in the radius with those in the calcaneus of the same patients. In the radius, images were acquired at an in-plane resolution of 156 μm and a slice thickness of 0.5 mm. Stereologic measures such as the apparent trabecular thickness (app. Tb.Th), fractional trabecular bone volume (app. BV/TV), trabecular spacing (app. Tb.Sp) and trabecular number (app. Tb.N) were derived from the images. Measures of app. Tb.Sp and app. Tb.N in the distal radius showed significant (p<0.05) differences between the two groups, as did hip BMD measures. However, radial trabecular BMD measures showed only a marginal difference (p= 0.05). Receiver operating curve analysis was used to determine the diagnostic efficacy of BMD, structural measures and a combination of the two. The area under the curve (AUC) for total hip BMD was 0.73, and for radial trabecular BMD was 0.69. AUC for most of the measures of trabecular bone structure at the distal radius was lower than for hip BMD measures; however, AUC for app. Tb.N at the radius was 0.69, comparable to trabecular BMD using pQCT. The AUC for combined BMD (hip) and structure measures was higher (0.87) when radius and calcaneus structure was included. Measures of trabecular architecture derived from MR images combined with BMD measures improve the discrimination between subjects with hip fractures and normal age-matched controls.

Accuracy of cortical and trabecular bone measurements with peripheral quantitative computed tomography (pQCT)

In order to assess the accuracy of peripheral QCT (Stratec XCT 960) we analysed scans of the European Forearm Phantom and another phantom consisting of K2HPO4 encased in aluminium tubes to simulate cortical walls. Additionally 14 cadaveric forearm specimen scans were compared to CT scans acquired on a GE9800Q. The accuracy for density assessment of the European Forearm Phantom was better than 3%. A small increase in density was observed with increasing thickness of the aluminium wall (10% for each mm). Density measurements within the wall were confounded by limited spatial resolution. For a thickness of less than 4 mm, the density within the wall was underestimated by up to 40%. The measurement of mineral content was not influenced by this effect and showed an accuracy error of less than 6%. The agreement of density measurements on the different CT systems was very strong (R2 > 0.96; RMSE < 6.2%). Our findings suggest that the Stratec pQCT scanner very accurately measures volumetric trabecular and total bone mineral densities at the distal radius while the assessment of cortical density is associated with considerable inaccuracies due to limited spatial resolution.

Quantitative Bone Mineral Assessment at the Forearm: A Review

Abstract: Bone mineral density and geometric properties of the human forearm can be measured to determine the amount of bone or bone loss at the scanning site and to predict the risk of forearm fractures. These forearm mesurements are also used to estimate bone mass at remote anatomical locations and thereby estimate the risk for spine, hip and other fractures. The peripheral location of the human forearm, with its relatively small amount of surrounding soft tissue, improves the accuracy and the precision of bone mass measurement and has made this site an early choice for the assessment of a subjects bone mineral status. Furthermore, the anatomy of the human radius enables the examination of both cortical and cancellous bone. This review describes the procedures for non-invasive bone assessment at peripheral sites including some of the more recently developed systems dedicated to assessment of the distal radius. The accuracy, precision and normative values they provide are presented. Responses to different forms of therapies as well as the ability to discriminate or predict osteoporotic fractures are also assessed. Low radiation dose, comfortable and fast handling, moderate cost, and a strong association with the risk of non-spine fractures, promote the use of forearm scanning as a widely applied screening procedure for the detection of generalised osteoporotic bone loss. However, a higher accuracy of fracture risk prediction at the spine or at the hip can be achieved by a direct bone density measurement at these sites. The monitoring of treatment at the distal forearm appears to require a longer follow-up time due to its decreased responsiveness compared with such highly trabecular load-bearing sites as the spine and the proximal femur.

Fractal analysis of radiographs: Assessment of trabecular bone structure and prediction of elastic modulus and strength

The purpose of this study was to determine whether fractal dimension of radiographs provide measures of trabecular bone structure which correlate with bone mineral density (BMD) and bone biomechanics, and whether these relationships depend on the technique used to calculate the fractal dimension. Eighty seven cubic specimen of human trabecular bone were obtained from the vertebrae and femur. The cubes were radiographed along all three orientations--superior-inferior (SI), medial-lateral (ML), and anterior-posterior (AP), digitized, corrected for background variations, and fractal based techniques were applied to quantify trabecular structure. Three different techniques namely, semivariance, surface area, and power spectral methods were used. The specimens were tested in compression along three orientations and the Youngs modulus (YM) was determined. Compressive strength was measured along the SI direction. Quantitative computed tomography was used to measure trabecular BMD. High-resolution magnetic-resonance images were used to obtain three-dimensional measures of trabecular architecture such as the apparent bone volume fraction, trabecular thickness, spacing, and number. The measures of trabecular structure computed in the different directions showed significant differences (p<0.05). The correlation between BMD, YM, strength, and the fractal dimension were direction and technique dependent. The trends of variation of the fractal dimension with BMD and biomechanical properties also depended on the technique and the range of resolutions over which the data was analyzed. The fractal dimension showed varying trends with bone mineral density changes, and these trends also depended on the range of frequencies over which the fractal dimension was measured. For example, using the power spectral method the fractal dimension increased with BMD when computed over a lower range of spatial frequencies and decreased for higher ranges. However, for the surface area technique the fractal dimension increased with increasing BMD. Fractal measures showed better correlation with trabecular spacing and number, compared to trabecular thickness. In a multivariate regression model inclusion of some of the fractal measures in addition to BMD improved the prediction of strength and elastic modulus. Thus, fractal based texture analysis of radiographs are technique dependent, but may be used to quantify trabecular structure and have a potentially valuable impact in the study of osteoporosis.

Trabecular Bone Mineral and Calculated Structure of Human Bone Specimens Scanned by Peripheral Quantitative Computed Tomography: Relation to Biomechanical Properties

The relationship of cortical bone mineral density (BMD), and geometry to bone strength has been well documented. In this study, we used peripheral quantitative computerized tomography (pQCT) to acquire trabecular BMD and high‐resolution images of trabeculae from specimens to determine their relationship with biomechanical properties. Fifty‐eight human cubic trabecular bone specimens, including 26 from the vertebral bodies, were scanned in water and air. Trabecular structure was quantitated using software developed with Advanced Visual Systems interfaced on a Sun/Sparc Workstation. BMD was also obtained using a whole‐body computerized tomography scanner (QCT). Nondestructive testing of the specimens was performed to assess their elastic modulus. QCT and pQCT measurements of BMD of specimens in water were strongly correlated (r2 = 0.95, p < 0.0001), with a slope (0.96) statistically not significantly different from 1. Strong correlations were found between pQCT measurements of specimens in water and in air, for BMD (r2 = 0.96, p < 0.0001), and for apparent trabecular structural parameters (r2 = 0.89–0.93, p < 0.0001). Correlations were moderate between BMD and apparent trabecular structural parameters (r2 = 0.37–0.64, p < 0.0001). Precision as coefficient of variation (CV) and standardized coefficient of variation (SCV) for these measurements was < 5%. For the vertebral specimens, the correlation was higher between elastic modulus and BMD (r2 = 0.76, p < 0.0001) than between elastic modulus and apparent trabecular structural parameters (r2 = 0.58–0.72, p < 0.0001), while the addition of apparent trabecular nodes and branches to BMD in a multivariate regression model significantly increased the correlation with the elastic modulus (r2 = 0.86, p < 0.01). Thus, pQCT can comparably and reproducibly measure trabecular bone mineral in water or air, and trabecular structure can be quantitated from pQCT images. The combination of volumetric BMD with trabecular structural parameters rather than either alone improves the prediction of biomechanical properties. Such a noninvasive approach may be useful for the preclinical study of osteoporosis.