Heterogeneity of trabecular bone structure in the calcaneus using magnetic resonance imaging

Treatment of Severe Avascular Necrosis of the Talus Using a Novel Keystone-Shaped 3D-Printed Titanium Truss Implant

BACKGROUND

Avascular necrosis (AVN) of the talus most commonly occurs secondary to trauma. Significant bone loss and collapse in severe talar AVN remains an operative challenge. Tibiotalocalcaneal arthrodesis (TTC) using femoral head allograft is at risk of collapse and subsidence. The use of a void-filling titanium truss can mitigate against this. This study describes the use of a novel keystone shaped 3D-printed titanium truss for treatment of severe talar AVN.

METHODS

Three patients with end-stage AVN of the talus were included. Each patient underwent a TTC arthrodesis with a custom-made, 3D-printed, keystone-shaped, truss implant in conjunction with a hindfoot intramedullary nail. Modified patient American Orthopaedic Foot & Ankle Society (AOFAS) scores were recorded at the preoperative, 6-month, 12-month, and annual postoperative timepoints.

RESULTS

All patients progressed to satisfactory radiological union by one year. Mean follow up time was 32 months (24-48 months). Mean preoperative modified AOFAS score was 5. There was progressive improvement in AOFAS scores from 6 months postoperatively. Mean modified AOFAS score improved from 28 at 6 months to 37 at 2 years postoperatively.

CONCLUSION

Custom-made 3D-printed titanium trusses provide promising outcomes for treating severe AVN of the talus. The "keystone" design is advantageous as it allows for bone stock preservation and conforms to the shape of the native calcaneum. All patients showed progressive improvements in outcomes at sequential time intervals postoperatively. The implant provides a strong mechanical structure resisting collapse and subsidence during the arthrodesis process.

LEVEL OF EVIDENCE

Level IV, retrospective case series.

Trabecular bone in the calcaneus of runners

Trabecular bone of the human calcaneus is subjected to extreme repetitive forces during endurance running and should adapt in response to this strain. To assess possible bone functional adaptation in the posterior region of the calcaneus, we recruited forefoot-striking runners (n = 6), rearfoot-striking runners (n = 6), and non-runners (n = 6), all males aged 20-41 for this institutionally approved study. Foot strike pattern was confirmed for each runner using a motion capture system. We obtained high resolution peripheral computed tomography scans of the posterior calcaneus for both runners and non-runners. No statistically significant differences were found between runners and nonrunners or forefoot strikers and rearfoot strikers. Mean trabecular thickness and mineral density were greatest in forefoot runners with strong effect sizes (<0.80). Trabecular thickness was positively correlated with weekly running distance (r2 = 0.417, p<0.05) and years running (r2 = 0.339, p<0.05) and negatively correlated with age at onset of running (r2 = 0.515, p<0.01) Trabecular thickness, mineral density and bone volume ratio of nonrunners were highly correlated with body mass (r2 = 0.824, p<0.05) and nonrunners were significantly heavier than runners (p<0.05). Adjusting for body mass revealed significantly thicker trabeculae in the posterior calcaneus of forefoot strikers, likely an artifact of greater running volume and earlier onset of running in this subgroup; thus, individuals with the greatest summative loading stimulus had, after body mass adjustment, the thickest trabeculae. Further study with larger sample sizes is necessary to elucidate the role of footstrike on calcaneal trabecular structure. To our knowledge, intraspecific body mass correlations with measures of trabecular robusticity have not been reported elsewhere. We hypothesize that early adoption of running and years of sustained moderate volume running stimulate bone modeling in trabeculae of the posterior calcaneus.

Imaging methods for bone mass evaluation during childhood and adolescence: an update

The objective of the work was to prepare an update on imaging methods for bone evaluation during childhood and adolescence. The text was based on original and review articles on imaging methods for clinical evaluation of bone mass in children and adolescents up to 20 years old. They were selected from BIREME and PUBMED by means of the following keywords: bone density; osteoporosis/diagnosis; densitometry; tomography; ultrasonography; magnetic resonance imaging; and radiogrammetry and published in Portuguese or English, in the last 10 years (2006-2016). The article was organized into topics with the description of peculiarities, advantages and disadvantages of each imaging method and their possible clinical applicability. Despite the emergence of new technologies, dual energy X-ray absorptiometry (DXA) remains the gold standard method for low bone mass diagnosis in all age groups. However, interpretation is complex in children and adolescents and demands skilled people. Quantitative computed tomography (QCT) [central QCT, peripheral QCT (pQCT) and high resolution-pQCT (HR-pQCT)] and magnetic resonance imaging (MRI) evaluate real bone density, but are not yet available for routine use. Quantitative bone ultrasound (QUS) shows good perspectives for its use in primary prevention actions. Automated radiogrammetry shows promise as a method able to flag individuals who might benefit from a complete bone assessment, but the clinical value of the measures still needs to be established.

Unusually High Calcaneal Speed of Sound Measurements in Children with Small Foot Size

The purpose of this clinical note is to describe the performance of the Lunar Achilles Insight device in assessing bone quality at the calcaneus in 142 children between the ages of 5 and 11 y accessing healthcare in Johannesburg, South Africa. We observed an asymmetric bimodal distribution in speed of sound (SOS). The minor mode consisted of unusually high SOS values (≥1625 m/s), which were primarily observed among children with foot size <19 cm and height <119 cm. Cortical regions of the bone may have been inadvertently included in the region of interest for smaller feet, causing unusually high SOS values. The unusually high SOS values indicate that the validity of SOS in this device, as it is currently used for measuring bone quality in young children, is questionable. Future studies using this device in young children should develop new methodology to account for smaller foot size.

Potential diagnostic role of the MRI-derived internal magnetic field gradient in calcaneus cancellous bone for evaluating postmenopausal osteoporosis at 3T

INTRODUCTION

Bone mineral density (BMD) result has a low predictive value on patients' risk for future fractures. Thus, new approaches for examining patients at risk for developing osteoporosis would be desirable. Magnetic resonance (MR) investigations in cancellous bone have been shown to yield useful quantitative information on both trabecular-bone microstructure and bone marrow composition. This work was undertaken to address the hypothesis that the effective internal magnetic field gradient (IMFG), a new MR parameter, discriminates between healthy, osteopenic and osteoporotic postmenopausal women, classified on the basis of bone mineral density (BMD) criteria. The work builds on preliminary results indicating that IMFG, measured in trabecular-bone pores and quantified by spin-echo decay and water diffusion MR near the bone-bone marrow interface depends on both the bone marrow water rate of diffusion and the magnetic susceptibility difference (ΔX) between water and bone.

MATERIALS AND METHODS

MR relaxometry, MR spectroscopy and diffusion-weighted MR imaging of the heel was performed in fifty-five women (mean age, 62.9±6.6years) at 3T. Moreover, in order to study the reproducibility of IMFG measurement, five young women (mean age 31.0±3.2years; age range, 28-36years) were scanned and rescanned. The study protocol was approved by the local Ethics Committee. Quantitative Computer Tomography (QCT) of the L1-L3 vertebral segments was performed to classify the postmenopausal women into three groups according to QCT BMD: healthy (n=8); osteopenic (n=25); and osteoporotic (n=22). In all subjects, BMD T-scores, marrow fat content (Mfc), T2*, apparent diffusion coefficient (ADC) and IMFG (estimated from the additional spin-echo decay due to diffusion of water in local magnetic field gradients), were assessed in the whole calcaneus as well as in three calcaneal subregions: subtalar, tuber calcaneus, and cavum calcaneus. Between-group comparisons to assess group differences and Pearson correlation analysis were performed. Short and long-term coefficients of variation (CVS and CVL, respectively) were evaluated in young subjects.

RESULTS

Reproducibility of the IMFG measurement was satisfactory. No significant difference was found in the IMFG measurement performed in both calcaneus and subtalar calcaneal region between the two separate sessions comprised of five young women. Mfc did not significantly differ between groups. The IMFG in the subtalar region was significantly different between all three groups (P<0.01), being greatest in healthy women, intermediate in those with osteopenia, and lowest in osteoporotic subjects. Conversely neither T2* nor ADC is able to discriminate healthy subjects from those with osteopenia and osteoporosis. Increased inter-trabecular space, as it typically occurs in patients with osteoporosis, modifies water diffusion, conferring higher ADC values, thereby lowering the IMFG.

CONCLUSION

The IMFG measured in the calcaneal subtalar region shows a high ability in identifying healthy subjects. The new quantitative MR method based on measurement of the IMFG may provide a new means for assessing patients with osteoporosis.

Cylinders or walls? A new computational model to estimate the MR transverse relaxation rate dependence on trabecular bone architecture

OBJECTIVE

Bone density is distributed in a complex network of interconnecting trabecular plates and rods that are interspersed with bone marrow. A computational model to assess the dependence of the relaxation rate on the geometry of bone can consider the distribution of bone material in the form of two components: cylinders and open walls (walls with gaps). We investigate whether the experimentally known dependence of the transverse relaxation rate on the trabecular bone structure can be usefully interpreted in terms of these two components.

