Assessment of anisotropic tissue elasticity of cortical bone from high-resolution, angular acoustic measurements

Quantitative Sound Speed Imaging of Cortical Bone and Soft Tissue: Results From Observational Data Sets

This work presents the first quantitative ultrasonic sound speed images of ex vivo limb cross-sections containing both soft tissue and bone using Full Waveform Inversion (FWI) with level set (LS) and travel time regularization. The estimated bulk sound speed of bone and soft tissue are within 10% and 1%, respectively, of ground truth estimates. The sound speed imagery shows muscle, connective tissue and bone features. Typically, ultrasound tomography (UST) using FWI is applied to imaging breast tissue properties (e.g. sound speed and density) that correlate with cancer. With further development, UST systems have the potential to deliver volumetric operator independent tissue property images of limbs with non-ionizing and portable hardware platforms. This work addresses the algorithmic challenges of imaging the sound speed of bone and soft tissue by combining FWI with LS regularization and travel time methods to recover soft tissue and bone sound speed with improved accuracy and reduced soft tissue artifacts when compared to conventional FWI. The value of leveraging LS and travel time methods is realized by evidence of improved bone geometry estimates as well as promising convergence properties and reduced risk of final model errors due to un-modeled shear wave propagation. Ex vivo bulk measurements of sound speed and MRI cross-sections validates the final inversion results.

Ultrasound Scattering in Cortical Bone

Recent advances in imaging of bone microstructure have led to a growing recognition of the role of cortical microstructure in osteoporosis. It is now accepted that the assessment of the microstructure of cortical porosity is essential to assess bone mechanical competence and predict fracture risk. Cortical porosity affects the propagation of ultrasound waves because pores act as ultrasound scatterers. Scattering by the porosity is an opportunity that should be leveraged to extract quantitative information about cortical microstructure. Scattering by the pores affects a number of ultrasound parameters that should be quantified, including attenuation, backscatter coefficient, ultrasound diffusivity, and their frequency dependence. Measuring these ultrasound parameters and developing models that describe their dependence upon parameters of cortical microstructure is the key to solve inverse problems that will allow the quantitative assessment of cortical porosity and ultimately will improve the non-invasive ultrasound-based evaluation of bone mechanical competence and fracture risk. In this chapter, we present recent advances in measuring and modeling those parameters in cortical bone.

Estimation of Cortical Bone Microstructure From Ultrasound Backscatter

Multichannel pulse-echo ultrasound using linear arrays and single-channel data acquisition systems opens new perspectives for the evaluation of cortical bone. In combination with spectral backscatter analysis, it can provide quantitative information about cortical microstructural properties. We present a numerical study, based on the finite-difference time-domain method, to estimate the backscatter cross section of randomly distributed circular pores in a bone matrix. A model that predicts the backscatter coefficient using arbitrary pore diameter distributions was derived. In an ex vivo study on 19 human tibia bones (six males, 13 females, 83.7 ± 8.4 years), multidirectional ultrasound backscatter measurements were performed using an ultrasound scanner equipped with a 6-MHz 128-element linear array with sweep motor control. A normalized depth-dependent spectral analysis was performed to derive backscatter and attenuation coefficients. Site-matched reference values of tissue acoustic impedance Z , cortical thickness (Ct.Th), pore density (Ct.Po.Dn), porosity (Ct.Po), and characteristic parameters of the pore diameter (Ct.Po.Dm) distribution were obtained from 100-MHz scanning-acoustic microscopy images. Proximal femur areal bone mineral density (aBMD), stiffness S , and ultimate force Fu from the same donors were available from a previous study. All pore structure and material properties could be predicted using linear combinations of backscatter parameters with a median to high accuracy (0.28 ≤ adjusted R² ≤ 0.59). The combination of cortical thickness and backscatter parameter provided similar or better prediction accuracies than aBMD. For the first time, a method for the noninvasive assessment of the pore diameter distribution in cortical bone by ultrasound is proposed. The combined assessment of cortical thickness, sound velocity, and pore size distribution in a mobile, nonionizing measurement system could have a major impact on preventing osteoporotic fractures.

Cortical thinning and accumulation of large cortical pores in the tibia reflect local structural deterioration of the femoral neck

INTRODUCTION

Cortical bone thinning and a rarefaction of the trabecular architecture represent possible causes of increased femoral neck (FN) fracture risk. Due to X-ray exposure limits, the bone microstructure is rarely measurable in the FN of subjects but can be assessed at the tibia. Here, we studied whether changes of the tibial cortical microstructure, which were previously reported to be associated with femur strength, are also associated with structural deteriorations of the femoral neck.

METHODS

The cortical and trabecular architectures in the FN of 19 humans were analyzed ex vivo on 3D microcomputed tomography images with 30.3 μm voxel size. Cortical thickness (Ct.Th_tibia), porosity (Ct.Po_tibia) and pore size distribution in the tibiae of the same subjects were measured using scanning acoustic microscopy (12 μm pixel size). Femur strength during sideways falls was simulated with homogenized voxel finite element models.

RESULTS

Femur strength was associated with the total (vBMD_tot; R² = 0.23, p < 0.01) and trabecular (vBMD_trab; R² = 0.26, p < 0.01) volumetric bone mineral density (vBMD), with the cortical thickness (Ct.Th_FN; R² = 0.29, p < 0.001) and with the trabecular bone volume fraction (Tb.BV/TV_FN; R² = 0.34, p < 0.001), separation (Tb.Sp_FN; R² = 0.25, p < 0.01) and number (Tb.N_FN; R² = 0.32, p < 0.001) of the femoral neck. Moreover, smaller Ct.Th_tibia was associated with smaller Ct.Th_FN (R² = 0.31, p < 0.05), lower Tb.BV/TV_FN (R² = 0.29, p < 0.05), higher Tb.Sp_FN (R² = 0.33, p < 0.05) and lower Tb.N_FN (R² = 0.42, p < 0.01). A higher prevalence of pores with diameter > 100 μm in tibial cortical bone (relCt.Po_100μm-tibia) indicated higher Tb.Sp_FN (R² = 0.36, p < 0.01) and lower Tb.N_FN (R² = 0.45, p < 0.01).

