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

The respective and dependent effects of scattering and bone matrix absorption on ultrasound attenuation in cortical bone

Cortical bone is characterized by a dense solid matrix permeated by fluid-filled pores. Ultrasound scattering has potential for the non-invasive evaluation of changes in bone porosity. However, there is an incomplete understanding of the impact of ultrasonic absorption in the solid matrix on ultrasound scattering. In this study, maps were derived from scanning acoustic microscopy images of human femur cross-sections. Finite-difference time domain ultrasound scatter simulations were conducted on these maps. Pore density, diameter distribution of the pores, and nominal absorption values in the solid and fluid matrices were controlled. Ultrasound pulses with a central frequency of 8.2 MHz were propagated, both in through-transmission and backscattering configurations. From these data, the scattering, bone matrix absorption, and attenuation extinction lengths were calculated. The results demonstrated that as absorption in the solid matrix was varied, the scattering, absorption, and attenuation extinction lengths were significantly impacted. It was shown that for lower values of absorption in the solid matrix (less than 2 dB mm-1), attenuation due to scattering dominates, whereas at higher values of absorption (more than 2 dB mm-1), attenuation due to absorption dominates. This will impact how ultrasound attenuation and scattering parameters can be used to extract quantitative information on bone microstructure.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy

BACKGROUND AND OBJECTIVE

Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure.

METHODS

hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases.

RESULTS

The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role.

CONCLUSIONS

hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.

Does the bone mineral density measured on a preoperative CT scan before total hip arthroplasty reflect the bone's mechanical properties?

INTRODUCTION

No method exists to quantify the bone quality and factors that will ensure osteointegration of total hip arthroplasty (THA) implants. A preoperative CT scan can be used to evaluate the bone mineral density (BMD) when planning a THA procedure. The aim of this study was to validate BMD measurement as a marker of bone quality based on a preoperative CT scan.

HYPOTHESIS

BMD reflects the bone's mechanical properties for the purposes of preoperative THA planning.

METHODS

Patients who underwent primary THA for hip osteoarthritis or dysplasia with cementless implants and 3D preoperative plan were enrolled prospectively. The cortical BMD was calculated on CT scans used in the preoperative planning process. During the surgical procedure, the femoral head and neck were collected. These bone samples were subsequently scanned with a calibrated micro-CT scanner. The BMD was derived from the micro-CT scan and used as input for a finite element model to determine the bone's mechanical properties. Correlations between BMD, apparent moduli of elasticity and porosity were calculated.

RESULTS

The values of cortical BMD measured on the micro-CT and CT scan were significantly correlated (cc=0.52). The mean angular cortical BMD measured with the micro-CT scan was 1472.33mg/cm³ (SD: 357.53mg/cm³, 980.64-2830.6mg/cm³). There was no significant correlation between cortical BMD and the various apparent moduli of elasticity, except for E_yy and G_zy. Cortical BMD and porosity were inversely correlated with a Spearman coefficient of -0.41 (CI₉₅: [-0.71; -0.02], p=0.03). There was also an inverse correlation between the apparent moduli of elasticity (independent of their orientation) and porosity (p<0.01).

DISCUSSION

BMD provides information about porosity, which is a major factor when evaluating the bone's mechanical properties before THA.

LEVEL OF EVIDENCE

IV.

Cortex or cancellous-which is early for the decrease of bone content for vertebral body in health?

OBJECTIVE

To obtain the cortex and cancellous parameter of the vertebral bone of healthy subjects using QCT. To explore which is earlier or faster for bone loss with age.

MATERIALS AND METHODS

733 physical examiners underwent chest low-dose CT examination were recruited, from April 1, 2021 to October 1, 2021. QCT sequence was used to obtain the bone mineral density of T12-L2 vertebral body without additional radiation. The mass and area of vertebral cortex and cancellous at the central level of L2 vertebral body were measured. The age -related characteristics of vertebral cortex and cancellous between male and female was analyzed and compared.

RESULTS

The vBMD of T12-L2 vertebral body decreased with age. Significant differences were found in volumetric bone mineral density (vBMD) of T12-L2 vertebral body. For female, significant differences were found in bone content involving cortical mass, cancellous mass, cortical area, cancellous area, cortical mass/cancellous mass and cortical area/cancellous area in different age groups, respectively. The cortical mass decreased with age in female. The cancellous mass of female increased and then decreased with peak at 31-40 y. The cortical area of female decreased gradually before 71 y. The cancellous area of female increased and then decreased with peak at 51-60 y. The values of mass ratio and area ratio in female showed a slowly downward trend with age. Significant differences of bone content between non-menopausal and menopausal women were found except the cancellous mass. For male, no significant differences were found in all parameters of bone content.

CONCLUSION

The changes of vertebral BMD, bone content of cortex and cancellous have different characteristics in different age. The change of cortex in female maybe earlier and faster than that of cancellous, especially in menopausal women.

