CT image analysis of the vertebral trabecular network in vivo

Regional Distribution of Stress on the Distal Femur in Advanced Osteoarthritis

Background

The aim of this study was to analyze regional differences in the stress distribution within the distal femur in advanced osteoarthritis (OA).

Methods

Distal femoral specimens with primary OA were obtained from 10 female donor cadavers (mean age, 65 years; range, 53–79). As controls, distal femurs without OA were obtained from 10 age- and sex-matched female cadavers (mean age, 67 years; range, 58–81). The articular surface of the distal femur was divided into anterior, middle, and posterior regions on each condyle. Mechanical properties and microstructure were assessed for each region with micro-computed tomography and finite element model analysis.

Results

The control group showed differences in stress distribution among 6 regions on the distal femur (P=0.037), but there was no regional difference in stress distribution among 6 regions on the distal femur in the advanced OA group (P=0.179).

Conclusions

Regional stress distribution in the distal femur was different between advanced OA and normal groups. There were no regional differences in stress distribution in the advanced OA group. Altered loading patterns, bone remodeling, and chemical composition will affect stress distribution.

Qualitative and quantitative ultrashort-TE MRI of cortical bone

Osteoporosis causes over 1.5 million fractures per year, costing about $15 billion annually in the USA. Current guidelines utilize bone mineral density (BMD) to assess fracture risk; however, BMD alone only accounts for 30-50% of fractures. The other two major components of bone, organic matrix and water, contribute significantly to bone mechanical properties, but cannot be assessed with conventional imaging techniques in spite of the fact that they make up about 57% of cortical bone by volume. Conventional clinical MRI usually detects signals from water in tissues without difficulty, but cannot detect the water bound to the organic matrix, or the free water in the microscopic pores of the Haversian and the lacunar-canalicular system of cortical bone, because of their very short apparent transverse relaxation times (T2 *). In recent years, a new class of sequences, ultrashort-TE (UTE) sequences, with nominal TEs of less than 100 µs, which are much shorter than the TEs available with conventional sequences, have received increasing interest. These sequences can detect water signals from within cortical bone and provide an opportunity to study disease of this tissue in a new way. This review summarizes the recent developments in qualitative UTE imaging (techniques and contrast mechanisms to produce bone images with high contrast) and quantitative UTE imaging (techniques to quantify the MR properties, including T1 , T2 * and the magnetization transfer ratio, and tissue properties, including bone perfusion, as well as total, bound and free water content) of cortical bone in vitro and in vivo. The limitations of the current techniques for clinical applications and future directions are also discussed.

The Founder's Lecture 2009: advances in imaging of osteoporosis and osteoarthritis

The objective of this review article is to provide an update on new developments in imaging of osteoporosis and osteoarthritis over the past three decades. A literature review is presented that summarizes the highlights in the development of bone mineral density measurements, bone structure imaging, and vertebral fracture assessment in osteoporosis as well as MR-based semiquantitative assessment of osteoarthritis and quantitative cartilage matrix imaging. This review focuses on techniques that have impacted patient management and therapeutic decision making or that potentially will affect patient care in the near future. Results of pertinent studies are presented and used for illustration. In summary, novel developments have significantly impacted imaging of osteoporosis and osteoarthritis over the past three decades.

Trabecular bone structure analysis in the osteoporotic spine using a clinical in vivo setup for 64-slice MDCT imaging: comparison to microCT imaging and microFE modeling

Assessment of trabecular microarchitecture may improve estimation of biomechanical strength, but visualization of trabecular bone structure in vivo is challenging. We tested the feasibility of assessing trabecular microarchitecture in the spine using multidetector CT (MDCT) on intact human cadavers in an experimental in vivo-like setup. BMD, bone structure (e.g., bone volume/total volume = BV/TV; trabecular thickness = Tb.Th; structure model index = SMI) and bone texture parameters were evaluated in 45 lumbar vertebral bodies using MDCT (mean in-plane pixel size, 274 microm(2); slice thickness, 500 microm). These measures were correlated with structure measures assessed with microCT at an isotropic spatial resolution of 16 microm and to microfinite element models (microFE) of apparent modulus and stiffness. MDCT-derived BMD and structure measures showed significant correlations to the density and structure obtained by microCT (BMD, R(2) = 0.86, p < 0.0001; BV/TV, R(2) = 0.64, p < 0.0001; Tb.Th, R(2) = 0.36, p < 0.01). When comparing microCT-derived measures with microFE models, the following correlations (p < 0.001) were found for apparent modulus and stiffness, respectively: BMD (R(2) = 0.58 and 0.66), BV/TV (R(2) = 0.44 and 0.58), and SMI (R(2) = 0.44 and 0.49). However, the overall highest correlation (p < 0.001) with microFE app. modulus (R(2) = 0.75) and stiffness (R(2) = 0.76) was achieved by the combination of QCT-derived BMD with the bone texture measure Minkowski Dimension. In summary, although still limited by its spatial resolution, trabecular bone structure assessment using MDCT is overall feasible. However, when comparing with microFE-derived bone properties, BMD is superior compared with single parameters for microarchitecture, and correlations further improve when combining with texture measures.

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

Objectives

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

Materials and methods

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

Results

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

Conclusion

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

Review of nonprimate, large animal models for osteoporosis research

Large animal models are required for preclinical prevention and intervention studies related to osteoporosis research. The challenging aspect of this requirement is that no single animal model exactly mimics the progression of this human-specific chronic condition. There are pros and cons associated with the skeletal, hormonal, and metabolic conditions of each species that influence their relevance and applicability to human physiology. Of all larger mammalian species, nonhuman primates (NHPs) are preeminent in terms of replicating important aspects of human physiology. However, NHPs are very expensive, putting them out of reach of the vast majority of researchers. Practical, cost-effective alternatives to NHPs are sought after among ungulate (porcine, caprine, and ovine) and canine species that are the focus of this review. The overriding caveat to using large lower-order species is to take the time in advance to understand and appreciate the limitations and strengths of each animal model. Under these circumstances, experiments can be strategically designed to optimize the potential of an animal to develop the cardinal features of postmenopausal bone loss and/or yield information of relevance to treatment.

Advanced CT bone imaging in osteoporosis

