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

Bone collagen tensile properties of the aging human proximal femur

Despite the dominant role of bone mass in osteoporotic fractures, aging bone tissue properties must be thoroughly understood to improve osteoporosis management. In this context, collagen content and integrity are considered important factors, although limited research has been conducted on the tensile behavior of demineralized compact bone in relation to its porosity and elastic properties in the native mineralized state. Therefore, this study aims (i) at examining the age-dependency of mineralized bone and collagen micromechanical properties; (ii) to test whether, and if so to which extent, collagen properties contribute to mineralized bone mechanical properties. Two cylindrical cortical bone samples from fresh frozen human anatomic donor material were extracted from 80 proximal diaphyseal sections from a cohort of 24 female and 19 male donors (57 to 96 years at death). One sample per section was tested in uniaxial tension under hydrated conditions. First, the native sample was tested elastically (0.25 % strain), and after demineralization, up to failure. Morphology and composition of the second specimen was assessed using micro-computed tomography, Raman spectroscopy, and gravimetric methods. Simple and multiple linear regression were employed to relate morphological, compositional, and mechanical variables with age and sex. Macro-tensile properties revealed that only elastic modulus of native samples was age dependent whereas apparent elastic modulus was sex dependent (p < 0.01). Compositional and morphological analysis detected a weak but significant age and sex dependency of relative mineral weight (r = -0.24, p < 0.05) and collagen disorder ratio (I∼1670/I∼1640, r = 0.25, p < 0.05) and a strong sex dependency of bone volume fraction while generally showing consistent results in mineral content assessment. Young's modulus of demineralized bone was significantly related to tissue mineral density and Young's modulus of native bone. The results indicate that mechanical properties of the organic phase, that include collagen and non-collagenous proteins, are independent of donor age. The observed reduction in relative mineral weight and corresponding overall stiffer response of the collagen network may be caused by a reduced number of mineral-collagen connections and a lack of extrafibrillar and intrafibrillar mineralization that induces a loss of waviness and a collagen fiber pre-stretch.

Bone Biomarkers Based on Magnetic Resonance Imaging

Magnetic resonance imaging (MRI) is increasingly used to evaluate the microstructural and compositional properties of bone. MRI-based biomarkers can characterize all major compartments of bone: organic, water, fat, and mineral components. However, with a short apparent spin-spin relaxation time (T2*), bone is invisible to conventional MRI sequences that use long echo times. To address this shortcoming, ultrashort echo time MRI sequences have been developed to provide direct imaging of bone and establish a set of MRI-based biomarkers sensitive to the structural and compositional changes of bone. This review article describes the MRI-based bone biomarkers representing total water, pore water, bound water, fat fraction, macromolecular fraction in the organic matrix, and surrogates for mineral density. MRI-based morphological bone imaging techniques are also briefly described.

3D T1rho sequences with FASE, UTE, and MAPSS acquisitions for knee evaluation

PURPOSE

For biochemical evaluation of soft tissues of the knee, T1rho magnetic resonance imaging (MRI) has been proposed. Purpose of this study was to compare three T1rho sequences based on fast advanced spin echo (FASE), ultrashort echo time (UTE), and magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots (MAPSS) acquisitions for the knee evaluation.

MATERIALS AND METHODS

We developed two T1rho sequences using 3D FASE or 3D radial UTE acquisitions. 3D MAPSS T1rho was provided by the manufacturer. Agarose phantoms with varying concentrations were imaged. Additionally, bilateral knees of asymptomatic subjects were imaged sagittally. T1rho values of the phantoms and 4 regions of interest (ROI) of the knees (i.e., anterior and posterior meniscus, femoral and tibial cartilage) were determined.

RESULTS

In phantoms, all T1rho values monotonically decreased with increasing agarose concentration. 3D MAPSS T1rho values of 51, 34, and 38 ms were found for 2, 3, and 4% agarose, respectively, similar to published values on another platform. In the knee, the raw images were detailed with good contrast. Cartilage and meniscus T1rho values varied with the pulse sequence, being the lowest in the 3D UTE T1rho sequence. Comparing different ROIs, menisci generally had lower T1rho values compared to cartilage, as expected in healthy knees.

CONCLUSION