CONCLUSION

Bone resorption and structural decline of the femoral neck may be identified in vivo by measuring cortical bone thickness and large pores in the tibia.

Estimation of Thickness and Speed of Sound in Cortical Bone Using Multifocus Pulse-Echo Ultrasound

Most bone loss during the development of osteoporosis occurs in cortical bone at the peripheral skeleton. Decreased cortical thickness (Ct.Th) and the prevalence of large pores at the tibia are associated with reduced bone strength at the hip. Ct.Th and cortical sound velocity, i.e., a surrogate marker for changes of cortical porosity (Ct.Po), are key biomarkers for the identification of patients at high fracture risk. In this study, we have developed a method using a conventional ultrasound array transducer to determine thickness (Ct.Th) and the compressional sound velocity propagating in the radial bone direction (Ct. ν₁₁ ) using a refraction-corrected multifocus imaging approach. The method was validated in-silico on porous bone plate models using a 2-D finite-difference time-domain method and ex vivo on plate-shaped plastic reference materials and on plate-shaped cortical bovine tibia samples. Plane-wave pulse-echo measurements provided reference values to assess precision and accuracy of our method. In-silico results revealed the necessity to account for inclination-dependent transmission losses at the bone surface. Moreover, the dependence of Ct. ν₁₁ on both porosity and pore density was observed. Ct.Th and Ct. ν₁₁ obtained ex vivo showed a high correlation ) with reference values. The ex-vivo accuracy and precision for Ct. ν₁₁ were 29.9 m/s and 0.94%, respectively, and those for Ct.Th were 0.04 mm and 1.09%, respectively. In conclusion, this numerical and experimental study demonstrates an accurate and precise estimation of Ct.Th and Ct. ν₁₁ . The developed multifocus technique may have high clinical potential to improve fracture risk prediction using noninvasive and nonionizing conventional ultrasound technology with image guidance.

Acoustic diffusion constant of cortical bone: Numerical simulation study of the effect of pore size and pore density on multiple scattering

While osteoporosis assessment has long focused on the characterization of trabecular bone, the cortical bone micro-structure also provides relevant information on bone strength. This numerical study takes advantage of ultrasound multiple scattering in cortical bone to investigate the effect of pore size and pore density on the acoustic diffusion constant. Finite-difference time-domain simulations were conducted in cortical microstructures that were derived from acoustic microscopy images of human proximal femur cross sections and modified by controlling the density (Ct.Po.Dn) ∈[5-25] pore/mm² and size (Ct.Po.Dm) ∈[30-100] μm of the pores. Gaussian pulses were transmitted through the medium and the backscattered signals were recorded to obtain the backscattered intensity. The incoherent contribution of the backscattered intensity was extracted to give access to the diffusion constant D. At 8 MHz, significant differences in the diffusion constant were observed in media with different porous micro-architectures. The diffusion constant was monotonously influenced by either pore diameter or pore density. An increase in pore size and pore density resulted in a decrease in the diffusion constant (D =285.9Ct.Po.Dm^-1.49, R²=0.989 , p=4.96×10^-5,RMSE=0.06; D=6.91Ct.Po.Dn^-1.01, R²=0.94, p=2.8×10^-3 , RMSE=0.09), suggesting the potential of the proposed technique for the characterization of the cortical microarchitecture.

Large cortical bone pores in the tibia are associated with proximal femur strength

Alterations of structure and density of cortical bone are associated with fragility fractures and can be assessed in vivo in humans at the tibia. Bone remodeling deficits in aging women have been recently linked to an increase in size of cortical pores. In this ex vivo study, we characterized the cortical microarchitecture of 19 tibiae from human donors (aged 69 to 94 years) to address, whether this can reflect impairments of the mechanical competence of the proximal femur, i.e., a major fracture site in osteoporosis. Scanning acoustic microscopy (12 μm pixel size) provided reference microstructural measurements at the left tibia, while the bone vBMD at this site was obtained using microcomputed tomography (microCT). The areal bone mineral density of both left and right femoral necks (aBMD_neck) was measured by dual‐energy X‐ray absorptiometry (DXA), while homogenized nonlinear finite element models based on high-resolution peripheral quantitative computed tomography provided hip stiffness and strength for one-legged standing and sideways falling loads. Hip strength was associated with aBMD_neck (r = 0.74 to 0.78), with tibial cortical thickness (r = 0.81) and with measurements of the tibial cross-sectional geometry (r = 0.48 to 0.73) of the same leg. Tibial vBMD was associated with hip strength only for standing loads (r = 0.59 to 0.65). Cortical porosity (Ct.Po) of the tibia was not associated with any of the femoral parameters. However, the proportion of Ct.Po attributable to large pores (diameter > 100 μm) was associated with hip strength in both standing (r = -0.61) and falling (r = 0.48) conditions. When added to aBMD_neck, the prevalence of large pores could explain up to 17% of the femur ultimate force. In conclusion, microstructural characteristics of the tibia reflect hip strength as well as femoral DXA, but it remains to be tested whether such properties can be measured in vivo.

Ex vivo cortical porosity and thickness predictions at the tibia using full-spectrum ultrasonic guided-wave analysis

The estimation of cortical thickness (Ct.Th) and porosity (Ct.Po) at the tibia using axial transmission ultrasound was successfully validated ex vivo against site-matched micro-computed tomography. The assessment of cortical parameters based on full-spectrum guided-wave analysis might improve the prediction of bone fractures in a cost-effective and radiation-free manner.

PURPOSE

Cortical thickness (Ct.Th) and porosity (Ct.Po) are key parameters for the identification of patients with fragile bones. The main objective of this ex vivo study was to validate the measurement of Ct.Po and Ct.Th at the tibia using a non-ionizing, low-cost, and portable 500-kHz ultrasound axial transmission system. Additional ultrasonic velocities and site-matched reference parameters were included in the study to broaden the analysis.

METHODS

