Sensitivity of qus parameters to controlled variations of bone strength assessed with a cellular model

A Combined Ultrasonic Backscatter Parameter for Bone Status Evaluation in Neonates

Metabolic bone disease (MBD) is one of the major complications of prematurity. Ultrasonic backscatter technique has the potential to be a portable and noninvasive method for early diagnosis of MBD. This study firstly applied CAS to neonates, which was defined as a linear combination of the apparent integrated backscatter coefficient (AIB) and spectral centroid shift (SCS). The objective was to evaluate the feasibility of ultrasonic backscatter technique for assessing neonatal bone health using AIB, SCS, and CAS. Ultrasonic backscatter measurements at 3.5 MHz, 5.0 MHz, and 7.5 MHz were performed on a total of 505 newborns within 48 hours after birth. The values of backscatter parameters were calculated and compared among gestational age groups. Correlations between backscatter parameters, gestational age, anthropometric indices, and biochemical markers were analyzed. The optimal predicting models for CAS were determined. The results showed term infants had lower SCS and higher AIB and CAS than preterm infants. Gestational age and anthropometric indices were negatively correlated with SCS (|r| = 0.45 – 0.57, P < 0.001), and positively correlated with AIB (|r| = 0.36 – 0.60, P < 0.001) and CAS (|r| = 0.56 – 0.69, P < 0.001). Biochemical markers yielded weak or nonsignificant correlations with backscatter parameters. CAS had relatively stronger correlations with the neonatal variables than AIB and SCS. At 3.5 MHz and 5.0 MHz, only gestational age (P < 0.001) independently contributed to the measurements of CAS, and could explain up to 40.5% – 44.3% of CAS variation. At 7.5 MHz, the combination of gestational age (P < 0.001), head circumference (P = 0.002), and serum calcium (P = 0.037) explained up to 40.3% of CAS variation. This study suggested ultrasonic backscatter technique was feasible to evaluate neonatal bone status. CAS was a promising parameter to provide more information about bone health than AIB or SCS alone.

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

Ultrasonic Backscatter Difference Measurement of Bone Health in Preterm and Term Newborns

Metabolic bone disease of prematurity remains a significant problem for preterm infants. Quantitative ultrasound (QUS) has potential as a non-invasive tool for assessing bone health of newborns. The aim of this study was to assess bone health in preterm and term newborns using ultrasonic backscatter difference measurement. This study analyzed a total of 493 neonates, including 239 full-term infants (gestational age [GA] >37 wk), 201 preterm I infants (GA: 32-37 wk) and 53 extreme preterm II infants (GA <32 wk). Ultrasonic backscatter measurements were performed on the calcaneus of infants at birth, and the normalized mean of the backscatter difference spectrum (nMBD) was calculated as an ultrasonic index of neonatal bone status. Simple and multiple linear regressions were performed to determine the association of ultrasonic nMBD with GA, anthropometric characteristics and biochemical markers. Statistically significant differences in GA, anthropometric characteristics (birth weight, birth length [BL], birth head circumference and body mass index [BMI]) and biochemical markers (alkaline phosphatase, serum calcium and serum phosphate) were observed among preterm and term infants. The nMBD for term infants (median = 3.72 dB/μs, interquartile range [IR] = 1.95 dB/μs) was significantly higher than that for preterm I infants (median = 1.95 dB/μs, IR = 3.12 dB/μs), which was, in turn, significantly higher than that for preterm II infants (median = 0.19 dB/μs, IR = 3.50 dB/μs). The nMBD yielded moderate correlations (ρ = 0.57-0.62, p < 0.001) with GA and anthropometric characteristics and weak correlations (|ρ| = 0.08-0.21, p < 0.001 or not significant) with biochemical markers. Multivariate regressions revealed that only BL (p = 0.002) and BMI (p = 0.032) yielded significantly independent contributions to the nMBD measurement, and combinations of BL and BMI could explain up to 42% of the variation of nMBD in newborn infants. This study found that ultrasonic backscatter difference measurement might be helpful in bone health evaluation in preterm and term newborns. The utility of ultrasonic backscatter measurement in diagnosis of metabolic bone disease in infants should be investigated further.