We have successfully developed and implemented the new T1rho sequences and validated them using agarose phantoms and volunteer knees. All sequences were optimized to be clinically feasible (~ 5 min or less) and yielded satisfactory image quality and T1rho values consistent with the literature.

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.

An Update in Qualitative Imaging of Bone Using Ultrashort Echo Time Magnetic Resonance

Bone is comprised of mineral, collagenous organic matrix, and water. X-ray-based techniques are the standard approach for bone evaluation in clinics, but they are unable to detect the organic matrix and water components in bone. Magnetic resonance imaging (MRI) is being used increasingly for bone evaluation. While MRI can non-invasively assess the proton pools in soft tissues, cortical bone typically appears as a signal void with clinical MR techniques because of its short T2*. New MRI techniques have been recently developed to image bone while avoiding the ionizing radiation present in x-ray-based methods. Qualitative bone imaging can be achieved using ultrashort echo time (UTE), single inversion recovery UTE (IR-UTE), dual-inversion recovery UTE (Dual-IR-UTE), double-inversion recovery UTE (Double-IR-UTE), and zero echo time (ZTE) sequences. The contrast mechanisms as well as the advantages and disadvantages of each technique are discussed.

Emerging quantitative MR imaging biomarkers in inflammatory arthritides

PURPOSE

To review quantitative magnetic resonance imaging (qMRI) methods for imaging inflammation in connective tissues and the skeleton in inflammatory arthritis. This review is designed for a broad audience including radiologists, imaging technologists, rheumatologists and other healthcare professionals.

METHODS

We discuss the use of qMRI for imaging skeletal inflammation from both technical and clinical perspectives. We consider how qMRI can be targeted to specific aspects of the pathological process in synovium, cartilage, bone, tendons and entheses. Evidence for the various techniques from studies of both adults and children with inflammatory arthritis is reviewed and critically appraised.

RESULTS

qMRI has the potential to objectively identify, characterize and quantify inflammation of the connective tissues and skeleton in both adult and pediatric patients. Measurements of tissue properties derived using qMRI methods can serve as imaging biomarkers, which are potentially more reproducible and informative than conventional MRI methods. Several qMRI methods are nearing transition into clinical practice and may inform diagnosis and treatment decisions, with the potential to improve patient outcomes.

CONCLUSIONS

qMRI enables specific assessment of inflammation in synovium, cartilage, bone, tendons and entheses, and can facilitate a more consistent, personalized approach to diagnosis, characterisation and monitoring of disease.

Cortical bone vessel identification and quantification on contrast-enhanced MR images

Background

Cortical bone porosity is a major determinant of bone strength. Despite the biomechanical importance of cortical bone porosity, the biological drivers of cortical porosity are unknown. The content of cortical pore space can indicate pore expansion mechanisms; both of the primary components of pore space, vessels and adipocytes, have been implicated in pore expansion. Dynamic contrast-enhanced MRI (DCE-MRI) is widely used in vessel detection in cardiovascular studies, but has not been applied to visualize vessels within cortical bone. In this study, we have developed a multimodal DCE-MRI and high resolution peripheral QCT (HR-pQCT) acquisition and image processing pipeline to detect vessel-filled cortical bone pores.

Methods

For this in vivo human study, 19 volunteers (10 males and 9 females; mean age =63±5) were recruited. Both distal and ultra-distal regions of the non-dominant tibia were imaged by HR-pQCT (82 µm nominal resolution) for bone structure segmentation and by 3T DCE-MRI (Gadavist; 9 min scan time; temporal resolution =30 sec; voxel size 230×230×500 µm³) for vessel visualization. The DCE-MRI was registered to the HR-pQCT volume and the voxels within the MRI cortical bone region were extracted. Features of the DCE data were calculated and voxels were categorized by a 2-stage hierarchical kmeans clustering algorithm to determine which voxels represent vessels. Vessel volume fraction (volume ratio of vessels to cortical bone), vessel density (average vessel count per cortical bone volume), and average vessel volume (mean volume of vessels) were calculated to quantify the status of vessel-filled pores in cortical bone. To examine spatial resolution and perform validation, a virtual phantom with 5 channel sizes and an applied pseudo enhancement curve was processed through the proposed image processing pipeline. Overlap volume ratio and Dice coefficient was calculated to measure the similarity between the detected vessel map and ground truth.