Guided waves were successfully measured ex vivo in 17 human tibiae using a novel 500-kHz bi-directional axial transmission probe. Theoretical dispersion curves of a transverse isotropic free plate model with invariant matrix stiffness were fitted to the experimental dispersion curves in order to estimate Ct.Th and Ct.Po. In addition, the velocities of the first arriving signal (υ_FAS) and A₀ mode (υ_A0) were measured. Reference Ct.Po, Ct.Th, and vBMD were obtained from site-matched micro-computed tomography. Scanning acoustic microscopy (SAM) provided the acoustic impedance of the axial cortical bone matrix.

RESULTS

The best predictions of Ct.Po (R² = 0.83, RMSE = 2.2%) and Ct.Th (R² = 0.92, RMSE = 0.2 mm, one outlier excluded) were obtained from the plate model. The second best predictors of Ct.Po and Ct.Th were vBMD (R² = 0.77, RMSE = 2.6%) and υ_A0 (R² = 0.28, RMSE = 0.67 mm), respectively.

CONCLUSIONS

Ct.Th and Ct.Po were accurately predicted at the human tibia ex vivo using a transverse isotropic free plate model with invariant matrix stiffness. The model-based predictions were not further enhanced when we accounted for variations in axial tissue stiffness as reflected by the acoustic impedance from SAM.

EXPLORING CORTICAL BONE DENSITY THROUGH THE ULTRASOUND INTEGRATED REFLECTION COEFFICIENT

Objective:

This work evaluates the relationship between ultrasonic reflection and bone density from fourteen cylindrical bovine cortical bone samples (3.0-cm thick).

Methods:

Twenty US reflection signals per sample were acquired along the bone surface (2.0-mm step). The Integrated Reflection Coefficient (IRC) from each signal was compared to Quantitative Computed Tomography (QCT).

Results:

Seven IRC and QCT curves presented Pearson's Correlation R-values above 0.5. For weak correlation curves, QCT and IRC showed similar trends in several segments.

Conclusion:

IRC was sensitive to bone density variation. Level of Evidence: Experimental Study, Investigating a Diagnostic Test.

BMD-based assessment of local porosity in human femoral cortical bone

Cortical pores are determinants of the elastic properties and of the ultimate strength of bone tissue. An increase of the overall cortical porosity (Ct.Po) as well as the local coalescence of large pores cause an impairment of the mechanical competence of bone. Therefore, Ct.Po represents a relevant target for identifying patients with high fracture risk. However, given their small size, the in vivo imaging of cortical pores remains challenging. The advent of modern high-resolution peripheral quantitative computed tomography (HR-pQCT) triggered new methods for the clinical assessment of Ct.Po at the peripheral skeleton, either by pore segmentation or by exploiting local bone mineral density (BMD). In this work, we compared BMD-based Ct.Po estimates with high-resolution reference values measured by scanning acoustic microscopy. A calibration rule to estimate local Ct.Po from BMD as assessed by HR-pQCT was derived experimentally. Within areas of interest smaller than 0.5 mm², our model was able to estimate the local Ct.Po with an error of 3.4%. The incorporation of the BMD inhomogeneity and of one parameter from the BMD distribution of the entire scan volume led to a relative reduction of the estimate error of 30%, if compared to an estimate based on the average BMD. When applied to the assessment of Ct.Po within entire cortical bone cross-sections, the proposed BMD-based method had better accuracy than measurements performed with a conventional threshold-based approach.

Multimodal correlative investigation of the interplaying micro-architecture, chemical composition and mechanical properties of human cortical bone tissue reveals predominant role of fibrillar organization in determining microelastic tissue properties

The mechanical competence of bone is crucially determined by its material composition and structural design. To investigate the interaction of the complex hierarchical architecture, the chemical composition and the resulting elastic properties of healthy femoral bone at the level of single bone lamellae and entire structural units, we combined polarized Raman spectroscopy (PRS), scanning acoustic microscopy (SAM) and synchrotron X-ray phase contrast nano tomography (SR-nanoCT). In line with earlier studies, mutual correlation analysis strongly suggested that the characteristic elastic modulations of bone lamellae within single units are the result of the twisting fibrillar orientation, rather than compositional variations, modulations of the mineral particle maturity, or mass density deviations. Furthermore, we show that predominant fibril orientations in entire tissue units can be rapidly assessed from Raman parameter maps. Coexisting twisted and oscillating fibril patterns were observed in all investigated tissue domains. Ultimately, our findings demonstrate in particular the potential of combined PRS and SAM measurements in providing multi-scalar analysis of correlated fundamental tissue properties. In future studies, the presented approach can be applied for non-destructive investigation of small pathologic samples from bone biopsies and a broad range of biological materials and tissues.

STATEMENT OF SIGNIFICANCE

Bone is a complex structured composite material consisting of collagen fibrils and mineral particles. Various studies have shown that not only composition, maturation, and packing of its components, but also their structural arrangement determine the mechanical performance of the tissue. However, prominent methodologies are usually not able to concurrently describe these factors on the micron scale and complementary tissue characterization remains challenging. In this study we combine X-ray nanoCT, polarized Raman imaging and scanning acoustic microscopy and propose a protocol for fast and easy assessment of predominant fibril orientations in bone. Based on our site-matched analysis of cortical bone, we conclude that the elastic modulations of bone lamellae are mainly determined by the fibril arrangement.

Differences in acoustic impedance of fresh and embedded human trabecular bone samples-Scanning acoustic microscopy and numerical evaluation