Non-invasive and/or non-destructive techniques can provide structural information about bone, beyond simple bone densitometry. While the latter provides important information about osteoporotic fracture risk, many studies indicate that BMD only partly explains bone strength. Quantitative assessment of macro- and microstructural features may improve our ability to estimate bone strength. Methods for quantitatively assessing macrostructure include (besides conventional radiographs) DXA and CT, particularly volumetric quantitative CT (vQCT). Methods for assessing microstructure of trabecular bone non-invasively and/or non-destructively include high-resolution CT (hrCT), microCT (μCT), high-resolution magnetic resonance (hrMR) and microMR (μMR). vQCT, hrCT and hrMR are generally applicable in vivo; μCT and μMR are principally applicable in vitro. Despite recent progress made with these advanced imaging techniques, certain issues remain. The important balances between spatial resolution and sampling size, or between signal-to-noise and radiation dose or acquisition time, need further consideration, as do the complexity and expense of the methods vs their availability and accessibility. Clinically, the challenges for bone imaging include balancing the advantages of simple bone densitometry vs the more complex architectural features of bone or the deeper research requirements vs the broader clinical needs. The biological differences between the peripheral appendicular skeleton and the central axial skeleton must be further addressed. Finally, the relative merits of these sophisticated imaging techniques must be weighed with respect to their applications as diagnostic procedures, requiring high accuracy or reliability, compared with their monitoring applications, requiring high precision or reproducibility.

Trabecular structure quantified with the MRI-based virtual bone biopsy in postmenopausal women contributes to vertebral deformity burden independent of areal vertebral BMD

In postmenopausal women with a wide range of vertebral deformities, MRI-based structural measures of topology and scale at the distal radius are shown to account for as much as 30% of vertebral deformity, independent of integral vertebral BMD.

INTRODUCTION

Trabecular bone architecture has been postulated to contribute to overall bone strength independent of vertebral BMD measured by DXA. However, there has thus far been only sparse in vivo evidence to support this hypothesis.

MATERIALS AND METHODS

Postmenopausal women, 60-80 yr of age, were screened by DXA, and those with T-scores at either the hip or spine falling within the range of -2.5 +/- 1.0 were studied with the MRI-based virtual bone biopsy, along with heel broadband ultrasound absorption and pQCT of the tibia. The data from 98 subjects meeting the enrollment criteria were subjected to microMRI at the distal tibia and radius, and measures of topology and scale of the trabecular bone network were computed. A spinal deformity index (SDI) was obtained from morphometric measurements in midline sagittal MR images of the thoracic and lumbar spine to evaluate associations between structure and deformity burden.

RESULTS

A number of structural indices obtained at the distal radius were correlated with the SDI. Among these were the topological surface density (a measure of trabecular plates) and trabecular bone volume fraction, which were inversely correlated with SDI (p < 0.0001). Combinations of two structural parameters accounted for up to 30% of the variation in SDI (p < 0.0001) independent of spinal BMD, which was not significantly correlated. pQCT trabecular BMD was also weakly associated, whereas broadband ultrasound absorption was not. No significant association between SDI and structural indices were found at the tibia.

CONCLUSIONS

Structural measures at the distal radius obtained in vivo by microMRI explained a significant portion of the variation in total spinal deformity burden in postmenopausal women independent of areal BMD.

Monitoring teriparatide-associated changes in vertebral microstructure by high-resolution CT in vivo: results from the EUROFORS study

We introduce a method for microstructural analysis of vertebral trabecular bone in vivo based on HRCT. When applied to monitor teriparatide treatment, changes in structural variables exceeded and were partially independent of changes in volumetric BMD.

INTRODUCTION

Monitoring of osteoporosis therapy based solely on bone densitometry is insufficient to assess anti-fracture efficacy. Assessing bone microstructure in vivo is therefore of importance. We studied whether it is possible to monitor effects of teriparatide on vertebral trabecular microstructure independent of BMD by high-resolution CT (HRCT).

MATERIALS AND METHODS

In a subset of 65 postmenopausal women with established osteoporosis who participated in the EUROFORS study, HRCT scans of T(12), quantitative CT of L(1)-L(3), and DXA of L(1)-L(4) were performed after 0, 6, and 12 mo of teriparatide treatment (20 microg/d). We compared BMD and 3D microstructural variables in three groups of women, based on prior antiresorptive treatment: treatment-naïve; pretreated; and pretreated women showing inadequate response to treatment.

RESULTS

We found statistically highly significant increases in most microstructural variables and BMD 6 mo after starting teriparatide. After 12 mo, apparent bone volume fraction (app. BV/TV) increased by 30.6 +/- 4.4% (SE), and apparent trabecular number (app. Tb.N.) increased by 19.0 +/- 3.2% compared with 6.4 +/- 0.7% for areal and 19.3 +/- 2.6% for volumetric BMD. The structural changes were partially independent of BMD as shown by a significantly larger standardized increase and a standardized long-term precision at least as good as DXA. Patients who had shown inadequate response to prior osteoporosis treatment did show improvements in BMD and structural measures comparable to treatment-naïve patients.

CONCLUSIONS

HRCT is a feasible method for longitudinal microstructural analysis of human vertebrae in vivo, offers information beyond BMD, and is sufficiently precise to show profound effects of teriparatide after 12 mo.

The long-term effects of ovariectomy on bone metabolism in sheep

Osteoporosis and associated fractures are major public health concerns, and as such require appropriate large animal models to further our understanding of this disease. Although sheep appear to be an ideal model with which to study bone loss caused by estrogen depletion, limited data are available concerning the long-term effect of ovariectomy on bone in sheep. The goal of the present study was to observe the ovariectomy-induced changes in bone mass, structure, and metabolism in sheep over a period of 18 months. Six ewes were ovariectomized (OVX) and compared to an age-matched control group by analyzing bone mineral density, trabecular structure, biochemical markers of bone formation and resorption, and plasma estrogen levels. Bone loss (13%, P < 0.01) occurred during the first 4 months after surgery, then stabilized and returned to pre-OVX levels for the remainder of the study. Trabecular architecture was also altered and tended toward osteopenia with recovery to baseline values. Markers of bone formation and resorption were elevated up to 6 months postovariectomy, after which time levels returned to baseline values. Although estradiol measurements demonstrated a clear decline following surgical ovariectomy, levels returned to normal after 6 months. Therefore, the detrimental effect of ovariectomy on sheep bone metabolism seems to be reversible, with normal bone parameters being reestablished within 6 months after surgery. These data seem to indicate that the sheep is not an appropriate model for human postmenopausal osteoporosis.

Protective Effects of Sodium Daidzein Sulfonate on Trabecular Bone in Ovariectomized Rats