Results

In the human study, mean vessel volume fraction was 2.2%±1.0%, mean vessel density was 0.68±0.27 vessel/mm³, and mean average vessel volume was 0.032±0.012 mm³/vessel. Signal intensity for detected vessel voxels increased during the scan, while signal for non-vessel voxels within pores did not enhance. In the validation phantom, channels with diameter 250 µm or greater were detected successfully, with volume ratio equal to 1 and Dice coefficient above 0.6. Both statistics decreased dramatically for channel sizes less than 250 µm.

Conclusions

We have a developed a multi-modal image acquisition and processing pipeline that successfully detects vessels within cortical bone pores. The performance of this technique degrades for vessel diameters below the in-plane spatial resolution of the DCE-MRI acquisition. This approach can be applied to investigate the biological systems associated with cortical pore expansion.

Computed tomography porosity and spherical indentation for determining cortical bone millimetre-scale mechanical properties

The cortex of the femoral neck is a key structural element of the human body, yet there is not a reliable metric for predicting the mechanical properties of the bone in this critical region. This study explored the use of a range of non-destructive metrics to measure femoral neck cortical bone stiffness at the millimetre length scale. A range of testing methods and imaging techniques were assessed for their ability to measure or predict the mechanical properties of cortical bone samples obtained from the femoral neck of hip replacement patients. Techniques that can potentially be applied in vivo to measure bone stiffness, including computed tomography (CT), bulk wave ultrasound (BWUS) and indentation, were compared against in vitro techniques, including compression testing, density measurements and resonant ultrasound spectroscopy. Porosity, as measured by micro-CT, correlated with femoral neck cortical bone's elastic modulus and ultimate compressive strength at the millimetre length scale. Large-tip spherical indentation also correlated with bone mechanical properties at this length scale but to a lesser extent. As the elastic mechanical properties of cortical bone correlated with porosity, we would recommend further development of technologies that can safely measure cortical porosity in vivo.

Three-dimensional ultrashort echo time imaging with tricomponent analysis for human cortical bone

PURPOSE

To investigate tricomponent analysis of human cortical bone using a multipeak fat signal model with 3D ultrashort TE Cones sequences on a clinical 3T scanner.

METHODS

Tricomponent fitting of bound water, pore water, and fat content using a multipeak fat spectra model was proposed for 3D ultrashort TE imaging of cortical bone. Three-dimensional ultrashort TE Cones acquisitions combined with tricomponent analysis were used to investigate bound and pore water T 2 ∗ and fractions, as well as fat T 2 ∗ and fraction in cortical bone. Feasibility studies were performed on 9 human cortical bone specimens with regions of interest selected from the endosteum to the periosteum in 4 circumferential regions. Microcomputed tomography studies were performed to measure bone porosity and bone mineral density for comparison and validation of the bound and pore water analyses.

RESULTS

The oscillation of the signal decay was well-fitted with the proposed tricomponent model. The sum of the pore water and fat fractions from tricomponent analysis showed a high correlation with microcomputed tomography porosity (R = 0.74, P < 0.01). Estimated bound-water fraction also demonstrated a high correlation with bone mineral density (R = 0.70, P < 0.01).

CONCLUSION

Tricomponent analysis significantly improves the estimation of bound-water and pore-water fractions in human cortical bone.

Manipulating the Amount and Structure of the Organic Matrix Affects the Water Compartments of Human Cortical Bone