Estimation of Thickness and Speed of Sound for Transverse Cortical Bone Imaging Using Phase Aberration Correction Methods: An In Silico and Ex Vivo Validation Study

Delay-and-sum (DAS) beamforming of backscattered echoes is used for conventional ultrasound imaging. Although DAS beamforming is well suited for imaging in soft tissues, refraction, scattering, and absorption, porous mineralized tissues cause phase aberrations of reflected echoes and subsequent image degradation. The recently developed refraction corrected multi-focus technique uses subsequent focusing of waves at variable depths, the tracking of travel times of waves reflected from outer and inner cortical bone interfaces, the estimation of the shift needed to focus from one interface to another to determine cortical thickness (Ct.Th), and the speed of sound propagating in a radial bone direction (Ct.ν11). The method was validated previously in silico and ex vivo on plate shaped samples. The aim of this study was to correct phase aberration caused by bone geometry (i.e., curvature and tilt with respect to the transducer array) and intracortical pores for the multi-focus approach. The phase aberration correction methods are based on time delay estimation via bone geometry differences to flat bone plates and via the autocorrelation and cross correlation of the reflected ultrasound waves from the endosteal bone interface. We evaluate the multi-focus approach by incorporating the phase aberration correction methods by numerical simulation and one experiment on a human tibia bone, and analyze the precision and accuracy of measuring Ct.Th and Ct.ν11. Site-matched reference values of the cortical thickness of the human tibia bone were obtained from high-resolution peripheral computed tomography. The phase aberration correction methods resulted in a more precise (coefficient of variation of 5.7%) and accurate (root mean square error of 6.3%) estimation of Ct.Th, and a more precise (9.8%) and accurate (3.4%) Ct.ν11 estimation, than without any phase aberration correction. The developed multi-focus method including phase aberration corrections provides local estimations of both cortical thickness and sound velocity and is proposed as a biomarker of cortical bone quality with high clinical potential for the prevention of osteoporotic fractures.

Pore-Size Distribution and Frequency-Dependent Attenuation in Human Cortical Tibia Bone Discriminate Fragility Fractures in Postmenopausal Women With Low Bone Mineral Density

Osteoporosis is a disorder of bone remodeling leading to reduced bone mass, structural deterioration, and increased bone fragility. The established diagnosis is based on the measurement of areal bone mineral density by dual‐energy X‐ray absorptiometry (DXA), which poorly captures individual bone loss and structural decay. Enlarged cortical pores in the tibia have been proposed to indicate structural deterioration and reduced bone strength in the hip. Here, we report for the first time the in vivo assessment of the cortical pore‐size distribution together with frequency‐dependent attenuation at the anteromedial tibia midshaft by means of a novel ultrasonic cortical backscatter (CortBS) technology. We hypothesized that the CortBS parameters are associated with the occurrence of fragility fractures in postmenopausal women (n = 55). The discrimination performance was compared with those of DXA and high‐resolution peripheral computed tomography (HR‐pQCT). The results suggest a superior discrimination performance of CortBS (area under the receiver operating characteristic curve [AUC]: 0.69 ≤ AUC ≤ 0.75) compared with DXA (0.54 ≤ AUC ≤ 0.55) and a similar performance compared with HR‐pQCT (0.66 ≤ AUC ≤ 0.73). CortBS is the first quantitative bone imaging modality that can quantify microstructural tissue deteriorations in cortical bone, which occur during normal aging and the development of osteoporosis. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Structure-function relationships of the human vertebral endplate

BACKGROUND

Although deformation and fracture of the vertebral endplate have been implicated in spinal conditions such as vertebral fracture and disc degeneration, few biomechanical studies of this structure are available. The goal of this study was to quantify the mechanical behavior of the vertebral endplate.

METHODS

Eight-five rectangular specimens were dissected from the superior and/or inferior central endplates of human lumbar spine segments L1 to L4. Micro-computed tomography (μCT) imaging, four-point-bend testing, and ashing were performed to quantify the apparent elastic modulus and yield stress (modulus and yield stress, respectively, of the porous vertebral endplate), tissue yield stress (yield stress of the tissue of the vertebral endplate, excluding pores), ultimate strain, fracture strain, bone volume fraction (BV/TV), bone mineral density (BMD), and various measures of tissue density and composition (tissue mineral density, ash fraction, and ash density). Regression was used to assess the dependence of mechanical properties on density and composition.

RESULTS

Wide variations in elastic and failure properties, and in density and tissue composition, were observed. BMD and BV/TV were good predictors of many of the apparent-level mechanical properties, including modulus, yield stress, and in the case of the inferior vertebral endplate, failure strains. Similar values of the mechanical properties were noted between superior and inferior vertebral endplates. In contrast to the dependence of apparent stiffness and strength on BMD and BV/TV, none of the mechanical properties depended on any of the tissue-level density measurements.

