Computationally-optimized bone mechanical modeling from high-resolution structural images

Computationally efficient dominant load-based local bone microstructure reconstruction method using topology optimization

The bone microstructure of the human proximal femur is clinically crucial for diagnosing skeletal pathologies, such as osteoporosis and bone metastases. The topology optimization-based bone microstructure method obtains these bone microstructures by converting low-resolution (LR) images into high-resolution images. However, this method is inherently computationally inefficient as it requires numerous finite elements, iterative analyses, and parallel computations. Therefore, this study proposes a novel topology optimization-based localised bone microstructure reconstruction method using the dominant load, which highly affects the selected region of interest (ROI), for efficient resolution enhancement. The load dependency of selected ROIs is quantified with a load dependency score. Then, the localised finite element model is constructed based on the local load estimation. Finally, the selected dominant load is applied as an input for the topology optimization-based bone microstructure reconstruction method. The reconstructed bone microstructure was similar to that of the conventional method. The localised finite element model applied by the dominant load effectively and accurately reconstructed the bone morphology and exhibited high computational efficiency. In conclusion, the dominant load-based approach can be used to construct a reasonable trabecular bone structure for ROI with high computational efficiency. The predictive performance of the proposed method was validated and showed promise for accurate trabecular bone structure prediction without additional radiation exposure.

Six-minute, in vivo MRI quantification of proximal femur trabecular bone 3D microstructure

BACKGROUND

Assessment of proximal femur trabecular bone microstructure in vivo by magnetic resonance imaging has recently been validated for acquiring information independent of bone mineral density in osteoporotic patients. However, the requisite signal-to-noise ratio (SNR) and resolution for interrogation of the trabecular microstructure at this anatomical location prolongs the scan duration and renders the imaging protocol clinically infeasible. Parallel imaging and compressed sensing (PICS) techniques can reduce the scan duration of the imaging protocol without substantially compromising image quality. The present work investigates the limits of acceleration for a commonly used PICS technique, ℓ1-ESPIRiT, for the purpose of quantifying measures of trabecular bone microarchitecture. Based on a desired error tolerance, a six-minute, prospectively accelerated variant of the imaging protocol was developed and assessed for intersession reproducibility and agreement with the longer reference scan.

PURPOSE

To investigate the limits of acceleration for MRI-based trabecular bone quantification by parallel imaging and compressed sensing reconstruction, and to develop a prototypical imaging protocol for assessing the proximal femur microstructure in a clinically practical scan time.

METHODS

Healthy participants (n = 11) were scanned by a 3D balanced steady-state free precession (bSSFP) sequence satisfying the Nyquist criterion with a scan duration of about 18 min. The raw data were retrospectively undersampled and reconstructed to mimic various acceleration factors ranging from 2 to 6. Trabecular volumes-of-interest in four major femoral regions (greater trochanter, intertrochanteric region, femoral neck, and femoral head) were analyzed and six relevant measures of trabecular bone microarchitecture (bone volume fraction, surface-to-curve ratio, erosion index, elastic modulus, trabecular thickness, plates-to-rods ratio) were obtained for images of all accelerations. To assess agreement, median percent error and intraclass correlation coefficients (ICCs) were computed using the fully-sampled data as reference. Based on this analysis, a prospectively 3-fold accelerated sequence with a duration of about 6 min was developed and the analysis was repeated.

RESULTS

A prospective acceleration factor of 3 demonstrated comparable performance in reproducibility and absolute agreement to the fully-sampled scan. The median CoV over all image-derived metrics was generally <6 % and ICCs >0.70. Also, measurements from prospectively 3-fold accelerated scans demonstrated in general median percent errors of <7 % and ICCs >0.70.

CONCLUSION

The present work proposes a method to make in vivo quantitative assessment of proximal femur trabecular microstructure with a clinically practical scan duration of about 6 min.

A Super-Resolution Diffusion Model for Recovering Bone Microstructure from CT Images

Purpose

To use a diffusion-based deep learning model to recover bone microstructure from low-resolution images of the proximal femur, a common site of traumatic osteoporotic fractures.

Materials and Methods