Being predictors of the mechanical properties of human cortical bone, bound and pore water measurements by magnetic resonance (MR) imaging are being developed for the clinical assessment of fracture risk. While pore water is a surrogate of cortical bone porosity, the determinants of bound water are unknown. Manipulation of organic matrix properties by oxidative deproteinization, thermal denaturation, or nonenzymatic glycation lowers bone toughness. Because bound water contributes to bone toughness, we hypothesized that each of these matrix manipulations affect bound water fraction (V_bw/V_bone). Immersing cadaveric bone samples in sodium hypochlorite (NaClO) for 96 hours did not affect tissue mineral density or cortical porosity, but rather decreased V_bw/V_bone and increased short‐T₂ pore water signals as determined by ¹H nuclear MR relaxometry (¹H NMR). Moreover, the post treatment V_bw/V_bone linearly correlated with the remaining weight fraction of the organic matrix. Heating bone samples at 110°C, 120°C, 130°C, and then 140°C (∼24 hours per temperature and rehydration for ∼24 hours before ¹H NMR analysis) did not affect V_bw/V_bone. After subsequently heating them at 200°C, V_bw/V_bone increased. Boiling bone samples followed by heating at 110°C, 120°C, and then 130°C in water under pressure (8 hours per temperature) had a similar effect on V_bw/V_bone. Raman spectroscopy analysis confirmed that the increase in V_bw/V_bone coincided with an increase in an Amide I subpeak ratio that is sensitive to changes in the helical structure of collagen I. Glycation of bone by ribose for 4 weeks, but not in glucose for 16 weeks, decreased V_bw/V_bone, although the effect was less pronounced than that of oxidative deproteinization or thermal denaturation. We propose that MR measurements of bound water reflect the amount of bone organic matrix and can be modulated by collagen I helicity and by sugar‐derived post translational modifications of the matrix. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

The Role of Matrix Composition in the Mechanical Behavior of Bone

PURPOSE OF REVIEW

While thinning of the cortices or trabeculae weakens bone, age-related changes in matrix composition also lower fracture resistance. This review summarizes how the organic matrix, mineral phase, and water compartments influence the mechanical behavior of bone, thereby identifying characteristics important to fracture risk.

RECENT FINDINGS

In the synthesis of the organic matrix, tropocollagen experiences various post-translational modifications that facilitate a highly organized fibril of collagen I with a preferred orientation giving bone extensibility and several toughening mechanisms. Being a ceramic, mineral is brittle but increases the strength of bone as its content within the organic matrix increases. With time, hydroxyapatite-like crystals experience carbonate substitutions, the consequence of which remains to be understood. Water participates in hydrogen bonding with organic matrix and in electrostatic attractions with mineral phase, thereby providing stability to collagen-mineral interface and ductility to bone. Clinical tools sensitive to age- and disease-related changes in matrix composition that the affect mechanical behavior of bone could potentially improve fracture risk assessment.

Bone quality changes associated with aging and disease: a review

Bone quality encompasses all the characteristics of bone that, in addition to density, contribute to its resistance to fracture. In this review, we consider changes in architecture, porosity, and composition, including collagen structure, mineral composition, and crystal size. These factors all are known to vary with tissue and animal ages, and health status. Bone morphology and presence of microcracks, which also contribute to bone quality, will not be discussed in this review. Correlations with mechanical performance for collagen cross-linking, crystallinity, and carbonate content are contrasted with mineral content. Age-dependent changes in humans and rodents are discussed in relation to rodent models of disease. Examples are osteoporosis, osteomalacia, osteogenesis imperfecta (OI), and osteopetrosis in both humans and animal models. Each of these conditions, along with aging, is associated with increased fracture risk for distinct reasons.

Quantitative in vivo assessment of bone microarchitecture in the human knee using HR-pQCT

OBJECTIVE

High-resolution peripheral quantitative computed tomography (HR-pQCT) is a novel imaging modality capable of visualizing bone microarchitecture in vivo at human peripheral sites such as the distal radius and distal tibia. This research has extended the technology to provide a non-invasive assessment of bone microarchitecture at the human knee by establishing new hardware, imaging protocols and data analysis.

DESIGN

A custom leg holder was developed to stabilize a human knee centrally within a second generation HR-pQCT field of view. Five participants with anterior cruciate ligament reconstructions had their knee joint imaged in a continuous scan of 6cm axially. The nominal isotropic voxel size was 60.7μm. Bone mineral density and microarchitecture were assessed within the weight-bearing regions of medial and lateral compartments of the knee at three depths from the weight-bearing articular bone surface, including both the cortical and trabecular bone regions.

RESULTS

Scan duration was approximately 18min per knee and produced 5GB of projection data and 10GB of reconstructed image data (2304×2304 image matrix, 1008 slices). Motion during the scan was minimized by the leg holder and was similar in magnitude as a scan of the distal tibia. Bone mineral density and microarchitectural parameters were assessed for 16 volumes of interest in the tibiofemoral joint.