MATERIALS AND METHODS

We established a computer model based on an elementary computational cell. The model includes a variable number of open walls and infinitely long cylinders as well as multiple geometric parameters. The transverse relaxation rate is computed as a function of these parameters. Within the model, increasing the trabecular spacing with a fixed trabecular radius is equivalent to thinning the trabeculae while maintaining constant spacing.

RESULTS

Increasing the number of cylinder and wall gap elements beyond their nearest neighbors does not change the transverse relaxation rate. Although the absolute contribution to the relaxation due to open walls is on average more important than that due to cylinders, the latter drops off rapidly. The change on transverse relaxation rate is larger for changing cylinder geometry than for changing wall geometry, as it can be seen from the effect on the relaxation rate when trabecular spacing is varied, compared to varying the size of wall gaps.

CONCLUSION

Our results provide strong evidence that trabecular thinning, which is associated with increasing age, decreases the relaxation rates. The effect of thinning plates and rods on the transverse relaxation can be understood in terms of simple cylinders and open walls. A reduction in the relaxation rate can be seen as an indication of thinning cylinders, corresponding to reduced bone stability and ultimately, osteoporosis.

Ultrasonic evaluation of bone quality in cadaver ilia

Several imaging modalities have traditionally been utilized to assess bone health. However, none of these standards is capable of providing a clear rendition or display of the damaged bone layers caused, for instance, by osteoporosis. This study examines the use of ultrasound for non-invasive monitoring of bone quality in bone samples with various degrees of porosity. A user-defined region of interest (ROI) in the iliac portion of extracted human cadaver coxal bones is monitored and quantified. Raster C-scan images of the ROI were acquired and compared to basic physical measurements, and to bone scans using dual energy X-ray absorptiometry (DXA). A quantitative measure of the superficial sub-surface composite matrix (ScM) content was analyzed using linear regression with all physical and DXA measures. The trend in the degree of percent bone loss (PBL) measured by ultrasound (US) was found to be closely paralleled with that measured by DXA (R(2) = 0.82, p < .0005). Also, the trend in which PBL (US) correlated with bone mineral density (BMD) (R(2) = 0.62, p < .01) was found to exhibit a similar behavior when the latter was compared to dry mass density (DmD) of the bone samples (R(2) = 0.63, p < .01). However, when PBL (DXA) was compared to DmD, it did reveal a better linearity (R(2) = 0.69, p < .005) than the one obtained when PBL (US) was compared with the same DmD (R(2) = 0.45, p < .05). A similar outcome was observed when PBL (US) was compared with percent porosity (R(2) = 0.51, p < .05), as opposed to the better linearity exhibited between PBL (DXA) and porosity (R(2) = 0.86, p < .0005). Despite these slight variations, further analyses on the statistical significance between these correlations suggest that ultrasound can be an effective imaging technique in assessing the degree of bone damage, and can be used to assess the structural integrity of bones.

Diagnosis of osteoarthritis and prognosis of tibial cartilage loss by quantification of tibia trabecular bone from MRI

A longitudinal study was used to investigate the quantification of osteoarthritis and prediction of tibial cartilage loss by analysis of the tibia trabecular bone from magnetic resonance images of knees. The Kellgren Lawrence (KL) grades were determined by radiologists and the levels of cartilage loss were assessed by a segmentation process. Aiming to quantify and potentially capture the structure of the trabecular bone anatomy, a machine learning approach used a set of texture features for training a classifier to recognize the trabecular bone of a knee with radiographic osteoarthritis. Using cross-validation, the bone structure marker was used to estimate for each knee both the probability of having radiographic osteoarthritis (KL >1) and the probability of rapid cartilage volume loss. The diagnostic ability reached a median area under the receiver-operator-characteristics curve of 0.92 (P < 0.0001), and the prognosis had odds ratio of 3.9 (95% confidence interval: 2.4-6.5). The medians of cartilage loss of the subjects classified as slow and rapid progressors were 1.1% and 4.9% per year, respectively. A preliminary radiological reading of the high and low risk knees put forward an hypothesis of which pathologies the bone marker could be capturing to define the prognosis of cartilage loss.

Accuracy of high-resolution in vivo micro magnetic resonance imaging for measurements of microstructural and mechanical properties of human distal tibial bone

Micro magnetic resonance imaging (µMRI) is an in vivo imaging method that permits 3D quantification of cortical and trabecular bone microstructure. µMR images can also be used for building microstructural finite element (µFE) models to assess bone stiffness, which highly correlates with bone's resistance to fractures. In order for µMRI-based microstructural and µFE analyses to become standard clinical tools for assessing bone quality, validation with a current gold standard, namely, high-resolution micro computed tomography (µCT), is required. Microstructural measurements of 25 human cadaveric distal tibias were performed for the registered µMR and µCT images, respectively. Next, whole bone stiffness, trabecular bone stiffness, and elastic moduli of cubic subvolumes of trabecular bone in both µMR and µCT images were determined by voxel-based µFE analysis. The bone volume fraction (BV/TV), trabecular number (Tb.N*), trabecular spacing (Tb.Sp*), cortical thickness (Ct.Th), and structure model index (SMI) based on µMRI showed strong correlations with µCT measurements (r(2) = 0.67 to 0.97), and bone surface-to-volume ratio (BS/BV), connectivity density (Conn.D), and degree of anisotropy (DA) had significant but moderate correlations (r(2) = 0.33 to 0.51). Each of these measurements also contributed to one or many of the µFE-predicted mechanical properties. However, model-independent trabecular thickness (Tb.Th*) based on µMRI had no correlation with the µCT measurement and did not contribute to any mechanical measurement. Furthermore, the whole bone and trabecular bone stiffness based on µMRI were highly correlated with those of µCT images (r(2) = 0.86 and 0.96), suggesting that µMRI-based µFE analyses can directly and accurately quantify whole bone mechanical competence. In contrast, the elastic moduli of the µMRI trabecular bone subvolume had significant but only moderate correlations with their gold standards (r(2) = 0.40 to 0.58). We conclude that most microstructural and mechanical properties of the distal tibia can be derived efficiently from µMR images and can provide additional information regarding bone quality.

Radiographic texture analysis of densitometric calcaneal images: relationship to clinical characteristics and to bone fragility

Osteoporotic fractures are related not only to bone mineral density (BMD) but also to bone structure or microarchitecture, which is not assessed routinely with currently available methods. We have developed radiographic texture analysis (RTA) for calcaneal images from a peripheral densitometer as an easy, noninvasive method for assessing bone structure. We conducted a cross-sectional study of the relationship between RTA and prevalent vertebral fractures (n = 148) among 900 subjects (ages 19 to 99 years, 94 males) referred for bone densitometry as part of their routine medical care. RTA features were derived from Fourier-based image analysis of the radiographic texture pattern (including root mean square, first moment, and power spectral analyses). RTA features were associated with age, weight, gender, and race, as well as glucocorticoid use. When controlling for clinical risk factors and BMD (or a summary measure calculated using FRAX algorithms), RTA features were significantly different for subjects with and without prevalent vertebral fractures [adjusted odds ratio (OR) = 1.5 per 1 standard deviation (SD) decrease in RTA feature beta, 95% confidence interval (CI) 1.2-1.8, p = .001]. Gender and use of pharmacologic therapy for osteoporosis did not significantly affect this association, suggesting that RTA can be applied to a wide range of densitometry patients. We conclude that RTA obtained using a portable instrument has a potential as a noninvasive method to enhance identification of patients at increased risk of osteoporotic fractures.

Reproducibility and sources of variability in radiographic texture analysis of densitometric calcaneal images

Radiographic texture analysis (RTA) is a computerized analysis of the spatial pattern of radiographic images used as a way of evaluating bone structure. We have shown that RTA performed on high-resolution heel images obtained using a portable densitometer differentiates subjects with and without osteoporotic fractures. In the present study, short-term precision of RTA was examined on densitometric heel images obtained from 33 subjects scanned 8 times each, with 3 observers placing a region of interest (ROI) 3 times on each image. The long-term precision was examined on images obtained from 10 subjects 3 times on each of 3 days separated by 1 week, with 2 observers placing an ROI on each image. The RTA features examined included the root mean square (RMS) variation, a measure of the contrast between the light and dark areas of the image, the first moment of the power spectrum, a measure of the spatial frequency of the trabecular pattern, and Minkowski fractal (MINK), a measure of roughness/smoothness of the trabecular pattern. The precision of the RTA features expressed as coefficient of variation ranged between the lowest of 0.5-0.7% for MINK and the highest of 14-16% for RMS. The short- and long-term precision was similar, and was not significantly influenced by repositioning and rescanning, or by ROI placement by the same or different observers. Significant sources of variability of RTA were the between-subject differences and differences between regions of the heel, but not differences due to repositioning, rescanning in the same position, or ROI placement by the same or different observers. We conclude that technical aspects of image acquisition and processing are adequate to allow further development of RTA of the densitometric images for clinical application as a method for noninvasive assessment of bone structure.