Training and testing data in this retrospective study consisted of high-resolution cadaveric micro-CT scans (n = 26), which served as ground truth. The images were downsampled prior to use for model training. The model was used to increase spatial resolution in these low-resolution images threefold, from 0.72 mm to 0.24 mm, sufficient to visualize bone microstructure. Model performance was validated using microstructural metrics and finite element simulation-derived stiffness of trabecular regions. Performance was also evaluated across a handful of image quality assessment metrics. Correlations between model performance and ground truth were assessed using intraclass correlation coefficients (ICCs) and Pearson correlation coefficients.

Results

Compared with popular deep learning baselines, the proposed model exhibited greater accuracy (mean ICC of proposed model, 0.92 vs ICC of next best method, 0.83) and lower bias (mean difference in means, 3.80% vs 10.00%, respectively) across the physiologic metrics. Two gradient-based image quality metrics strongly correlated with accuracy across structural and mechanical criteria (r > 0.89).

Conclusion

The proposed method may enable accurate measurements of bone structure and strength with a radiation dose on par with current clinical imaging protocols, improving the viability of clinical CT for assessing bone health.Keywords: CT, Image Postprocessing, Skeletal-Appendicular, Long Bones, Radiation Effects, Quantification, Prognosis, Semisupervised Learning Online supplemental material is available for this article. © RSNA, 2023.

MRI-based mechanical competence assessment of bone using micro finite element analysis (micro-FEA): Review

Areal bone mineral density (aBMD) from dual-energy x-ray absorptiometry (DEXA) and volumetric bone mineral density (vBMD) have demonstrated limited capabilities in the evaluation of bone mechanical competence and prediction of bone fracture. Predicting the macroscopic mechanical behavior of the bone structure has been challenging because of the heterogeneous and anisotropic nature of bone, such as the dependencies on loading direction, anatomical location, and sample dimensions. Magnetic resonance imaging (MRI) has been introduced as a promising modality that can be coupled with finite element analysis (FEA) for the assessment of bone mechanical competence. This review article describes studies investigating MRI-based micro-FEA as a potential non-invasive method to predict bone mechanical competence and facilitate bone fracture risk estimation without exposure to ionizing radiation. Specifically, the steps, applications, and future potential of FEA using indirect and direct bone imaging are discussed.

Heterogeneous material mapping methods for patient-specific finite element models of pelvic trabecular bone: A convergence study

Patient-specific finite element (FE) models of bone require the assignment of heterogeneous material properties extracted from the subject's computed tomography (CT) images. Though node-based (NB) and element-based (EB) material mapping methods (MMMs) have been proposed, the sensitivity and convergence of FE models to MMM for varying mesh sizes are not well understood. In this work, CT-derived and synthetic bone material data were used to evaluate the effect of MMM on results from FE analyses. Pelvic trabecular bone data was extracted from CT images of six subjects, while synthetic data were created to resemble trabecular bone properties. The numerical convergence of FE bone models using different MMMs were evaluated for strain energy, von-Mises stress, and strain. NB and EB MMMs both demonstrated good convergence regarding total strain energy, with the EB method having a slight edge over the NB. However, at the local level (e.g., maximum stress and strain), FE results were sensitive to the field type, MMM, and the FE mesh size. The EB method exhibited superior performance in finer meshes relative to the voxel size. The NB method converged better than did the EB method for coarser meshes. These findings may lead to higher-fidelity patient-specific FE bone models.

Effect of Low-Intensity Vibration on Bone Strength, Microstructure, and Adiposity in Pre-Osteoporotic Postmenopausal Women: A Randomized Placebo-Controlled Trial

