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

Menopausal hormone therapy for the management of osteoporosis

Postmenopausal osteoporosis is a frequent clinical condition which affects nearly 1 in 3 women. Estrogen deficiency leads to rapid bone loss which is maximal within the first 2-3 years after the menopause transition and can be prevented by menopause hormone therapy (MHT). Not only, MHT prevents bone loss and the degradation of the bone microarchitecture but it significantly reduces the risk of fracture at all bone sites by 20-40%. It is the only anti-osteoporotic therapy that has a proven efficacy regardless of basal level of risk, even in low-risk women for fracture. Following the publication of the WHI results, use of MHT has considerably declined due to safety concerns which raise the question as to whether it might still be used in the prevention of osteoporosis. Over the last years, subsequent re-analyses of the WHI and further trials have challenged the initial conclusions of the WHI. It is now clearer that the individual benefit-risk balance of MHT is dependent on the individual risk profile in each woman as well as whether estrogen is opposed or unopposed, the type of estrogens and progestogens or doses and routes of administration. It must be also reminded that to date osteoporosis is a chronic disease that cannot be cured. The choice of the 1st treatment option should thus always be made in the context of a more comprehensive long-term strategy. This is particular true in early postmenopausal women found to be at low/moderate risk of fragility fracture over the first 10 years after menopause but who may have a much greater lifetime risk. In the absence of contraindication, use of MHT should be considered as a 1st option for the maintenance of bone health in those women where specific bone active medications are not warranted. Subsequent reassessment of the individual benefit-risk balance of MHT is thereafter recommended, with the possibility of switching to another osteoporosis treatment if the balance is not considered as favourable as at the beginning of the menopause for women still at high risk of fracture.

[Menopause, menopause hormone therapy and osteoporosis. Postmenopausal women management: CNGOF and GEMVi clinical practice guidelines]

Postmenopausal osteoporosis is a frequent clinical condition, which affects nearly 1 in 3 women. Estrogen deficiency leads to rapid bone loss, which is maximal within the first years after the menopause transition and can be prevented by menopause hormone therapy (MHT). Assessment of the individual risk of osteoporosis is primarily based on the measurement of bone mineral density (BMD) at the spine and femur by DXA. Clinical risk factors (CRFs) for fractures taken either alone or in combination in the FRAX score were shown not to reliably predict fractures and/or osteoporosis (as defined by a T-score<-2.5) in early postmenopausal women. If DXA measurement is indicated in all women with CRFs for fractures, it can be proposed on a case-by-case basis, when knowledge of BMD is likely to condition the management of women at the beginning of menopause, particularly the benefit-risk balance of MHT. MHT prevents both bone loss and degradation of the bone microarchitecture in early menopause. It significantly reduces the risk of fracture at all bone sites by 20 to 40% regardless of basal level of risk with an estrogen-dependent dose-effect. Given the inter-individual variability in bone response, individual monitoring of the bone effect of MHT is warranted when prescribed for the prevention of osteoporosis. This monitoring is based on repeated measurement of lumbar and femoral BMD (on the same DXA measurement system) after 2years of MHT, the response criterion being no significant bone loss. Discontinuation of treatment is associated with a resumption of transient bone loss although there is a large variability in the rate of bone loss among women. Basically, there is a return to the level of fracture risk comparable to that of in untreated woman of the same age within 2 to 5years. Therefore, when MHT is prescribed for the prevention of osteoporosis in women with an increased risk at the beginning of menopause, measurement of BMD is recommended when MHT is stopped in order to consider further management of the risk of fracture whenever necessary (with possibly another anti-osteoporotic treatment).

Survey of MRI Usefulness for the Clinical Assessment of Bone Microstructure

Bone microarchitecture has been shown to provide useful information regarding the evaluation of skeleton quality with an added value to areal bone mineral density, which can be used for the diagnosis of several bone diseases. Bone mineral density estimated from dual-energy X-ray absorptiometry (DXA) has shown to be a limited tool to identify patients' risk stratification and therapy delivery. Magnetic resonance imaging (MRI) has been proposed as another technique to assess bone quality and fracture risk by evaluating the bone structure and microarchitecture. To date, MRI is the only completely non-invasive and non-ionizing imaging modality that can assess both cortical and trabecular bone in vivo. In this review article, we reported a survey regarding the clinically relevant information MRI could provide for the assessment of the inner trabecular morphology of different bone segments. The last section will be devoted to the upcoming MRI applications (MR spectroscopy and chemical shift encoding MRI, solid state MRI and quantitative susceptibility mapping), which could provide additional biomarkers for the assessment of bone microarchitecture.

Deep Learning Based High-Resolution Reconstruction of Trabecular Bone Microstructures from Low-Resolution CT Scans using GAN-CIRCLE

Osteoporosis is a common age-related disease characterized by reduced bone density and increased fracture-risk. Microstructural quality of trabecular bone (Tb), commonly found at axial skeletal sites and at the end of long bones, is an important determinant of bone-strength and fracture-risk. High-resolution emerging CT scanners enable in vivo measurement of Tb microstructures at peripheral sites. However, resolution-dependence of microstructural measures and wide resolution-discrepancies among various CT scanners together with rapid upgrades in technology warrant data harmonization in CT-based cross-sectional and longitudinal bone studies. This paper presents a deep learning-based method for high-resolution reconstruction of Tb microstructures from low-resolution CT scans using GAN-CIRCLE. A network was developed and evaluated using post-registered ankle CT scans of nineteen volunteers on both low- and high-resolution CT scanners. 9,000 matching pairs of low- and high-resolution patches of size 64×64 were randomly harvested from ten volunteers for training and validation. Another 5,000 matching pairs of patches from nine other volunteers were used for evaluation. Quantitative comparison shows that predicted high-resolution scans have significantly improved structural similarity index (p < 0.01) with true high-resolution scans as compared to the same metric for low-resolution data. Different Tb microstructural measures such as thickness, spacing, and network area density are also computed from low- and predicted high-resolution images, and compared with the values derived from true high-resolution scans. Thickness and network area measures from predicted images showed higher agreement with true high-resolution CT (CCC = [0.95, 0.91]) derived values than the same measures from low-resolution images (CCC = [0.72, 0.88]).

CT-Based Characterization of Transverse and Longitudinal Trabeculae and Its Applications

Osteoporosis is a common age-related disease characterized by reduced bone mineral density (BMD), micro-structural deterioration, and enhanced fracture-risk. Although, BMD is clinically used to define osteoporosis, there are compelling evidences that bone micro-structural properties are strong determinants of bone strength and fracture-risk. Reliable measures of effective trabecular bone (Tb) micro-structural features are of paramount clinical significance. Tb consists of transverse and longitudinal micro-structures, and there is a hypothesis that transverse trabeculae improve bone strength by arresting buckling of longitudinal trabeculae. In this paper, we present an emerging clinical CT-based new method for characterizing transverse and longitudinal trabeculae, validate the method, and examine its application in human studies. Specifically, we examine repeat CT scan reproducibility, and evaluate the relationships of these measures with gender and body size using human CT data from the Iowa Bone Development Study (IBDS) (n = 99; 49 female). Based on a cadaveric ankle study (n = 12), both transverse and longitudinal Tb measures are found reproducible (ICC > 0.94). It was observed in the IBDS human data that males have significantly higher trabecular bone measures than females for both inner (p < 0.05) and outer (p < 0.01) regions of interest (ROIs). For weight, Spearman correlations ranged 0.43-0.48 for inner ROI measures and 0.50-0.52 for outer ROI measures for females versus 0.30-0.34 and 0.23-0.25 for males. Correlation with height was lower (0.36-0.39), but still mostly significant for females. No association of trabecular measures with height was found for males.

Segmentation of the Proximal Femur from MR Images using Deep Convolutional Neural Networks