Trabecular bone samples are traditionally embedded and polished for scanning acoustic microscopy (SAM). The effect of sample processing, including dehydration, on the acoustic impedance of bone is unknown. In this study, acoustic impedance of human trabecular bone samples (n = 8) was experimentally assessed before (fresh) and after embedding using SAM and two-dimensional (2-D) finite-difference time domain simulations. Fresh samples were polished with sandpapers of different grit (P1000, P2500, and P4000). Experimental results indicated that acoustic impedance of samples increased significantly after embedding [mean values 3.7 MRayl (fresh), 6.1 MRayl (embedded), p < 0.001]. After polishing with different papers, no significant changes in acoustic impedance were found, even though higher mean values were detected after polishing with finer (P2500 and P4000) papers. A linear correlation (r = 0.854, p < 0.05) was found between the acoustic impedance values of embedded and fresh bone samples polished using P2500 SiC paper. In numerical simulations dehydration increased the acoustic impedance of trabecular bone (38%), whereas changes in surface roughness of bone had a minor effect on the acoustic impedance (-1.56%/0.1 μm). Thereby, the numerical simulations corroborated the experimental findings. In conclusion, acoustic impedance measurement of fresh trabecular bone is possible and may provide realistic material values similar to those of living bone.

Skeletal maturation substantially affects elastic tissue properties in the endosteal and periosteal regions of loaded mice tibiae

Although it is well known that the bone adapts to changes in the mechanical environment by forming and resorbing the bone matrix, little is known about the influence of mechanical loading on tissue material properties of the pre-existing and newly formed bone. In this study, we analyzed the newly formed and pre-existing tissue after two weeks of controlled in vivo axial compressive loading in tibia of young (10 week-old) and adult (26 week-old) female mice and compared to the control contralateral limb, by means of scanning acoustic microscopy. Additionally, we used quantitative backscattered electron imaging to determine the bone mineral density distribution within the newly formed and pre-existing bone of young mice. No significant differences were found in tissue stiffness or mineral density in the pre-existing bone tissue as a result of external loading. In the endosteal region, 10 and 26 week loaded animals showed a 9% reduction in bone tissue stiffness compared to control animals. An increase of 200% in the mineral apposition rate in this region was observed in both age groups. In the periosteal region, the reduction in bone tissue stiffness and the increase in bone mineral apposition rate as a result of loading were two times higher in the 10 compared to the 26 week old animals. These data suggest that, during growth and skeletal maturation, the response of bone to mechanical loading is a deposition of new bone matrix, where the tissue amount but not its mineral or elastic properties are influenced by animal age.

Distribution of mesoscale elastic properties and mass density in the human femoral shaft

Cortical bone properties are determined by tissue composition and structure at several hierarchical length scales. In this study, the spatial distribution of micro- and mesoscale elastic properties within a human femoral shaft has been investigated. Microscale tissue degree of mineralization (DMB), cortical vascular porosity Ct.Po and the average transverse isotropic stiffness tensor C(Micro) of cylindrical-shaped samples (diameter: 4.4 mm, N = 56) were obtained from cortical regions between 20 and 85% of the total femur length and around the periphery (anterior, medial, posterior and lateral quadrants) by means of synchrotron radiation µCT (SRµCT) and 50-MHz scanning acoustic microscopy (SAM). Within each cylinder, the volumetric bone mineral density (vBMD) and the mesoscale stiffness tensor C(Meso) were derived using a numerical homogenization approach. Moreover, microelastic maps of the axial elastic coefficient c33 measured by SAM at distinct cross-sectional locations along the femur were used to construct a 3-D multiscale elastic model of the femoral shaft. Variations of vBMD (6.1%) were much lower than the variations of mesoscale elastic coefficients (11.1-21.3%). The variation of DMB was only a minor predictor for variations of the mesoscale elastic properties (0.05 ≤ R(2) ≤ 0.34). Instead, variations of the mesoscale elastic properties could be explained by variations of cortical porosity and microscale elastic properties. These data were suitable inputs for numerical evaluations and may help to unravel the relations between structure and composition on the elastic function in cortical bone.

Development of new criteria for cortical bone histomorphometry in femoral neck: intra- and inter-observer reproducibility

Histomorphometry is commonly applied to study bone remodeling. Histological definitions of cortical bone boundaries have not been consistent. In this study, new criteria for specific definition of the transitional zone between the cortical and cancellous bone in the femoral neck were developed. The intra- and inter-observer reproducibility of this method was determined by quantitative histomorphometry and areal overlapping analysis. The undecalcified histological sections of femoral neck specimens (n = 6; from men aged 17-59 years) were processed and scanned to acquire histological images of complete bone sections. Specific criteria were applied to define histological boundaries. "Absolute cortex area" consisted of pure cortical bone tissue only, and was defined mainly based on the size of composite canals and their distance to an additional "guide" boundary (so-called "preliminary cortex boundary," the clear demarcation line of density between compact cortex and sparse trabeculae). Endocortical bone area was defined by recognizing characteristic endocortical structures adjacent to the preliminary cortical boundary. The present results suggested moderate to high reproducibility for low-magnification parameters (e.g., cortical bone area). The coefficient of variation (CV %) ranged from 0.02 to 5.61 in the intra-observer study and from 0.09 to 16.41 in the inter-observer study. However, the intra-observer reproducibility of some high-magnification parameters (e.g., osteoid perimeter/endocortical perimeter) was lower (CV %, 0.33-87.9). The overlapping of three histological areas in repeated analyses revealed highest intra- and inter-observer reproducibility for the absolute cortex area. This study provides specific criteria for the definition of histological boundaries for femoral neck bone specimens, which may aid more precise cortical bone histomorphometry.

Ultrasound to assess bone quality

Bone quality is determined by a variety of compositional, micro- and ultrastructural properties of the mineralized tissue matrix. In contrast to X-ray-based methods, the interaction of acoustic waves with bone tissue carries information about elastic and structural properties of the tissue. Quantitative ultrasound (QUS) methods represent powerful alternatives to ionizing x-ray based assessment of fracture risk. New in vivo applicable methods permit measurements of fracture-relevant properties, [eg, cortical thickness and stiffness at fragile anatomic regions (eg, the distal radius and the proximal femur)]. Experimentally, resonance ultrasound spectroscopy and acoustic microscopy can be used to assess the mesoscale stiffness tensor and elastic maps of the tissue matrix at microscale resolution, respectively. QUS methods, thus, currently represent the most promising approach for noninvasive assessment of components of fragility beyond bone mass and bone microstructure providing prospects for improved assessment of fracture risk.

Modeling of femoral neck cortical bone for the numerical simulation of ultrasound propagation