There has been evidence that cyclical mechanical stimulation may be osteogenic, thus providing opportunities for nonpharmacological treatment of degenerative bone disease. Here, we applied this technology to a cohort of postmenopausal women with varying bone mineral density (BMD) T-scores at the total hip (-0.524 ± 0.843) and spine (-0.795 ± 1.03) to examine the response to intervention after 1 year of daily treatment with 10 minutes of vibration therapy in a randomized double-blinded trial. The device operates either in an active mode (30 Hz and 0.3 g) or placebo. Primary endpoints were changes in bone stiffness at the distal tibia and marrow adiposity of the vertebrae, based on 3 Tesla high-resolution MRI and spectroscopic imaging, respectively. Secondary outcome variables included distal tibial trabecular microstructural parameters and vertebral deformity determined by MRI, volumetric and areal bone densities derived using peripheral quantitative computed tomography (pQCT) of the tibia, and dual-energy X-ray absorptiometry (DXA)-based BMD of the hip and spine. Device adherence was 83% in the active group (n = 42) and 86% in the placebo group (n = 38) and did not differ between groups (p = .7). The mean 12-month changes in tibial stiffness in the treatment group and placebo group were +1.31 ± 6.05% and -2.55 ± 3.90%, respectively (group difference 3.86%, p = .0096). In the active group, marrow fat fraction significantly decreased after 12 months of intervention (p = .0003), whereas no significant change was observed in the placebo group (p = .7; group difference -1.59%, p = .029). Mean differences of the changes in trabecular bone volume fraction (p = .048) and erosion index (p = .044) were also significant, as was pQCT-derived trabecular volumetric BMD (vBMD; p = .016) at the tibia. The data are commensurate with the hypothesis that vibration therapy is protective against loss in mechanical strength and, further, that the intervention minimizes the shift from the osteoblastic to the adipocytic lineage of mesenchymal stem cells. © 2020 American Society for Bone and Mineral Research (ASBMR).

Quantitative Magnetic Resonance Imaging of Cortical and Trabecular Bone

Bone is a composite material consisting of mineral, organic matrix, and water. Water in bone can be categorized as bound water (BW), which is bound to bone mineral and organic matrix, or as pore water (PW), which resides in Haversian canals as well as in lacunae and canaliculi. Bone is generally classified into two types: cortical bone and trabecular bone. Cortical bone is much denser than trabecular bone that is surrounded by marrow and fat. Magnetic resonance (MR) imaging has been increasingly used for noninvasive assessment of both cortical bone and trabecular bone. Bone typically appears as a signal void with conventional MR sequences because of its short T2*. Ultrashort echo time (UTE) sequences with echo times 100 to 1,000 times shorter than those of conventional sequences allow direct imaging of BW and PW in bone. This article summarizes several quantitative MR techniques recently developed for bone evaluation. Specifically, we discuss the use of UTE and adiabatic inversion recovery prepared UTE sequences to quantify BW and PW, UTE magnetization transfer sequences to quantify collagen backbone protons, UTE quantitative susceptibility mapping sequences to assess bone mineral, and conventional sequences for high-resolution imaging of PW as well as the evaluation of trabecular bone architecture.

MRI-based assessment of proximal femur strength compared to mechanical testing

Half of the women who sustain a hip fracture would not qualify for osteoporosis treatment based on current DXA-estimated bone mineral density criteria. Therefore, a better approach is needed to determine if an individual is at risk of hip fracture from a fall. The objective of this study was to determine the association between radiation-free MRI-derived bone strength and strain simulations compared to results from direct mechanical testing of cadaveric femora. Imaging was conducted on a 3-Tesla MRI scanner using two sequences: one balanced steady-state free precession sequence with 300 μm isotropic voxel size and one spoiled gradient echo with anisotropic voxel size of 234 × 234 × 1500 μm. Femora were dissected free of soft-tissue and 4350-ohm strain-gauges were securely applied to surfaces at the femoral shaft, inferior neck, greater trochanter, and superior neck. Cadavers were mechanically tested with a hydraulic universal test frame to simulate loading in a sideways fall orientation. Sideways fall forces were simulated on MRI-based finite element meshes and bone stiffness, failure force, and force for plastic deformation were computed. Simulated bone strength metrics from the 300 μm isotropic sequence showed strong agreement with experimentally obtained values of bone strength, with stiffness (r = 0.88, p = 0.0002), plastic deformation point (r = 0.89, p < 0.0001), and failure force (r = 0.92, p < 0.0001). The anisotropic sequence showed similar trends for stiffness, plastic deformation point, and failure force (r = 0.68, 0.70, 0.84; p = 0.02, 0.01, 0.0006, respectively). Surface strain-gauge measurements showed moderate to strong agreement with simulated magnitude strain values at the greater trochanter, superior neck, and inferior neck (r = -0.97, -0.86, 0.80; p ≤0.0001, 0.003, 0.03, respectively). The findings from this study support the use of MRI-based FE analysis of the hip to reliably predict the mechanical competence of the human femur in clinical settings.

A multi-imaging modality study of bone density, bone structure and the muscle - bone unit in end-stage renal disease