Magnetic resonance imaging (MRI) has been proposed as a complimentary method to measure bone quality and assess fracture risk. However, manual segmentation of MR images of bone is time-consuming, limiting the use of MRI measurements in the clinical practice. The purpose of this paper is to present an automatic proximal femur segmentation method that is based on deep convolutional neural networks (CNNs). This study had institutional review board approval and written informed consent was obtained from all subjects. A dataset of volumetric structural MR images of the proximal femur from 86 subjects were manually-segmented by an expert. We performed experiments by training two different CNN architectures with multiple number of initial feature maps, layers and dilation rates, and tested their segmentation performance against the gold standard of manual segmentations using four-fold cross-validation. Automatic segmentation of the proximal femur using CNNs achieved a high dice similarity score of 0.95 ± 0.02 with precision = 0.95 ± 0.02, and recall = 0.95 ± 0.03. The high segmentation accuracy provided by CNNs has the potential to help bring the use of structural MRI measurements of bone quality into clinical practice for management of osteoporosis.

Age-related mechanical strength evolution of trabecular bone under fatigue damage for both genders: Fracture risk evaluation

Bone tissue is a living composite material, providing mechanical and homeostatic functions, and able to constantly adapt its microstructure to changes in long term loading. This adaptation is conducted by a physiological process, known as "bone remodeling". This latter is manifested by interactions between osteoclasts and osteoblasts, and can be influenced by many local factors, via effects on bone cell differentiation and proliferation. In the current work, age and gender effects on damage rate evolution, throughout life, have been investigated using a mechanobiological finite element modeling. To achieve the aim, a mathematical model has been developed, coupling both cell activities and mechanical behavior of trabecular bone, under cyclic loadings. A series of computational simulations (ABAQUS/UMAT) has been performed on a 3D human proximal femur, allowing to investigate the effects of mechanical and biological parameters on mechanical strength of trabecular bone, in order to evaluate the fracture risk resulting from fatigue damage. The obtained results revealed that mechanical stimulus amplitude affects bone resorption and formation rates, and indicated that age and gender are major factors in bone response to the applied loadings.

Effect of the different implant protocols on peri-implant bone: a 6-month prospective study in beagle dogs

OBJECTIVE

To compare the effects of different implant placement and loading protocols on the marginal bone loss (MBL) in beagles by intraoral radiography.

METHODS AND MATERIALS

61 dental implants were inserted on 9 beagle dogs at bilateral lower posteriors according to 8 different protocols: immediate implant placement and immediate loading for 3 months (IIP + IL3) or 6 months (IIP + IL6) and unloading (IIP + UL), immediate implant placement and delayed loading for 3 months (IIP + DL3) or 6 months (IIP + DL6), delayed implant placement and immediate loading for 3 months (DIP + IL3) or delayed loading for 3 months (DIP + DL3) and unloading (DIP + UL). Intraoral radiography was performed to analyze the MBL during each surgery, before and after the implant placement and at 3-month intervals after the procedure.

RESULTS

In total, 57 samples were included. There was less MBL (p<0.05) in the IIP + IL3 group (1.22 ± 0.63 mm) compared to the DIP + IL3 group (1.89 ± 0.9 mm). The longer the loading time, the more bone loss appeared in the IIP + IL group; however, the results were reversed in the IIP + DL group. The MBL during the latter 3-month period was dramatically decreased compared to the former 3-month period in the IIP + DL3 group (p<0.05).

CONCLUSIONS

The IIP + IL group seems superior to the DL protocol and the MBL changed significantly during the first three months and thereafter became stable.

MRI assessment of bone structure and microarchitecture

Osteoporosis is a disease of weak bone and increased fracture risk caused by low bone mass and microarchitectural deterioration of bone tissue. The standard-of-care test used to diagnose osteoporosis, dual-energy x-ray absorptiometry (DXA) estimation of areal bone mineral density (BMD), has limitations as a tool to identify patients at risk for fracture and as a tool to monitor therapy response. Magnetic resonance imaging (MRI) assessment of bone structure and microarchitecture has been proposed as another method to assess bone quality and fracture risk in vivo. MRI is advantageous because it is noninvasive, does not require ionizing radiation, and can evaluate both cortical and trabecular bone. In this review article, we summarize and discuss research progress on MRI of bone structure and microarchitecture over the last decade, focusing on in vivo translational studies. Single-center, in vivo studies have provided some evidence for the added value of MRI as a biomarker of fracture risk or treatment response. Larger, prospective, multicenter studies are needed in the future to validate the results of these initial translational studies.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 5 J. MAGN. RESON. IMAGING 2017;46:323-337.

Changes in Nonosteoporotic Bone Density and Subsequent Fractures in Women

OBJECTIVES

Osteopenia is considerably more common than osteoporosis and accounts for most of the fracture burden in women older than 50 years. It is uncertain when to initiate treatment in osteopenia. We sought to determine in women with osteopenia what effect transitioning to lower categories had on subsequent fracturing.

METHODS

We surveyed 1150 women from office-based practices who had initial normal or osteopenic bone mineral densities (BMDs) and who were retested after 5.75 years. We classified categories related to baseline T scores as follows: normal (>-1.0), mild osteopenia (-1.0 to -1.49), moderate osteopenia (-1.5 to -1.99), and severe osteopenia (-2.0 to -2.49). We determined during a 9.6-year follow-up period the fracture occurrence in those who maintained their initial category status or transitioned into lower categories.

RESULTS

Transitioning to lower categories was not significantly different among baseline osteopenic categories but significantly more than normal baseline BMDs. Total fractures, individuals fracturing, and major fractures were significantly more, with baseline T scores of ≤-1.5 (<0.001). Although only 10.2% transitioned to osteoporosis, 90.5% of these transitions occurred with baseline T scores ≤-1.5 and accounted for significantly more fractures than baseline T scores of >-1.5.

CONCLUSIONS

Most subsequent fractures and transitions to osteoporosis occurred with baseline T scores ≤-1.5. Clinical risk factors need to be used to determine at what T score threshold treatment would be cost effective.

Digital Topology and Geometry in Medical Imaging: A Survey

Digital topology and geometry refers to the use of topologic and geometric properties and features for images defined in digital grids. Such methods have been widely used in many medical imaging applications, including image segmentation, visualization, manipulation, interpolation, registration, surface-tracking, object representation, correction, quantitative morphometry etc. Digital topology and geometry play important roles in medical imaging research by enriching the scope of target outcomes and by adding strong theoretical foundations with enhanced stability, fidelity, and efficiency. This paper presents a comprehensive yet compact survey on results, principles, and insights of methods related to digital topology and geometry with strong emphasis on understanding their roles in various medical imaging applications. Specifically, this paper reviews methods related to distance analysis and path propagation, connectivity, surface-tracking, image segmentation, boundary and centerline detection, topology preservation and local topological properties, skeletonization, and object representation, correction, and quantitative morphometry. A common thread among the topics reviewed in this paper is that their theory and algorithms use the principle of digital path connectivity, path propagation, and neighborhood analysis.

Magnetic resonance imaging assessed cortical porosity is highly correlated with μCT porosity

Cortical bone is typically regarded as "MR invisible" with conventional clinical magnetic resonance imaging (MRI) pulse sequences. However, recent studies have demonstrated that free water in the microscopic pores of cortical bone has a short T2* but a relatively long T2, and may be detectable with conventional clinical spin echo (SE) or fast spin echo (FSE) sequences. In this study we describe the use of a conventional two-dimensional (2D) FSE sequence to assess cortical bone microstructure and measure cortical porosity using a clinical 3T scanner. Twelve cadaveric human cortical bone samples were studied with MRI and microcomputed tomography (μCT) (downsampled to the same spatial resolution). Preliminary results show that FSE-determined porosity is highly correlated (R(2)=0.83; P<0.0001) with μCT porosity. Bland-Altman analysis suggested a good agreement between FSE and μCT with tight limit of agreement at around 3%. There is also a small bias of -2% for the FSE data, which suggested that the FSE approach slightly underestimated μCT porosity. The results demonstrate that cortical porosity can be directly assessed using conventional clinical FSE sequences. The clinical feasibility of this approach was also demonstrated on six healthy volunteers using 2D FSE sequences as well as 2D ultrashort echo time (UTE) sequences with a minimal echo time (TE) of 8μs, which provide high contrast imaging of cortical bone in vivo.