CONCLUSIONS

This is a new non-invasive in vivo assessment tool for bone microarchitecture in the human knee that provides an opportunity to gain insight into normal, injured and surgically reconstructed human knee bone architecture in cross-sectional or longitudinal studies.

Multiscale and multimodality computed tomography for cortical bone analysis

In clinical studies, high resolution peripheral quantitative computed tomography (HR-pQCT) is used to separately evaluate cortical bone and trabecular bone with an isotropic voxel of 82 µm³, and typical cortical parameters are cortical density (D.comp), thickness (Ct.Th), and porosity (Ct.Po). In vitro, micro-computed tomography (micro-CT) is used to explore the internal cortical bone micro-structure with isotropic voxels and high resolution synchrotron radiation (SR); micro-CT is considered the 'gold standard'. In 16 tibias and 8 femurs, HR-pQCT measurements were compared to conventional micro-CT measurements. To test modality effects, conventional micro-CT measurements were compared to SR micro-CT measurements at 7.5 µm³; SR micro-CT measurements were also tested at different voxel sizes for the femurs, specifically, 7.5 µm³ versus 2.8 µm³. D.comp (r  =  -0.88, p  <  10^-3) was the parameter best correlated with porosity (Po.V/TV). The correlation was not affected by the removal of pores under 130 µm. Ct.Th was also significantly highly correlated (r  =  -0.89 p  <  10^-3), while Ct.Po was correlated with its counterpart Po.V/TV (r  =  0.74, p  <  10^-3). From SR micro-CT and conventional micro-CT at 7.5 µm³ in matching areas, Po.V/TV and pore diameter were underestimated in conventional micro-CT with mean  ±  standard deviation (SD) biases of  -2.5  ±  1.9% and  -0.08  ±  0.08 mm, respectively. In contrast, pore number (Po.N) and pore separation (Po.Sp) were overestimated with mean  ±  SD biases of  +0.03  ±  0.04 mm^-1 and  +0.02  ±  0.04 mm, respectively. The results from the tibia and femur were similar when the results of SR micro-CT at 7.5 µm³ and 2.8 µm³ were compared. Po.V/TV, specific surface of pores (Po.S/Po.V), and Po.N were underestimated with mean biases of  -1.7  ±  0.9%, -4.6  ±  4.4 mm^-1, and  -0.26  ±  0.15 mm^-1, respectively. In contrast, pore spacing was overestimated at 7.5 µm³ compared to 2.8 µm³ with mean biases of 0.05  ±  0.03 mm. Cortical bone measurements from HR-pQCT images provided consistent results compared to those obtained using conventional micro-CT at the distal tibia. D.comp was highly correlated to Po.V/TV because it considers both the micro-porosity (Haversian systems) and macro-porosity (resorption lacunae) of cortical bone. The complexity of canal organization, (including shape, connectivity, and surface) are not fully considered in conventional micro-CT in relation to beam hardening and cone beam reconstruction artifacts. With the exception of Po.V/TV measurements, morphological and topological measurements depend on the characteristics of the x-ray beam, and to a lesser extent, on image resolution.

Cortical Bone Porosity: What Is It, Why Is It Important, and How Can We Detect It?

There is growing recognition of the role of micro-architecture in osteoporotic bone loss and fragility. This trend has been driven by advances in imaging technology, which have enabled a transition from measures of mass to micro-architecture. Imaging trabecular bone has been a key research focus, but advances in resolution have also enabled the detection of cortical bone micro-architecture, particularly the network of vascular canals, commonly referred to as 'cortical porosity.' This review aims to provide an overview of what this level of porosity is, why it is important, and how it can be characterized by imaging. Moving beyond a 'trabeculocentric' view of bone loss holds the potential to improve diagnosis and monitoring of interventions. Furthermore, cortical porosity is intimately linked to the remodeling process, which underpins bone loss, and thus a larger potential exists to improve our fundamental understanding of bone health through imaging of both humans and animal models.

Fast volumetric imaging of bound and pore water in cortical bone using three-dimensional ultrashort-TE (UTE) and inversion recovery UTE sequences