Quantitative ultrasound assessment of the cortical compartment of the femur neck (FN) is investigated with the goal of achieving enhanced fracture risk prediction. Measurements at the FN are influenced by bone size, shape and material properties. The work described here was aimed at determining which FN material properties have a significant impact on ultrasound propagation around 0.5 MHz and assessing the relevancy of different models. A methodology for the modeling of ultrasound propagation in the FN, with a focus on the modeling of bone elastic properties based on scanning acoustic microscopy data, is introduced. It is found that the first-arriving ultrasound signal measured in through-transmission at the FN is not influenced by trabecular bone properties or by the heterogeneities of the cortical bone mineralized matrix. In contrast, the signal is sensitive to variations in cortical porosity, which can, to a certain extent, be accounted for by effective properties calculated with the Mori-Tanaka method.

Alterations of mass density and 3D osteocyte lacunar properties in bisphosphonate-related osteonecrotic human jaw bone, a synchrotron µCT study

Osteonecrosis of the jaw, in association with bisphosphonates (BRONJ) used for treating osteoporosis or cancer, is a severe and most often irreversible side effect whose underlying pathophysiological mechanisms remain largely unknown. Osteocytes are involved in bone remodeling and mineralization where they orchestrate the delicate equilibrium between osteoclast and osteoblast activity and through the active process called osteocytic osteolysis. Here, we hypothesized that (i) changes of the mineralized tissue matrix play a substantial role in the pathogenesis of BRONJ, and (ii) the osteocyte lacunar morphology is altered in BRONJ. Synchrotron µCT with phase contrast is an appropriate tool for assessing both the 3D morphology of the osteocyte lacunae and the bone matrix mass density. Here, we used this technique to investigate the mass density distribution and 3D osteocyte lacunar properties at the sub-micrometer scale in human bone samples from the jaw, femur and tibia. First, we compared healthy human jaw bone to human tibia and femur in order to assess the specific differences and address potential explanations of why the jaw bone is exclusively targeted by the necrosis as a side effect of BP treatment. Second, we investigated the differences between BRONJ and control jaw bone samples to detect potential differences which could aid an improved understanding of the course of BRONJ. We found that the apparent mass density of jaw bone was significantly smaller compared to that of tibia, consistent with a higher bone turnover in the jaw bone. The variance of the lacunar volume distribution was significantly different depending on the anatomical site. The comparison between BRONJ and control jaw specimens revealed no significant increase in mineralization after BP. We found a significant decrease in osteocyte-lacunar density in the BRONJ group compared to the control jaw. Interestingly, the osteocyte-lacunar volume distribution was not altered after BP treatment.

Influence of porosity, pore size, and cortical thickness on the propagation of ultrasonic waves guided through the femoral neck cortex: a simulation study

The femoral neck is a common fracture site in elderly people. The cortical shell is thought to be the major contributor to the mechanical competence of the femoral neck, but its microstructural parameters are not sufficiently accessible under in vivo conditions with current X-ray-based methods. To systematically investigate the influences of pore size, porosity, and thickness of the femoral neck cortex on the propagation of ultrasound, we developed 96 different bone models (combining 6 different pore sizes with 4 different porosities and 4 different thicknesses) and simulated the ultrasound propagation using a finite-difference time-domain algorithm. The simulated single-element emitter and receiver array consisting of 16 elements (8 inferior and 8 superior) were placed at anterior and posterior sides of the bone, respectively (transverse transmission). From each simulation, we analyzed the waveform collected by each of the inferior receiver elements for the one with the shortest time of flight. The first arriving signal of this waveform, which is associated with the wave traveling through the cortical shell, was then evaluated for its three different waveform characteristics (TOF: time point of the first point of inflection of the received signal, Δt: difference between the time point at which the signal first crosses the zero baseline and TOF, and A: amplitude of the first extreme of the first arriving signal). From the analyses of these waveform characteristics, we were able to develop multivariate models to predict pore size, porosity, and cortical thickness, corresponding to the 96 different bone models, with remaining errors in the range of 50 μm for pore size, 1.5% for porosity, and 0.17 mm for cortical thickness.

Multiscale, converging defects of macro-porosity, microstructure and matrix mineralization impact long bone fragility in NF1

Bone fragility due to osteopenia, osteoporosis or debilitating focal skeletal dysplasias is a frequent observation in the Mendelian disease Neurofibromatosis type 1 (NF1). To determine the mechanisms underlying bone fragility in NF1 we analyzed two conditional mouse models, Nf1Prx1 (limb knock-out) and Nf1Col1 (osteoblast specific knock-out), as well as cortical bone samples from individuals with NF1. We examined mouse bone tissue with micro-computed tomography, qualitative and quantitative histology, mechanical tensile analysis, small-angle X-ray scattering (SAXS), energy dispersive X-ray spectroscopy (EDX), and scanning acoustic microscopy (SAM). In cortical bone of Nf1Prx1 mice we detected ectopic blood vessels that were associated with diaphyseal mineralization defects. Defective mineral binding in the proximity of blood vessels was most likely due to impaired bone collagen formation, as these areas were completely devoid of acidic matrix proteins and contained thin collagen fibers. Additionally, we found significantly reduced mechanical strength of the bone material, which was partially caused by increased osteocyte volume. Consistent with these observations, bone samples from individuals with NF1 and tibial dysplasia showed increased osteocyte lacuna volume. Reduced mechanical properties were associated with diminished matrix stiffness, as determined by SAM. In line with these observations, bone tissue from individuals with NF1 and tibial dysplasia showed heterogeneous mineralization and reduced collagen fiber thickness and packaging. Collectively, the data indicate that bone fragility in NF1 tibial dysplasia is partly due to an increased osteocyte-related micro-porosity, hypomineralization, a generalized defect of organic matrix formation, exacerbated in the regions of tensional and bending force integration, and finally persistence of ectopic blood vessels associated with localized macro-porotic bone lesions.

The early phases of bone healing can be differentiated in a rat osteotomy model by focused transverse-transmission ultrasound