Quantitative ultrasound methods to assess bone mineral status in children: technical characteristics, performance, and clinical application

Measurement of bone mineral status may be a useful tool in identifying the children who could be exposed to an increased risk of osteoporosis in adulthood. Dual energy x-ray absorptiometry and peripheral quantitative computed tomography may be used to this purpose, but the exposure to ionizing radiation is a limiting factor for preventive studies in large populations of children. In the last years, quantitative ultrasound (QUS) methods have been developed to assess bone mineral status in some peripheral skeletal sites such as calcaneus, phalanges of the hand, and tibia. QUS techniques are safe, easy to use, radiation-free, and devices are portable, so that they are particularly indicated to assess bone mineral status in children. This review will concentrate on the main methodological principles of ultrasounds and the QUS variables derived from their application to bone tissue, technical differences and performance of QUS methods, factors influencing QUS measurements, normative data and results obtained in children with disturbances of growth or affected by disorders of bone and mineral metabolism, including the assessment of fracture risk, and comparison among QUS, dual energy x-ray absorptiometry, and peripheral quantitative computed tomography methods.

Direct measurement of trabecular bone anisotropy using directional fractal dimension and principal axes of inertia

OBJECTIVE

Precise in vivo measurement of the trabecular bone's mechanical properties is very important for endosseous dental implant treatment and design in clinical practice. The fractal structure of trabecular bone shows directional anisotropy of the architecture, as is shown in most biological fractals. To analyze the anisotropy of the trabecular bone, the fractal geometry technique was applied to 2-dimensional plain radiographs.

STUDY DESIGN

The power spectrum was used to calculate the fractal dimensions (FD) of the trabecular bone. The FDs calculated as a function of orientation yielded the fractal information reflecting the spatial characteristics of the trabecular bone in each direction. A polar plot of directional FDs was defined as an ellipse of inertia. The principal loading direction in a local region of the trabecular bone was determined from the minimum moment of inertia for the ellipse of FDs. The anisotropy was calculated directly as the ratio of the 2 principal moments of inertia from the ellipse.

RESULTS

The anisotropies were measured for radiographs from the angle and incisor region of 21 human mandibles based on the principal axes of inertia and the best-fitting ellipse. The anisotropy of the angle region was significantly greater than that of the incisor region of the mandibles.

CONCLUSION

The method using directional FDs as determined by the principal axis of inertia measures the anisotropy directly, using 2-dimensional plain radiographs. It can quantify the anisotropy of trabecular bone in vivo. The investigation can be applied to the analysis of the relationships between in vivo 2-dimensional parameters and 3-dimensional mechanical properties, which enables us to predict the bone mechanical properties such as strength in vivo in various regions of the mandible.

Comparison of trabecular bone anisotropies based on fractal dimensions and mean intercept length determined by principal axes of inertia

The mechanical quality of trabecular bone depends on both its stiffness and its strength characteristics, which can be predicted indirectly by the combination of bone volume fraction and architectural anisotropy. To analyze the directional anisotropy of the trabecular bone, we applied the fractal geometry technique to plain radiographs. The anisotropy of the bone was quantified from an ellipse, based on the directional fractal dimensions (FD), by the principal axes of inertia. The anisotropies based on the FD were compared with those determined using the common method of mean intercept length (MIL). The directional FD gave the fractal information obtained from a projection along the MIL orientation. For this reason, the spatial variations associated with the bone length in any direction were manifested in a related frequency band of the power spectrum determined along the direction. The directional FD and MIL plots were highly correlated, although they originated from quite different geometries. Of the angle, premolar, and incisor regions of the human mandible, the anisotropies calculated using both FD and MIL showed the highest correlation in the trabecular bone of the angle region. The method using directional FDs as determined by the principal axis of inertia measures the anisotropy directly, using two-dimensional plain radiographs. This kind of method will be a useful to provide better estimates of bone quality in vivo compared with the density measurements alone, especially for the indirect diagnosis of jawbone quality in dental clinics.

[The calcaneus as the site of manifestation for osteoporosis-associated fractures: age- and sex-specific changes in calcaneal morphology correlate with the incidence and severity of intra-articular calcaneal fractures]

BACKGROUND

While it is recognized that trauma energy at the time of injury is an important factor in the pathogenesis and severity of calcaneal fractures, the possible role of changes in calcaneal microarchitecture remains largely undefined. The purpose of this study was to determine whether the calcaneal bone structure changes with age and to address if local bone mass is of clinical relevance in respect to the occurrence and complexity of calcaneal fractures.

MATERIAL AND METHODS

The radiographic and clinical data of 182 patients with intra-articular calcaneal fractures were analyzed to provide correlative clinical evidence for a relation between local bone mass and fractures of the calcaneus. To measure bone mass, 60 calcanei were harvested from 30 age- and gender-matched patients at autopsy.

RESULTS

The average age at the time of fracture was higher in females (46.0+/-18.3 years) than in males (39.9+/-13.9 years). Furthermore, the relative frequency of fractures during aging shifted from males to females and the frequency of compound fractures was higher in females (65%) than in males (48%). The calcaneal bone mass was significantly reduced by 19% in older females (female symbol 20-40 years: 292 mg/cm(3); female symbol 61-80 years: 237 mg/cm(3); p<0.05).

CONCLUSION

The calcaneus displayed age- and gender-related changes in its microarchitecture that are known to reduce the biomechanical stability of trabecular bone. These results suggest that bone mass and structure are risk factors in respect to the occurrence and severity of calcaneal fractures.

The 3D-based scaling index algorithm: a new structure measure to analyze trabecular bone architecture in high-resolution MR images in vivo

INTRODUCTION

The purpose of this study was to obtain different structure measures as the three-dimensional (3D)-based scaling index method (SIM) and standard two-dimensional (2D) bone histomorphometric parameters from high-resolution (HR) magnetic resonance (MR) images of the distal radius and to compare these parameters with bone mineral density (BMD) in their diagnostic performance to differentiate postmenopausal patients with and without vertebral fractures.

METHODS

Axial HR-MR images of the distal radius were obtained at 1.5 T in 40 postmenopausal women (17 with osteoporotic spine fractures and 23 controls). Trabecular microarchitecture analysis was performed using the new structure measure mP(alpha), derived from the SIM, as well as standard morphological 2D parameters. BMD of the spine was obtained using quantitative computed tomography (QCT). Receiver operating characteristic (ROC) analyses were used to determine diagnostic performance in differentiating both groups. Results were validated by bootstrapping techniques.

RESULTS

Significant differences between both patient groups were obtained using mP(alpha), 2D parameters, and spine BMD (p<0.05). In comparison with the 2D texture parameters [area under the curve (AUC) up to 0.67], diagnostic performance was significantly higher for mP(alpha)(AUC=0.85; p<0.05). There was a trend for a higher AUC value for mP(alpha) compared with BMD of the spine (AUC=0.71; p=0.81).

CONCLUSION

mP(alpha) yielded a robust measure of trabecular bone microarchitecture for HR-MR images of the radius, which significantly improved the diagnostic performance in differentiating postmenopausal women with and without osteoporotic spine fractures compared with standard 2D bone histomorphometric parameters. This 3D characterization of trabecular microarchitecture may provide a new approach to better assess the strength of human cancellous bone using HR-MR image data.

Radiographic texture analysis of densitometer-generated calcaneus images differentiates postmenopausal women with and without fractures

INTRODUCTION

Bone fragility is determined by bone mass, measured as bone mineral density (BMD), and by trabecular structure, which cannot be easily measured using currently available noninvasive methods. In previous studies, radiographic texture analysis (RTA) performed on the radiographic images of the spine, proximal femur, and os calcis differentiated subjects with and without osteoporotic fractures. The present cross-sectional study was undertaken to determine whether such differentiation could also be made using high-resolution os calcis images obtained on a peripheral densitometer.

METHODS

In 170 postmenopausal women (42 with and 128 without prevalent vertebral fractures) who had no secondary causes of osteoporosis and were not receiving treatment for osteoporosis, BMD of the lumbar spine, proximal femur, and os calcis was measured using dual energy x-ray absorptiometry. Vertebral fractures were diagnosed on densitometric spine images. RTA, including Fourier-based and fractal analyses, was performed on densitometric images of os calcis.

RESULTS

BMD at all three sites and all texture features was significantly different in subjects with and without fractures, with the most significant differences observed for the femoral neck and total hip measurements and for the RTA feature Minkowski fractal (p<0.001). In univariate logistic regression analysis, Minkowski fractal predicted the presence of vertebral fractures as well as femoral neck BMD (p<0.001). In multivariate logistic regression analysis, both femoral neck BMD and Minkowski fractal yielded significant predictive effects (p=0.001), and when age was added to the model, the effect of RTA remained significant (p=0.002), suggesting that RTA reflects an aspect of bone fragility that is not captured by age or BMD. Finally, when RTA was compared in 42 fracture patients and 42 nonfracture patients matched for age and BMD, the RTA features were significantly different between the groups (p=0.003 to p=0.04), although BMD and age were not.