The ovariectomized (OVX) rat, as an established animal model of human osteoporosis, was adopted in the present experiment to study the protective effects of sodium daidzein sulfonate (SDS) on trabecular bone. Six-month-old Sprague-Dawley rats were sham-operated or ovariectomized. Five days later, the OVX rats were randomly assigned to one of three experimental groups and treated for 90 days with vehicle, 17beta-estradiol (E(2)) or SDS. Compared with OVX rats, SDS administration (15 mg/kg) prevented OVX-induced decrease in lumbar vertebral and femoral bone mineral density (BMD), and significantly increased bone mechanical strength parameters, including ultimate stress and elastic modulus. In the OVX group, the structure of trabecular plate in the femoral head was absorbed and became progressively thinner or was removed completely, accompanied by enlargement of marrow cavities and amalgamation of two or more marrow cavities. Administration of SDS and E(2 )prevented the change of trabecular bone microarchitecture induced by OVX, increasing the trabecular bone area and trabecular thickness, while decreasing the trabecular separation. These results indicate that SDS administration prevents OVX-induced decrease in BMD and bone mechanical strength, and has a moderate protective effect on the microarchitecture of trabecular bone in aged Sprague-Dawley rats.

Neue Techniken in der Osteoporosediagnostik

In this review article current developments and applications in quantitative osteoporosis imaging are presented. Developments in the field of DXA include geometrical parameters of the proximal femur such as the "hip axis length" and new ROIs to determine BMD. Advances in QCT are new volumetric techniques to quantify BMD at the lumbar spine and the proximal femur. In addition techniques to determine BMD in standard contrast-enhanced abdominal computed tomography studies are described. Currently with the new bone quality concept in full bloom techniques to quantify trabecular bone architecture as new surrogates of bone strength are of increasing significance. Spatial high-resolution techniques such as magnetic resonance imaging and new computed tomography techniques have shown their potential in assessing trabecular bone structure. In addition ultrasound is considered a low-cost technique to explore bone quality.

Clinical utility of microarchitecture measurements of trabecular bone

Osteoporosis is a metabolic disorder that manifests changes in bone density and structure accompanied by an increased susceptibility to fractures. Recent studies have demonstrated the potential contributions of trabecular bone microarchitecture in the assessment of the therapeutic efficacy of emerging treatments, and also in the assessment of fracture risk. The main goal of this paper is to emphasize the clinical implementation of bone microarchitecture measurements. Thus, this paper provides an overview of the main imaging modalities for depicting trabecular bone microarchitecture and a corresponding description of common computed structural bone parameters. The imaging modalities presented to characterize the complex three-dimensional trabecular bone network include micro-CT, quantitative CT, and magnetic resonance imaging. Two-dimensional analyses of radiographic patterns are also discussed. Results demonstrating the ability to distinguish between different populations based on trabecular bone microarchitecture in longitudinal studies are also presented for the various imaging modalities.

Advanced Imaging Assessment of Bone Quality

Noninvasive and/or nondestructive techniques can provide structural information about bone, beyond simple bone densitometry. While the latter provides important information about osteoporotic fracture risk, many studies indicate that bone mineral density (BMD) only partly explains bone strength. Quantitative assessment of macrostructural characteristics, such as geometry, and microstructural features, such as relative trabecular volume, trabecular spacing, and connectivity, may improve our ability to estimate bone strength. Methods for quantitatively assessing macrostructure include (besides conventional radiographs) dual X ray absorptiometry (DXA) and computed tomography (CT), particularly volumetric quantitative computed tomography (vQCT). Methods for assessing microstructure of trabecular bone noninvasively and/or nondestructively include high-resolution computed tomography (hrCT), microcomputed tomography (micro-CT), high-resolution magnetic resonance (hrMR), and micromagnetic resonance (micro-MR). vQCT, hrCT, and hrMR are generally applicable in vivo; micro-CT and micro-MR are principally applicable in vitro. Despite progress, problems remain. The important balances between spatial resolution and sampling size, or between signal-to-noise and radiation dose or acquisition time, need further consideration, as do the complexity and expense of the methods versus their availability and accessibility. Clinically, the challenges for bone imaging include balancing the advantages of simple bone densitometry versus the more complex architectural features of bone, or the deeper research requirements versus the broader clinical needs. The biological differences between the peripheral appendicular skeleton and the central axial skeleton must be further addressed. Finally, the relative merits of these sophisticated imaging techniques must be weighed with respect to their applications as diagnostic procedures, requiring high accuracy or reliability, versus their monitoring applications, requiring high precision or reproducibility.

[Vertebral trabecular network analysis on CT images]

PURPOSE

Comparing texture analysis, density measurement and visual quantification of trabecular network on spine CT images, to better evaluate bone architecture in osteoporosis.

METHOD AND MATERIALS

Seventeen patients, aged 19 to 84 years, were included. One patient presented osteoporotic fractures. High resolution computed tomographic (HR-CT) images of the third lumbar vertebra were acquired using a Somatom 4 plus CT (Siemens) in a strict axial orientation with FOV of 12 cm and slice thickness of 1 mm. The size of the Region Of Interest was 1,6 cm(2). Three analyses were performed on this ROI: Density (in Hounsfield Unity), texture analysis (run length) and features inspired from bone histomorphometry (Bone Volume/Tissue Volume).

RESULTS

Density measurement, run length methods and BV/TV provided consistent results with regards to age. Indeed density, run length and BV/TV results were lower for older patients with more advanced bone trabeculra alterations.

CONCLUSION

Only BV/TV and run length parameters seemed to show additional information on trabecular network architecture. The contribution of these two measurements to diagnose and classify osteoporosis will be the goal of a clinical study.

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.

Trabecular Bone Structure Obtained From Multislice Spiral Computed Tomography of the Calcaneus Predicts Osteoporotic Vertebral Deformities

PURPOSE

To compare multislice computed tomography (MSCT)-derived parameters of the trabecular bone structure of the calcaneus with bone mineral density (BMD) in their ability to differentiate between donors with and without osteoporotic fractures of the spine and to optimize CT scan protocols.

METHODS

Forty-two postmortem calcanei (81.2 +/- 10 years) were imaged with a 16-detector row MSCT system using 4 different scan protocols varying spatial resolution (12-24 lp/cm) and radiation dose. Structural parameters of trabecular bone were derived from these images, and BMDs of the calcanei were determined using dual x-ray absorptiometry. Vertebral deformities of the spine were radiographically classified using the Spinal Fracture Index. Diagnostic performance in differentiation between donors with and without vertebral fractures was assessed using receiver operating characteristic (ROC) analysis.

RESULTS