A μCT-based investigation of the influence of tissue modulus variation, anisotropy and inhomogeneity on ultrasound propagation in trabecular bone

Ultrasound propagation is widely used in the diagnosis of osteoporosis by providing information on bone mechanical quality. When it loses calcium, the tissue properties will first decrease. However, limited research about the influence of tissue properties on ultrasound propagation have been done due to the cumbersome experiment. The goal of this study was to explore the relationships between tissue modulus (Es) and speed of sound (SOS) through numerical simulations, and to study the influence of Es on the acoustical behavior in characterizing the local structural anisotropy and inhomogeneity. In this work, three-dimensional finite element (FE) simulations were performed on a cubic high-resolution (15μm) bovine trabecular bone sample (4×4×4mm(3), BV/TV=0.18) mapped from micro-computed tomography. Ultrasound excitations of 50kHz, 500kHz and 2MHz were applied in three orthogonal axes and the first arriving signal (FAS) was collected to quantify wave velocity. In this study, a strong power law relationship between Es and SOS was measured with estimated exponential index β=2.08-3.44 for proximal-distal (PD), anterior-posterior (AP) and medial-lateral (ML), respectively (all R(2)>0.95). For various Es, a positive dispersion of sound speed with respect to sound frequency was observed and the velocity dispersion magnitude (VDM) was measured. Also, with Es=15GPa in three orientations, the SOS in PD axis is 2009±120m/s, faster than that of AP (1762±106m/s) and ML (1798±132m/s) (f=2MHz) directions. Besides, the standard deviation of SOS increases with the sound frequency and the Es in all directions except for that at 50kHz. For the mechanical properties, the apparent modulus with certain Es was highest in the longitudinal direction compared with the transverse directions. It indicates that the tissue modulus combining with anisotropy and inhomogeneity has great influence on ultrasound propagation. Simulation results agree well with theoretical and experimental results.

Ultrasound Speed of Sound Measurements in Trabecular Bone Using the Echographic Response of a Metallic Pin

Bone quality is an important parameter in spine surgery, but its clinical assessment remains difficult. The aim of the work described here was to demonstrate in vitro the feasibility of employing quantitative ultrasound to retrieve bone mechanical properties using an echographic technique taking advantage of the presence of a metallic pin inserted in bone tissue. A metallic pin was inserted in bone tissue perpendicular to the transducer axis. The echographic response of the bone sample was determined, and the echo of the pin inserted in bone tissue and water were compared to determine speed of sound, which was compared with bone volume fraction. A 2-D finite-element model was developed to assess the effect of positioning errors. There was a significant correlation between speed of sound and bone volume fraction (R(2) = 0.6). The numerical results indicate the relative robustness of the measurement method, which could be useful to estimate bone quality intra-operatively.

Numerical Analysis of Ultrasound Backscattered Waves in Cancellous Bone Using a Finite-Difference Time-Domain Method: Isolation of the Backscattered Waves From Various Ranges of Bone Depths

Using a finite-difference time-domain method, ultrasound backscattered waves inside cancellous bone were numerically analyzed to investigate the backscatter mechanism. Two bone models with different thicknesses were modeled with artificial absorbing layers positioned at the back surfaces of the model, and an ultrasound pulse wave was transmitted toward the front surface. By calculating the difference between the simulated waveforms obtained using the two bone models, the backscattered waves from a limited range of depths in cancellous bone could be isolated. The results showed that the fast and slow longitudinal waves, which have previously been observed only in the ultrasound waveform transmitted through the bone, could be distinguished in the backscattered waveform from a deeper bone depth when transmitting the ultrasound wave parallel to the main orientation of the trabecular network. The amplitudes of the fast and slow backscattered waves were more closely correlated with the bone porosity [R2 = 0.84 and 0.66 (p < 0.001), respectively] than the amplitude of the whole (nonisolated) backscattered waves [R2 = 0.48 (p < 0.001)]. In conclusion, the nonisolated backscattered waves could be regarded as the superposition of the fast and slow waves reflected from various bone depths, returning at different times.