CONCLUSION

This study suggests that RTA of densitometer-generated calcaneus images provides an estimate of bone fragility independent of and complementary to BMD measurement and age.

Bone microarchitecture of the calcaneus and its changes in aging: a histomorphometric analysis of 60 human specimens

Bone structure and quality are an important parameter in the propensity of bone to fracture. Although the calcaneus is used as diagnostic reference site for osteoporosis by ultrasound, its structure has never been analyzed in detail. The purpose of this study was therefore to histomorphometrically analyze the trabecular microarchitecture of the calcaneus, and to determine whether the calcaneal bone structure is changing with age. Sixty complete human calcanei were harvested from thirty age- and gender-matched patients at autopsy. Each of the three different age groups (group I: 20 to 40, group II: 41 to 60, group III: 61 to 80 years of age) was represented by 20 specimens. The specimens were subjected to radiographic, microCT, and histologic analysis. Bone structure and bone mass of the calcaneus were quantified for three different regions of interest: the anterior ROI, the superior ROI (the subtalar region under the posterior facet), and the posterior ROI. An iliac crest biopsy was obtained from all patients to exclude any metabolic bone disease. Histomorphometric analysis revealed significant differences in bone volume within the calcaneus with highest values in the superior ROI: age group I: 31.3% (27.8-34.8%); II: 25.5% (22.1-28.9%); III: 18.9% (16.6-21.2%) and lowest bone volumes in the anterior ROI; age group I: 6.2% (4.8-7.6%); II: 3.6% (2.1-5.1%); III: 3.9% (2.9-4.9%). There was a significant age-related decrease in bone volume (BV/TV) in aging. Interestingly, this bone loss was most prominent in the superior ROI, with a 39% decrease in BV/TV between age group I and III (p < 0.001). Qualitative and structural analysis of trabecular number, thickness, and spacing demonstrated that the bone loss in the thalamic portion of the calcaneus was due to the transition of plate-like trabecular elements into a rod-like structure. In conclusion, our study demonstrated that the calcaneus displayed age-related changes in its microarchitecture that are known to reduce the biomechanical stability of trabecular bone, and that the age-related bone loss was most prominent in the region adjacent to the posterior facet (superior ROI). These results suggest that bone mass and structure are risk factors in respect to the occurrence and severity of calcaneal fractures, and indicate that calcaneal fractures are at least in part osteoporotic fractures.

Trabecular bone volume fraction measurements of a large number of subjects using a compact MRI

Trabecular bone volume fraction (TBVF) and speed of sound (SOS) were measured for the right calcanei of 416 female volunteers. The TBVF was measured with a compact MRI developed in our laboratory. The SOS was measured with a commercial quantitative ultrasound system. It was observed that the correlation coefficient between TBVF and SOS and that between TBVF and age varied depending on the location of region of interest (ROI) in the calcaneus. As a result, an optimum circular ROI with a diameter of 20 mm was determined so that the correlation coefficients were maximized. In the optimum ROI, transverse relaxation time (T2) of the bone marrow protons of the calcaneus was found to be concentrated in a narrow range over the subjects. This result suggested that a 50% scan time reduction in the TBVF measurements could be made by skipping the T2 correction procedure.

Fractal analysis of bone architecture at distal radius

Bone strength depends on bone quality (architecture, turnover, damage accumulation, and mineralization) as well as bone mass. In this study, human bone architecture was analyzed using fractal image analysis, and the clinical relevance of this method was evaluated. The subjects were 12 healthy female controls and 16 female patients suspected of having osteoporosis (age range, 22-70 years; mean age, 49.1 years). High-resolution CT images of the distal radius were acquired and analyzed using a peripheral quantitative computed tomography (pQCT) system. On the same day, bone mineral densities of the lumbar spine (L-BMD), proximal femur (F-BMD), and distal radius (R-BMD) were measured by dual-energy X-ray absorptiometry (DXA). We examined the correlation between the fractal dimension and six bone mass indices. Subjects diagnosed with osteopenia or osteoporosis were divided into two groups (with and without vertebral fracture), and we compared measured values between these two groups. The fractal dimension correlated most closely with L-BMD (r=0.744). The coefficient of correlation between the fractal dimension and L-BMD was very similar to the coefficient of correlation between L-BMD and F-BMD (r=0.783) and the coefficient of correlation between L-BMD and R-BMD (r=0.742). The fractal dimension was the only measured value that differed significantly between both the osteopenic and the osteoporotic subjects with and without vertebral fracture. The present results suggest that the fractal dimension of the distal radius can be reliably used as a bone strength index that reflects bone architecture as well as bone mass.

Multi-detector row CT imaging of vertebral microstructure for evaluation of fracture risk

We applied MDCT for in vivo evaluation of the microarchitecture of human vertebrae. Microstructure parameters, such as structure model index, Euler's number, and bone volume fraction, revealed higher relative risk for prevalent vertebral fracture than did BMD obtained by DXA. Thus, microstructure analysis by MDCT, together with simultaneously obtained volumetric BMD values, is useful for clinical assessment of fracture risk.

INTRODUCTION

BMD measurement by DXA alone has limitations in predicting fracture, and methods for clinical assessment of bone quality, such as microstructure, are awaited. This study was undertaken to examine the applicability of multidetector row CT (MDCT) for in vivo evaluation of trabecular microstructure.

MATERIALS AND METHODS

Optimal conditions for MDCT scanning were determined at a spatial resolution of 250 x 250 x 500 mum, using muCT data of excised human vertebra specimens as a reference. We analyzed the trabecular microstructure of the vertebrae of 82 postmenopausal women (55-76 years old), including 39 women with and 43 without a recent vertebral fracture.

RESULTS

Microstructure indices obtained by MDCT scanning revealed higher relative risk for prevalent vertebral fracture (OR: 16.0 for structure model index, 13.6 for bone volume fraction, and 13.1 for Euler's number) than did spinal BMD obtained by DXA (OR: 4.8). MDCT could also provide volumetric BMD data, which had higher diagnostic value (OR: 12.7) than did DXA.

CONCLUSION

Vertebral microarchitecture can be visualized by MDCT, and microstructure parameters obtained by MDCT, together with volumetric BMD, provided better diagnostic performance for assessing fracture risk than DXA measurement.

Anisotropy changes in post-menopausal osteoporosis: characterization by a new index applied to trabecular bone radiographic images

Bone intrinsic strength is conditioned by several factors, including material property and trabecular micro-architecture. Bone mineral density (BMD) is a good surrogate for material property. Architectural anisotropy is of special interest in mechanics-architecture relations and characterizes the degree of directional organization of a material. We have developed anisotropy indices from the Fast Fourier Transform (FFT) on bone radiographs. We have validated these indices in a cross-sectional uni-center case-control study including 39 postmenopausal women with vertebral fracture and 70 age-matched control cases. BMD was measured at the lumbar spine and femoral neck. A fractal analysis of texture was performed on calcaneus radiographs at three regions of interest (ROIs), and the result was expressed as the H parameter (fractal dimension =H-2). The anisotropy evaluation was based on the FFT spectrum of these three ROIs extracted on calcaneus radiographs. On the FFT spectrum, we have measured the spreading angle of the longitudinal trabeculae called the dispersion longitudinal index (DLI) and the spreading angle of the transversal trabeculae called the dispersion transversal index (DTI). From the measured parameters, an anisotropy index was derived, and the degree of anisotropy (DA) calculated with DLI and DTI. We have compared the results from the vertebral fracture cases and control cases. The best distinction was obtained for the largest ROI located in the great tuberosity of the calcaneus for all parameters ( P <10(-4)).( )The DA parameter showed a higher value in vertebral fracture cases (1.746+/-0.169) than in control cases (1.548+/-0.136); P <10(-4), and the difference persisted after removal of the subjects with hormonal replacement therapy. The analysis of the receiver operating characteristics (ROC) has shown the best results with DA and Hmean: areas under curves (AUCs) respectively of 0.765 and 0.683, while AUCs associated to LS-BMD and FN-BMD were 0.614 and 0.591 lower, respectively. We determined the odds ratios (OR) by uni- and multivariate analysis. Crude ORs were respectively 3.91 (95% CI: 2.22-6.87) and 3.08 (95% CI: 1.72-5.52) for DA and Hmean. Crude ORs were respectively 1.71 (95% CI: 1.15-2.56) and 1.56 (95% CI: 1.05-2.31) for LS-BMD and FN-BMD. All ORs were statistically significant, and those associated to Hmean and anisotropy indices were higher than those of BMD measurements. From a multivariate analysis including anisotropy indices, Hmean, age and FN-BMD, the remaining significant ORs were respectively 6.33 (95% CI: 2.80-14.30) and 3.08 (95% CI: 1.48-6.37) for DA and Hmean. These data have shown that anisotropy indices on calcaneus radiographs can distinguish vertebral fracture cases from control cases. This analysis provides complementary information concerning the BMD and fractal parameter. These data suggest that we can improve the fracture risk evaluation by adding information related to the directional organization of trabecular bone derived from the FFT spectrum on conventional radiographic images.

Gender differences in trabecular bone architecture of the distal radius assessed with magnetic resonance imaging and implications for mechanical competence

High-resolution magnetic resonance imaging (hrMRI) has recently made it possible to evaluate trabecular bone structure in vivo. Despite obvious gender differences in fracture incidence at the distal radius, little is known about gender differences in trabecular bone microarchitecture and its relationship to the structural strength of the forearm. The aim of this study was to determine trabecular bone structure in the distal radius of elderly women and men and its correlation with failure loads of the distal radius as determined in a fall configuration. Specifically, we tested the hypotheses that structural indices differ between women and men and that they offer information that is independent from BMD for predicting structural strength. Intact right arms were obtained from 73 formalin-fixed cadavers (age 80+/-11 years, 43 women, 30 men). Trabecular structural indices (apparent bone volume fraction [app. BV/TV], trabecular number [app. Tb.N], trabecular separation [app. Tb.Sp], trabecular thickness [app. Tb.Th] and fractal dimension [Frac.Dim]) were assessed in the distal metaphysis, using hrMRI with 156 microm in-plane resolution and proprietary digital image analysis, while BMD was measured with dual X-ray absorptiometry (DXA). Women displayed significantly lower BMD (-29.8%, p <0.001), app. BV/TV (-8.2%, p <0.05) and app. Tb.Th (-10.2%, p <0.001) than men, whereas app. Tb.N, app. Tb.Sp. and fractal dimension did not differ significantly. Structural parameters differed between normal and osteopenic women (BV/TV: -11%, p <0.01; Tb.Th: -8%, p <0.001) and between normal and osteoporotic women BV/TV: -21%, p <0.001; Tb.Th: -16%, p <0.001). App. BV/TV, app. Tb.Th and fractal dimension provided information independent from BMD in the prediction of radial failure loads in multiple regression models. These findings imply that it should be of clinical interest to monitor both bone mass and trabecular microstructure for predicting osteoporotic fracture risk.