There were significant case-control differences for many of the structural parameters measured (P < 0.05). The highest ROC values were found for apparent trabecular thickness using the high-resolution and high-dose protocols. Statistically significant correlations were found between most structure parameters and BMD (up to r = 0.85, P < 0.01).

CONCLUSION

Structural parameters of trabecular bone as obtained from high-resolution MSCT images of the calcaneus can be used to differentiate between donors with and without osteoporotic vertebral fractures, using a high-resolution and high-dose CT protocol.

Femoral neck trabecular microstructure in ovariectomized ewes treated with calcitonin: MRI microscopic evaluation

Ovariectomy induces deterioration of the trabecular structure in the femoral neck of ewes, as depicted by MR microscopic imaging. This structural deterioration is prevented by salmon calcitonin treatment.

INTRODUCTION

This study evaluated the trabecular (Tb) microarchitecture of an ovariectomy (OVX)-induced osteoporotic model in ewes and determined the effects of salmon calcitonin (sCT), an osteoclast inhibitor, on the Tb structure. This is the first report of OVX-induced changes in the Tb structure in the femoral neck in the ewes and effect of sCT on the microarchitecture.

MATERIALS AND METHODS

Ewes (5-8 years old, n = 28) were equally allocated into sham (Sham), OVX injected with vehicle, or OVX injected with sCT at 50 or 100 IU, three injections per week. They were killed 6 months after OVX. The femoral neck was examined with an MR imager at 9.4 T in axial, coronal, and sagittal planes. An internal calibration procedure as a means of standardizing image analysis was used to adjust the segmentation threshold. Data from all three planes were averaged.

RESULTS AND CONCLUSIONS

Compared with Sham, OVX induced significant changes (p < 0.0125) in the MRI-derived femoral neck Tb structure: Tb bone volume fraction (BV/TV), -18%; Tb number, -20%; Tb separation, +23%; number of free ends, +28%; number of nodes, -39%; number of Tb branches, -23%; mean length of Tb branches, -19%. Compared with OVX, treatment of sCT at 100 IU significantly improved all the Tb structural parameters to the Sham level (p < 0.0001 approximately p = 0.0281), whereas 50 IU significantly increased the Tb number and the mean length of the Tb branches. BV/TV explained 74% of the variation of compressive stress of the trabecular cylinder cores of the femoral neck. Combining all structural parameters in a multivariate regression analysis significantly improved the explanation to 84%, and adding BMD further improved the predictive ability of the model to 92%. We conclude that OVX induces deterioration of the MRI-derived Tb microstructure in the femoral neck of ewes. sCT treatment prevents OVX-induced changes. The femoral neck microarchitecture significantly correlates with its biomechanical properties. Combining microstructural parameters with BMD further improves the prediction of bone biomechanical properties. The effects of sCT on OVX ewes may help explain reduced fracture risk in postmenopausal osteoporotic women treated with sCT.

Trabecular bone structure of the distal radius, the calcaneus, and the spine: which site predicts fracture status of the spine best?

RATIONALE AND OBJECTIVES

To compare trabecular bone structure measures obtained in magnetic resonance images of the distal radius and the calcaneus as well as computed tomographic images of the spine versus bone mineral density (BMD) of the spine and the calcaneus in the prediction of osteoporotic spine fracture status.

MATERIAL AND METHODS

High-resolution magnetic resonance images of the calcaneus and the distal radius and thin-section computed tomographic images of thoracic and lumbar vertebrae were obtained from 74 cadavers. Structure analysis was performed using parameters analogous to standard histomorphometry. BMD of the spine was determined by using quantitative computed tomography and of the calcaneus by using dual x-ray absorptiometry. Spine radiographs of these cadavers were assessed concerning vertebral deformities.

RESULTS

The diagnostic performance in differentiating fracture and nonfracture subjects was highest for structure parameters in the spine and slightly lower for these parameters in the distal radius and for BMD of the spine.

CONCLUSION

In this study structure parameters in the spine were best suited to predict the osteoporotic fracture status of the spine.

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.

Effects of h-PTH on cancellous bone mass, connectivity, and bone strength in ovariectomized rats with and without sciatic-neurectomy

The purpose of this study was to determine whether h-PTH (1-34) treatment would recover cancellous bone connectivity and bone strength in ovariectomized (OVX) or ovariectomized and sciatic-neurectomized (OVX+NX) rats. Seven-month-old female Wistar rats were treated with h-PTH or vehicle (6.0 microg/kg, six times a week, subcutaneously) for four weeks beginning 4, 8, or 12 weeks after OVX or OVX+NX. These were compared to age-matched baseline and sham-operated groups. Right tibiae were used for bone histomorphometry and node-strut analysis, and left tibiae were used for mechanical testing. The bone formation rates in the OVX and OVX+NX rats treated with h-PTH were significantly higher than those in their baseline controls. h-PTH treatment increased the node numbers and failure energies in the OVX rats, compared to their baseline controls, at all time points. However, in the OVX+NX rats, the effects of h-PTH treatment on the node number and failure energy were observed only at four weeks after surgery, but not at eight weeks or 12 weeks after surgery. These results suggest that the lowest limit, at which trabecular connectivity and bone strength are able to be restored by h-PTH, occurred between four and eight weeks in OVX+NX rats, but not in OVX rats. h-PTH cannot recover trabecular connectivity and bone strength in advanced osteopenia.

Analysis of the influence of image resolution on the discriminating power of trabecular bone architectural parameters

In the study the influence of image resolution on the discriminating power of different parameters quantifying the architecture of trabecular bone is investigated. High-resolution images of 200-microm-thick sections of young and old trabecular bone are collected. Different architectural parameters are evaluated and shown to be statistically different in both groups of images. Then the resolution of the images is artificially degraded to the level comparable with resolution achievable in vivo. It is shown that although the errors of evaluation are quite large for low-resolution images, the statistical difference present in the original data is still observed for parameters that depend only on the global characteristics of trabecular structure or marrow space, or depend only linearly on the number and area of disconnected parts of marrow space or trabecular structure. The parameters that fulfill such conditions are Euler number, mean area of marrow cavities, star volume and trabecular spacing. The statistical difference disappears for parameters that depend on the number and area of disconnected parts of marrow space or trabecular structure in a nonlinear manner, that is, marrow and trabecular disconnection probabilities, two-point distance and second moment of marrow cavities area distribution. It is shown that the mechanism that leads to the losing of discriminating power does not depend crucially on the estimation error or noise level. This mechanism, which involves an interplay between changes of mean values and standard deviations, is the manifestation of partial volume effect amplified with image segmentation.