End stage renal disease (ESRD) is associated with sarcopenia and skeletal fragility. The objectives of this cross-sectional study were to (1) characterize body composition, bone mineral density (BMD) and bone structure in hemodialysis patients compared with controls, (2) assess whether DXA areal BMD (aBMD) correlates with peripheral quantitative CT (pQCT) measures of volumetric BMD (vBMD), cortical dimensions and MRI measures of trabecular microarchitecture, and (3) determine the magnitude of bone deficits in ESRD after adjustment for muscle mass. Thirty ESRD participants, ages 25 to 64 years, were compared with 403 controls for DXA and pQCT outcomes and 104 controls for MRI outcomes; results were expressed as race- and sex- specific Z-scores relative to age. DXA appendicular lean mass index (ALMI kg/m2) and total hip, femoral neck, ultradistal and 1/3rd radius aBMD were significantly lower in ESRD, vs. controls (all p < 0.01). pQCT trabecular vBMD (p < 0.01), cortical vBMD (p < 0.001) and cortical thickness (due to a greater endosteal circumference, p < 0.02) and MRI measures of trabecular number, trabecular thickness, and whole bone stiffness were lower (all p < 0.01) in ESRD, vs. controls. ALMI was positively associated with total hip, femoral neck, ultradistal radius and 1/3rd radius aBMD and with tibia cortical thickness (R = 0.46 to 0.64). Adjustment for ALMI significantly attenuated bone deficits at these sites: e.g. mean femoral neck aBMD was 0.79 SD lower in ESRD, compared with controls and this was attenuated to 0.33 with adjustment for ALMI. In multivariate models within the dialysis participants, pQCT trabecular vBMD and cortical area Z-scores were significant and independently (all p < 0.02) associated with DXA femoral neck, total hip, and ultradistal radius aBMD Z-scores. Cortical vBMD (p = 0.01) and cortical area (p < 0.001) Z-scores were significantly and independently associated with 1/3rd radius areal aBMD Z-scores (R2 = 0.62). These data demonstrate that DXA aBMD captures deficits in trabecular and cortical vBMD and cortical area. The strong associations with ALMI, as an index of skeletal muscle, highlight the importance of considering the role of sarcopenia in skeletal fragility in patients with ESRD.

The reproducibility of measuring trabecular bone parameters using a commercially available high-resolution magnetic resonance imaging approach: A pilot study

Bone imaging is currently the best non-invasive way to assess changes to bone associated with aging or chronic disease. However, common imaging techniques such as dual energy x-ray absorptiometry are associated with limitations. Magnetic resonance imaging (MRI) is a radiation-free technique that can measure bone microarchitecture. However, published MRI bone assessment protocols use specialized MRI coils and sequences and therefore have limited transferability across institutions. We developed a protocol on a Siemens 3 Tesla MRI machine, using a commercially available coil (Siemens 15 CH knee coil), and manufacturer supplied sequences to acquire images at the tibia. We tested the reproducibility of the FSE and the GE Axial sequences and hypothesized that both would generate reproducible trabecular bone parameters. Eight healthy adults (age 25.5 ± 5.4 years) completed three measurements of each MRI sequence at the tibia. Each of the images was processed for 8 different bone parameters (such as volumetric bone volume fraction). We computed the coefficient of variation (CV) and intraclass correlation coefficients (ICC) to assess reproducibility and reliability. Both sequences resulted in trabecular parameters that were reproducible (CV <5% for most) and reliable (ICC >80% for all). Our study is one of the first to report that a commercially available MRI protocol can result in reproducible data, and is significant as MRI may be an accessible method to measure bone microarchitecture in clinical or research environments. This technique requires further testing, including validation and evaluation in other populations.

•

Trabecular bone is difficult to measure using commercial MRI techniques

•

Reproducibility of a MRI protocol measuring trabecular bone was assessed

•

Tibia trabecular bone was reproducible using a knee coil and a FSE Axial sequence

•

Tibia trabecular bone was reproducible using a knee coil and a GE Axial sequence

Accuracy of MRI-based finite element assessment of distal tibia compared to mechanical testing