A new anisotropy index on trabecular bone radiographic images using the fast Fourier transform

BACKGROUND

The degree of anisotropy (DA) on radiographs is related to bone structure, we present a new index to assess DA.

METHODS

In a region of interest from calcaneus radiographs, we applied a Fast Fourier Transform (FFT). All the FFT spectra involve the horizontal and vertical components corresponding respectively to longitudinal and transversal trabeculae. By visual inspection, we measured the spreading angles: Dispersion Longitudinal Index (DLI) and Dispersion Transverse Index (DTI) and calculated DA = 180/(DLI+DTI). To test the reliability of DA assessment, we synthesized images simulating radiological projections of periodic structures with elements more or less disoriented.

RESULTS

Firstly, we tested synthetic images which comprised a large variety of structures from highly anisotropic structure to the almost isotropic, DA was ranging from 1.3 to 3.8 respectively. The analysis of the FFT spectra was performed by two observers, the Coefficients of Variation were 1.5% and 3.1 % for intra-and inter-observer reproducibility, respectively. In 22 post-menopausal women with osteoporotic fracture cases and 44 age-matched controls, DA values were respectively 1.87 +/- 0.15 versus 1.72 +/- 0.18 (p = 0.001). From the ROC analysis, the Area Under Curve (AUC) were respectively 0.65, 0.62, 0.64, 0.77 for lumbar spine, femoral neck, total femoral BMD and DA.

CONCLUSION

The highest DA values in fracture cases suggest that the structure is more anisotropic in osteoporosis due to preferential deletion of trabeculae in some directions.

Quantification of trabecular bone structure using magnetic resonance imaging at 3 Tesla--calibration studies using microcomputed tomography as a standard of reference

The purpose of this study is to use high-resolution magnetic resonance (MR) imaging at 3 Tesla (3T) to quantify trabecular bone structure in vitro using femoral head specimens, and to correlate the calculated structure measures with those that were determined using microcomputed tomography (microCT), the standard of reference. Fifteen cylindrical cores were obtained from fresh femoral heads after total hip arthroplasty. MR images were obtained at 3T using a transmit-receive wrist coil. High-resolution coronal images were acquired using a modified three-dimensional (3D) fast-gradient echo sequence. From these data sets two-dimensional (2D) structural parameters analogous to bone histomorphometry were derived by using both mean intercept length (MIL) methods based on the plate model and the more recent model-assumption free 3D distance-transformation (DT) methods. The parameters measured by the 2D plate model-based MIL method and the DT method included apparent (App). BV/TV (bone volume/total volume), App. Tb.Th (trabecular thickness), App. Tb.Sp (trabecular separation), and App. Tb.N (trabecular number). Identical regions of interest were analyzed in the MR images and the microCT data sets, and similar structure measures were derived. The means and standard deviations of the parameters over all slices were calculated and MR-derived measures were correlated with those derived from the microCT data sets using linear regression analyses. Structure measures were overestimated with MRI, for example, the mean App. BV/TV was 0.45 for MRI and 0.20 for microT, and the slope of the graph was 1.45. App. Tb.Th was overestimated by a factor of 1.9, whereas App. Tb.Sp was underestimated; Tb.N showed the smallest effect. Correlations between the individual parameters were excellent (App. BV/TV, r2 = 0.82; App. Tb.Sp, r2 = 0.84; App. Tb.N, r2 = 0.81), except for App.Tb.Th (r2 = 0.67). The results of this study show that trabecular bone structure measures may be obtained using 3T MR imaging. These measures, although higher than the standard of reference, show a highly significant correlation with true structure measures obtained by microCT.

Comparison of high-resolution MRI, optical microscopy and SEM for quantitation of trabecular architecture in the rat femur

Magnetic resonance imaging (MRI) has been used to analyze trabecular bone architecture in femur heads taken from adult Wistar rats. The aim of this study was to validate the use of MRI in assessing trabecular structure and morphology by comparing standard measures of bone morphology in the rat femur obtained from high resolution MRI with those obtained by conventional optical microscopy and by scanning electron microscopy (SEM). MR images were obtained on a Bruker 4.7 T micro-imaging system using a three-dimensional spin echo sequence with spatial resolution of 23 microm in-plane and a slice thickness of 39 microm. Optical images were obtained by de-calcifying the bone in EDTA and then sectioning 5-microm-thick slices. SEM images were obtained from bone embedded in epoxy resin with surface preparation by diamond polishing. Values of standard bone morphological parameters were compared and correlation coefficients between the MRI and the optical- and SEM-derived measures of morphology were calculated. Partial volume effects in MRI were minimized in this study by the use of very thin slices, yielding better agreement with optical- and SEM-derived measures of trabecular bone morphology than have been obtained in previous studies. Correlations between the MRI and optical data were significantly lower than those between the MRI and SEM data. Effects of de-calcification were also investigated. The results indicate that comparison of MRI with thin (de-calcified) optical images may be inherently flawed due to the destructive de-calcification and sectioning process used to prepare samples for the optical imaging.

Discrimination of osteoporotic patients with quantitative ultrasound using imaging or non-imaging device

OBJECTIVES

Quantitative ultrasound (QUS) has emerged as a new tool in the assessment of fracture risk. The aim of this study was to compare the clinical utility of QUS parameters measured using imaging and non-imaging devices in the discrimination of osteoporotic patients.

METHODS

QUS (Broadband Ultrasound Attenuation, BUA dB/MHz, and Speed of Sound, SOS m/s) were measured and then statistical analyses were performed.

RESULTS

The 106 women included were 65 +/- 8 years aged. Using DXA, T score was < or = -2.5 at either lumbar spine or hip in 59% of patients, and 25% had osteoporotic fractures. QUS results were different among devices, and these differences were highly dependent on the measured value. There was a similar effect of age and duration of menopause for all parameters. To obtain 90% of sensitivity for the diagnosis of osteoporosis, the thresholds were 50.80 and 71.70 dB/MHz for BUA and 1544.80 and 1551.50 m/s for SOS, using imaging and non-imaging devices, respectively. Belonging to the highest tertile of QUS had a negative predictive value for osteoporosis ranging from 59% to 65%. In the lowest tertile of QUS, the proportion of osteoporotic women was between 73% and 80%. All QUS parameters, except BUA measured with the non-imaging device, were able to discriminate post-menopausal women with fractures after adjustment for age and hip BMD.

CONCLUSIONS

Our data suggest that an imaging system improves the utility of BUA measurement, but not SOS, for post-menopausal osteoporosis assessment.

Relationship between computed tomographic image analysis and histomorphometry for microarchitectural characterization of human calcaneus