Effects of combined treatment of insulin and human parathyroid hormone(1-34) on cancellous bone mass and structure in streptozotocin-induced diabetic rats

The purpose of this study was to test the hypothesis that combined treatment with insulin and human parathyroid hormone (hPTH) is more effective than treatment with insulin or hPTH alone in improving cancellous bone mass, connectivity, and strength in insulin-dependent diabetic rats. Diabetes was induced by intraperitoneal injection of streptozotocin (STZ) in 7-month-old female Wistar rats. The diabetic rats received insulin, hPTH, insulin and hPTH, or hPTH vehicle for 4 weeks, starting 8 weeks after STZ injection. They were compared with baseline controls and normal controls that received STZ alone and STZ vehicle alone, respectively. The rats' proximal right tibias were processed to serve as undecalcified Villanueva-stained bone sections for histomorphometry. Changes in trabecular connectivity were determined through node-strut analysis. The decreased cancellous bone volume (BV/TV) and bone formation in diabetic rats improved in all the drug-treated groups compared with baseline controls. Furthermore, recovery of BV/TV was greater in rats that received the combination of insulin and hPTH than in those that received insulin or hPTH alone. In node-strut analysis, the node-related parameter (N.Nd/TV) significantly increased in rats that received the combination of insulin and hPTH, but did not increase in those that received insulin or hPTH alone. In addition to these results, the combination treatment significantly increased bone mineral density of the femur and bone strength in the femoral metaphysis compared with treatment with insulin or hPTH alone. These results indicate that the doses of insulin and hPTH employed in the combination treatment were more effective in improving not only bone mass but also trabecular connectivity and bone strength than treatment with insulin or hPTH alone in insulin-dependent diabetic rats.

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.

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.

An investigation of thoracic and lumbar cancellous vertebral architecture using power-spectral analysis of plain radiographs

The internal architecture of the vertebral bodies spanning the levels T1 to L5 in seven male columns was studied using mammographic-resolution radiographs of 2.5-mm-thick planar parasagittal slices. The overlapping radiographic shadows of vertebral trabeculae combined in the image to form a series of 'elements', broadly representative of the cancellous structure. The orientations and sizes of these elements were analysed by applying the Fast Fourier transform (FFT) to the digitized radiographic images. Elements aligned in the 'vertical' orientation, along the long axis of the column, were the most prominent for all vertebral levels. The relative prominence of horizontal to vertical elements was generally constant along the column below T5. In contrast, the relative prominence of oblique to vertical elements declined in the cranio-caudal direction, particularly in individuals aged > or = 60 years. The ratio of 'large' (x > 0.3 mm) to 'small' (0.15 mm < or = x < or = 0.3 mm) elements was unchanged cranio-caudally in specimens < 60 years. However, in individuals > or = 60 years, large elements increased in relative prominence in the caudal direction. These results suggest that a basic orthogonal pattern of trabeculae is found along the male human spine, regardless of differences in vertebral body size. Power-spectral analysis is shown to yield information summarizing the predominant orientations and sizes of radiographically rendered architectural elements of vertebral cancellous bone, to define the effects of ageing on architecture, and to identify broad structural differences between vertebral levels in the adult male spine.

Surface area overestimation within three-dimensional digital images and its consequence for skeletal dosimetry

The most recent methods for trabecular bone dosimetry are based on Monte Carlo transport simulations within three-dimensional (3D) images of real human bone samples. Nuclear magnetic resonance and micro-computed tomography have been commonly used as imaging tools for studying trabecular microstructure. In order to evaluate the accuracy of these techniques for radiation dosimetry, a previous study was conducted that showed an overestimate in the absorbed fraction of energy for low-energy electrons emitted within the marrow space and irradiating the bone trabeculae. This problem was found to be related to an overestimate of the surface area of the true bone-marrow interface within the 3D digital images, and was identified as the surface-area effect. The goal of the present study is to better understand how this surface-area effect occurs in the case of single spheres representing individual marrow cavities within trabecular bone. First, a theoretical study was conducted which showed that voxelization of the spherical marrow cavity results in a 50% overestimation of the spherical surface area. Moreover, this overestimation cannot be reduced through a reduction in the voxel size (e.g., improved image resolution). Second, a series of single-sphere marrow cavity models was created with electron sources simulated within the sphere (marrow source) and outside the sphere (bone trabeculae source). The series of single-sphere models was then voxelized to represent 3D digital images of varying resolution. Transport calculations were made for both marrow and bone electron sources within these simulated images. The study showed that for low-energy electrons (<100 keV), the 50% overestimate of the bone-marrow interface surface area can lead to a 50% overestimate of the cross-absorbed fraction. It is concluded that while improved resolution will not reduce the surface area effects found within 3D image-based transport models, a tenfold improvement in current image resolution would compensate the associated errors in cross-region absorbed fractions for low-energy electron sources. Alternatively, other methods of defining the bone-marrow interface, such as with a polygonal isosurface, would provide improvements in dosimetry without the need for drastic reductions in image voxel size.

Segmentation of cancellous bone from high-resolution computed tomography images: influence on trabecular bone measurements

The quantification of cancellous bone network from computed tomography (CT) images requires a segmentation step which is crucial and difficult because of the partial volume effect in CT images. In this paper, we present and evaluate a new approach for segmenting cancellous bone network from high-resolution CT (HRCT) slices. The idea is first to detect a skeleton from the crest lines of the structure and then to thicken it to extract the whole bone structure by satisfying local neighborhood constraints. The segmentation requires the adjustment of relative and not absolute parameters like most methods. We quantified the influence of these parameters on architectural measurements. Results were first validated by using a physical phantom and then examined on a series of 12 HRCT images of human lumbar vertebra of different ages. We demonstrated that the choice of segmentation parameters yielded important variability on architectural measurements (up to 20%), but less variability than a more commonly used approach. This stresses the importance of settle on the segmentation parameters for once, which is possible with the proposed method.

Fractal analysis of lumbar vertebral cancellous bone architecture