CONCLUSION

The dependence of many of the mechanical properties of the vertebral endplate on BV/TV and BMD suggests possibilities for noninvasive assessment of how this region of the spine behaves during habitual and injurious loading. Further study of the nonmineral components of the endplate tissue is required to understand how the composition of this tissue may influence the overall mechanical behavior of the vertebral endplate.

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.

Design and properties of biomimetic irregular scaffolds for bone tissue engineering

The treatment of sizeable segmental bone defects remains a challenge encountered by surgeons. In addition to bone transplantation, porous scaffolds have become a common option. Although the mechanical and biological properties of porous scaffold have recently been the subject of intense research, pore irregularity as a critical characteristic has been poorly explored. Therefore, this study aimed to design an irregular biomimetic scaffold for use in bone tissue engineering applications. The irregular scaffold was based on the Voronoi tessellation method for similarity with the primary histomorphological indexes of bone (porosity, trabecular thickness, cortical bone thickness, and surface to volume ratio). Moreover, a new gradient method was adopted, in which porosity was maintained constant, and the strut diameter was changed to generate a gradient in the irregular scaffold. The permeability and stress concentration characteristics of the irregular scaffold were compared against three conventional scaffolds (the octet, body-centered cubic, pillar body-centered cubic). The results illustrated that the microstructure of the irregular scaffold could be controlled similarly to that of the cortical/cancellous bone unit. Simultaneously, a broad range of permeability was identified for the irregular scaffold, and gradient irregular scaffolds performed better in terms of both permeability and stress distribution than regular scaffolds. This study describes a novel method for the design of irregular scaffolds, which have good controllability and excellent permeability.

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.

Femur strength predictions by nonlinear homogenized voxel finite element models reflect the microarchitecture of the femoral neck

In the human femoral neck, the contribution of the cortical and trabecular architecture to mechanical strength is known to depend on the load direction. In this work, we investigate if QCT-derived homogenized voxel finite element (hvFE) simulations of varying hip loading conditions can be used to study the architecture of the femoral neck. The strength of 19 pairs of human femora was measured ex vivo using nonlinear hvFE models derived from high-resolution peripheral QCT scans (voxel size: 30.3 µm). Standing and side-backwards falling loads were modeled. Quasi-static mechanical tests were performed on 20 bones for comparison. Associations of femur strength with volumetric bone mineral density (vBMD) or microstructural parameters of the femoral neck obtained from high-resolution QCT were compared between mechanical tests and simulations and between standing and falling loads. Proximal femur strength predictions by hvFE models were positively associated with the vBMD of the femoral neck (R² > 0.61, p < 0.001), as well as with its cortical thickness (R² > 0.27, p < 0.001), trabecular bone volume fraction (R² = 0.42, p < 0.001) and with the first two principal components of the femoral neck architecture (R² > 0.38, p < 0.001). Associations between femur strength and femoral neck microarchitecture were stronger for one-legged standing than for side-backwards falling. For both loading directions, associations between structural parameters and femur strength from hvFE models were in good agreement with those from mechanical tests. This suggests that hvFE models can reflect the load-direction-specific contribution of the femoral neck microarchitecture to femur strength.

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.

Cortical porosity assessment in the distal radius: A comparison of HR-pQCT measures with Synchrotron-Radiation micro-CT-based measures

OBJECTIVE

To determine the agreement between cortical porosity derived from high resolution peripheral quantitative computed tomography (HR-pQCT) (via standard threshold, mean density and density inhomogeneity methods) and synchrotron radiation micro-CT (SR-μCT) derived porosity at the distal radius.

METHODS

We scanned 10 cadaveric radii (mean donor age: 79, SD 11 years) at the standard distal region using HR-pQCT and SR-μCT at voxel sizes of 82 μm and 17.7 μm, respectively. Common cortical regions were delineated for each specimen in both imaging modalities. HR-pQCT images were analyzed for cortical porosity using the following methods: Standard threshold, mean density, and density inhomogeneity (via recommended and optimized equations). We assessed agreement in porosity measures between HR-pQCT methods and SR-μCT by reporting predicted variance from linear regression and mean bias with limits of agreement (LOA).

RESULTS

The standard threshold and mean density methods predicted 85% and 89% of variance and indicated underestimation (mean bias -9.1%, LOA -15.9% to -2.2%) and overestimation (10.4%, 4.6% to 16.2%) of porosity, respectively. The density inhomogeneity method with recommended equation predicted 89% of variance and mean bias of 14.9% (-4.3 to 34.2) with systematic over-estimation of porosity in more porous specimens. The density inhomogeneity method with optimized equation predicted 91% of variance without bias (0.0%, -5.3 to 5.2).

CONCLUSION

HR-pQCT imaged porosity assessed with the density inhomogeneity method with optimized equation indicated the best agreement with SR-μCT derived porosity.