We report the three-dimensional ultrashort-TE (3D UTE) and adiabatic inversion recovery UTE (IR-UTE) sequences employing a radial trajectory with conical view ordering for bi-component T2 * analysis of bound water (T2 *(BW) ) and pore water (T2 *(PW) ) in cortical bone. An interleaved dual-echo 3D UTE acquisition scheme was developed for fast bi-component analysis of bound and pore water in cortical bone. A 3D IR-UTE acquisition scheme employing multiple spokes per IR was developed for bound water imaging. Two-dimensional UTE (2D UTE) and IR-UTE sequences were employed for comparison. The sequences were applied to bovine bone samples (n = 6) and volunteers (n = 6) using a 3-T scanner. Bi-component fitting of 3D UTE images of bovine samples showed a mean T2 *(BW) of 0.26 ± 0.04 ms and T2 *(PW) of 4.16 ± 0.35 ms, with fractions of 21.5 ± 3.6% and 78.5 ± 3.6%, respectively. The 3D IR-UTE signal showed a single-component decay with a mean T2 *(BW) of 0.29 ± 0.05 ms, suggesting selective imaging of bound water. Similar results were achieved with the 2D UTE and IR-UTE sequences. Bi-component fitting of 3D UTE images of the tibial midshafts of healthy volunteers showed a mean T2 *(BW) of 0.32 ± 0.08 ms and T2 *(PW) of 5.78 ± 1.24 ms, with fractions of 34.2 ± 7.4% and 65.8 ± 7.4%, respectively. Single-component fitting of 3D IR-UTE images showed a mean T2 *(BW) of 0.35 ± 0.09 ms. The 3D UTE and 3D IR-UTE techniques allow fast volumetric mapping of bound and pore water in cortical bone. Copyright © 2016 John Wiley & Sons, Ltd.

Separation of collagen-bound and porous bone water transverse relaxation in mice: proposal of a multi-step approach

The separation and quantification of collagen-bound water (CBW) and pore water (PW) components of the cortical bone signal are important because of their different contribution to bone mechanical properties. Ultrashort TE (UTE) imaging can be used to exploit the transverse relaxation from CBW and PW, allowing their quantification. We tested, for the first time, the feasibility of UTE measurements in mice for the separation and quantification of the transverse relaxation of CBW and PW in vivo using three different approaches for T2 * determination. UTE sequences were acquired at 4.7 T in six mice with 10 different TEs (50-5000 μs). The transverse relaxation time T2 * of CBW (T2 *cbw ) and PW (T2 *pw ) and the CBW fraction (bwf) were computed using a mono-exponential (i), a standard bi-exponential (ii) and a new multi-step bi-exponential (iii) approach. Regions of interest were drawn at multiple levels of the femur and vertebral body cortical bone for each mouse. The sum of the normalized squared residuals (Res) and the homogeneity of variance were tested to compare the different methods. In the femur, approach (i) yielded mean T2 * ± standard deviation (SD) of 657 ± 234 μs. With approach (ii), T2 *cbw , T2 *pw and bwf were 464 ± 153 μs, 15 777 ± 10 864 μs and 57.6 ± 9.9%, respectively. For approach (iii), T2 *cbw , T2 *pw and bwf were 387 ± 108 μs, 7534 ± 2765 μs and 42.5 ± 6.2%, respectively. Similar values were obtained from vertebral bodies. Res with approach (ii) was lower than with the two other approaches (p < 0.007), but T2 *pw and bwf variance was lower with approach (iii) than with approach (ii) (p < 0.048). We demonstrated that the separation and quantification of cortical bone water components with UTE sequences is feasible in vivo in mouse models. The direct bi-exponential approach exhibited the best approximation to the measured signal curve with the lowest residuals; however, the newly proposed multi-step algorithm resulted in substantially lower variability of the computed parameters. Copyright © 2016 John Wiley & Sons, Ltd.

Separation of collagen-bound and porous bone-water longitudinal relaxation in mice using a segmented inversion recovery zero-echo-time sequence

PURPOSE

Cortical bone mechanical properties are related to the collagen-bound water (CBW) and pore water (PW) components of cortical bone. The study evaluates the feasibility of zero-echo-time imaging in mice in vivo for longitudinal relaxation time (T1) measurements in cortical bone and separation of CBW and PW components.

METHODS