NUMERICAL SIMULATIONS OF CHANGE IN TRABECULAR STRUCTURE DUE TO BONE REMODELING UNDER ULTRASOUND PROPAGATION

Bone remodeling is defined as the coupling of bone formation and resorption on the bone surface. Numerical simulations of the remodeling in cancellous bone were performed to reproduce the change in the trabecular structure. Assuming that the formation/resorption in cancellous bone could be generated on the trabecular surface, where the local stress under the mechanical load was larger/smaller than the averaged stress on the surrounding surface, voxel trabecular elements in a numerical model of bovine cancellous bone were added/removed. An ultrasound continuous wave in the frequency range 0.1–1.0 MHz was applied as the mechanical load, and then, the local stress was analyzed using a finite-difference time-domain (FDTD) method. Using the remodeling simulations, both changes in the trabecular structure could be reproduced with decreasing and increasing porosity. In changes, the trabecular elements and the pore spaces became strongly oriented in the direction of ultrasound propagation. In addition, the remodeling simulations indicated that both bone formation and resorption lessened as the frequency increased.

Anatomical distribution of the degree of mineralization of bone tissue in human femoral neck: impact on biomechanical properties

Osteoporotic hip fractures represent a major public health problem associated with high human and economic costs. The anatomical variation of the tissue mineral density (TMD) and of the elastic constants in femoral neck cortical bone specimens is an important determinant of bone fragility. The purpose of this study was to show that a Synchrotron radiation microcomputed tomography system coupled with a multiscale biomechanical model allows the determination of the 3-D anatomical dependence of TMD and of the elastic constants (i.e. the mechanical properties of an anisotropic material) in human femoral neck. Bone specimens from the inferior femoral neck were obtained from 18 patients undergoing standard hemiarthroplasty. The specimens were imaged using 3-D synchrotron micro-computed tomography with a voxel size of 10.13 μm, leading to the determination of the anatomical distributions of porosity and TMD. The elastic properties of bone tissue were computed using a multiscale model. The model uses the experimental data obtained at the scale of several micrometers to estimate the components of the elastic tensor of bone at the scale of the organ. Statistical analysis (ANOVA) revealed a significant effect of the radial position on porosity and TMD and a significant effect of axial position on TMD only. Porosity was found to increase in the radial direction moving from the periosteum inwards (p<10(-5)). At any given distance from the periosteum, porosity does not vary noticeably along the bone axis. TMD was found to be significantly higher (p<10(-5)) in the periosteal region than in other bone locations and decreases from the periosteal to the endosteal region with an average slope of 10.05 g.cm(-3).m(-1), the decrease being faster in the porous part of the samples (average slope equal of 30.04 g.cm(-3).m(-1)) than in dense cortical bone. TMD was found to decrease from the distal to the proximal part of the femur neck (average slope of 6.5 g.cm(-3).m(-1)). Considering TMD variations in the radial direction induces weak changes of bone properties compared to constant TMD. TMD variations in the axial direction are responsible for a significant variation of elastic constants. These results demonstrate that the anatomical variations of TMD affect the bone elastic properties, which could be explained by the complex stress field in bone affecting bone remodeling. TMD spatial variations should be taken into account to properly describe the spatial heterogeneity of elastic coefficients of bone tissue at the organ scale.

Modeling and analysis of multiple scattering of acoustic waves in complex media: application to the trabecular bone