The present study aimed to characterize the relationships between several variables reflecting bone microarchitecture assessed by both computed tomographic (CT) image analysis and histomorphometry (conventional CT system) at the calcaneus. A total of 24 cadaveric specimens were studied. The mean age at death was 78 +/- 10 years (range, 53-93 years). A total of 15 sagittal sections (1 mm in width and spaced 2 mm apart) were selected for CT analysis; 6 undecalcified sections (7 microm) were analyzed for histomorphometry. The histomorphometric analysis was performed on a Leica Quantimet Q570 image analyzer. Features measured by both methods were: bone volume/tissue volume (BV/TV), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp), trabecular number (Tb.N), interconnectivity index (ICI), number of nodes (N Nd), number of terminus (N Tm), node-to-node strut count (NNS), node-to-terminus strut count (NTS), terminus-to-terminus strut count (TTS), marrow space star volume (SV), Euler number (EN), and fractal dimension (FD). The coefficient of correlations' values (simple linear regression) between histomorphometry and CT image analysis varied according to the parameters selected. R values were high for BV/TV, Tb.N, and Tb.Sp (range, 0.69-0.90; P < 0.01). R values were less significant for some variables also obtained from the binary image: SV (0.5, P < 0.05) and EN (0.43, P < 0.05). Finally R values were also significant for (two) variables obtained from skeletonized images, i.e., N Nd (0.4, P < 0.05) and N Tm (0.61, P < 0.01). Other correlations were not statistically significant. Moreover, for some variables the relationships between the two methods (CT analysis and histomorphometry) seemed best-described by using nonlinear models. For example, a logarithmic model was more appropriate for SV (r = 0.71, P < 0.01), N Nd (r = 0.52, P < 0.01). Finally the relationship between apparent (App) N Tm and N Tm was most satisfying when using an exponential model (r = 0.64, P < 0.01). In conclusion, trabecular bone structure measures determined on CT images show highly significant correlations with those determined using histomorphometry. The level of correlation varies according to the type of method used for characterizing bone structure, however, and the strongest correlations were found for the most basic features (Parfitt's parameters). Finally, for some variables, nonlinear models seem more appropriate.

Differences between the right and the left foot in calcaneal quantitative ultrasound measurements

A total of 621 women, aged 30-80 years, who were grouped according to having single (group 1) or duplicate measurements of their both calcanei with quantitative ultrasound (QUS) on the same day (group 2) or on a different day than the first measurement (group 3) were evaluated for differences between the right and left foot. Despite similar mean values of QUS indices on both sides, individual percentage differences were found varying from 7.3 to 9.5% in the quantitative ultrasound index (QUI), from 11.1 to 12.5% in broadband ultrasound attenuation (BUA), from 0.62 to 0.86% in speed of sound (SOS) and from 8.9 to 10.9% in estimated heel bone mineral density as measured using the Sahara Clinical Bone Sonometer in three groups of subjects. The percentage of subjects with a proportional difference exceeding the coefficient of variation of duplicate measurements of the same heel was the highest for BUA, varying from 63 to 76.7%, and ranged between 43.1 and 76.7% in other QUS indices. We conclude that there is a real inter-individual difference between the right and left foot in QUS parameters, whether measured once or twice or on different occasions. We recommend measuring both sides when using QUS to avoid misleading implications regarding a subject's bone status.

Singular spectrum analysis applied to backscattered ultrasound signals from in vitro human cancellous bone specimens

Mean scatterer spacing (MSS) holds particular promise for the detection of changes in quasiperiodic tissue microstructures such as may occur during development of disease in the liver, spleen, or bones. Many techniques that may be applied for MSS estimation (temporal and spectral autocorrelation, power spectrum and cepstrum, higher order statistics, and quadratic transformation) characterize signals that contain a mixture of periodic and nonperiodic contributions. In contrast, singular spectrum analysis (SSA), a method usually applied in nonlinear dynamics, first identifies components of signals corresponding to periodic structures and, second, identifies dominant periodicity. Thus, SSA may better separate periodic structures from nonperiodic structures and noise. Using an ultrasound echo simulation model, we previously demonstrated SSA's potential to identify MSS of structures in quasiperiodic scattering media. The current work aims to observe the behavior of MSS estimation by SSA using ultrasound measurements in phantom materials (two parallel, nylon-line phantoms and four foam phantoms of different densities). The SSA was able to estimate not only the nylon-line distances but also nylon-line thickness. The method also was sensitive to the average pore-size differences of the four sponges. The algorithms then were applied to characterize human cancellous bone microarchitectures. Using 1-MHz center-frequency, radio-frequency ultrasound signals, MSS was measured in 24 in vitro bone samples and ranged from 1.0 to 1.7 mm. The SSA MSS estimates correlate significantly to MSS measured independently from synchrotron microtomography, r2 = 0.68. Thus, application of SSA to backscattered ultrasound signals seems to be useful for providing information linked to tissue microarchitecture that is not evident from clinical images.

Development of a compact MRI system for trabecular bone volume fraction measurements

A compact MRI system for measuring trabecular bone volume fraction (TBVF) of the calcaneus was developed with the use of a 0.21 T permanent magnet and portable MRI console. The entire system weighed < 600 kg and was installed in a 2 m x 2 m space. Two cross-sectional spin-echo images of a heel acquired with external reference phantoms (total measurement time = 5 min) were used to quantify the TBVF of the calcaneus. The linearity and reproducibility of the measurements were evaluated by means of proton density-adjusted phantoms. Comparative measurements with quantitative ultrasound (QUS) in groups of healthy female volunteers showed a relatively high positive correlation (R(2) = 0.4539, 0.2693) between TBVF and the speed of sound (SOS). These results demonstrate the potential of this new system for measuring bone density. Magn Reson Med 52:440-444, 2004.

Osteoporosis imaging

Because osteoporotic fractures may be prevented, diagnostic techniques are essential in the assessment of osteoporosis. Conventional radiographs of the spine are not suited for diagnosing early osteoporosis, but they show fractures that may have no clinical symptoms. The radiologist should be aware of the enormous significance of these fractures for future osteoporotic fractures. Bone mass measurements are standard techniques in the diagnosis of osteoporosis, which are the basis of the WHO definition of osteoporosis. In this article the authors presented these standard techniques and newer diagnostic techniques that provide insights in the structure of trabecular bone.

Bone homogeneity factor: an advanced tool for the assessment of osteoporotic bone structure in high-resolution magnetic resonance images

RATIONALE AND OBJECTIVES

Osteoporosis is characterized by low bone mass and inferior structural competence. In this study we introduce the bone homogeneity factor (BHF) as a quantitative measurement of bone structure, which could be equally important as bone mineral density.

METHODS

BHF represents an advanced texture analysis tool based on the spatial autocorrelation function calculated in 9 different directions. These calculations were performed on high-resolution magnetic resonance images of the calcaneus at 3.0 T and compared with dual-energy x-ray absorptiometry measurements of the femoral neck.

RESULTS

The quality and resolution of the high-resolution magnetic resonance images is sufficient for reliably calculating BHF. The mean BHF of the control group (n = 5, mean BHF = 525,0) with normal bone is significantly (P = 0.009, Mann-Whitney U test) higher than in the osteoporotic group (n = 7, mean BHF = 137,8). The BHF correlates with the DXA measurements of the femoral neck (correlation coefficient = 0.75).

CONCLUSIONS

By calculating the BHF, it was possible to distinguish between osteoporotic and nonosteoporotic bone structure. Hence, BHF could be a possible candidate for noninvasive assessment of osteoporotic bone structure giving additional information to routinely used bone mineral densitometry.

Noninvasive assessment of bone structure

Bone fragility is determined by bone mass and trabecular structure. While bone mass can be readily measured as bone density, bone trabecular structure cannot be easily assessed by currently available methods. The realization of the importance of bone structure in determining fracture risk has led to the development of several imaging modalities aimed at evaluating the contribution of bone quality to its biomechanical strength and fragility. High-resolution magnetic resonance imaging and computed tomography have limited spatial resolution and high cost but have a potential to generate true three-dimensional images of trabecular structure in vivo. Bone radiographs subjected to various forms of texture analysis have higher resolution and lower cost but provide only a two-dimensional representation of bone structure. Both two- and three-dimensional methods have been shown to predict biomechanical strength in vitro and to differentiate between subjects with and without fractures in vivo. Therefore, all of these methods deserve closer evaluation and also need further technical improvements before they can be considered for use in clinical practice.

Trabecular bone structure of the calcaneus: preliminary in vivo MR imaging assessment in men with osteoporosis

PURPOSE

To use magnetic resonance (MR) imaging to evaluate potential differences in bone structure between men with and men without osteoporosis.

MATERIALS AND METHODS

Sagittal MR images of the calcaneus were obtained in 50 men (26 patients with osteoporosis and 24 age-matched healthy control subjects). Osteoporosis was defined as a low bone mineral density (at least 2.5 SDs below the normal value for young adults at either the lumbar spine or proximal femur) as measured with dual-energy x-ray absorptiometry. Seventeen patients had a history of osteoporotic fractures. For each participant, 10 consecutive sagittal three-dimensional gradient-echo MR sections were analyzed by using a rectangular region of interest. Twenty structural measurements were obtained from these images. Additionally, density measurements at the calcaneus were obtained in 46 participants. The significance of differences between the two groups was calculated by using the unpaired Student t test. The odds ratios for fracture per 1 SD decrease in the control group were calculated with logistic regression analysis. Adjustment for participant weight and height was performed if necessary.