Zero-echo-time data were acquired at 4.7 Tesla in six mice with 14 different inversion times (0-2,600 ms). Region-of-interest analysis was performed at level of femur diaphysis. The T1 of cortical bone and of CBW (T1cbw) and PW (T1pw) as well as the CBW fraction (cbwf) was computed using a mono-exponential and a bi-exponential fitting approach, respectively. The sum of the squared residuals (Res) to the fit was provided for both approaches.

RESULTS

For the mono-exponential model, mean T1 ± standard deviation (SD) was 1,057 ± 160 ms. The bi-exponential approach provided a reliable separation of two different bone-water components, with a mean T1cbw of 213 ± 95 ms, T1pw of 2,152 ± 894 ms, and cbwf of 7.4 ± 2.7 %. Lower Res was obtained with bi-exponential approach (P < 0.001), and Res mean values ± SD were 0.016 ± 0.007 (bi-exponential) and 0.033 ± 0.016 (mono-exponential).

CONCLUSION

Zero-echo-time imaging allows for longitudinal relaxation measurements of cortical bone in vivo in mice models, with a reliable separation of PW and CBW components using a bi-exponential curve fitting approach. Magn Reson Med 77:1909-1915, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Evaluation of bound and pore water in cortical bone using ultrashort-TE MRI

Bone water exists in different states with the majority bound to the organic matrix and to mineral, and a smaller fraction in 'free' form in the pores of cortical bone. In this study, we aimed to develop and evaluate ultrashort-TE (UTE) MRI techniques for the assessment of T2*, T1 and concentration of collagen-bound and pore water in cortical bone using a 3-T clinical whole-body scanner. UTE MRI, together with an isotope study using tritiated and distilled water (THO-H2O) exchange, as well as gravimetric analysis, were performed on ten sectioned bovine bone samples. In addition, 32 human cortical bone samples were prepared for comparison between the pore water concentration measured with UTE MRI and the cortical porosity derived from micro-computed tomography (μCT). A short T2* of 0.27 ± 0.03 ms and T1 of 116 ± 6 ms were observed for collagen-bound water in bovine bone. A longer T2* of 1.84 ± 0.52 ms and T1 of 527 ± 28 ms were observed for pore water in bovine bone. UTE MRI measurements showed a pore water concentration of 4.7-5.3% by volume and collagen-bound water concentration of 15.7-17.9% in bovine bone. THO-H2O exchange studies showed a pore water concentration of 5.9 ± 0.6% and collagen-bound water concentration of 18.1 ± 2.1% in bovine bone. Gravimetric analysis showed a pore water concentration of 6.3 ± 0.8% and collagen-bound water concentration of 19.2 ± 3.6% in bovine bone. A mineral water concentration of 9.5 ± 0.6% was derived in bovine bone with the THO-H2O exchange study. UTE-measured pore water concentration is highly correlated (R(2) = 0.72, p < 0.0001) with μCT porosity in the human cortical bone study. Both bovine and human bone studies suggest that UTE sequences could reliably measure collagen-bound and pore water concentration in cortical bone using a clinical scanner.

The Role of Water Compartments in the Material Properties of Cortical Bone

Comprising ~20% of the volume, water is a key determinant of the mechanical behavior of cortical bone. It essentially exists in two general compartments: within pores and bound to the matrix. The amount of pore water-residing in the vascular-lacunar-canalicular space-primarily reflects intracortical porosity (i.e., open spaces within the matrix largely due to Haversian canals and resorption sites) and as such is inversely proportional to most mechanical properties of bone. Movement of water according to pressure gradients generated during dynamic loading likely confers hydraulic stiffening to the bone as well. Nonetheless, bound water is a primary contributor to the mechanical behavior of bone in that it is responsible for giving collagen the ability to confer ductility or plasticity to bone (i.e., allows deformation to continue once permanent damage begins to form in the matrix) and decreases with age along with fracture resistance. Thus, dehydration by air-drying or by solvents with less hydrogen bonding capacity causes bone to become brittle, but interestingly, it also increases stiffness and strength across the hierarchical levels of organization. Despite the importance of matrix hydration to fracture resistance, little is known about why bound water decreases with age in hydrated human bone. Using (1)H nuclear magnetic resonance (NMR), both bound and pore water concentrations in bone can be measured ex vivo because the proton relaxation times differ between the two water compartments, giving rise to two distinct signals. There are also emerging techniques to measure bound and pore water in vivo with magnetic resonance imaging (MRI). The NMR/MRI-derived bound water concentration is positively correlated with both the strength and toughness of hydrated bone and may become a useful clinical marker of fracture risk.