Osteoporosis is characterized by bone mineral density (BMD) decreasing and spongy bone rearrangement with consequent loss of elasticity and increased bone fragility. Quantitative computed tomography (QCT) quantifies bone mineral content but does not describe spongy architecture. Analysis of trabecular pattern may provide additional information to evaluate osteoporosis. The aim of this study was to determine whether the fractal analysis of the microradiography of lumbar vertebrae provides a reliable assessment of bone texture, which correlates with the BMD. The lumbar segment of the spine was removed from 22 cadavers with no history of back pain and examined with standard x-ray, traditional tomography, and quantitative computed tomography to measure BMD. The fractal dimension, which quantifies the image fractal complexity, was calculated on microradiographs of axial sections of the fourth lumbar vertebra to determine its characteristic spongy network. The relationship between the values of the BMD and those of the fractal dimension was evaluated by linear regression and a statistically significant correlation (R = 0.96) was found. These findings suggest that the application of fractal analysis to radiological analyses can provide valuable information on the trabecular pattern of vertebrae. Thus, fractal dimensions of trabecular bone structure should be considered as a supplement to BMD evaluation in the assessment of osteoporosis.

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.

A comparison of the texture of computed tomography and projection radiography images of vertebral trabecular bone using fractal signature and lacunarity

The structural integrity of trabecular bone is an important factor characterizing the biomechanical strength of the vertebra, and is determined by the connectivity of the bone network and the trabeculation pattern. These can be assessed using texture measures such as the fractal signature and lacunarity from a high resolution projection radiograph. Using central sections of lumbar vertebrae we compared the results obtained from high-resolution transverse projection images with those obtained from spatially registered low-resolution images from a conventional clinical CT scanner to determine whether clinical CT data can provide useful structural information. Provided the power spectra of the CT images are corrected for image system blurring, the resulting fractal signature is similar for both modalities. Although the CT images are blurred relative to the projection images, with a consequent reduction in lacunarity, the estimated trabecular separation obtained from the lacunarity plots is similar for both modalities. This suggests that these texture measures contain essential information on trabecular microarchitecture, which is present even in low resolution CT images. Such quantitative texture measurements from CT or MRI images are potentially useful in monitoring bone strength and predicting future fracture risk.

Inhomogeneity of human vertebral cancellous bone: systematic density and structure patterns inside the vertebral body

In the spine, cancellous bone quality is usually assessed for the whole vertebral body in a transverse central slice. Correct identification and assessment of the weakest parts of the cancellous bone may lead to better prediction of fracture risk. The density and structural parameters were systematically investigated inside the thoracic (T-9), thoracolumbar (T12-L1), and lumbar (L-4) vertebral bodies of nine subjects. On both sides of the median sagittal plane, anterior and posterior 8.2 mm vertical cores were harvested in the thoracic vertebra. In the thoracolumbar and lumbar vertebrae, external samples were also cored. Peripheral quantitative computed tomographic (pQCT) density analysis of the 136 cores was performed at four different levels, from the lower to the upper endplate. The relatively thin slice thickness (300 microm) and small pixel size (70 microm x 70 microm) was considered sufficient to investigate the structural parameters on the four transverse slices and in the sagittal and coronal planes (total of 816 images). Using a constant threshold a binary image was generated and the morphometric data were extracted. The binary image was further skeletonized and classical strut analysis was performed. Cancellous bone density was 20% higher in the posterior cores than in the anterior and external cores. Moreover, clear vertical inhomogeneity was noted because the lowest half of the vertebral body presented lower density than the upper half (differences ranging from 25% to 15%). All structural parameters were strongly dependent on the location of the measurement. Structural differences between anterior, posterior, and external areas were mild and followed the density patterns. On the other hand, vertical inhomogeneity of the structural parameters was important. For example, in the thoracolumbar and lumbar vertebrae, the numbers of nodes or node-to-node struts were almost twofold higher in the inferior half than in the superior half (p < 0.01), whereas trabecular thickness and number of free-ends presented a center/close-to-endplate structural pattern, with central trabeculae being 15% thicker (p < 0.05) and presenting 30% fewer free-ends (p < 0.01) than the close-to-endplate ones. Variability of density and structural parameters was high and a substantial part of this variability could be explained by the place inside the vertebral body where the measurement was made. The weak part was not in the center of the body but in its upper half where the lower density did not seem to be compensated by a higher structural architecture. Further clinical investigation could enhance fracture prediction by tracking and focusing on the weakest part of the vertebral body.

Texture analysis of X-ray radiographs is a more reliable descriptor of bone loss than mineral content in a rat model of localized disuse induced by the Clostridium botulinum toxin

Botulism is a generalized paralyzing disease caused by the toxin of Clostridium botulinum (BTX). The toxin acts 3-4 days after injection by blocking the release of acetylcholine to the muscle. Six Wistar rats received a 2-U injection of BTX in the right quadriceps. Six rats were similarly injected with saline and were used as control. Paralysis of the quadriceps was obtained 4-5 days after the injection. Animals were killed 4 weeks after the BTX injection. The bone mineral content (BMC) was measured by dual-energy X-ray absorptiometry (DXA) on the femur and tibia. No side-to-side difference was observed for BMC on the whole tibia and femur in the BTX group. When subregions were selected in the bones, a significant decrease in BMC was obtained on the proximal tibia (-17.4 +/- 2.5%, p < 0.02). No significant difference could be observed on the proximal or distal femur, nor on the diaphyseal shafts. Numeric X-rays were done and a region of interest was transferred to an image analyzer. The texture of the trabecular bone was analyzed by the run length and fractal methods (skyscrapers and blanket). Significant differences were obtained on the proximal tibia for all methods except with fractal skyscrapers. On the distal femur, significant differences were obtained with the run length method, and the skyscrapers and the blanket method in the vertical direction. No differences were obtained with any method on the tibia and femur from control animals. Bone is a highly anisotropic material and its architecture at the microscopic level is conditioned by strains. The trabecular pattern differs in the proximal tibia than in the distal femural. Depending on the trabecular anisotropy, the algorithms can be more or less pertinent. BTX induced a significant bone loss on the bony subparts that are directly influenced by disuse. Texture analysis of X-ray images can reveal differences that were not evidenced by naked eyes. However, a combination of several methods appears necessary to appreciate the bone loss.

Predictive value of bone mineral density and morphology determined by peripheral quantitative computed tomography for cancellous bone strength of the proximal femur