RESULTS

Thirteen of 20 structural parameters, especially connectivity parameters, showed significant differences between control subjects and patients (P <.05). Differences between the two groups were more significant (P <.001) for apparent bone marrow skeleton length, apparent node count, apparent node-to-node strut count, and apparent terminus-to-terminus strut count. Odds ratios for 11 of 13 structural parameters but not for calcaneus density were significant (P <.05). After adjustment for calcaneus density, these parameters were still significant predictors of osteoporotic fracture.

CONCLUSION

Structural measurements derived from MR images of the calcaneus may be used in vivo to characterize trabecular bone architecture in men with osteoporosis.

Local 3D scaling properties for the analysis of trabecular bone extracted from high-resolution magnetic resonance imaging of human trabecular bone: comparison with bone mineral density in the prediction of biomechanical strength in vitro

RATIONALE AND OBJECTIVES

A novel, nonlinear morphologic measure [DeltaP(alpha)] based on local 3D scaling properties was applied to high-resolution magnetic resonance images (HR-MRI) of human trabecular bone to predict biomechanical strength in vitro.

METHODS

We extracted DeltaP(alpha) and traditional morphologic parameters (apparent trabecular volume fraction, apparent trabecular separation) from HR-MR images of 32 femoral and 13 spinal bone specimens. Furthermore, bone mineral density (BMD) and maximum compressive strength (MCS) were determined. The morphologic measures were compared with BMD in predicting the biomechanical strength.

RESULTS

In the vertebral (femoral) specimens, R2 for MCS versus DeltaP(alpha) was 0.87 (0.61) (P < 0.001). Correlation between BMD and MCS was 0.53 (P = 0.05) (0.79 [P < 0.001]) for the vertebral (femoral) specimens. For the femoral specimens, prediction of MCS could be improved further by combining BMD and morphologic parameters by multiple regression (R2 = 0.88).

CONCLUSIONS

Morphologic measures extracted from HR-MRI considering local 3D-scaling properties can be used to predict biomechanical properties of bone in vitro. They are superior to 2-dimensional standard linear morphometric measures and, depending on the anatomic location, more reliably predict bone strength as measured by MCS than does BMD.

Quantitative magnetic resonance imaging in the calcaneus and femur of women with varying degrees of osteopenia and vertebral deformity status

Quantitative magnetic resonance imaging (QMRI) allows measurement of two parameters that are related to the integrity of the trabecular bone: R2*, the rate constant of the free induction signal, and trabecular bone volume fraction (BVF), the counterpart of apparent density. In this work, R2* and BVF were measured in 68 women (mean age, 58.2 +/- 9.5 years) of varying spinal bone mineral density (BMD) T scores (mean, -1.37 +/- 1.54) and vertebral fracture status on a commercial 1.5 T whole-body imager using customized image acquisition and processing techniques. Twenty-five of the patients had vertebral fractures, characterized by the total cumulative deformity burden exceeding 200%. R2* was measured in the calcaneus and proximal femur and BVF could be measured in the calcaneus only. On a pixel-by-pixel basis, calcaneal R2* and BVF within each subject were highly positively correlated (r2 = 0.61 +/- 0.11) but the correlation of region-of-interest (ROI) means for different calcaneal sites among patients was weaker (r2 = 0.34; p < 0.0001). The strongest discriminator of vertebral deformity was R2* of the calcaneus, which was lower in the fracture group, consistent with lower trabecular density. Among the calcaneal sites examined, the subtalar region, a location characterized by dense nearly horizontal trabeculae that transmit the stresses imparted by body weight from the tibia to the heel, best discriminated the two groups (p = 0.0001), with 77% diagnostic accuracy as determined from the area under the receiver operating characteristic (ROC) curve (compared with 66% for vertebral BMD). The cavum calcanei, an anterior site of low trabecular density, and the tuber calcanei (the location ordinarily used for ultrasound measurements) also had significantly reduced R2* in the fracture group (p < 0.005 and p = 0.01, respectively). The R2av*, computed as the average of all pixels in the calcaneus, was a strong discriminator as well (p < 0.005). On the other hand, calcaneal BVF was only marginally discriminating (p = 0.05). Among the BMD sites examined, the lumbar spine (average L1-L4) was significant (p = 0.005, 66% diagnostic accuracy), as was the femoral neck (p = 0.01). The data suggest the calcaneus to be suited as a surrogate site to assess vertebral osteoporosis and that R2* is sensitive to alterations in bone quality not captured by density.

Local differences in the trabecular bone structure of the proximal femur depicted with high-spatial-resolution MR imaging and multisection CT

RATIONALE AND OBJECTIVES

The authors performed this study to investigate structural variations in the trabecular bone of the proximal femur at high-resolution magnetic resonance (MR) imaging and high-resolution multisection computed tomography (CT).

MATERIALS AND METHODS

Bone mineral density (BMD) was measured in 36 proximal human femur specimens by using dual x-ray absorptiometry. High-resolution MR imaging was performed at 1.5 T with an in-plane spatial resolution of 0.195 x 0.195 mm and a section thickness of 0.3 and 0.9 mm. Multisection CT was performed with an ultra-high-resolution protocol; images were obtained with an in-plane spatial resolution of 0.25 mm and a section thickness of 1 mm. In a subset of these specimens, micro CT was performed with an isotropic spatial resolution of 30 microm. Identical regions of interest (ROIs) were used to analyze images obtained with MR imaging, multisection CT, and micro CT. Trabecular bone structural parameters were obtained, and the parameters from the individual imaging modalities and BMD were correlated.

RESULTS

Significant differences concerning the trabecular microarchitecture between the individual ROIs were demonstrated with multisection CT and MR imaging. A number of the correlations between structural parameters derived with multisection CT, MR imaging, micro CT, and BMD measurements were significant. For MR imaging, threshold technique and section thickness had an effect on structural parameters.

CONCLUSION

Structural parameters obtained in the proximal femur with multisection CT and high-resolution MR imaging show regional differences. These techniques may be useful for depicting the trabecular architecture in the diagnosis of osteoporosis.

Three-dimensional-line skeleton graph analysis of high-resolution magnetic resonance images: a validation study from 34-microm-resolution microcomputed tomography

The resolution achievable in vivo by magnetic resonance imaging (MRI) techniques is not sufficient to depict precisely individual trabeculae and, thus, does not permit the quantification of the "true" trabecular bone morphology and topology. Nevertheless, the characterization of the "apparent" trabecular bone network derived from high-resolution MR images (MRIs) and their potential to provide information in addition to bone mineral density (BMD) alone has been established in studies of osteoporosis. The aim of this work was to show the ability of the three-dimensional-line skeleton graph analysis (3D-LSGA) to characterize high-resolution MRIs of trabecular bone structure. Fifteen trabecular bone samples of the distal radius were imaged using the high-resolution MRI (156 x 156 x 300 microm3) and microcomputed tomography (microCT; 34 x 34 x 34 microm3). After thresholding, the 3D skeleton graph of each binary image was obtained. To remove the assimilated-noise branches of the skeleton graph and smooth this skeleton graph before it was analyzed, we defined a smoothing length criterion (l(c)), such that all "termini" branches having a length lower than l(c) were removed. Local topological and morphological LSGA measurements were performed from MRIs and microCT images of the same samples. The correlations between these two sets of measurements were dependent on the smoothing criterion l(c), reaching R2 = 0.85 for topological measurements and R2 = 0.57-0.64 for morphological measurements. 3D-LSGA technique could be applied to in vivo high-resolution MRIs of trabecular bone structure, giving an indirect characterization of the microtrabecular bone network.

Influence of region of interest and bone size on calcaneal BMD: implications for the accuracy of quantitative ultrasound assessments at the calcaneus