High-resolution MRI-derived finite element analysis (FEA) has been used in translational research to estimate the mechanical competence of human bone. However, this method has yet to be validated adequately under in vivo imaging spatial resolution or signal-to-noise conditions. We therefore compared MRI-based metrics of bone strength to those obtained from direct, mechanical testing. The study was conducted on tibiae from 17 human donors (12 males and five females, aged 33 to 88years) with no medical history of conditions affecting bone mineral homeostasis. A 25mm segment from each distal tibia underwent MR imaging in a clinical 3-Tesla scanner using a fast large-angle spin-echo (FLASE) sequence at 0.137mm×0.137mm×0.410mm voxel size, in accordance with in vivo scanning protocol. The resulting high-resolution MR images were processed and used to generate bone volume fraction maps, which served as input for the micro-level FEA model. Simulated compression was applied to compute stiffness, yield strength, ultimate strength, modulus of resilience, and toughness, which were then compared to metrics obtained from mechanical testing. Moderate to strong positive correlations were found between computationally and experimentally derived values of stiffness (R²=0.77, p<0.0001), yield strength (R²=0.38, p=0.0082), ultimate strength (R²=0.40, p=0.0067), and resilience (R²=0.46, p=0.0026), but only a weak, albeit significant, correlation was found for toughness (R²=0.26, p=0.036). Furthermore, experimentally derived yield strength and ultimate strength were moderately correlated with MRI-derived stiffness (R²=0.48, p=0.0022 and R²=0.58, p=0.0004, respectively). These results suggest that high-resolution MRI-based finite element (FE) models are effective in assessing mechanical parameters of distal skeletal extremities.

The Efficacy of Low-intensity Vibration to Improve Bone Health in Patients with End-stage Renal Disease Is Highly Dependent on Compliance and Muscle Response

RATIONAL AND OBJECTIVES

Low intensity vibration (LIV) may represent a nondrug strategy to mitigate bone deficits in patients with end-stage renal disease.

MATERIALS AND METHODS

Thirty end-stage renal patients on maintenance hemodialysis were randomized to stand for 20 minutes each day on either an active or placebo LIV device. Analysis at baseline and completion of 6-month intervention included magnetic resonance imaging (tibia and fibula stiffness; trabecular thickness, number, separation, bone volume fraction, plate-to-rod ratio; and cortical bone porosity), dual-energy X-ray absorptiometry (hip and spine bone mineral density [BMD]), and peripheral quantitative computed tomography (tibia trabecular and cortical BMD; calf muscle cross-sectional area).

RESULTS

Intention-to-treat analysis did not show any significant changes in outcomes associated with LIV. Subjects using the active device and with greater than the median adherence (70%) demonstrated an increase in distal tibia stiffness (5.3%), trabecular number (1.7%), BMD (2.3%), and plate-to-rod ratio (6.5%), and a decrease in trabecular separation (-1.8%). Changes in calf muscle cross-sectional area were associated with changes in distal tibia stiffness (R = 0.85), trabecular bone volume/total volume (R = 0.91), number (R = 0.92), and separation (R = -0.94) in the active group but not in the placebo group. Baseline parathyroid hormone levels were positively associated with increased cortical bone porosity over the 6-month study period in the placebo group (R = 0.55) but not in the active group (R = 0.01). No changes were observed in the nondistal tibia locations for either group except a decrease in hip BMD in the placebo group (-1.7%).

CONCLUSION

Outcomes and adherence thresholds identified from this pilot study could guide future longitudinal studies involving vibration therapy.

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.

Patient-specific Hip Fracture Strength Assessment with Microstructural MR Imaging-based Finite Element Modeling

Purpose To describe a nonlinear finite element analysis method by using magnetic resonance (MR) images for the assessment of the mechanical competence of the hip and to demonstrate the reproducibility of the tool. Materials and Methods This prospective study received institutional review board approval and fully complied with HIPAA regulations for patient data. Written informed consent was obtained from all subjects. A nonlinear finite element analysis method was developed to estimate mechanical parameters that relate to hip fracture resistance by using MR images. Twenty-three women (mean age ± standard deviation, 61.7 years ± 13.8) were recruited from a single osteoporosis center. To thoroughly assess the reproducibility of the finite element method, three separate analyses were performed: a test-retest reproducibility analysis, where each of the first 13 subjects underwent MR imaging on three separate occasions to determine longitudinal variability, and an intra- and interoperator reproducibility analysis, where a single examination was performed in each of the next 10 subjects and four operators independently performed the analysis two times in each of the subjects. Reproducibility of parameters that reflect fracture resistance was assessed by using the intraclass correlation coefficient and the coefficient of variation. Results For test-retest reproducibility analysis and inter- and intraoperator analyses for proximal femur stiffness, yield strain, yield load, ultimate strain, ultimate load, resilience, and toughness in both stance and sideways-fall loading configurations each had an individual median coefficient of variation of less than 10%. Additionally, all measures had an intraclass correlation coefficient higher than 0.99. Conclusion This experiment demonstrates that the finite element analysis model can consistently and reliably provide fracture risk information on correctly segmented bone images. ^© RSNA, 2016 Online supplemental material is available for this article.