Here we describe the use of a 5-MHz focused transmission system to image the bone repair region and to distinguish the early healing phases in a rat osteotomy (OT) model. Twelve-month-old female rats underwent a 2-mm OT. After 6 wk of consolidation, 2-D projection images of time-of-flight, speed of sound, and ultrasound attenuation were measured in vitro. The tissue types in the OT gap region were assessed by site-matched histology sections and micro-computed tomography (μCT). In the cases investigated, OT gap regions containing fibrous tissue (group A) were found to have similar properties compared with adjacent muscle tissue, whereas regions filled with cartilage and mineralized callus tissues (group B) differed significantly. Analysis of variance revealed that the healing group had a stronger effect on acoustic parameters (F < 35) than on μCT-based parameters (F < 22). This pilot study reports the feasibility of transverse transmission quantitative ultrasound in assessment of the onset of cartilage formation during callus formation.

Longitudinal elastic properties and porosity of cortical bone tissue vary with age in human proximal femur

Tissue level structural and mechanical properties are important determinants of bone strength. As an individual ages, microstructural changes occur in bone, e.g., trabeculae and cortex become thinner and porosity increases. However, it is not known how the elastic properties of bone change during aging. Bone tissue may lose its elasticity and become more brittle and prone to fractures as it ages. In the present study the age-dependent variation in the spatial distributions of microstructural and microelastic properties of the human femoral neck and shaft were evaluated by using acoustic microscopy. Although these properties may not be directly measured in vivo, there is a major interest to investigate their relationships with the linear elastic measurements obtained by diagnostic ultrasound at the most severe fracture sites, e.g., the femoral neck. However, before the validity of novel in vivo techniques can be established, it is essential to understand the age-dependent variation in tissue elastic properties and porosity at different skeletal sites. A total of 42 transverse cross-sectional bone samples were obtained from the femoral neck (Fn) and proximal femoral shaft (Ps) of 21 men (mean±SD age 47.1±17.8, range 17-82years). Samples were quantitatively imaged using a scanning acoustic microscope (SAM) equipped with a 50MHz ultrasound transducer. Distributions of the elastic coefficient (c33) of cortical (Ct) and trabecular (Tr) tissues and microstructure of cortex (cortical thickness Ct.Th and porosity Ct.Po) were determined. Variations in c33 were observed with respect to tissue type (c33Tr<c33Ct), location (c33(Ct.Ps)=37.7GPa>c33(Ct.Fn)=35.3GPa>c33(Tr.Ps)=33.8GPa>c33(Tr.Fn)=31.9GPa), and cadaver age (R(2)=0.28-0.46, p<0.05). Regional variations in porosity were found in the neck (superior 13.1%; inferior 6.1%; anterior 10.1%; posterior 8.6%) and in the shaft (medial 9.5%; lateral 7.7%; anterior 8.6%; posterior 12.0%). In conclusion, significant variations in elastic coefficients were detected between femoral neck and shaft as well as between the quadrants of the cross-sections of neck and shaft. Moreover, an age-related increase in cortical porosity and a stiffening of the bone tissue were observed. These findings may explain in part the increase in susceptibility to suffer low energy fractures during aging and highlight the potential of ultrasound in clinical osteoporosis diagnostics.

Microfibril orientation dominates the microelastic properties of human bone tissue at the lamellar length scale

The elastic properties of bone tissue determine the biomechanical behavior of bone at the organ level. It is now widely accepted that the nanoscale structure of bone plays an important role to determine the elastic properties at the tissue level. Hence, in addition to the mineral density, the structure and organization of the mineral nanoparticles and of the collagen microfibrils appear as potential key factors governing the elasticity. Many studies exist on the role of the organization of collagen microfibril and mineral nanocrystals in strongly remodeled bone. However, there is no direct experimental proof to support the theoretical calculations. Here, we provide such evidence through a novel approach combining several high resolution imaging techniques: scanning acoustic microscopy, quantitative scanning small-Angle X-ray scattering imaging and synchrotron radiation computed microtomography. We find that the periodic modulations of elasticity across osteonal bone are essentially determined by the orientation of the mineral nanoparticles and to a lesser extent only by the particle size and density. Based on the strong correlation between the orientation of the mineral nanoparticles and the collagen molecules, we conclude that the microfibril orientation is the main determinant of the observed undulations of microelastic properties in regions of constant mineralization in osteonal lamellar bone. This multimodal approach could be applied to a much broader range of fibrous biological materials for the purpose of biomimetic technologies.

Using the gradient of human cortical bone properties to determine age-related bone changes via ultrasonic guided waves

Bone fragility depends not only on bone mass but also on bone quality (structure and material). To accurately evaluate fracture risk or propose therapeutic treatment, clinicians need a criterion, which reflects the determinants of bone strength: geometry, structure and material. In human long bone, the changes due to aging, accentuated by osteoporosis are often revealed through the trabecularization of cortical bone, i.e., increased porosity of endosteal bone inducing a thinning of the cortex. Consequently, the intracortical porosity gradient corresponding to the spatial variation in porosity across the cortical thickness is representative of loss of mass, changes in geometry (thinning) and variations in structure (porosity). This article examines the gradient of material properties and its age-related evolution as a relevant parameter to assess bone geometry, structure and material. By applying a homogenization process, cortical bone can be considered as an anisotropic functionally graded material with variations in material properties. A semi-analytical method based on the sextic Stroh formalism is proposed to solve the wave equation in an anisotropic functionally graded waveguide for two geometries, a plate and a tube, without using a multilayered model to represent the structure. This method provides an analytical solution called the matricant and explicitly expressed under the Peano series expansion form. Our findings indicate that ultrasonic guided waves are sensitive to the age-related evolution of realistic gradients in human bone properties across the cortical thickness and have their place in a multimodal clinical protocol.

Relative contributions of porosity and mineralized matrix properties to the bulk axial ultrasonic wave velocity in human cortical bone