There is considerable technological diversity among quantitative ultrasound (QUS) devices used to assess osteoporosis. Because the distance between the transducer and the footplate remains constant, the location of the calcaneus measured will vary with foot size. This study was designed to quantify the variation in bone mineral density (BMD) between a manufacturer's region of interest (ROI_M), which is fixed relative to the footplate, and an anatomical region of interest (ROI_A), which is defined as 20% of calcaneal length. The effect of foot length and width on QUS variables measured using two Food and Drug Administration cleared QUS devices, the Sahara (Hologic) and the Achilles+ (Lunar) was assessed. 26 healthy subjects (12 male and 14 female), aged 22-54 years (35.6+/-10 years) and with foot lengths of 21.5 cm to 29.7 cm (25.1+/-2.3 cm) were recruited. QUS assessments were performed at the right calcaneus. In addition, a Hologic 4500 densitometer was used to measure the BMD of the calcaneus in the ROI_M and ROI_A. The sizes of the ROIs were approximated to the sizes of the transducers of the Sahara and Achilles+ devices. The results showed a significant difference in BMD between the two ROI locations for the Sahara device (BMD 0.642+/-0.135 g cm(-2) vs 0.616+/-0.114 g cm(-2), p=0.014), but no significant difference was found in BMD between the two locations for the Achilles device (BMD 0.661+/-0.120 g cm(-2) vs 0.662+/-0.123 g cm(-2), p=0.818). At the ROI_A, there was a significant difference in BMD between the two QUS devices (p<0.001). The correlation between QUS variables and BMD was slightly higher for the ROI_M (r=0.68-0.79, since this is site-matched) than the ROI_A (r=0.59-0.70) for the Achilles device, while for the Sahara device the correlations were r=0.35-0.40 and r=0.51-0.54, respectively. The smaller ROI of the Sahara device resulted in more than 50% of the subjects having BMD differences of greater than 5% between the ROI_A and the ROI_M, compared with only 20% of the subjects on the Achilles device. ROIs containing cortical bone edge and other soft tissues were found in 58% of cases for the Achilles device and 46% of cases for the Sahara device. The greatest differences occurred in very small and very large feet. Calcaneal length correlated significantly with Sahara speed of sound (SOS), and heel width correlated significantly with Achilles SOS. Heel width also correlated significantly with Sahara broadband ultrasound attenuation (BUA) but not Achilles+ BUA. These results suggest that variation in ROI and bone size might affect the accuracy of QUS measurements, since the calcaneus is heterogeneous both in terms of its external geometry and its internal structure and density.

Future methods in the assessment of bone mass and structure

There have been major advances in the diagnosis of osteoporosis over the last few decades not only in the definitions that are now used but also in the technology that is available. The future will see further development of the techniques currently in common clinical use, such us dual energy X-ray absorptiometry and quantitative ultrasound. In addition new techniques for assessing bone structure, including MRI and fractal analysis of X-rays, may add significantly to our understanding of the pathophysiology of osteoporosis and to the prediction of fracture risk.

Fractal analysis of radiographic trabecular bone texture and bone mineral density: two complementary parameters related to osteoporotic fractures

Trabecular bone microarchitecture and bone mineral density (BMD) are two main factors related to osteoporotic fractures. Currently, however, microarchitecture is not evaluated. We have developed and validated a trabecular bone texture analysis from radiographic images. The objective was to determine if the fractal analysis of texture was able to distinguish osteoporotic fracture groups from control groups, either in vertebrae, hip, or wrist fractures, and to determine if this indicator and BMD were independent and complementary. In this cross-sectional unicenter case-control population study in postmenopausal women, 107 fracture cases were enrolled and age-matched with 197 control cases. This population comprised 40 vertebral fractures (with 70 controls), 30 hip fractures (55 controls), and 37 wrist fractures (62 controls). Hip and lumbar spine BMD were measured by double-energy X-ray absorptiometry. Fractal analysis of texture was performed on calcaneus radiographs and the result was expressed as the H parameter (H = 2-fractal dimension). The H parameter showed a lower value (0.679 +/- 0.053 SD) in fracture cases versus control cases (0.696 +/- 0.030; p = 0.007), the statistical significance persisting after adjustment for age and for lumbar spine (LS) or hip BMD. This result was confirmed in vertebral fractures (p = 0.0001) and hip fractures (p = 0.003) but not wrist fractures (p = 0.07). We determined the threshold between high and low H values and then the odds ratios (OR) of fracture for low H for BMD < or = -2.5 SD in T score and for the combinations of both parameters. The OR of fracture for low H was 1.6 (95% CI, 1.1-2.6). For LS BMD < or = -2.5 SD the OR of 6.1 (3.4-10.8) shifted to 9.0 (4.0-20.4) when we added low H and for hip BMD it shifted from 5.6 (3.3-9.4) to 8.1 (4.0-16.8). In vertebral, hip, and wrist fracture cases the results were also significant. These data have shown that the fractal analysis of texture on calcaneus radiographs can distinguish osteoporotic fracture groups from control groups. This analysis and BMD provide independent and complementary information. These data suggest that we can improve the fracture risk evaluation by adding information related to microarchitecture, derived from analysis of conventional radiographic images.

Does the trabecular bone structure depicted by high-resolution MRI of the calcaneus reflect the true bone structure?

RATIONALE AND OBJECTIVES

The purpose of this study was to compare trabecular bone structure parameters assessed with high-resolution magnetic resonance imaging (HR-MRI) with those determined in specimen sections.

METHODS

High-resolution MR images were obtained for 30 calcaneus specimens with a three-dimensional, T1-weighted spin-echo sequence (spatial in-plane resolution 0.195 mm, slice thicknesses of 0.3 and 0.9 mm). Thirty-eight sections were obtained from the specimens, and contact radiography was performed. In the corresponding sections, structural parameters analogous to bone histomorphometry were determined.

RESULTS

Significant correlations between MRI-derived structural parameters and those derived from macro pathological sections were found: r values of up to 0.75 were obtained (P < 0.01). The highest correlations were found for apparent bone volume/total volume and trabecular thickness. Image thresholding techniques showed a significant impact on these correlations (P < 0.01). The thinner MR sections were less susceptible to the different thresholding algorithms.

CONCLUSIONS

Trabecular bone structure depicted by HR-MR images is significantly correlated with that shown in macro sections (P < 0.01); however, a number of limitations have to be considered, including the substantial impact of thresholding techniques and slice thickness.

In vitro measurement of the frequency-dependent attenuation in cancellous bone between 0.2 and 2 MHz

Our goal was to evaluate the frequency dependence of the ultrasonic attenuation coefficient in cancellous bone. Estimates were obtained in immersion, using a substitution method in the through-transmit mode, by scanning 14 human bone specimens (calcaneus). Measurements were performed with three pairs of focused transducers with a center frequency of 0.5, 1.0, and 2.25 MHz, respectively in order to cover an extended frequency bandwidth (0.2-1.7 MHz). When the experimental attenuation coefficient values were modeled with a nonlinear power fit alpha(f)=alpha0 +alpha(I)f(n), the attenuation coefficient was found to increase as f(1.09+/-0.3) over the measurement bandwidth. However, a substantial variation of the exponent n (0.4-2.2) within specimens and also between specimens was observed. The acoustical parameters were compared to bone mineral density. A highly significant relationship was noted between alpha1 and BMD (r2= 0.75, p< 10(-4)). No correlation was found between n and BMD. Several attenuation mechanisms are discussed as well as the potential impact these results may have in in vivo quantitative measurements.

Quantitative ultrasound in postmenopausal osteoporosis

Ultrasound has been proposed as a low-cost, radiation-free method for osteoporosis assessment in postmenopausal women. Large prospective studies have shown that ultrasound parameters can be used for fracture risk estimate in this population, providing that adequate quality control is performed. The places of both ultrasound and the current gold standard method for bone assessment, dual energy x-ray absorptiometry, are still to be determined. Further studies are needed on the diagnosis of osteoporosis using ultrasound, because current diagnostic thresholds, designed by the World Health Organization, do not apply to this-new technology. Monitoring of skeletal changes and treatment effects by ultrasound cannot be recommended.

Quantitative ultrasound of the calcaneus with parametric imaging: correlation with bone mineral density at different sites and with anthropometric data in menopausal women

OBJECTIVE

To prospectively study the relationship of quantitative ultrasound of the calcaneus with anthromopometric variables and with bone mineral density (BMD) assessed at the level of the calcaneus as well as at other sites.

METHOD

Osteosonography of the non-dominant calcaneus was performed in 135 menopausal women, using a DTU-one device with parametric imaging. Broadband ultrasound attenuation (BUA) and speed of sound (SOS) were assessed. BMD of the calcaneus (BMDcal) was measured using dual energy X-ray absorptiometry (DXA), in a subregion matched with the region of interest for osteosonography. BMD of the lumbar trabecular bone was measured using quantitative computed tomography (BMD QCT) while the non-dominant hip was studied using DXA, which provided the total bone mineral density (BMDhip) and that of the Ward triangle (BMDWard).

RESULTS

The Pearson correlation coefficients between BUA, SOS and the various measurements of BMD ranged from 0.305 (SOS versus BMDhip) to 0.717 (BUA versus BMDcal). BMD QCT and BMDWard were found to depend on age, but not on weight or height, while BUA, SOS, BMDcal, BMDhip were unrelated to age, but correlated with weight (SOS, BMDhip) or with weight and height (BUA, BMDcal). In a multiple stepwise regression analysis, age was a significant predictor for BMD QCT, BMD hip and BMDWard; BMD QCT, BMDWard and BMDhip admitted BUA as sole predictor, while BMDcal was significantly related to both BUA and SOS.

CONCLUSION

BUA and SOS of the calcaneus, assessed in 135 menopausal women using a parametric imaging device, reflected BMDcal, measured with DXA at a matched region of interest, and did not decline significantly with age.