Effects of testosterone and growth hormone on the structural and mechanical properties of bone by micro-MRI in the distal tibia of men with hypopituitarism

CONTEXT

Severe deficiencies of testosterone (T) and GH are associated with low bone mineral density (BMD) and increased fracture risk. Replacement of T in hypogonadal men improves several bone parameters. Replacement of GH in GH-deficient men improves BMD.

OBJECTIVE

Our objective was to determine whether T and GH treatment together improves the structural and mechanical parameters of bone more than T alone in men with hypopituitarism.

DESIGN AND SUBJECTS

This randomized, prospective, 2-year study included 32 men with severe deficiencies of T and GH due to panhypopituitarism.

INTERVENTION

Subjects were randomized to receive T alone (n = 15) or T and GH (n = 17) for 2 years.

MAIN OUTCOME MEASURES

We evaluated magnetic resonance microimaging-derived structural (bone volume fraction [BVF] and trabecular thickness) and mechanical (axial stiffness [AS], a measure of bone strength) properties of the distal tibia at baseline and after 1 and 2 years of treatment.

RESULTS

Treatment with T and GH did not affect BVF, thickness, or AS differently from T alone. T treatment in all subjects for 2 years increased trabecular BVF by 9.6% (P < .0001), trabecular thickness by 2.6% (P < .001), and trabecular AS by 9.8% (P < .001). In contrast, testosterone treatment in all subjects significantly increased cortical thickness by 2.4% (P < .01) but decreased cortical BVF by -4.7% (P < .01) and cortical AS by -6.9% (P < .01).

CONCLUSION

Combined T and GH treatment of men with hypopituitarism for 2 years did not improve the measured structural or mechanical parameters of the distal tibia more than T alone. However, testosterone significantly increased the structural and mechanical properties of trabecular bone but decreased most of these properties of cortical bone, illustrating the potential importance of assessing trabecular and cortical bone separately in future studies of the effect of testosterone on bone.

Vertebral deformities and fractures are associated with MRI and pQCT measures obtained at the distal tibia and radius of postmenopausal women

SUMMARY

We investigated the association of postmenopausal vertebral deformities and fractures with bone parameters derived from distal extremities using MRI and pQCT. Distal extremity measures showed variable degrees of association with vertebral deformities and fractures, highlighting the systemic nature of postmenopausal bone loss.

INTRODUCTION

Prevalent vertebral deformities and fractures are known to predict incident further fractures. However, the association of distal extremity measures and vertebral deformities in postmenopausal women has not been fully established.

METHODS

This study involved 98 postmenopausal women (age range 60-88 years, mean 70 years) with DXA BMD T-scores at either the hip or spine in the range of -1.5 to -3.5. Wedge, biconcavity, and crush deformities were computed on the basis of spine MRI. Vertebral fractures were assessed using Eastell's criterion. Distal tibia and radius stiffness was computed using MRI-based finite element analysis. BMD at the distal extremities were obtained using pQCT.

RESULTS

Several distal extremity MRI and pQCT measures showed negative association with vertebral deformity on the basis of single parameter correlation (r up to 0.67) and two-parameter regression (r up to 0.76) models involving MRI stiffness and pQCT BMD. Subjects who had at least one prevalent vertebral fracture showed decreased MRI stiffness (up to 17.9 %) and pQCT density (up to 34.2 %) at the distal extremities compared to the non-fracture group. DXA lumbar spine BMD T-score was not associated with vertebral deformities.

CONCLUSIONS

The association between vertebral deformities and distal extremity measures supports the notion of postmenopausal osteoporosis as a systemic phenomenon.

Effects of immediate and delayed loading on peri-implant trabecular structures: a cone beam CT evaluation

PURPOSE

To develop a method for characterizing trabecular bone microarchitecture using cone beam computed tomography (CBCT) and to evaluate trabecular bone changes after rehabilitation using immediate versus delayed implant protocols.

MATERIALS AND METHODS

Six mongrel dogs randomly received 27 titanium implants in the maxillary incisor or mandibular premolar areas, following one of four protocols: (1) normal extraction socket healing; (2) immediate implant placement and immediate loading; (3) delayed implant placement and delayed loading; (4) delayed implant placement and immediate loading. The animals were euthanized at 8 weeks, and block biopsies were scanned using high resolution CBCT. Standard bone structural variables were assessed in coronal, middle, and apical levels.

RESULTS

Coronal and middle regions had more compact, more platelike, and thicker trabeculae. Protocols (2), (3), and (4) had significantly higher values (p < 0.001) than protocol (1) for bone surface density, bone surface volume ratio, and connectivity density, while significantly lower values (p < 0.001) were found for trabecular separation and fractal dimension. However, protocols (2), (3), and (4) did not show significantly different bone remodeling.

CONCLUSIONS

Compared with normal extraction healing, the implant protocols have an improved bone structural integration. Results do not suggest a different bone remodeling pattern when a delayed versus an immediate implant protocol is used.

Magnetic resonance of calcified tissues

MRI of the human body is largely made possible by the favorable relaxation properties of protons of water and triacyl glycerides prevalent in soft tissues. Hard tissues--key among them bone--are generally less amenable to measurement with in vivo MR imaging techniques, not so much as a result of the lower proton density but rather due to the extremely short life-times of the proton signal in water bound to solid-like entities, typically collagen, or being trapped in micro-pores. Either mechanism can enhance T2 relaxation by up to three orders of magnitude relative to their soft-tissue counterparts. Detection of these protons requires solid-state techniques that have emerged in recent years and that promise to add a new dimension to the study of hard tissues. Alternative approaches to probe calcified tissues exploit their characteristic magnetic properties. Bone, teeth and extra-osseous calcium-containing biomaterials are unique in that they are more diamagnetic than all other tissues and thus yield information indirectly by virtue of the induced magnetic fields present in their vicinity. Progress has also been made in methods allowing very high-resolution structural imaging of trabecular and cortical bone relying on detection of the surrounding soft-tissues. This brief review, much of it drawn from work conducted in the author's laboratory, seeks to highlight opportunities with focus on early-stage developments for image-based assessment of structure, function, physiology and mechanics of calcified tissues in humans via liquid and solid-state approaches, including proton, deuteron and phosphorus NMR and MRI.

Cortical and trabecular bone structure analysis at the distal radius-prediction of biomechanical strength by DXA and MRI