Velocity of ultrasound waves has proved to be a useful indicator of bone biomechanical competence. A detailed understanding of the dependence of ultrasound parameters such as velocity on bone characteristics is a key to the development of bone quantitative ultrasound (QUS). The objective of this study is to investigate the relative contributions of porosity and mineralized matrix properties to the bulk compressional wave velocity (BCV) along the long bone axis. Cross-sectional slabs from the diaphysis of four human femurs were included in the study. Seven regions of interest (ROIs) were selected in each slab. BCV was measured in through-transmission at 5 MHz. Impedance of the mineralized matrix (Z(m)) and porosity (Por) were obtained from 50 MHz scanning acoustic microscopy. Por and Z(m) had comparable effects on BCV along the bone axis (R=-0.57 and R=0.72, respectively).

Principal stiffness orientation and degree of anisotropy of human osteons based on nanoindentation in three distinct planes

Haversian systems or ‘osteons’ are cylindrical structures, formed by bone lamellae, that make up the major part of human cortical bone. Despite their discovery centuries ago in 1691 by Clopton Havers, their mechanical properties are still poorly understood.

The objective of this study is a detailed identification of the anisotropic elastic properties of the secondary osteon in the lamella plane. Additionally, the principal material orientation with respect to the osteon is assessed. Therefore a new nanoindentation method was developed which allows the measurement of indentation data in three distinct planes on a single osteon.

All investigated osteons appeared to be anisotropic with a preferred stiffness alignment along the axial direction with a small average helical winding around the osteon axis. The mean degree of anisotropy was 1.75 ± 0.36 and the mean helix angle was 10.3^°±0.8°.

These findings oppose two well established views of compact bone microstructure: first, the generally clear axial stiffness orientation contradicts a regular ‘twisted plywood’ collagen fibril orientation pattern in lamellar bone that would lead to a more isotropic behavior. Second, the class of transverse osteons were not observed from the mechanical point of view.

Spatial-temporal mapping of bone structural and elastic properties in a sheep model following osteotomy

The course of bone healing in animal models is conventionally monitored by morphologic approaches, which do not allow the determination of the material properties of the tissues involved. Mechanical characterization techniques are either dedicated to the macroscopic evaluation of the entire organ or to the microscopic evaluation of the tissue matrix. The latter provides insight to regionally specific alterations at the tissue level in the course of healing. In this study, quantitative scanning acoustic microscopy was used at 50 MHz to investigate microstructural and elastic alterations of mineralized callus and cortical tissue after transverse osteotomy in sheep tibiae. Analyses were performed after 2, 3, 6 and 9 weeks of consolidation with stabilization by either a rigid or a semi-rigid external fixator. Increased stiffness and decreased porosity were observed in the callus tissue over the course of the healing process, which was dependent on the fixator type. In the adjacent cortical tissue, stiffness decreased during the first 3 weeks. Cortical porosity increased over time but the time-dependence was different between the two fixator types. The changes of stiffness of cortical and callus tissues were measured with respect to the distance to the periosteal cortex-callus boundary. Stiffness of cortex and callus tissue smoothly decreased as a function of the distance from the inner cortical region. The data obtained in this study can help to understand the processes involved in tissue maturation during endogenous bone healing.

Ultrasound velocity and attenuation of porcine soft tissues with respect to structure and composition: I. Muscle

Ultrasound velocity and attenuation of soft tissues have been widely investigated. However, few studies completely covered considerable variations of both, structure and composition. The aim of this study was to collect acoustic reference data of porcine Longissimus muscle and associate them with compositional traits. In addition, measurements were conducted on fresh, formalin fixed, and frozen-thawed samples to evaluate the effect of processing on ultrasound parameters and comparisons with earlier investigations. Measurement conditions (temperature and fibre orientation) were realised close to hanging carcasses conditions. Sound velocity ranged from 1617 ± 6 to 1622 ± 5 ms(-1), while attenuation mostly ranged from 1.0 ± 0.3 to 1.2 ± 0.3 dB MHz(-1)cm(-1). Only formalin fixed samples showed significantly higher attenuation (2.2 ± 0.6 dB MHz(-1)cm(-1)). Highest correlations have been observed between intramuscular fat and attenuation (up to r = .7). The obtained results are anticipated to improve ultrasound based estimation of the intramuscular fat of pig muscle on intact carcasses.

Density imaging using a multiple-frequency DBIM approach

Current inverse scattering methods for quantitative density imaging have limitations that keep them from practical experimental implementations. In this work, an improved approach, termed the multiple-frequency distorted Born iterative method (MF-DBIM) algorithm, was developed for imaging density variations. The MF-DBIM approach consists of inverting the wave equation by solving for a single function that depends on both sound speed and density variations at multiple frequencies. Density information was isolated by using a linear combination of the reconstructed single-frequency profiles. Reconstructions of targets using MF-DBIM from simulated data were compared with reconstructions using methods currently available in the literature, i.e., the dual-frequency DBIM (DF-DBIM) and T-matrix approaches. Useful density reconstructions, i.e., root mean square errors (RMSEs) less than 30%, were obtained with MF-DBIM even with 2% Gaussian noise in the simulated data and using frequency ranges spanning less than an order of magnitude. Therefore, the MFDBIM approach outperformed both the DF-DBIM method (which has problems converging with noise even an order of magnitude smaller) and the T-matrix method (which requires a ka factor close to unity to achieve convergence). However, the convergence of all the density imaging algorithms was compromised when imaging targets with object functions exhibiting high spatial frequency content.

Adaptive remodeling of trabecular bone core cultured in 3-D bioreactor providing cyclic loading: an acoustic microscopy study