Estimation of Local Bone Loads for the Volume of Interest

Computational bone remodeling simulations have recently received significant attention with the aid of state-of-the-art high-resolution imaging modalities. They have been performed using localized finite element (FE) models rather than full FE models due to the excessive computational costs of full FE models. However, these localized bone remodeling simulations remain to be investigated in more depth. In particular, applying simplified loading conditions (e.g., uniform and unidirectional loads) to localized FE models have a severe limitation in a reliable subject-specific assessment. In order to effectively determine the physiological local bone loads for the volume of interest (VOI), this paper proposes a novel method of estimating the local loads when the global musculoskeletal loads are given. The proposed method is verified for the three VOI in a proximal femur in terms of force equilibrium, displacement field, and strain energy density (SED) distribution. The effect of the global load deviation on the local load estimation is also investigated by perturbing a hip joint contact force (HCF) in the femoral head. Deviation in force magnitude exhibits the greatest absolute changes in a SED distribution due to its own greatest deviation, whereas angular deviation perpendicular to a HCF provides the greatest relative change. With further in vivo force measurements and high-resolution clinical imaging modalities, the proposed method will contribute to the development of reliable patient-specific localized FE models, which can provide enhanced computational efficiency for iterative computing processes such as bone remodeling simulations.

An in silico parametric model of vertebrae trabecular bone based on density and microstructural parameters to assess risk of fracture in osteoporosis

Osteoporosis is a progressive bone disease characterized by deterioration in the quantity and quality of bone, leading to inferior mechanical properties and an increased risk of fracture. Current assessment of osteoporosis is typically based on bone densitometry tools such as Quantitative Computed Tomography (QCT) and Dual Energy X-ray absorptiometry (DEXA). These assessment modalities mainly rely on estimating the bone mineral density (BMD). Hence present densitometry tools describe only the deterioration of the quantity of bone associated with the disease and not the affected morphology or microstructural changes, resulting in potential incomplete assessment, many undetected patients, and unexplained fractures. In this study, an in-silico parametric model of vertebral trabecular bone incorporating both material and microstructural parameters was developed towards the accurate assessment of osteoporosis and the consequent risk of bone fracture. The model confirms that the mechanical properties such as strength and stiffness of vertebral trabecular tissue are highly influenced by material properties as well as morphology characteristics such as connectivity, which reflects the quality of connected inter-trabecular parts. The FE cellular solid model presented here provides a holistic approach that incorporates both material and microstructural elements associated with the degenerative process, and hence has the potential to provide clinical practitioners and researchers with more accurate assessment method for the degenerative changes leading to inferior mechanical properties and increased fracture risk associated with age and/or disease such as Osteoporosis.

Novel assessment tools for osteoporosis diagnosis and treatment

This review describes new technologies for the diagnosis and treatment, including fracture risk prediction, of postmenopausal osteoporosis. Four promising technologies and their potential for clinical translation and basic science studies are discussed. These include reference point indentation (RPI), Fourier transform infrared spectroscopy (FTIR), Raman spectroscopy, and magnetic resonance imaging (MRI). While each modality exploits different physical principles, the commonality is that none of them require use of ionizing radiation. To provide context for the new developments, brief summaries are provided for the current state of biomarker assays, fracture risk assessment (FRAX), and other fracture risk prediction algorithms and quantitative ultrasound (QUS) measurements.

Registration-based autofocusing technique for automatic correction of motion artifacts in time-series studies of high-resolution bone MRI

PURPOSE

To develop a registration-based autofocusing (RAF) motion correction technique for high-resolution trabecular bone (TB) imaging and to evaluate its performance on in vivo MR data.

MATERIALS AND METHODS