Bi-component T2 * analysis of bound and pore bone water fractions fails at high field strengths

Osteoporosis involves the degradation of the bone's trabecular architecture, cortical thinning and enlargement of cortical pores. Increased cortical porosity is a major cause of the decreased strength of osteoporotic bone. The majority of cortical pores, however, are below the resolution limit of MRI. Recent work has shown that porosity can be evaluated by MRI-based quantification of bone water. Bi-exponential T2 * fitting and adiabatic inversion preparation are the two most common methods purported to distinguish bound and pore water in order to quantify matrix density and porosity. To assess the viability of T2 * bi-component analysis as a method for the quantification of bound and pore water fractions, we applied this method to human cortical bone at 1.5, 3, 7 and 9.4 T, and validated the resulting pool fractions against micro-computed tomography-derived porosity and gravimetrically determined bone densities. We also investigated alternative methods: two-dimensional T1 -T2 * bi-component fitting by incorporation of saturation recovery, one- and two-dimensional fitting of Carr-Purcell-Meiboom-Gill (CPMG) echo amplitudes, and deuterium inversion recovery. The short-T2 * pool fraction was moderately correlated with porosity (R(2)  = 0.70) and matrix density (R(2)  = 0.63) at 1.5 T, but the strengths of these associations were found to diminish rapidly as the field strength increased, falling below R(2)  = 0.5 at 3 T. The addition of the T1 dimension to bi-component analysis only slightly improved the strengths of these correlations. T2 *-based bi-component analysis should therefore be used with caution. The performance of deuterium inversion recovery at 9.4 T was also poor (R(2)  = 0.50 vs porosity and R(2)  = 0.46 vs matrix density). The CPMG-derived short-T2 fraction at 9.4 T, however, was highly correlated with porosity (R(2)  = 0.87) and matrix density (R(2)  = 0.88), confirming the utility of this method for independent validation of bone water pools.

Identifying Novel Clinical Surrogates to Assess Human Bone Fracture Toughness

Fracture risk does not solely depend on strength but also on fracture toughness; ie, the ability of bone material to resist crack initiation and propagation. Because resistance to crack growth largely depends on bone properties at the tissue level, including collagen characteristics, current X-ray based assessment tools may not be suitable to identify age-related, disease-related, or treatment-related changes in fracture toughness. To identify useful clinical surrogates that could improve the assessment of fracture resistance, we investigated the potential of (1)H nuclear magnetic resonance spectroscopy (NMR) and reference point indentation (RPI) to explain age-related variance in fracture toughness. Harvested from cadaveric femurs (62 human donors), single-edge notched beam (SENB) specimens of cortical bone underwent fracture toughness testing (R-curve method). NMR-derived bound water showed the strongest correlation with fracture toughness properties (r = 0.63 for crack initiation, r = 0.35 for crack growth, and r = 0.45 for overall fracture toughness; p < 0.01). Multivariate analyses indicated that the age-related decrease in different fracture toughness properties were best explained by a combination of NMR properties including pore water and RPI-derived tissue stiffness with age as a significant covariate (adjusted R(2)  = 53.3%, 23.9%, and 35.2% for crack initiation, crack growth, and overall toughness, respectively; p < 0.001). These findings reflect the existence of many contributors to fracture toughness and emphasize the utility of a multimodal assessment of fracture resistance. Exploring the mechanistic origin of fracture toughness, glycation-mediated nonenzymatic collagen crosslinks and intracortical porosity are possible determinants of bone fracture toughness and could explain the sensitivity of NMR to changes in fracture toughness. Assuming fracture toughness is clinically important to the ability of bone to resist fracture, our results suggest that improvements in fracture risk assessment could potentially be achieved by accounting for water distribution (quantitative ultrashort echo time magnetic resonance imaging) and by a local measure of tissue resistance to indentation, RPI.