Scanning acoustic microscopy (SAM) provides high-resolution mapping of acoustic impedance related to tissue stiffness. This study investigates changes in tissue acoustic impedance resulting from mechanical loading in trabecular bone cores cultured in 3-D bioreactor. Trabecular bone cores were extracted from bovine sternum (n = 15) and ulna metaphysis (n = 15). From each bone, the samples were divided in three groups. The basal control (BC) group was fixed post-extraction, the control (C) and loaded (L) groups were maintained as viable in a controlled culture-loading cell over three weeks. Samples of L group underwent a dynamic compressive strain, whereas C samples were left free from loading. After three weeks, L and C samples were embedded in polymethylmethacrylate and all samples were explored with a 200-MHz SAM. For each specimen, the acoustic impedance distribution was obtained over flat and polished section of bone blocks prepared parallel to the loading axis. Our results showed that in basal controls, the acoustic impedance varied with bone anatomical location and was 15% higher in weight-bearing ulna compared with nonweight-bearing sternum. The comparison between loaded and nonloaded groups showed that sternum-only exhibited significant change in acoustic impedance (L vs. C sternum: +9%). This result suggests that when the applied load is comparable with the stress naturally experienced by a weight-bearing bone (ulna), the tissue material properties (manifested by acoustic impedance) remained unchanged. In conclusion, SAM is a potentially relevant tool for the assessment of subtle changes in intrinsic microelastic properties of bone induced by adaptive remodeling process in response to mechanical loading.

Microelastic imaging of bone

Several high-frequency ultrasound techniques have been developed during the last decade with the intention of assessing elastic properties of bone at the tissue level. The basic measurement principles can be divided into: 1) measurement of the compressional wave velocity in thin tissue sections; 2) measurement of surface acoustic wave velocities in thick sections; and 3) derivation of the acoustic impedance from the confocal reflection amplitude in thick sections. In this paper, the 3 principles are described with example measurements given in the frequency range from 50 MHz to 1.2 GHz. The measurements were made with 2 microscopes operating in the pulse-echo mode, either with frequencies up to 200 MHz and time-resolved detection or between 100 MHz and 2 GHz and amplitude detection. The methods are compared and their application potentials and limitations are discussed with respect to the hierarchical structure of cortical bone. Mapping of the confocal reflection amplitude has superior capabilities for deriving quantitative elastic and structural parameters in the heterogeneous bone material. Even at low frequencies (50 MHz), the mineralized tissue matrix can be separated from the larger pores (Haversian canals), and the elastic coefficient in the probing direction can be measured in 2 dimensions. Depending on the type of sample surface preparation (flat or cylindrically shaped), local distribution of a single elastic coefficient or the average transverse isotropic stiffness tensor can be derived. With frequencies in the GHz range, the lamellar bone structure can be analyzed. However, at one GHz, the acoustic wavelength is still one order of magnitude larger than the individual mineralized collagen fibrils. Although the thickness of a lamellar unit can easily be assessed from the acoustic image, the derivation of the anisotropic elastic properties of the mineralized collagen fibrils as well as the detailed structure of a lamella can only be accomplished with further model assumptions.

Spatial distribution of anisotropic acoustic impedance assessed by time-resolved 50-MHz scanning acoustic microscopy and its relation to porosity in human cortical bone

We used quantitative scanning acoustic microscopy (SAM) to assess tissue acoustic impedance and microstructure of cortical bone of human radii with the aim to provide data on regional distribution of acoustic impedance along the circumferential and across the radial directions in the entire cross-section of the radius diaphysis as well as to determine the range of impedance values in transverse (perpendicular to bone axis) and longitudinal (parallel to bone axis) cross-sections. Several microstructural features related to cortical porosity were analyzed in order to determine whether these features differ in different parts of the cortex and to assess the relationship between the microstructure and tissue acoustic impedance. Fifteen fresh bone specimens (human radius) were investigated using a SAM (center frequency of 50 MHz and -6 dB lateral resolution of approximately 23 microm). The sample acoustic impedance was obtained by means of a calibration curve correlating the reflected signal amplitude of reference materials with their corresponding well-known acoustic impedance. Tissue acoustic impedance and microstructural features were derived from the morphometric analysis of the segmented impedance images. A higher porosity was found in the inner cortical layer (mean+/-SD=8.9+/-2.3%) compared to the peripheral layer (2.7+/-1.5%) (paired t-test, p<10(-5)). ANOVA showed that most of the variance can be explained by the regional effect across the radial direction with a minor contribution due to between-sample variability. Similar to porosity, the number and diameter of pores were greater in the inner layer. In contrast to porosity, ANOVA showed that impedance variability can mostly be explained by between-specimen variability. Two-way ANOVA revealed that after compensation for the between-sample variability the variation in acoustic impedance across the radial direction was much larger than that along the circumferential direction. In addition to the significant difference between the inner cortical layer (8.25+/-0.4 Mrayl) and peripheral layer (8.0+/-0.5 Mrayl) (unilateral paired t-test, p<10(-4)), the values in the anterior region (8.2+/-0.5 Mrayl) were found to be significantly higher than those of the posterior region (7.9+/-0.6 Mrayl). Impedance mean value of longitudinal sections was lower than mean value measured in transverse cross-sections, resulting in an impedance acoustic anisotropy ratio of 1.17+/-0.03 in the inner cortical layer and 1.19+/-0.02 in the peripheral layer. SAM is a valuable tool to provide data on the spatial distribution of microstructural and microelastic bone properties that is useful to improve our understanding of the impact of bone microstructure on tissue material properties.

Prediction of biomechanical stability after callus distraction by high resolution scanning acoustic microscopy

Accurate clinical prediction of the resistance to fracture after callus distraction requires a detailed understanding of structural and elastic properties of the newly formed bone. We investigated 26 sheep that underwent middiaphyseal callus distraction at a rate of 0.5 mm every 12 h for 30 d using a standard unilateral fixator system. The sample population included four groups undergoing different treatments to improve bone healing, including bone grafting and the local application of growth factors. All animals were sacrificed eight weeks after the end of distraction. The fracture forces of the lengthened tibia and the contralateral control tibia from each animal were evaluated by biomechanical (four-point bending) testing. The microstructure and anisotropic acoustic impedance distributions were assessed by quantitative 50-MHz scanning acoustic microscopy. The relationships between resistance to fracture, structural properties and acoustic impedance of the newly formed callus tissue and adjacent cortical tissue were investigated. A significant linear multivariate regression model was developed that predicts the fracture force with a high accuracy (RMSE = 248 N, R(2) = 0.86, p < 0.0001).