The technique combines serial registration with a previously developed motion correction technique - autofocusing - for automatic correction of subject movement degradation of MR images acquired in longitudinal studies. The method was tested on in vivo images of the distal radius to measure improvements in serial reproducibility of parameters in 12 women (ages 50-75 years), and to compare with the navigator echo-based correction and autofocusing. Furthermore, the technique's ability to optimize the sensitivity to detect simulated bone loss was ascertained.

RESULTS

The new technique yielded superior reproducibility of image-derived structural and mechanical parameters. Average coefficient of variation across all parameters improved by 12.5%, 27.0%, 33.5%, and 37.0%, respectively, following correction by navigator echoes, autofocusing, and the RAF technique (without and with correction for rotational motion); average intra-class correlation coefficient increased by 1.2%, 2.2%, 2.8%, and 3.2%, respectively. Furthermore, simulated bone loss (5%) was well recovered independent of the choice of reference image (4.71% or 4.86% with respect to using either the original or the image subjected to bone loss) in the time series.

CONCLUSION

The data suggest that our technique simultaneously corrects for intra-scan motion corruption while improving inter-scan registration. Furthermore, the technique is not biased by small changes in bone architecture between time-points.

Finite element analysis applied to 3-T MR imaging of proximal femur microarchitecture: lower bone strength in patients with fragility fractures compared with control subjects

PURPOSE

To determine the feasibility of using finite element analysis applied to 3-T magnetic resonance (MR) images of proximal femur microarchitecture for detection of lower bone strength in subjects with fragility fractures compared with control subjects without fractures.

MATERIALS AND METHODS

This prospective study was institutional review board approved and HIPAA compliant. Written informed consent was obtained. Postmenopausal women with (n = 22) and without (n = 22) fragility fractures were matched for age and body mass index. All subjects underwent standard dual-energy x-ray absorptiometry. Images of proximal femur microarchitecture were obtained by using a high-spatial-resolution three-dimensional fast low-angle shot sequence at 3 T. Finite element analysis was applied to compute elastic modulus as a measure of strength in the femoral head and neck, Ward triangle, greater trochanter, and intertrochanteric region. The Mann-Whitney test was used to compare bone mineral density T scores and elastic moduli between the groups. The relationship (R(2)) between elastic moduli and bone mineral density T scores was assessed.

RESULTS

Patients with fractures showed lower elastic modulus than did control subjects in all proximal femur regions (femoral head, 8.51-8.73 GPa vs 9.32-9.67 GPa; P = .04; femoral neck, 3.11-3.72 GPa vs 4.39-4.82 GPa; P = .04; Ward triangle, 1.85-2.21 GPa vs 3.98-4.13 GPa; P = .04; intertrochanteric region, 1.62-2.18 GPa vs 3.86-4.47 GPa; P = .006-.007; greater trochanter, 0.65-1.21 GPa vs 1.96-2.62 GPa; P = .01-.02), but no differences in bone mineral density T scores. There were weak relationships between elastic moduli and bone mineral density T scores in patients with fractures (R(2) = 0.25-0.31, P = .02-.04), but not in control subjects. CONCLUSION Finite element analysis applied to high-spatial-resolution 3-T MR images of proximal femur microarchitecture can allow detection of lower elastic modulus, a marker of bone strength, in subjects with fragility fractures compared with control subjects. MR assessment of proximal femur strength may provide information about bone quality that is not provided by dual-energy x-ray absorptiometry.

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.

Potential of in vivo MRI-based nonlinear finite-element analysis for the assessment of trabecular bone post-yield properties

PURPOSE

Bone strength is the key factor impacting fracture risk. Assessment of bone strength from high-resolution (HR) images have largely relied on linear micro-finite element analysis (μFEA) even though failure always occurs beyond the yield point, which is outside the linear regime. Nonlinear μFEA may therefore be more informative in predicting failure behavior. However, existing nonlinear models applied to trabecular bone (TB) have largely been confined to micro-computed tomography (μCT) and, more recently, HR peripheral quantitative computed tomography (HR-pQCT) images, and typically have ignored evaluation of the post-yield behavior. The primary purpose of this work was threefold: (1) to provide an improved algorithm and program to assess TB yield as well as post-yield properties; (2) to explore the potential benefits of nonlinear μFEA beyond its linear counterpart; and (3) to assess the feasibility and practicality of performing nonlinear analysis on desktop computers on the basis of micro-magnetic resonance (μMR) images obtained in vivo in patients.