Peripheral quantitative computed tomography (pQCT) is an established diagnostic method for assessment of bone mineral density in the diagnosis of osteoporosis. However, the capacity of structural parameters of cancellous bone measured by high-resolution computed tomography remains to be explored. In 33 patients, bone mineral density (BMD) of the proximal femur was measured in vitro by pQCT using cylindrical biopsies from the intertrochanteric region harvested before the implantation of an artificial hip joint. By digital image analysis of CT scans, parameters derived from histomorphometry describing the microarchitecture of cancellous bone were measured. The biopsies were also loaded to failure by an uniaxial compression test to determine the biomechanical parameters, Young's modulus, strength, and maximum energy absorption (E(max)). Strong correlations were found for BMD vs. mechanical parameters (r = 0.73 for Young's modulus, r = 0.82 for strength, and r = 0.79 for E(max); p < 0.001, n = 29). The morphological parameters, bone volume per trabecular volume (BV/TV), apparent trabecular thickness (app.Tb.Th), apparent trabecular separation (app.Tb.Sp), and trabecular number (Tb.N), correlated significantly with all mechanical parameters. The combination of morphological parameters with BMD in a multivariate regression model led to an overall, but only moderate, increase in R(2) in all cases. Our data confirm the high predictive value of BMD for the mechanical competence of cancellous bone of the intertrochanteric region. However, quantification of cancellous bone structure by image analysis of CT scans may provide additional qualitative information for the analysis of bone strength.

Changes in calcaneal trabecular bone structure after heart transplantation: an MR imaging study

PURPOSE

To use high-spatial-resolution magnetic resonance (MR) imaging to analyze the trabecular bone structure of the calcaneus in patients before and after heart transplantation and to compare this technique with bone mineral density (BMD) measurement in predicting therapy-induced bone loss and vertebral fracture status.

MATERIALS AND METHODS

High-spatial-resolution 1.5-T MR imaging of the calcaneus was performed in 40 men 11-120 months after heart transplantation, in 11 men before heart transplantation, and in 10 age-matched male volunteers. Sagittal and transverse T1-weighted spin-echo images with a voxel size of 0.195 x 0.195 x 1.000 mm were obtained, and structure measurements analogous to bone histomorphometric values were calculated. In addition, the BMD of the lumbar spine was determined in the transplant recipients pre- and postoperatively by using quantitative computed tomography, and vertebral fracture status was assessed.

RESULTS

Significant differences in structure and BMD measurements were found between patients before and after heart transplantation (P <. 05). In 17 (42%) of 40 transplant recipients, vertebral fractures were found. Although structure measurements were significantly different between patients with and those without fractures (P <.05), BMDs were not. Correlations between time after transplantation and some structure measurements were moderately significant (P <. 05), but such correlations with BMD measurements were not.

CONCLUSION

MR imaging-derived structure measurements in the calcaneus are useful for monitoring bone changes after heart transplantation and assessing vertebral fracture status.

In vivo assessment of trabecular bone structure using fractal analysis of distal radius radiographs

Our purpose in this study was (i) to measure trabecular bone structure using fractal analysis of distal radius radiographs in subjects with and without osteoporotic hip fractures, and (ii) to compare these measures with bone mineral density (BMD) as well as with measures of trabecular bone structure derived from high resolution magnetic resonance (MR) images. Distal radius radiographs were obtained using semi-industrial films (55 kVp, 400 mAs) in 30 postmenopausal patients, who had suffered osteoporotic hip fractures (74.8+/-8.2 years) in the last 24 months and 27 postmenopausal age-matched (74.6+/-6.6 yr) normal volunteers. Radiographs were digitized at 50 microm. A Fourier power spectrum-based fractal dimension (FD) characterizing the trabecular pattern was measured in a region of interest proximal to the joint line. The fractal dimension was calculated over two spatial frequency (f) ranges: FD1 was calculated over 0.5<log(f)<l.0, FD2 over the higher range 1.0<log(f)<1.5. Trabecular BMD in the radius was obtained using peripheral quantitative computed tomography (pQCT) (Stratec GmbH, Germany). In addition BMD of the proximal femur was determined using dual x-ray absorptiometry (DXA) (QDR 2000, Hologic, MA). In a subset of patients (16 controls and 18 with hip fractures), high resolution MR imaging of the distal radius (spatial resolution of 156 x 156 x 500 microm) was used to obtain measures analogous to bone histomorphometry. There were significant differences (p<0.05) between the fracture and nonfracture groups in the total femur BMD (13%), trabecular BMD in the distal radius (4%), and the fractal dimension in the radiographs (FD2) (3%). The correlations between FD2 and the total femur BMD as well as trabecular bone BMD in the distal radius were -0.48 (p<0.006) and -0.22 (p<0.33); respectively; FD1 increased with BMD and showed lower correlations. FD2 showed good correlations with App. Tb.N (-0.71) and App. Tb.Sp (0.69) (p<0.01), moderate correlation with App BV/TV (-0.53) (p<0.05), and no significant correlation with App. Tb.Th. The correlations between structural measures and FD1 showed the inverse trend and were typically lower. The odds ratios for a hip fracture were 2.44 for total femur BMD, 1.5 for trabecular BMD (radius), and 1.5 for FD2, respectively. In summary, the fractal measures derived from radiographs of the radius show differences between subjects with and without hip fractures, the predictive power of measures in the distal radius are comparable to radial trabecular BMD but lower than that of total hip BMD.

Does high-resolution computed tomography image analysis of the distal radius provide information independent of bone mass?

This study aimed to investigate the usefulness of computed tomography (CT) image analysis of the distal radius for comparing two groups of postmenopausal women matched for age and bone mineral density at both the lumbar spine and femoral neck. The first one consisted of 16 women with at least one vertebral fracture and the second consisted also of 16 women without disease affecting bone mass or bone metabolism. Eight slices were selected in each patient: four consecutive coronal slices and four consecutive axial slices. Bone texture analysis was performed using structural methods leading to the measurement of 24 features. Most of the structural variables derived from histomorphometric parameters and were measured after segmentation from a binary or a skeletonized image. Nine variables were significantly different between the two groups on axial slices: valley number, valley surface area, apparent bone volume/tissue volume (BV/TV), apparent trabecular separation, apparent trabecular number, trabecular bone pattern factor, trabecular skeletal length, node count, and node-to-node strut count. Also four variables were significantly different between osteoporotic women and controls on coronal slices: apparent BV/TV, trabecular partition, node-to-node strut count, and terminus-to-terminus strut count. In conclusion this study suggests that bone texture analysis could yield additional data on bone mass for explaining bone strength and therefore could be used for improving the prediction of fracture risk.

In vivo comparison between computed tomography and magnetic resonance image analysis of the distal radius in the assessment of osteoporosis