The integral equations that describe scattering in the media with step-rise changing parameters have been numerically solved for the trabecular bone model. The model consists of several hundred discrete randomly distributed elements. The spectral distribution of scattering coefficients in subsequent orders of scattering has been presented. Calculations were carried on for the ultrasonic frequency ranging from 0.5 to 3 MHz. Evaluation of the contribution of the first, second, and higher scattering orders to total scattering of the ultrasounds in trabecular bone was done. Contrary to the approaches that use the μCT images of trabecular structure to modeling of the ultrasonic wave propagation condition, the 3D numerical model consisting of cylindrical elements mimicking the spatial matrix of trabeculae, was applied. The scattering, due to interconnections between thick trabeculae, usually neglected in trabecular bone models, has been included in calculations when the structure backscatter was evaluated. Influence of the absorption in subsequent orders of scattering is also addressed. Results show that up to 1.5 MHz, the influence of higher scattering orders on the total scattered field characteristic can be neglected while for the higher frequencies, the relatively high amplitude interference peaks in higher scattering orders clearly occur.

Numerical investigation of ultrasound refraction caused by oblique orientation of trabecular network in cancellous bone

Ultrasound propagation through cancellous bone can be greatly affected by the trabecular structure. In the present study, the ultrasound propagation for the oblique orientation of the trabecular network was numerically investigated using 3-D finite-difference time-domain (FDTD) simulations. The models of cancellous bone were reconstructed from X-ray microcomputed tomographic (μCT) images of a bovine bone. Cancellous bone models with various orientations of the trabecular network were realized by cutting the μCT images rotated from 0 to 90°. Ultrasound waveforms propagating through these cancellous bone models were simulated while changing the receiving position. The refraction of the ultrasound wave for the oblique angle of the main orientation was investigated on the basis of the variation in the arrival time and peak amplitude. As the propagation direction approached the direction parallel to the main orientation, the arrival time of the first peak became less and the peak amplitude became smaller. This means that the wave of the first peak, which corresponded to a fast wave, propagated in the direction perpendicular to the main orientation. In addition, a strong correlation between the first-peak amplitude and the arrival time was observed in the porosity range of 0.68 to 0.85, in which the slope of the amplitude with respect to time increased linearly with porosity.

Determination of the heterogeneous anisotropic elastic properties of human femoral bone: from nanoscopic to organ scale

Cortical bone is a multiscale composite material. Its elastic properties are anisotropic and heterogeneous across its cross-section, due to endosteal bone resorption which might affect bone strength. The aim of this paper was to describe a homogenization method leading to the estimation of the variation of the elastic coefficients across the bone cross-section and along the bone longitudinal axis. The method uses the spatial variations of bone porosity and of the degree of mineralization of the bone matrix (DMB) obtained from the analysis of 3-D synchrotron micro-computed tomography images. For all three scales considered (the foam (100 nm), the ultrastructure (5 microm) and the mesoscale (500 microm)), the elastic coefficients were determined using the Eshelby's inclusion problem. DMB values were used at the scale of the foam. Collagen was introduced at the scale of the ultrastructure and bone porosity was introduced at the mesoscale. The pores were considered as parallel cylinders oriented along the bone axis. Each elastic coefficient was computed for different regions of interest, allowing an estimation of its variations across the bone cross-section and along the bone longitudinal axis. The method was applied to a human femoral neck bone specimen, which is a site of osteoporotic fracture. The computed elastic coefficients for cortical bone were in good agreement with experimental results, but some discrepancies were obtained in the endosteal part (trabecular bone). These results highlight the importance of accounting for the heterogeneity of cortical bone properties across bone cross-section and along bone longitudinal axis.

2-D finite difference time domain model of ultrasound reflection from normal and osteoarthritic human articular cartilage surface