METHODS

A method for nonlinear μFE modeling of TB yield as well as post-yield behavior has been designed where material nonlinearity is captured by adjusting the tissue modulus iteratively according to the tissue-level effective strain obtained from linear analysis using a computationally optimized algorithm. The software allows for images at in vivo μMRI resolution as input with retention of grayscale information. Associations between axial stiffness estimated from linear analysis and yield as well as post-yield parameters from nonlinear analysis were investigated from in vivo μMR images of the distal tibia (N = 20; ages: 58-84) and radius (N = 20; ages: 50-75).

RESULTS

All simulations were completed in 1 h or less for 61 strain levels using a desktop computer (dual quad-core Xeon 3.16 GHz CPUs equipped with 40 GB of RAM). Although yield stress and ultimate stress correlated strongly (R(2) > 0.95, p < 0.001) with axial stiffness, toughness correlated moderately at the distal tibia (R(2) = 0.81, p < 0.001) and only weakly at the distal radius (R(2) = 0.34, p = 0.007). Further, toughness was found to vary by up to 16% for bone of very similar axial stiffness (<2%).

CONCLUSIONS

The work demonstrates the practicality of nonlinear μFE simulations at in vivo μMRI resolution, as well as its potential for providing additional information beyond that obtainable from linear analysis. The data suggest that a direct assessment of toughness may provide information not captured by stiffness.

Assessment of trabecular bone yield and post-yield behavior from high-resolution MRI-based nonlinear finite element analysis at the distal radius of premenopausal and postmenopausal women susceptible to osteoporosis

RATIONALE AND OBJECTIVES

To assess the performance of a nonlinear microfinite element model on predicting trabecular bone yield and post-yield behavior based on high-resolution in vivo magnetic resonance images via the serial reproducibility.

MATERIALS AND METHODS

The nonlinear model captures material nonlinearity by iteratively adjusting tissue-level modulus based on tissue-level effective strain. It enables simulations of trabecular bone yield and post-yield behavior from micro magnetic resonance images at in vivo resolution by solving a series of nonlinear systems via an iterative algorithm on a desktop computer. Measures of mechanical competence (yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness) were estimated at the distal radius of premenopausal and postmenopausal women (N = 20, age range 50-75) in whom osteoporotic fractures typically occur. Each subject underwent three scans (20.2 ± 14.5 days). Serial reproducibility was evaluated via coefficient of variation (CV) and intraclass correlation coefficient (ICC).

RESULTS

Nonlinear simulations were completed in an average of 14 minutes per three-dimensional image data set involving analysis of 61 strain levels. The predicted yield strain/strength, ultimate strain/strength, modulus of resilience, and toughness had a mean value of 0.78%, 3.09 MPa, 1.35%, 3.48 MPa, 14.30 kPa, and 32.66 kPa, respectively, covering a substantial range by a factor of up to 4. Intraclass correlation coefficient ranged from 0.986 to 0.994 (average 0.991); CV ranged from 1.01% to 5.62% (average 3.6%), with yield strain and toughness having the lowest and highest CV values, respectively.

CONCLUSIONS

The data suggest that the yield and post-yield parameters have adequate reproducibility to evaluate treatment effects in interventional studies within short follow-up periods.

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.

Bioreactor design for perfusion-based, highly-vascularized organ regeneration

Bioartificial or laboratory-grown organs is a growing field centered on developing replacement organs and tissues to restore body function and providing a potential solution to the shortage of donor organs for transplantation. With the entry of engineered planar tissues, such as bladder and trachea, into clinical studies, an increasing focus is being given to designing complex, three-dimensional solid organs. As tissues become larger, thicker and more complex, the vascular network becomes crucial for supplying nutrients and maintaining viability and growth of the neo-organ. Perfusion decellularization, the process of removing cells from an entire organ, leaves the matrix of the vascular network intact. Organ engineering requires a delicate process of decellularization, sterilization, reseeding with appropriate cells, and organ maturation and stimulation to ensure optimal development. The design of bioreactors to facilitate this sequence of events has been refined to the extent that some bioartificial organs grown in these systems have been transplanted into recipient animals with sustained, though limited, function. This review focuses on the state-of-art in bioreactor development for perfusion-based bioartificial organs and highlights specific design components in need of further refinement.