The purpose of this study was to investigate whether the combination of dual-energy X-ray absorptiometry (DXA)-based bone mass and magnetic resonance imaging (MRI)-based cortical and trabecular structural measures improves the prediction of radial bone strength. Thirty-eight left forearms were harvested from formalin-fixed human cadavers. Bone mineral content (BMC) and bone mineral density (BMD) of the distal radius were measured using DXA. Cortical and trabecular structural measures of the distal radius were computed in high-resolution 1.5T MR images. Cortical measures included average cortical thickness and cross-sectional area. Trabecular measures included morphometric and texture parameters. The forearms were biomechanically tested in a fall simulation to measure absolute radial bone strength (failure load). Relative radial bone strength was determined by dividing radial failure loads by age, body mass index, radius length, and average radius cross-sectional area, respectively. DXA derived BMC and BMD showed statistically significant (p < 0.05) correlations with absolute and relative radial bone strength (r ≤ 0.78). Correlation coefficients for cortical and trabecular structural measures with absolute and relative radial bone strength amounted up to r = 0.59 and r = 0.74, respectively, (p < 0.05). In combination with DXA-based bone mass, trabecular but not, cortical structural measures, added in multiple regression models significant (p < 0.05) information in predicting absolute and relative radial bone strength (up to R adj = 0.88). Thus, a combination of DXA-based bone mass and MRI-based trabecular structural measures most accurately predicted absolute and relative radial bone strength, whereas structural measures of the cortex did not provide significant additional information in combination with DXA.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1-L3 and total hip, QCT of L1-L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8-9.0], Tb.BMD L1-L3: 3.95 [1.8-8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29-5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r(2)=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

Osteoprotective effect of hormone therapy on bone microarchitecture before impaired bone mineral density in ovariectomized rats

OBJECTIVE

We aimed to determine the effect of hormone replacement therapy on bone microarchitecture in ovariectomized rats.

MATERIAL AND METHODS

In the Animal Ethics Committee approved-study, the effect of treatment with 17 β-estradiol 50 μg/kg and medroxyprogesterone 2.5 mg/kg on bone architecture and bone mineral density in rats versus ovariectomized control rats over the course of 20 days were evaluated. Femoral and lumbar bone mineral density levels and morphometric measurements were performed.

RESULTS

There were no significant differences in the femoral and lumbar bone mineral density levels between the groups. In the intact control group, the trabecular structures were significantly superior to those in the other groups. Additionally, the osteoblast count was significantly higher while the osteoclast count was significantly lower than in all other groups. Two parameters reflecting trabecular bone microarchitecture, which include the trabecular count and the trabecular area, demonstrated significant improvement in the hormone replacement group when compared to the ovariectomized control group. In the hormone replacement groups, the osteoblast count was significantly higher while the osteoclast count was significantly lower than in the ovariectomized control group.

CONCLUSION

We suggest that offering estrogen alone or in combination with progestogen can be a beneficial approach in preventing early postmenopausal bone loss regardless of bone mineral density.

Evaluation of trabecular microarchitecture in nonosteoporotic postmenopausal women with and without fracture

This study compared microscopic magnetic resonance imaging (µMRI) parameters of trabecular microarchitecture between postmenopausal women with and without fracture who have normal or osteopenic bone mineral density (BMD) on dual-energy X-ray absorptiometry (DXA). It included 36 postmenopausal white women 50 years of age and older with normal or osteopenic BMD (T-scores better than -2.5 at the lumbar spine, proximal femur, and one-third radius on DXA). Eighteen women had a history of low-energy fracture, whereas 18 women had no history of fracture and served as an age, race, and ultradistal radius BMD-matched control group. A three-dimensional fast large-angle spin-echo (FLASE) sequence with 137 µm × 137 µm × 400 µm resolution was performed through the nondominant wrist of all 36 women using the same 1.5T scanner. The high-resolution images were used to measure trabecular bone volume fraction, trabecular thickness, surface-to-curve ratio, and erosion index. Wilcoxon signed-rank tests were used to compare differences in BMD and µMRI parameters between postmenopausal women with and without fracture. Post-menopausal women with fracture had significantly lower (p < 0.05) trabecular bone volume fraction and surface-to-curve ratio and significantly higher (p < 0.05) erosion index than postmenopausal women without fracture. There was no significant difference between postmenopausal women with and without fracture in trabecular thickness (p = 0.80) and BMD of the spine (p = 0.21), proximal femur (p = 0.19), one-third radius (p = 0.47), and ultradistal radius (p = 0.90). Postmenopausal women with normal or osteopenic BMD who had a history of low-energy fracture had significantly different (p < 0.05) µMRI parameters than an age, race, and ultradistal radius BMD-matched control group of postmenopausal women with no history of fracture. Our study suggests that µMRI can be used to identify individuals without a DXA-based diagnosis of osteoporosis who have impaired trabecular microarchitecture and thus a heretofore-unappreciated elevated fracture risk.

Trabecular architecture and vertebral fragility in osteoporosis

Osteoporosis heightens vertebral fragility owing to the biomechanical effects of diminished bone structure and composition. These biomechanical effects are only partially explained by loss in bone mass, so additional factors that are independent of bone mass are also thought to play an important role in vertebral fragility. Recent advances in imaging equipment, imaging-processing methods, and computational capacity allow researchers to quantify trabecular architecture in the vertebra at the level of the individual trabecular elements and to derive biomechanics-based measures of architecture that are independent of bone mass and density. These advances have shed light on the role of architecture in vertebral fragility. In addition to the adverse biomechanical consequences associated with trabecular thinning and loss of connectivity, a reduction in the number of vertical trabecular plates appears to be particularly harmful to vertebral strength. In the clinic, detailed architecture analysis is primarily applied to peripheral sites such as the distal radius and tibia. Analysis of trabecular architecture at these peripheral sites has shown mixed results for discriminating between patients with and without a vertebral fracture independent of bone mass, but has the potential to provide unique insight into the effects of therapeutic treatments. Overall, it does appear that trabecular architecture has an independent role on vertebral strength. Additional research is required to determine how and where architecture should be measured in vivo and whether assessment of trabecular architecture in a clinical setting improves prospective fracture risk assessment for the vertebra.

Assessing fracture risk using gradient boosting machine (GBM) models

Advanced bone imaging with quantitative computed tomography (QCT) has had limited success in significantly improving fracture prediction beyond standard areal bone mineral density (aBMD) measurements. Thus, we examined whether a machine learning paradigm, gradient boosting machine (GBM) modeling, which can incorporate diverse measurements of bone density and geometry from central QCT imaging and of bone microstructure from high-resolution peripheral QCT imaging, can improve fracture prediction. We studied two cohorts of postmenopausal women: 105 with and 99 without distal forearm fractures (Distal Forearm Cohort) and 40 with at least one grade 2 or 3 vertebral deformity and 78 with no vertebral fracture (Vertebral Cohort). Within each cohort, individual bone density, structure, or strength variables had areas under receiver operating characteristic curves (AUCs) ranging from 0.50 to 0.84 (median 0.61) for discriminating women with and without fracture. Using all possible variables in the GBM model, the AUCs were close to 1.0. Fracture predictions in the Vertebral Cohort using the GBM models built with the Distal Forearm Cohort had AUCs of 0.82-0.95, while predictions in the Distal Forearm Cohort using models built with the Vertebral Cohort had AUCs of 0.80-0.83. Attempts at capturing a comparable parametric model using the top variables from the Distal Forearm Cohort resulted in resulted in an AUC of 0.81. Relatively high AUCs for differing fracture types suggest that an underlying fracture propensity is being captured by this modeling approach. More complex modeling, such as with GBM, creates stronger fracture predictions and may allow deeper insights into information provided by advanced bone imaging techniques.

Micro-MR imaging-based computational biomechanics demonstrates reduction in cortical and trabecular bone strength after renal transplantation

PURPOSE

To examine the ability of three-dimensional micro-magnetic resonance (MR) imaging-based computational biomechanics to detect mechanical alterations in trabecular bone and cortical bone in the distal tibia of incident renal transplant recipients 6 months after renal transplantation and compare them with bone mineral density (BMD) outcomes.

MATERIALS AND METHODS

The study was approved by the institutional review board and complied with HIPAA guidelines. Written informed consent was obtained from all subjects. Micro-MR imaging of distal tibial metaphysis was performed within 2 weeks after renal transplantation (baseline) and 6 months later in 49 participants (24 female; median age, 44 years; range, 19-61 years) with a clinical 1.5-T whole-body imager using a modified three-dimensional fast large-angle spin-echo pulse sequence. Micro-finite-element models for cortical bone, trabecular bone, and whole-bone section were generated from each image by delineating the endosteal and periosteal boundaries. Mechanical parameters (stiffness and failure load) were estimated with simulated uniaxial compression tests on the micro-finite-element models. Structural parameters (trabecular bone volume fraction [BV/TV, bone volume to total volume ratio], trabecular thickness [TbTh], and cortical thickness [CtTh]) were computed from micro-MR images. Total hip and spine areal BMD were determined with dual-energy x-ray absorptiometry (DXA). Parameters obtained at the follow-up were compared with the baseline values by using parametric or nonparametric tests depending on the normality of data.

RESULTS

All mechanical parameters were significantly lower at 6 months compared with baseline. Decreases in cortical bone, trabecular bone, and whole-bone stiffness were 3.7% (P = .03), 4.9% (P = .03), and 4.3% (P = .003), respectively. Decreases in cortical bone, trabecular bone, and whole-bone failure strength were 7.6% (P = .0003), 6.0% (P = .004), and 5.6% (P = .0004), respectively. Conventional structural measures, BV/TV, TbTh, and CtTh, did not change significantly. Spine BMD decreased by 2.9% (P < .0001), while hip BMD did not change significantly at DXA.

CONCLUSION

MR imaging-based micro-finite-element analysis suggests that stiffness and failure strength of the distal tibia decrease over a 6-month interval after renal transplantation.

On the significance of motion degradation in high-resolution 3D μMRI of trabecular bone

RATIONALE AND OBJECTIVES

Subtle subject movement during high-resolution three-dimensional micro-magnetic resonance imaging of trabecular bone (TB) causes blurring, thereby rendering the data unreliable for quantitative analysis. In this work, the effects of translational and rotational motion displacements were evaluated qualitatively and quantitatively.

MATERIALS AND METHODS

In experiment 1, motion was induced by applying various simulated and previously observed in vivo trajectories as phase shifts to k-space or rotation angles to k-space segments of a virtually motion-free data set. In experiment 2, images that were visually free of motion artifacts from two groups of 10 healthy individuals, differing in age, were selected to probe the effects of motion on TB parameters. In both experiments, images were rated for motion severity, and the scores were compared to a focus criterion, the normalized gradient squared.

RESULTS

Strong correlations were observed between the motion quality scores and the corresponding normalized gradient squared values (R(2) = 0.52-0.64, P < .01). The results from experiment 1 demonstrated consistently lower image quality and alterations in structural parameters of 9% to 45% with increased amplitude of displacements. In experiment 2, the significant differences in structural parameter group means of the motion-free images were lost upon motion degradation. Autofocusing, a postprocessing correction method, partially recovered the sharpness of the original motion-free images in 13 of 20 subjects.

CONCLUSIONS

Quantitative TB structural measures are highly sensitive to subtle motion-induced degradation, which adversely affects precision and statistical power. The results underscore the influence of subject movement in high-resolution three-dimensional micro-magnetic resonance imaging and its correction for TB structure analysis.

Performance of the MRI-based virtual bone biopsy in the distal radius: serial reproducibility and reliability of structural and mechanical parameters in women representative of osteoporosis study populations

Serial reproducibility and reliability critically determine sensitivity to detect changes in response to intervention and provide a basis for sample size estimates. Here, we evaluated the performance of the MRI-based virtual bone biopsy in terms of 26 structural and mechanical parameters in the distal radius of 20 women in the age range of 50 to 75 years (mean=62.0 years, S.D.=8.1 years), representative of typical study populations in drug intervention trials and fracture studies. Subjects were examined three times at average intervals of 20.2 days (S.D.=14.5 days) by MRI at 1.5 T field strength at a voxel size of 137×137×410 μm(3). Methods involved prospective and retrospective 3D image registration and auto-focus motion correction. Analyses were performed from a central 5×5×5 mm(3) cuboid subvolume and trabecular volume consisting of a 13 mm axial slab encompassing the entire medullary cavity. Whole-volume axial stiffness and sub-regional Young's and shear moduli were computed by finite-element analysis. Whole-volume-derived aggregate mean coefficient of variation of all structural parameters was 4.4% (range 1.8% to 7.7%) and 4.0% for axial stiffness; corresponding data in the subvolume were 6.5% (range 1.6% to 13.0%) for structural, and 5.5% (range 4.6% to 6.5%) for mechanical parameters. Aggregate ICC was 0.976 (range 0.947 to 0.986) and 0.992 for whole-volume-derived structural parameters and axial stiffness, and 0.946 (range 0.752 to 0.991) and 0.974 (range 0.965 to 0.978) for subvolume-derived structural and mechanical parameters, respectively. The strongest predictors of whole-volume axial stiffness were BV/TV, junction density, skeleton density and Tb.N (R(2) 0.79-0.87). The same parameters were also highly predictive of sub-regional axial modulus (R(2) 0.88-0.91). The data suggest that the method is suited for longitudinal assessment of the response to therapy. The underlying technology is portable and should be compatible with all general-purpose MRI scanners, which is appealing considering the very large installed base of this modality.

Performance of μMRI-Based virtual bone biopsy for structural and mechanical analysis at the distal tibia at 7T field strength

PURPOSE

To assess the performance of a 3D fast spin echo (FSE) pulse sequence utilizing out-of-slab cancellation through phase alternation and micro-magnetic resonance imaging (μMRI)-based virtual bone biopsy processing methods to probe the serial reproducibility and sensitivity of structural and mechanical parameters of the distal tibia at 7.0T.

MATERIALS AND METHODS

The distal tibia of five healthy subjects was imaged at three timepoints with a 3D FSE sequence at 137 × 137 × 410 μm(3) voxel size. Follow-up images were retrospectively 3D registered to baseline images. Coefficients of variation (CV) and intraclass correlation coefficients (ICCs) for measures of scale and topology of the whole tibial trabecular bone (TB) cross-section as well as finite-element-derived Young's and shear moduli of central cuboidal TB subvolumes (8 × 8 × 5 mm(3) ) were evaluated as measures of reproducibility and reliability. Four additional cubic TB subregions (anterior, medial, lateral, and posterior) of similar dimensions were extracted and analyzed to determine associations between whole cross-section and subregional structural parameters.

RESULTS

The mean signal-to-noise ratio (SNR) over the 15 image acquisitions was 27.5 ± 2.1. Retrospective registration yielded an average common analysis volume of 67% across the three exams per subject. Reproducibility (mean CV = 3.6%; range, 1.5%-5%) and reliability (ICCs, 0.95-0.99) of all parameters permitted parameter-based discrimination of the five subjects in spite of the narrow age range (26-36 years) covered. Parameters characterizing topology were better able to distinguish two individuals who demonstrated similar values for scalar measurements (≈ 34% difference, P < 0.001). Whole-section axial stiffness encompassing the cortex was superior at distinguishing two individuals relative to its central subregional TB counterpart (≈ 8% difference; P < 0.05). Interregion comparisons showed that although all parameters were correlated (mean R(2) = 0.78; range 0.57-0.99), the strongest associations observed were those for the erosion index (mean R(2) = 0.95, P ≤ 0.01).

CONCLUSION

The reproducibility and structural and mechanical parameter-based discriminative ability achieved in five healthy subjects suggests that 7T-derived μMRI of TB can be applied towards serial patient studies of osteoporosis and may enable earlier detection of disease or treatment-based effects.

Recent progress in bone imaging for osteoporosis research

Advances in bone imaging techniques have provided tools for analyzing bone structure at the macro-, micro- and nano-level. Quantitative assessment of macrostructure can be achieved using dual X-ray absorptiometry (DXA) and quantitative computed tomography (QCT), particularly volumetric quantitative CT (vQCT). In vivo quantitative techniques for assessing the microstructure of trabecular bone non-invasively and non-destructively include high-resolution CT (HR-CT) and high-resolution magnetic resonance (HR-MR). Compared with MR imaging, CT-based techniques have the advantage of directly visualizing the bone in the axial skeleton, with high spatial resolution, but the disadvantage of delivering a considerable radiation dose. Micro-CT (μCT), which provides a higher resolution of the microstructure and is principally applicable in vitro, has undergone technological advances such that it is now able to elucidate the physiological skeletal change mechanisms associated with aging and determine the effects of therapeutic intervention on the bone microstructure. In particular, synchrotron μCT (SR-CT) provides a more detailed view of trabecular structure at the nano-level. For the assessment of hip geometry, DXA-based hip structure analysis (HSA) and CT-based HSA have been developed. DXA-based HSA is a convenient tool for analyzing biomechanical properties and for assuming cross-sectional hip geometry based on two-dimensional (2D) data, whereas CT-based HSA provides these parameters three-dimensionally in robust relationship with biomechanical properties, at the cost of greater radiation exposure and the lengthy time required for the analytical procedure. Further progress in bone imaging technology is promising to bring new aspects of bone structure in relation to bone strength to light, and to establish a means for analyzing bone structural properties in the everyday clinical setting.

Influence of vertical trabeculae on the compressive strength of the human vertebra

Vertebral strength, a key etiologic factor of osteoporotic fracture, may be affected by the relative amount of vertically oriented trabeculae. To better understand this issue, we performed experimental compression testing, high-resolution micro-computed tomography (µCT), and micro-finite-element analysis on 16 elderly human thoracic ninth (T(9)) whole vertebral bodies (ages 77.5 ± 10.1 years). Individual trabeculae segmentation of the µCT images was used to classify the trabeculae by their orientation. We found that the bone volume fraction (BV/TV) of just the vertical trabeculae accounted for substantially more of the observed variation in measured vertebral strength than did the bone volume fraction of all trabeculae (r(2) = 0.83 versus 0.59, p < .005). The bone volume fraction of the oblique or horizontal trabeculae was not associated with vertebral strength. Finite-element analysis indicated that removal of the cortical shell did not appreciably alter these trends; it also revealed that the major load paths occur through parallel columns of vertically oriented bone. Taken together, these findings suggest that variation in vertebral strength across individuals is due primarily to variations in the bone volume fraction of vertical trabeculae. The vertical tissue fraction, a new bone quality parameter that we introduced to reflect these findings, was both a significant predictor of vertebral strength alone (r(2) = 0.81) and after accounting for variations in total bone volume fraction in multiple regression (total R(2) = 0.93). We conclude that the vertical tissue fraction is a potentially powerful microarchitectural determinant of vertebral strength.

Interrelationships between 3-T-MRI-derived cortical and trabecular bone structure parameters and quantitative-computed-tomography-derivedbone mineral density

Recently, 3-T magnetic resonance imaging (MRI) has been introduced for bone imaging. Through higher signal-to-noise ratios, as compared to 1.5-T MRI, it promises to be a more powerful tool for the assessment of cortical and trabecular bone measures. The goal of our study was to compare MRI-derived cortical and trabecular bone measures to quantitative computed tomography (QCT)-derived bone mineral density (BMD). Using 3-T MRI in 51 postmenopausal women, apparent (app.) measures of bone volume/total volume, trabecular number (Tb.N), trabecular thickness (Tb.Th) and trabecular separation were derived at the distal radius, distal tibia and calcaneus. Cortical thickness (Ct.Th) was calculated at the distal radius and distal tibia. These measures were compared to QCT-derived BMD of the spine, hip and radius. Significant correlations ((*)P<.05; (**)P<.001; (***)P<.0001) were found between spine BMD- and MRI-derived Ct.Th (r(radius)=.55, (*)P<.05; r(tibia)=.67, (***)P<.0001) and app. Tb.N (r(radius)=.33, (*)P<.05; r(tibia)=.35, (*)P<.05) at the radius and tibia. Furthermore, within the first 10 mm at the radius, an inverse correlation for Ct.Th and app. BV/TV (r(6mm)=-.56, P<.001; r(10mm)=-.36, P<.05) and app. Tb.Th (r(6mm)=-.54, P<.001; r(10mm)=-.41, P<.05) was found.

Bone density, geometry, microstructure, and stiffness: Relationships between peripheral and central skeletal sites assessed by DXA, HR-pQCT, and cQCT in premenopausal women

High-resolution peripheral quantitative computed tomography (HR-pQCT) is a new in vivo imaging technique for assessing 3D microstructure of cortical and trabecular bone at the distal radius and tibia. No studies have investigated the extent to which measurements of the peripheral skeleton by HR-pQCT reflect those of the spine and hip, where the most serious fractures occur. To address this research question, we performed dual-energy X-ray absorptiometry (DXA), central QCT (cQCT), HR-pQCT, and image-based finite-element analyses on 69 premenopausal women to evaluate relationships among cortical and trabecular bone density, geometry, microstructure, and stiffness of the lumbar spine, proximal femur, distal radius, and distal tibia. Significant correlations were found between the stiffness of the two peripheral sites (r = 0.86), two central sites (r = 0.49), and between the peripheral and central skeletal sites (r = 0.56-0.70). These associations were explained in part by significant correlations in areal bone mineral density (aBMD), volumetric bone mineral density (vBMD), and cross-sectional area (CSA) between the multiple skeletal sites. For the prediction of proximal femoral stiffness, vBMD (r = 0.75) and stiffness (r = 0.69) of the distal tibia by HR-pQCT were comparable with direct measurements of the proximal femur: aBMD of the hip by DXA (r = 0.70) and vBMD of the hip by cQCT (r = 0.64). For the prediction of vertebral stiffness, trabecular vBMD (r = 0.58) and stiffness (r = 0.70) of distal radius by HR-pQCT were comparable with direct measurements of lumbar spine: aBMD by DXA (r = 0.78) and vBMD by cQCT (r = 0.67). Our results suggest that bone density and microstructural and mechanical properties measured by HR-pQCT of the distal radius and tibia reflect the mechanical competence of the central skeleton.

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.

Mechanical implications of estrogen supplementation in early postmenopausal women

Whereas the structural implications of drug intervention are well established, there are few data on the possible mechanical consequences of treatment. In this work we examined the changes in elastic and shear moduli (EM and SM) in a region of trabecular bone in the distal radius and distal tibia of early postmenopausal women on the basis of MRI-based micro-finite-element (microFE) analysis. Whole-section axial stiffness (AS) encompassing both trabecular and cortical compartments was evaluated as well. The study was conducted on previously acquired high-resolution images at the two anatomic sites. Images were processed to yield a 3D voxel array of bone-volume fraction (BVF), which was converted to a microFE model of hexahedral elements in which tissue modulus was set proportional to voxel BVF. The study comprised 65 early postmenopausal women (age range 45 to 55 years), of whom 32 had chosen estrogen supplementation (estradiol group); the remainder had not (control group). Subjects had been scanned at baseline and 12 and 24 months thereafter. At the distal tibia, EM and SM were reduced by 2.9% to 5.5% in the control group (p < .05 to <.005), but there was no change in the estradiol subjects. AS decreased 3.9% (4.0%) in controls (p < .005) and increased by 5.8% (6.2%) in estradiol group subjects (p < .05) at 12 (24) months. At the distal radius, EM and SM changes from baseline were not significant, but at both time points AS was increased in estradiol group subjects and decreased in controls (p < .005 to <.05), albeit by a smaller margin than at the tibia. EM and SM were strongly correlated with BV/TV (r(2) = 0.44 to 0.92) as well as with topologic parameters expressing the ratio of plates to rods (r(2) = 0.45 to 0.82), jointly explaining up to 96% of the variation in the mechanical parameters. Finally, baseline AS was strongly correlated between the two anatomic sites (r(2) = 0.58), suggesting that intersubject variations in the bone's mechanical competence follows similar mechanisms. In conclusion, the results demonstrate that micro-MRI-based microFE models are suited for the study of the mechanical implications of antiresorptive treatment. The data further highlight the anabolic effect of short-term estrogen supplementation.

Changes in trabecular microarchitecture in postmenopausal women on bisphosphonate therapy

PURPOSE

In addition to bone mineral density (BMD), trabecular microstructure contributes to skeletal strength. Our goal was to examine changes in trabecular microstructure in women on therapy.

MATERIALS AND METHODS

We followed 10 postmenopausal women receiving a bisphosphonate, risedronate (35 mg once weekly), over 12 months and examined trabecular microarchitecture with high resolution wrist MR images (hr-MRI). MRI parameters included bone volume/total volume (BV/TV), surface density (representing plates), curve density (representing rods), surface-to-curve ratio and erosion index (depicting deterioration). We assessed BMD of the spine, hip and radius and markers of bone turnover.

RESULTS

Women had been receiving bisphosphonate therapy for 43+/-9 months (mean+/-SD) prior to the first MRI. Indices of hr-MRI demonstrated improvement in surface-to-curve ratio (13.0%) and a decrease in erosion index (12.1%) consistent with less deterioration (both p<0.05). BMD of the spine, hip and radius and markers of bone turnover remained stable. Parameters of hr-MRI were associated with 1/3 distal radius BMD (correlation coefficient 0.71 to 0.86, p<0.05).

DISCUSSION

We conclude that hr-MRI of the radius demonstrates improvements in trabecular microstructure not appreciated by conventional BMD and provides additional information on parameters that contribute to structural integrity in patients on antiresorptive therapy.

Advances in osteoporosis imaging

In the assessment of osteoporosis, the measurement of bone mineral density (BMD(a)) obtained from dual energy X-ray absorptiometry (DXA; g/cm(2)) is the most widely used parameter. However, bone strength and fracture risk are also influenced by parameters of bone quality such as micro-architecture and tissue properties. This article reviews the radiological techniques currently available for imaging and quantifying bone structure, as well as advanced techniques to image bone quality. With the recent developments in magnetic resonance (MR) techniques, including the availability of clinical 3T scanners, and advances in computed tomography (CT) technology (e.g. clinical Micro-CT), in-vivo imaging of the trabecular bone architecture is becoming more feasible. Several in-vitro studies have demonstrated that bone architecture, measured by MR or CT, was a BMD-independent determinant of bone strength. In-vivo studies showed that patients with, and without, osteoporotic fractures could better be separated with parameters of bone architecture than with BMD. Parameters of trabecular architecture were more sensitive to treatment effects than BMD. Besides the 3D tomographic techniques, projection radiography has been used in the peripheral skeleton as an additional tool to better predict fracture risk than BMD alone. The quantification of the trabecular architecture included parameters of scale, shape, anisotropy and connectivity. Finite element analyses required highest resolution, but best predicted the biomechanical properties of the bone. MR diffusion and perfusion imaging and MR spectroscopy may provide measures of bone quality beyond trabecular micro-architecture.

Relationship between ultrasonic parameters and apparent trabecular bone elastic modulus: a numerical approach

The physical principles underlying quantitative ultrasound (QUS) measurements in trabecular bone are not fully understood. The translation of QUS results into bone strength remains elusive. However, ultrasound being mechanical waves, it is likely to assess apparent bone elasticity. The aim of this study is to derive the sensitivity of QUS parameters to variations of apparent bone elasticity, a surrogate for strength. The geometry of 34 human trabecular bone samples cut in the great trochanter was reconstructed using 3-D synchrotron micro-computed tomography. Finite-difference time-domain simulations coupled to 3-D micro-structural models were performed in the three perpendicular directions for each sample and each direction. A voxel-based micro-finite element linear analysis was employed to compute the apparent Young's modulus (E) of each sample for each direction. For the antero-posterior direction, the predictive power of speed of sound and normalized broadband ultrasonic attenuation to assess E was equal to 0.9 and 0.87, respectively, which is better than what is obtained using bone density alone or coupled with micro-architectural parameters and of the same order of what can be achieved with the fabric tensor approach. When the direction of testing is parallel to the main trabecular orientation, the predictive power of QUS parameters decreases and the fabric tensor approach always gives the best results. This decrease can be explained by the presence of two longitudinal wave modes. Our results, which were obtained using two distinct simulation tools applied on the same set of samples, highlight the potential of QUS techniques to assess bone strength.

Geodesic topological analysis of trabecular bone microarchitecture from high-spatial resolution magnetic resonance images

In vivo assessment of trabecular bone microarchitecture could improve the prediction of fracture risk and the efficacy of osteoporosis treatment and prevention. Geodesic topological analysis (GTA) is introduced as a novel technique to quantify the trabecular bone microarchitecture from high-spatial resolution magnetic resonance (MR) images. Trabecular bone parameters that quantify the scale, topology, and anisotropy of the trabecular bone network in terms of its junctions are the result of GTA. The reproducibility of GTA was tested with in vivo images of human distal tibiae and radii (n = 6) at 1.5 Tesla; and its ability to discriminate between subjects with and without vertebral fracture was assessed with ex vivo images of human calcanei at 1.5 and 3.0 Tesla (n = 30). GTA parameters yielded an average reproducibility of 4.8%, and their individual areas under the curve (AUC) of the receiver operating characteristic curve analysis for fracture discrimination performed better at 3.0 than at 1.5 Tesla reaching values of up to 0.78 (p < 0.001). Logistic regression analysis demonstrated that fracture discrimination was improved by combining GTA parameters, and that GTA combined with bone mineral density (BMD) allow for better discrimination than BMD alone (AUC = 0.95; p < 0.001). Results indicate that GTA can substantially contribute in studies of osteoporosis involving imaging of the trabecular bone microarchitecture.

Implications of resolution and noise for in vivo micro-MRI of trabecular bone

Osteoporotic bone loss is accompanied by impaired structural integrity of the trabecular network, leading to a decrease in the overall mechanical properties of the bone. The development of the "virtual bone biopsy" (VBB), a method combining magnetic resonance microimaging (microMRI) and digital image processing techniques, has previously been shown to quantify topology and scale of human trabecular bone noninvasively. The aim of this work was to determine the extent to which structural parameters derived from images acquired in the limited spatial resolution regime of in vivo imaging are sensitive to resolution and noise and further, whether under these conditions, a small amount of bone loss and its associated structural manifestations can be detected. Toward these goals 3D models of trabecular bone representing multiple anatomic locations were generated on the basis of microCT images of human cadaveric bone cores. These images were binarized and the resulting data arrays representing pure bone (proton density=0) and pure marrow (proton density=255) subjected to simulated MR imaging by Cartesian sampling of k space, yielding, after 3D Fourier reconstruction, voxel sizes currently achievable in vivo. Subsequently, realistic levels of Gaussian noise were superimposed on the complex data and magnitude images were computed. The resulting images were subsequently VBB processed for a range of signal-to-noise ratio (SNR) values and image voxel sizes. For comparison of the predicted behavior to in vivo data, images from a recent patient study were evaluated as well. Systematic changes of the derived structural parameters changing progressively with decreasing SNR were noted, and it is shown that the errors are correctable using simple linear transformations, thereby allowing the data to be normalized. The predicted dependence of the structural parameters on SNR also closely parallel those observed in vivo. Finally, in order to assess the sensitivity of the VBB processing algorithms to detect bone loss during disease progression or regression in response to treatment, the high-resolution specimen data were subjected to 5% bone loss either by homogeneous or heterogeneous erosion and microMR images simulated at in vivo resolution and SNR. At typical in vivo SNR (SNR=12) and effective image resolution (160 microm isotropic and 137 x 137 X 410 microm3), VBB algorithms were able to detect the structural implications of a 5% loss in bone volume fraction with high statistical significance.

Retrospective 3D registration of trabecular bone MR images for longitudinal studies

PURPOSE

To evaluate an automatic 3D registration algorithm for serial high-resolution images of trabecular bone (TB) in studies designed to evaluate the response of the trabecular architecture to intervention or disease progression.

MATERIALS AND METHODS

An efficient algorithm for registering high-resolution 3D images of TB is presented. The procedure identifies the six parameters of rigid displacement between two scans performed at different timepoints. By assuming a relatively small through-plane rotation, considerable time is saved by combining the results of a collection of regional 2D registrations throughout the TB region of interest (ROI). The algorithm was applied to 26 pairs of MR images acquired 6 months apart. Reproducibility of local TB structural parameters (plate, rod, and junction density) computed in manually selected regions were compared between baseline and registered follow-up images.

RESULTS

All 26 registrations were completed successfully in less than 30 seconds per image pair. The resampled follow-up images agreed with baseline to around one pixel throughout the volume at 137 x 137 x 410 microm(3) image resolution. Structural parameters in each region correlated well from baseline to follow-up with intraclass correlation coefficients ranging between 85%-97% for TB plate density. Interregional variations in the parameters were large as compared with intraregion reproducibility.

CONCLUSION

The proposed algorithm was successful in automatically registering baseline and follow-up TB images in a translational study, and may be useful in regional analyses in longitudinal MR studies of TB architecture.

Trabecular bone structure analysis in the limited spatial resolution regime of in vivo MRI

RATIONALE AND OBJECTIVES

To develop a method for processing and visualization of trabecular bone networks on the basis of magnetic resonance (MR) images acquired in the limited spatial resolution regime of in vivo imaging at which trabecular thickness is comparable to voxel size.

MATERIALS AND METHODS

A sequence of processing steps for analyzing the topologic structure of trabecular bone networks is presented and evaluated using three types of datasets: images of synthetic structures with various levels of superimposed Gaussian noise, micro-computed tomographic images of human trabecular bone downsampled to in vivo resolution, and in vivo micro-MR images from a prior longitudinal study investigating the structural implications of testosterone treatment of hypogonadal men. The simulated images were analyzed at a voxel size of 150 microm(3), the clinical MR image data had been acquired with 137 x 137 x 410 microm(3) voxel size. The technique is a modification to the virtual bone biopsy processing chain that involves a sinc convolution step immediately preceding binarization, and employs the Manzanera-Bernard thinning algorithm for obtaining the three-dimensional skeleton before topologic classification. The detectability of plate and rod bone elements was also analyzed theoretically.

RESULTS

As compared with previously published techniques, the approach produced a more accurate bone skeleton in the micro-computed tomographic and simulation experiments, with clear improvement in preservation of rod and plate elements. Simulations suggest that rods are detectable down to a diameter of approximately 50% of the MR image voxel length, whereas plates can be detected at thicknesses of 20% or more of voxel length. For in vivo studies, it was shown that the method could recover the treatment response in terms of the ensuing topologic changes in patients undergoing antiresorptive treatment.

CONCLUSIONS

The algorithm for processing of in vivo micro-MR images of trabecular bone is superior to prior approaches in preserving the topology of the network in the presence of noise.

Inadequate responders to osteoporosis treatment: proposal for an operational definition

The concept of inadequate response to osteoporosis treatment is not clear. In the literature several criteria have been used. We propose an operational definition of an inadequate responder based on the changes observed in bone mineral density and incident fractures while on therapy.

INTRODUCTION

Fractures may occur in compliant patients even while on active treatment. These cases have been defined as inadequate responders (IR).

METHODS

We reviewed the basis for this concept and propose an operational definition for IR.

RESULTS

Good compliance and adequate calcium and vitamin D supplementation are the first requirement. The second requirement is a treatment period of at least 1 year, since before that time treatment may not have been fully effective. Fractures are the gold standard for measuring efficacy and changes in bone density and turnover markers may be surrogates. We propose classifying patient response as: Inadequate--incident fracture and a decrease in BMD greater than a significant change (Trend Assessment Margin or TAM); Possibly inadequate--incident fracture or a decrease in BMD greater than a significant change (TAM); and Appropriate--no fracture and no decrease in BMD greater than a significant change (TAM). Additional criteria (biochemical markers, bone quality parameters) may be taken into account.

CONCLUSION

A wide consensus on the IR concept is required given its clinical, regulatory, and reimbursement implications.

In vivo magnetic resonance detects rapid remodeling changes in the topology of the trabecular bone network after menopause and the protective effect of estradiol

INTRODUCTION

Estrogen depletion after menopause is accompanied by bone loss and architectural deterioration of trabecular bone. The hypothesis underlying this work is that the microMRI-based virtual bone biopsy can capture the temporal changes of scale and topology of the trabecular network and that estrogen supplementation preserves the integrity of the trabecular network.

MATERIALS AND METHODS

Subjects studied were early postmenopausal women, 45-55 yr of age (N = 65), of whom 32 were on estrogen (estradiol group), and the remainder were not (control group). Early menopause was defined by amenorrhea for 6-24 mo and elevated serum follicle-stimulating hormone (FSH) concentration. The subjects were evaluated with three imaging modalities at baseline and 12 and 24 mo to determine the temporal changes in trabecular and cortical architecture and density. microMRI of the distal radius and tibia was performed at 137 x 137 x 410-microm(3) voxel size. The resulting bone volume fraction maps were Fourier interpolated to a final voxel size of 45.7 x 45.7 x 136.7 microm(3), binarized, skeletonized, and subjected to 3D digital topological analysis (DTA). Skeletonization converts trabecular rods to curves and plates to surfaces. Parameters quantifying scale included BV/TV, whereas DTA parameters included the volume densities of curves (C) and surface (S)-type voxels, as well as composite parameters: the surface/curve ratio (S/C), and erosion index (EI, ratio of the sum of parameters expected to increase with osteoclastic resorption divided by the sum of those expected to decrease). For comparison, pQCT of the same peripheral locations was conducted, and trabecular density and cortical structural parameters were measured. Areal BMD of the lumbar vertebrae and hip was also measured.

RESULTS

Substantial changes in trabecular architecture of the distal tibia, in particular as they relate to topology of the network, were detected after 12 mo in the control group. S/C decreased 5.6% (p < 0.0005), and EI increased 7.1% (p < 0.0005). Most curve- and profile-type voxels (representative of trabecular struts), increased significantly (p < 0.001). Curve and profile edges resulting from disconnection of rod-like trabeculae increased by 9.8% and 5.1% (p = 0.0001 and <0.001, respectively). Similarly, DXA BMD in the spine and hip decreased 2.6% and 1.3% (p < 0.0001 and <0.005, respectively), and pQCT cortical area decreased 3.6% (p = 0.0001). However, neither trabecular density nor BV/TV changed. Furthermore, none of the parameters measured in the estradiol group were significantly different after 12 mo. Substantial differences in the mean changes from baseline between the estradiol treatment and control groups, in particular after 24 mo, were observed, with relative group differences as large as 13% (S/C, p = 0.005), and the relative changes in the two groups had the opposite sign for most parameters. The observed temporal alterations in architecture are consistent with remodeling changes that involve gradual conversion of plate-like to rod-like trabecular bone along with disconnection of trabecular elements, even in the absence of a net loss of trabecular bone. The high-resolution 3D rendered images provide direct evidence of the above remodeling changes in individual subjects. The radius structural data indicated similar trends but offered no definitive conclusions.

CONCLUSIONS

The short-term temporal changes in trabecular architecture after menopause, and the protective effects of estradiol ensuring maintenance of a more plate-like TB architecture, reported here, have not previously been observed in vivo. This work suggests that MRI-based in vivo micromorphometry of trabecular bone has promise as a tool for monitoring osteoporosis treatment.