In a prospective case-control cross-sectional study, we investigated the usefulness of both computed tomography (CT) and magnetic resonance (MR) image analysis of the distal radius for distinguishing a small sample of fractured osteoporotic women from control women regardless of bone mineral density. The study population included 12 subjects who were divided into two groups according to their bone status. The first group consisted of six women with at least one vertebral fracture occurring in the absence of high-energy trauma, and the second group consisted of six women without disease affecting bone mass or bone metabolism. Cross-sectional and coronal slices were obtained from both CT and MR systems. For CT images, the slice thickness was 1000 jim and the plane resolution was approx 200 jim x 200 jim. MR images were obtained from a 1.5-T imager with a two-dimensional spin-echo Ti-weighted sequence leading to a slice thickness of 2000 jim and a plane resolution of 195 jim x 195 jim. Bone texture analysis was performed using fractal and structural methods leading to the measurement of 23 features. Most of the structural variables derived from histomorphometric parameters and were measured after segmentation from a binary or a skeletonized image. Bone densitometry was measured by dual-energy X-ray absorptiometry both at the lumbar spine and the nondominant femoral neck. On cross-sectional slices, 12 variables, mainly obtained from structural analysis, were significantly different between the two groups for CT images (p < 0.05) against two variables only for MR images (p < 0.05). The number of variables statistically different between the two groups was significantly higher for CT images compared with MR images (p = 0.003). In the same way, odds ratios for fracture per 1 control group standard deviation decrease were significant for 10 variables on CT images, whereas, in contrast, none of the variables measured on MRI images led to significant odds ratios. The results obtained for the two methods on coronal slices were poorer without a difference between either CT or MR images in terms of discrimination between fracture and nonfracture subjects. In conclusion, this study suggests that bone texture analysis obtained from CT compared with MRI offers a best discrimination between controls and osteoporotic patients, probably the result to the spatial resolution. which is better for CT than for MR images.

Slice simulation from a model of the parenchymous vascularization to evaluate texture features: work in progress

RATIONALE AND OBJECTIVES

To demonstrate the usefulness of a model of the parenchymous vascularization to evaluate texture analysis methods.

METHODS

Slices with thickness varying from 1 to 4 mm were reformatted from a 3D vascular model corresponding to either normal tissue perfusion or local hypervascularization. Parameters of statistical methods were measured on 16128x128 regions of interest, and mean values and standard deviation were calculated. For each parameter, the performances (discrimination power and stability) were evaluated.

RESULTS

Among 11 calculated statistical parameters, three (homogeneity, entropy, mean of gradients) were found to have a good discriminating power to differentiate normal perfusion from hypervascularization, but only the gradient mean was found to have a good stability with respect to the thickness. Five parameters (run percentage, run length distribution, long run emphasis, contrast, and gray level distribution) were found to have intermediate results. In the remaining three, curtosis and correlation was found to have little discrimination power, skewness none.

CONCLUSION

This 3D vascular model, which allows the generation of various examples of vascular textures, is a powerful tool to assess the performance of texture analysis methods. This improves our knowledge of the methods and should contribute to their a priori choice when designing clinical studies.

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.

Morphometric texture analysis of spinal trabecular bone structure assessed using orthogonal radiographic projections

The measurement of bone microstructure as well as bone mineral density may improve the estimation of bone strength. Cubic specimens (N = 26, 12 mm X 12 mm X 12 mm) of human cadaver vertebrae were cut along three orthogonal anatomic orientations, i.e., superior-inferior (SI), medial-lateral (ML), and anterior-posterior (AP). Contact radiographs of the bone cubes along all three orientations were obtained and then digitized by a laser scanner with pixel size of 50 microns x 50 microns. The specimens were tested in compression along the 3 orthogonal orientations and the Young's modulus (YM) was calculated for each direction. Quantitative computed tomography (QCT) was used to obtain a measure of trabecular bone mineral density (BMD). Global gray level thresholding and local thresholding algorithms were used to extract the trabecular bone network. Apparent trabecular bone fraction (ABV/TV), mean intercept length (I.TH), mean intercept separation (I.SP), and number of nodes (N.ND) were measured from the extracted trabecular network. Fractal dimension (Fr.D) of the trabecular bone texture was also measured. Paired t-tests showed that the mean values of each texture parameter (except ABV/TV) and of YM along the SI direction were significantly different (p < 0.05) from those along the ML and AP direction. However, the mean values along the ML and AP directions were not significantly different. Multivariate regression of YM as a function of the texture parameters and BMD showed that without adjusting for the effect of BMD, YM was significantly explained by all the texture parameters (R2 = 0.2-0.6). When BMD was included in the regression, although the variations in YM of ML, AP, and SI orientations could be explained by BMD alone, some of the texture parameters did improve the overall prediction of the biomechanical properties, while, some parameters such as ABV/TV and Fr.D in the ML orientation showed a more significant overall effect in explaining mechanical strength than did BMD. In conclusion, trabecular texture parameters correlated significantly with BMD and YM. Trabecular texture parameters from projectional radiographs reflect the anisotropy of trabecular structure. Quantitative radiographic assessment of trabecular structure using fine-detail radiography can potentially improve the estimation of bone strength.

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 female (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 micron and in plane resolution of 195 x 195 micron. 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.

Instrumental diagnosis of osteoporosis

Considerable progress in the development of methods for assessing the skeleton now makes it possible to detect osteoporosis non-invasively and early. There is a variety of techniques available at present: single-photon (SPA) and single X-ray absorptiometry (SXA), dual-photon (DPA) and dual X-ray absorptiometry (DXA), quantitative computed tomography (QCT), radiographic absorptiometry (RA), and quantitative ultrasound (QUS), and their development has certainly been driven by the need to overcome the inherent shortcomings of plain radiography for this purpose. Both SPA and SXA methods make a quantitative assessment of the bone mineral content (BMC) or density (BMD) at peripheral sites of the skeleton possible. Single energy measurements are not possible at sites with variable soft tissue thickness and composition, i.e., the axial skeleton. For these purposes, DPA and DXA techniques were introduced. The main advantages of an X-ray system over a radionuclide system are shortened examination time, greater accuracy and precision limited to higher resolution, and removal of errors due to source decay correction. Low radiation dose, availability, capacity to evaluate multiple sites, and ease of use have made DXA the most widely used technique for measuring bone mineral density. QCT can determine the true volumetric density of trabecular or cortical bone in three dimensions at any skeletal site. Recently developed new computer-assisted methods have improved RA precision, thus providing a simple and inexpensive technique for screening of bone mineral status of large populations. QUS was reported to provide information regarding the structural characteristics of bone, which may be relevant to the appearance of osteoporotic fractures; indeed, some studies suggest a relationship between QUS and bone strength beyond that which can be explained by BMD. Recent experimental studies suggested that magnetic resonance might also constitute a promising tool for assessing osteoporosis.