Quantitative high-frequency ultrasonic evaluation of articular cartilage has shown a potential for the diagnosis of osteoarthritis, where the roughness of the surface, collagen and proteoglycan contents, and the density and mechanical properties of cartilage change concurrently. Experimentally, these factors are difficult to investigate individually and thus a numerical model is needed. The present study is the first one to use finite difference time domain modeling of pulse-echo measurements of articular cartilage. Ultrasound reflection from the surface was investigated with varying surface roughness, material parameters (Young's modulus, density, longitudinal, and transversal velocities) and inclination of the samples. The 2-D simulation results were compared with the results from experimental measurements of the same samples in an identical geometry. Both the roughness and the material parameters contributed significantly to the ultrasound reflection. The angular dependence of the ultrasound reflection was strong for a smooth cartilage surface but disappeared for the samples with a rougher surface. These results support the findings of previous experimental studies and indicate that ultrasound detects changes in the cartilage that are characteristic of osteoarthritis. In the present study there are differences between the results of the simulations and the experimental measurements. However, the systematic patterns in the experimental behavior are correctly reproduced by the model. In the future, our goal is to develop more realistic acoustic models incorporating inhomogeneity and anisotropy of the cartilage.

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

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

Numerical analysis of variability in ultrasound propagation properties induced by trabecular microstructure in cancellous bone

The manner by which the trabecular microstructure affects the propagation of ultrasound waves through cancellous bone is numerically investigated by finite difference time-domain (FDTD) simulation. Sixteen 3-D numerical models of 6.45 x 6.45 x 6.45 mm with a voxel size of 64.5 microm are reconstructed using a 3-D microcomputed tomographic (microCT) image taken from a bovine cancellous bone specimen of approximately 20 x 20 x 9 mm. All cancellous bone models have an oriented trabecular structure, and their trabecular elements are gradually eroded to increase the porosity using an image processing technique. Three erosion procedures are presented to realize various changes in the trabecular microstructure with increasing porosity. FDTD simulations of the ultrasound pulse waves propagating through the cancellous bone models at each eroded step are performed in 2 cases of the propagations parallel and perpendicular to the major trabecular orientation. The propagation properties of the wave amplitudes and propagation speeds are derived as a function of the porosity, and their variability due to the trabecular microstructure is revealed. To elucidate an effect of the microstructure, the mean intercept length (MIL), which is a microstructural parameter, is introduced, and the correlations of the propagation properties with the MILs of the trabecular elements and pore spaces are investigated.

Finite element simulation of ultrasound propagation in bone for quantitative ultrasound toward the diagnosis of osteoporosis

Osteoporosis is a serious bone disease which leads to the increased risk of bone fractures. For prevention and therapy, early detection of osteoporosis is critical. In general, for diagnosis of osteoporosis, dual-energy X-ray absoptiometry (DXA) or densitometry is most commonly used. However DXA exhibits some disadvantages such as ionizing radiation, relatively expensive cost, and limited information on mineralization and geometry of the bone. As an alternative method of DXA, quantitative ultrasound (QUS) is being investigated. In contrast to DXA, QUS is non-ionizing and relatively inexpensive. It can also provide some bone-related parameters (e.g., quantitative measurements including speed of sound and frequency-dependent attenuation). However the estimation of these parameters is difficult and few analytical solutions exist due to the complex behavior of ultrasound propagation in bone. As an alternative to the analytical methods, in most attempts, finite difference time domain (FDTD) method is used for simulation of ultrasound propagation in bone with a limited capability of modeling complex geometries of the bone. Finite element method (FEM) is a better solution since it can handle the complex geometry, but has been rarely applied due to its computational complexity. In this work, we propose an approach of FEM-based simulation of ultrasound propagation in bone. To validate our approach, we have tested simulated and real bone models from micro-CT using the index of speed-of-sound. Our results achieve an average of 97.54% in the computational accuracy.

A new quality of bone ultrasound research

Quantitative ultrasound (QUS) methods have strong power to predict osteoporotic fractures, but they are also very relevant for the assessment of bone quality. A representative sample of recent studies addressing these topics can be found in this special issue. Further pursuit of these methods will establish micro-QUS imaging methods as tools for measuring specific aspects of bone quality. Once this is achieved, we will be able to link such data to the clinical QUS methods used in vivo to determine which aspects of bone quality cause QUS to be a predictor of fracture risk that is independent of bone mineral density (BMD). Potentially this could lead to the development of a new generation of QUS devices for improved and expanded clinical assessment. Good quality of basic science work will thus lead to good quality of clinical patient examinations on the basis of a more detailed assessment of bone quality.