Ultrasonic characterization of human trabecular bone microstructure

Frequency dependence of the ultrasonic power reflected from the water-tissue interface of human cancellous bone in vitro

Numerous studies have performed in vitro ultrasonic measurements of cancellous bone in water to develop techniques for ultrasonic bone assessment. Because cancellous bone is a highly porous medium, ultrasonic reflections at the water-bone interface may be frequency dependent. The goal of this study was to investigate the effect of porosity on the frequency dependence of the reflected power. Ultrasonic measurements were performed in a water tank at room temperature on 15 specimens of cancellous bone prepared from the proximal end of 9 human femurs using single element, broadband transducers with center frequencies of 3.5, 5, 7.5, and 10 MHz. Power spectra of pulses reflected from the water-specimen interface were corrected for the frequency response of the measurement system to obtain the reflected power in decibels R^dB(f). To suppress random phase cancellation effects, R^dB(f) was averaged over multiple sites on multiple specimens. A frequency dependence of R^dB(f) was observed in the 2.6-10 MHz range. The frequency dependence was moderate, with a maximum change of less than 6 dB over the entire frequency range. R^dB(f) was greatest for low porosity specimens. The frequency averaged intensity reflection coefficient ranged from 7.4 × 10^-4 to 7.8 × 10^-3 for high and low porosity specimen groups, respectively.

In Vivo Comparison of Backscatter Techniques for Ultrasonic Bone Assessment at the Femoral Neck

Ultrasonic techniques are being developed to detect changes in cancellous bone caused by osteoporosis. The goal of this study was to test the relative in vivo performance of eight backscatter parameters developed over the last several years for ultrasonic bone assessment: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), frequency intercept of apparent backscatter (FIAB), normalized mean of the backscatter difference (nMBD), normalized slope of the backscatter difference (nSBD), normalized intercept of the backscatter difference (nIBD), normalized backscatter amplitude ratio (nBAR) and backscatter amplitude decay constant (BADC). Backscatter measurements were performed on the left and right femoral necks of 80 adult volunteers (age = 25 ± 11 y) using an imaging system equipped with a convex array transducer. For comparison, additional ultrasonic measurements were performed at the left and right heel using a commercially available heel-bone ultrasonometer that measured the stiffness index. Six of the eight backscatter parameters (all but nSBD and nIBD) exhibited similar and highly significant (p < 0.000001) left-right correlations (0.51 ≤ R ≤ 0.68), indicating sensitivity to naturally occurring variations in bone tissue. Left-right correlations for the stiffness index measured at the heel (R = 0.75) were not significantly better than those produced by AIB, FSAB and FIAB. The short-term precisions of AIB, nMBD, nBAR and BADC (7.8%-11.7%) were comparable to that of the stiffness index measured with the heel-bone ultrasonometer (7.5%).

Diagnosis of Bone Mineral Density Based on Backscattering Resonance Phenomenon Using Coregistered Functional Laser Photoacoustic and Ultrasonic Probes

Dual-energy X-ray absorptiometry (DXA) machines based on bone mineral density (BMD) represent the gold standard for osteoporosis diagnosis and assessment of fracture risk, but bone strength and toughness are strongly correlated with bone collagen content (CC). Early detection of osteoporosis combined with BMD and CC will provide improved predictability for avoiding fracture risk. The backscattering resonance (BR) phenomenon is present in both ultrasound (US) and photoacoustic (PA) signal transmissions through bone, and the peak frequencies of BR can be changed with BM and CC. This phenomenon can be explained by the formation of standing waves within the pores. Simulations were then conducted for the same bone µCT images and the resulting resonance frequencies were found to match those predicted using the standing wave hypothesis. Experiments were performed on the same bone sample using an 808 nm wavelength laser as the PA source and 3.5 MHz ultrasonic transducer as the US source. The backscattering resonance effect was observed in the transmitted waves. These results verify our hypothesis that the backscattering resonance phenomenon is present in both US and PA signal transmissions and can be explained using the standing waves model, which will provide a suitable method for the early detection of osteoporosis.

Ultrasonic Bone Assessment: Ability of Apparent Backscatter Techniques to Detect Changes in the Microstructure of Human Cancellous Bone

Ultrasonic backscatter techniques may offer a useful approach for detecting changes in bone caused by osteoporosis. The goal of this study was to investigate how bone mineral density (BMD) and the microstructure of human cancellous bone affect three ultrasonic backscatter parameters that have been identified as potentially useful for ultrasonic bone assessment purposes: the apparent integrated backscatter (AIB), the frequency slope of apparent backscatter (FSAB), and the frequency intercept of apparent backscatter (FIAB). Ultrasonic measurements were performed with a 3.5-MHz broadband transducer on 54 specimens of human cancellous bone prepared from the proximal femur. Microstructural parameters and BMD were measured using X-ray microcomputed tomography (micro-CT). Relationships between AIB, FSAB, FIAB, and the micro-CT parameters were investigated using univariate and multivariate statistical analysis techniques. Moderate-to-strong univariate correlations were observed between the backscatter parameters and microstructure and BMD in many cases. The partial correlation analysis indicated that the backscatter parameters are dependent on microstructure independently of BMD in some cases. Multiple stepwise linear regression analysis used to generate multivariate models found that microstructure was a significant predictor of the backscatter parameters in most cases.

The feasibility study of the transmission mode photoacoustic measurement of human calcaneus bone in vivo

The photoacoustic (PA) technique is uniquely positioned for biomedical applications primarily due to its ability to visualize optical absorption contrast in deep tissue at ultrasound resolution. In this work, via both three-dimensional (3D) numerical simulations and in vivo experiments on human subjects, we investigated the possibility of PA measurement of human calcaneus bones in vivo in a non-invasive manner, as well as its feasibility to differentiate osteoporosis patients from normal subjects. The results from the simulations and the experiments both demonstrated that, when one side of the heel is illuminated by laser with light fluence under the ANSI safety limit, the PA signal generated in the human calcaneus bone can be detected by an ultrasonic transducer at the other side of the heel (i.e. transmission mode). Quantitative power spectral analyses of the calcaneus bone PA signals were also conducted, demonstrating that the microarchitectural changes in calcaneus bone due to osteoporosis can be detected, as reflected by enhanced high frequency components in detected PA bone signal. Further statistical analysis of the experimental results from 10 osteoporosis patients and 10 healthy volunteers showed that the weighted frequency as a quantified PA spectral parameter can differentiate the two subject groups with statistical significance.

Three-dimensional-printed replica models of bone for experimentally decoupling trabecular bone properties contribution to ultrasound propagation parameters

A detailed investigation of the relationship between ultrasonic (US) properties and trabecular bone microstructure is difficult because of the great variability in the bone loss process. The aim of this work was twofold. First, to verify by compressive tests that the three-dimensional (3D)-printer is able to produce precisely and repeatedly "bone replica models" of different size and density. Following, replicas of the original specimens with two different polymers and thinned trabeculae models were used to investigate US properties (speed of sound, SOS, and backscatter coefficient), aiming to deconvolute the influence of material properties on ultrasound characteristics. The results revealed that matrix material properties influence only the magnitude of the backscatter coefficient, whereas the characteristic undulated patterns are related to the trabecular structure. Simulation of perforation and thinning of cancellous bone, associated with bone loss, showed that SOS and mechanical properties were reduced perfectly linearly with apparent density when structure deteriorated. The 3D-printed bone replicas have the potential to enable systematic investigations of the influence of structure on both acoustical and mechanical properties and evaluate changes caused by bone loss. The development of replicas from materials with properties close to those of bone will permit quantitative conclusions for trabecular bone.

Ultrasonic Bone Assessment Using the Backscatter Amplitude Decay Constant

Ultrasonic backscatter techniques are being developed to detect changes in bone caused by osteoporosis. The present study introduces a new technique that measures the exponential decay in the amplitude of the backscatter signal quantified by a parameter called the backscatter amplitude decay constant (BADC). Measurements were performed on 54 specimens of cancellous bone from 14 human femurs using a 3.5-MHz transducer. Six methods were tested to determine BADC. The recommended method measures the time slope of the natural log of the rectified signal. Measured values of BADC ranged from approximately 0.1 μs^-1 to 0.6 μs^-1. Moderate to strong correlations (Spearman's ρ >0.7) were found between BADC and the density and microstructural characteristics of the specimens determined using X-ray microcomputed tomography. The results of this study suggest that BADC may be able to detect changes in the density and microstructure of cancellous bone caused by osteoporosis and other diseases.

Relationships of the ultrasonic backscatter measurements with the bone mineral density and the microarchitectural parameters in bovine trabecular bone in vitro

Relationships of the backscatter coefficient (BC), the apparent integrated backscatter (AIB), and the integrated reflection coefficient (IRC) with the bone mineral density (BMD) and the microarchitectural parameters were investigated in 28 bovine femoral trabecular bone samples. The BC was highly correlated with the BMD and the microarchitectural parameters (R = -0.66 to 0.71). In contrast, the AIB and the IRC exhibited high correlations with the BMD and the bone volume fraction (R = -0.68 to 0.77) and relatively lower correlations with the remaining microarchitectural parameters (R = -0.62 to 0.60). The multiple regression models yielded the adjusted squared correlation coefficients of 0.54-0.76.

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.

The Ability of Ultrasonic Backscatter Parametric Imaging to Characterize Bovine Trabecular Bone

The ultrasonic backscatter technique holds the promise of characterizing bone density and microstructure. This paper conducts ultrasonic backscatter parametric imaging based on measurements of apparent integrated backscatter (AIB), spectral centroid shift (SCS), frequency slope of apparent backscatter (FSAB), and frequency intercept of apparent backscatter (FIAB) for representing trabecular bone mass and microstructure. We scanned 33 bovine trabecular bone samples using a 7.5 MHz focused transducer in a 20 mm × 20 mm region of interest (ROI) with a step interval of 0.05 mm. Images based on the ultrasonic backscatter parameters (i.e., AIB, SCS, FSAB, and FIAB) were constructed to compare with photographic images of the specimens as well as two-dimensional (2D) μ-CT images from approximately the same depth and location of the specimen. Similar structures and trabecular alignments can be observed among these images. Statistical analyses demonstrated that the means and standard deviations of the ultrasonic backscatter parameters exhibited significant correlations with bone density (|R| = 0.45-0.78, p < 0.01) and bone microstructure (|R| = 0.44-0.87, p < 0.001). Some bovine trabecular bone microstructure parameters were independently associated with the ultrasonic backscatter parameters (ΔR² = 4.18%-44.45%, p < 0.05) after adjustment for bone apparent density (BAD). The results show that ultrasonic backscatter parametric imaging can provide a direct view of the trabecular microstructure and can reflect information about the density and microstructure of trabecular bone.

Relationships of Ultrasonic Backscatter With Bone Densities and Microstructure in Bovine Cancellous Bone

This study was designed to investigate the associations among ultrasonic backscatter, bone densities, and microstructure in bovine cancellous bone. Ultrasonic backscatter measurements were performed on 33 bovine cancellous bone specimens with a 2.25-MHz transducer. Ultrasonic apparent backscatter parameters ("apparent" means not compensating for ultrasonic attenuation and diffraction) were calculated with optimal signals of interest. The results showed that ultrasonic backscatter was significantly related to bone densities and microstructure ( R² = 0.17 -0.88 and ). After adjusting the correlations by bone mineral density (BMD), the bone apparent density (BAD) and some trabecular structural features still contributed significantly to the adjusted correlations, with moderate additional variance explained ( ∆R² = 9.7 % at best). Multiple linear regressions revealed that both BAD and trabecular structure contributed significantly and independently to the prediction of ultrasound backscatter (adjusted R² = 0.75 -0.89 and ), explaining an additional 14% of the variance at most, compared with that of BMD measurements alone. The results proved that ultrasonic backscatter was primarily determined by BAD, not BMD, but the combination of bone structure and densities could achieve encouragingly better performances (89% of the variance explained at best) in predicting backscatter properties. This study demonstrated that ultrasonic apparent backscatter might provide additional density and structural features unrelated to current BMD measurement. Therefore, we suggest that ultrasonic backscatter measurement could play a more important role in cancellous bone evaluation.

Ultrasonic backscatter difference measurements of cancellous bone from the human femur: Relation to bone mineral density and microstructure

Ultrasonic backscatter techniques are being developed to detect changes in cancellous bone caused by osteoporosis. One technique, called the backscatter difference technique, measures the power difference between two portions of a backscatter signal. The goal of the present study is to investigate how bone mineral density (BMD) and the microstructure of human cancellous bone influence four backscatter difference parameters: the normalized mean of the backscatter difference (nMBD) spectrum, the normalized slope of the backscatter difference spectrum, the normalized intercept of the backscatter difference spectrum, and the normalized backscatter amplitude ratio (nBAR). Ultrasonic measurements were performed with a 3.5 MHz broadband transducer on 54 specimens of human cancellous bone from the proximal femur. Volumetric BMD and the microstructural characteristics of the specimens were measured using x-ray micro-computed tomography. Of the four ultrasonic parameters studied, nMBD and nBAR demonstrated the strongest univariate correlations with density and microstructure. Multivariate analyses indicated that nMBD and nBAR depended on trabecular separation and possibly other microstructural characteristics of the specimens independently of BMD. These findings suggest that nMBD and nBAR may be sensitive to changes in the density and microstructure of bone caused by osteoporosis.

Characterization of a polymer, open-cell rigid foam that simulates the ultrasonic properties of cancellous bone

Materials that simulate the ultrasonic properties of tissues are used widely for clinical and research purposes. However, relatively few materials are known to simulate the ultrasonic properties of cancellous bone. The goal of the present study was to investigate the suitability of using a polymer, open-cell rigid foam (OCRF) produced by Sawbones^®. Measurements were performed on OCRF specimens with four different densities. Ultrasonic speed of sound and normalized broadband ultrasonic attenuation were measured with a 0.5 MHz transducer. Three backscatter parameters were measured with a 5 MHz transducer: apparent integrated backscatter, frequency slope of apparent backscatter, and normalized mean of the backscatter difference. X-ray micro-computed tomography was used to measure the microstructural characteristics of the OCRF specimens. The trabecular thickness and relative bone volume of the OCRF specimens were similar to those of human cancellous bone, but the trabecular separation was greater. In most cases, the ultrasonic properties of the OCRF specimens were similar to values reported in the literature for cancellous bone, including dependence on density. In addition, the OCRF specimens exhibited an ultrasonic anisotropy similar to that reported for cancellous bone.

Relationships among ultrasonic and mechanical properties of cancellous bone in human calcaneus in vitro

Clinical bone sonometers applied at the calcaneus measure broadband ultrasound attenuation and speed of sound. However, the relation of ultrasound measurements to bone strength is not well-characterized. Addressing this issue, we assessed the extent to which ultrasonic measurements convey in vitro mechanical properties in 25 human calcaneal cancellous bone specimens (approximately 2×4×2cm). Normalized broadband ultrasound attenuation, speed of sound, and broadband ultrasound backscatter were measured with 500kHz transducers. To assess mechanical properties, non-linear finite element analysis, based on micro-computed tomography images (34-micron cubic voxel), was used to estimate apparent elastic modulus, overall specimen stiffness, and apparent yield stress, with models typically having approximately 25-30 million elements. We found that ultrasound parameters were correlated with mechanical properties with R=0.70-0.82 (p<0.001). Multiple regression analysis indicated that ultrasound measurements provide additional information regarding mechanical properties beyond that provided by bone quantity alone (p≤0.05). Adding ultrasound variables to linear regression models based on bone quantity improved adjusted squared correlation coefficients from 0.65 to 0.77 (stiffness), 0.76 to 0.81 (apparent modulus), and 0.67 to 0.73 (yield stress). These results indicate that ultrasound can provide complementary (to bone quantity) information regarding mechanical behavior of cancellous bone.

Porosity predicted from ultrasound backscatter using multivariate analysis can improve accuracy of cortical bone thickness assessment

A rapidly growing area of interest in quantitative ultrasound assessment of bone is to determine cortical bone porosity from ultrasound backscatter. Current backscatter analyses are based on numerical simulations, while there are no published reports of successful experimental measurements. In this study, multivariate analysis is applied to ultrasound reflections and backscatter to predict cortical bone porosity. The porosity is then applied to estimate cortical bone radial speed of sound (SOS) and thickness using ultrasound backscatter signals obtained at 2.25 and 5 MHz center frequencies from cortical bone samples (n = 43) extracted from femoral diaphyses. The study shows that the partial least squares regression technique could be employed to successfully predict (R²= 0.71-0.73) cortical porosity. It is found that this multivariate approach can reduce uncertainty in pulse-echo assessment of cortical bone thickness from 0.220 to 0.045 mm when porosity based radial SOS was applied, instead of a constant value from literature. Upon further validation, accurate estimation of cortical bone porosity and thickness may be applied as a financially viable option for fracture risk assessment of individuals.

Solid volume fraction estimation of bone:marrow replica models using ultrasound transit time spectroscopy

The acceptance of broadband ultrasound attenuation (BUA) for the assessment of osteoporosis suffers from a limited understanding of both ultrasound wave propagation through cancellous bone and its exact dependence upon the material and structural properties. It has recently been proposed that ultrasound wave propagation in cancellous bone may be described by a concept of parallel sonic rays; the transit time of each ray defined by the proportion of bone and marrow propagated. A Transit Time Spectrum (TTS) describes the proportion of sonic rays having a particular transit time, effectively describing the lateral inhomogeneity of transit times over the surface aperture of the receive ultrasound transducer. The aim of this study was to test the hypothesis that the solid volume fraction (SVF) of simplified bone:marrow replica models may be reliably estimated from the corresponding ultrasound transit time spectrum. Transit time spectra were derived via digital deconvolution of the experimentally measured input and output ultrasonic signals, and compared to predicted TTS based on the parallel sonic ray concept, demonstrating agreement in both position and amplitude of spectral peaks. Solid volume fraction was calculated from the TTS; agreement between true (geometric calculation) with predicted (computer simulation) and experimentally-derived values were R(2)=99.9% and R(2)=97.3% respectively. It is therefore envisaged that ultrasound transit time spectroscopy (UTTS) offers the potential to reliably estimate bone mineral density and hence the established T-score parameter for clinical osteoporosis assessment.

Ultrasonic backscatter from cancellous bone: the apparent backscatter transfer function

Ultrasonic backscatter techniques are being developed to detect changes in cancellous bone caused by osteoporosis. Many techniques are based on measurements of the apparent backscatter transfer function (ABTF), which represents the backscattered power from bone corrected for the frequency response of the measurement system. The ABTF is determined from a portion of the backscatter signal selected by an analysis gate of width τw delayed by an amount τd from the start of the signal. The goal of this study was to characterize the ABTF for a wide range of gate delays (1 μs ≤ τd ≤ 6 μs) and gate widths (1 μs ≤ τw ≤ 6 μs). Measurements were performed on 29 specimens of human cancellous bone in the frequency range 1.5 to 6.0 MHz using a broadband 5-MHz transducer. The ABTF was found to be an approximately linear function of frequency for most choices of τd and τw. Changes in τd and τw caused the frequency-averaged ABTF [quantified by apparent integrated backscatter (AIB)] and the frequency dependence of the ABTF [quantified by frequency slope of apparent backscatter (FSAB)] to change by as much as 24.6 dB and 6.7 dB/MHz, respectively. τd strongly influenced the measured values of AIB and FSAB and the correlation of AIB with bone density (-0.95 ≤ R ≤ +0.68). The correlation of FSAB with bone density was influenced less strongly by τd (-0.97 ≤ R ≤ -0.87). τw had a weaker influence than τd on the measured values of AIB and FSAB and the correlation of these parameters with bone density.

Effect of intervening tissues on ultrasonic backscatter measurements of bone: An in vitro study

Ultrasonic backscatter techniques are being developed to diagnose osteoporosis. Tissues that lie between the transducer and the ultrasonically interrogated region of bone may produce errors in backscatter measurements. The goal of this study is to investigate the effects of intervening tissues on ultrasonic backscatter measurements of bone. Measurements were performed on 24 cube shaped specimens of human cancellous bone using a 5 MHz transducer. Measurements were repeated after adding a 1 mm thick plate of cortical bone to simulate the bone cortex and a 3 cm thick phantom to simulate soft tissue at the hip. Signals were analyzed to determine three apparent backscatter parameters (apparent integrated backscatter, frequency slope of apparent backscatter, and frequency intercept of apparent backscatter) and three backscatter difference parameters [normalized mean backscatter difference (nMBD), normalized slope of the backscatter difference, and normalized intercept of the backscatter difference]. The apparent backscatter parameters were impacted significantly by the presence of intervening tissues. In contrast, the backscatter difference parameters were not affected by intervening tissues. However, only one backscatter difference parameter, nMBD, demonstrated a strong correlation with bone mineral density. Thus, among the six parameters tested, nMBD may be the best choice for in vivo backscatter measurements of bone when intervening tissues are present.

Photoacoustic and ultrasound imaging of cancellous bone tissue

We used ultrasound (US) and photoacoustic (PA) imaging modalities to characterize cattle trabecular bones. The PA signals were generated with an 805-nm continuous wave laser used for optimally deep optical penetration depth. The detector for both modalities was a 2.25-MHz US transducer with a lateral resolution of ~1 mm at its focal point. Using a lateral pixel size much larger than the size of the trabeculae, raster scanning generated PA images related to the averaged values of the optical and thermoelastic properties, as well as density measurements in the focal volume. US backscatter yielded images related to mechanical properties and density in the focal volume. The depth of interest was selected by time-gating the signals for both modalities. The raster scanned PA and US images were compared with microcomputed tomography (μCT) images averaged over the same volume to generate similar spatial resolution as US and PA. The comparison revealed correlations between PA and US modalities with the mineral volume fraction of the bone tissue. Various features and properties of these modalities such as detectable depth, resolution, and sensitivity are discussed.

Interstitial ultrasound ablation of vertebral and paraspinal tumours: parametric and patient-specific simulations

PURPOSE

Theoretical parametric and patient-specific models are applied to assess the feasibility of interstitial ultrasound ablation of tumours in and near the spine and to identify potential treatment delivery strategies.

METHODS

3D patient-specific finite element models (n = 11) of interstitial ultrasound ablation of tumours associated with the spine were generated. Gaseous nerve insulation and various applicator configurations, frequencies (3 and 7 MHz), placement trajectories, and tumour locations were simulated. Parametric studies with multilayered models investigated the impacts of tumour attenuation, tumour dimension, and the thickness of bone insulating critical structures. Temperature and thermal dose were calculated to define ablation (>240 equivalent minutes at 43 °C (EM43 °C)) and safety margins (<45 °C and <6 EM43 °C), and to determine performance and required delivery parameters.

RESULTS

Osteolytic tumours (≤44 mm) encapsulated by bone could be successfully ablated with 7 MHz interstitial ultrasound (8.1-16.6 W/cm(2), 120-5900 J, 0.4-15 min). Ablation of tumours (94.6-100% volumetric) 0-14.5 mm from the spinal canal was achieved within 3-15 min without damaging critical nerves. 3 MHz devices provided faster ablation (390 versus 930 s) of an 18 mm diameter osteoblastic (high bone content) volume than 7 MHz devices. Critical anatomy in proximity to the tumour could be protected by selection of appropriate applicator configurations, active sectors, and applied power schemas, and through gaseous insulation. Preferential ultrasound absorption at bone surfaces facilitated faster, more effective ablations in osteolytic tumours and provided isolation of ablative energies and temperatures.

CONCLUSIONS

Parametric and patient-specific studies demonstrated the feasibility and potential advantages of interstitial ultrasound ablation treatment of paraspinal and osteolytic vertebral tumours.

Ultrasound method applied to characterize healthy femoral diaphysis of Wistar rats in vivo

A simple experimental protocol applying a quantitative ultrasound (QUS) pulse-echo technique was used to measure the acoustic parameters of healthy femoral diaphyses of Wistar rats in vivo. Five quantitative parameters [apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), time slope of apparent backscatter (TSAB), integrated reflection coefficient (IRC), and frequency slope of integrated reflection (FSIR)] were calculated using the echoes from cortical and trabecular bone in the femurs of 14 Wistar rats. Signal acquisition was performed three times in each rat, with the ultrasound signal acquired along the femur's central region from three positions 1 mm apart from each other. The parameters estimated for the three positions were averaged to represent the femur diaphysis. The results showed that AIB, FSAB, TSAB, and IRC values were statistically similar, but the FSIR values from Experiments 1 and 3 were different. Furthermore, Pearson's correlation coefficient showed, in general, strong correlations among the parameters. The proposed protocol and calculated parameters demonstrated the potential to characterize the femur diaphysis of rats in vivo. The results are relevant because rats have a bone structure very similar to humans, and thus are an important step toward preclinical trials and subsequent application of QUS in humans.

The application of backscattered ultrasound and photoacoustic signals for assessment of bone collagen and mineral contents

BACKGROUND

This study examines the backscattered ultrasound (US) and back-propagating photoacoustic (PA) signals from trabecular bones, and their variations with reduction in bone minerals and collagen content. While the collagen status is directly related to the strength of the bone, diagnosis of its condition using US remains a challenge.

METHODS

For both PA and US methods, coded-excitation signals and matched filtering were utilized to provide high sensitivity of the detected signal. The optical source was a 805-nm CW laser and signals were detected employing a 2.2-MHz ultrasonic transducer. Bone decalcification and decollagenization were induced with mild ethylenediaminetetraacetic acid (EDTA) and sodium hypochlorite solutions, respectively.

RESULTS

The PA and US signals were measured on cattle bones, and apparent integrated backscatter/back-propagating (AIB) parameters were compared before and after demineralization and decollagenization.

CONCLUSIONS

The results show that both PA and US are sensitive to mineral changes. In addition, PA is also sensitive to changes in the collagen content of the bone, but US is not significantly sensitive to these changes.

Importance of Material Properties and Porosity of Bone on Mechanical Response of Articular Cartilage in Human Knee Joint—A Two-Dimensional Finite Element Study

Mechanical behavior of bone is determined by the structure and intrinsic, local material properties of the tissue. However, previously presented knee joint models for evaluation of stresses and strains in joints generally consider bones as rigid bodies or linearly elastic solid materials. The aim of this study was to estimate how different structural and mechanical properties of bone affect the mechanical response of articular cartilage within a knee joint. Based on a cadaver knee joint, a two-dimensional (2D) finite element (FE) model of a knee joint including bone, cartilage, and meniscus geometries was constructed. Six different computational models with varying properties for cortical, trabecular, and subchondral bone were created, while the biphasic fibril-reinforced properties of cartilage and menisci were kept unaltered. The simplest model included rigid bones, while the most complex model included specific mechanical properties for different bone structures and anatomically accurate trabecular structure. Models with different porosities of trabecular bone were also constructed. All models were exposed to axial loading of 1.9 times body weight within 0.2 s (mimicking typical maximum knee joint forces during gait) while free varus–valgus rotation was allowed and all other rotations and translations were fixed. As compared to results obtained with the rigid bone model, stresses, strains, and pore pressures observed in cartilage decreased depending on the implemented properties of trabecular bone. Greatest changes in these parameters (up to −51% in maximum principal stresses) were observed when the lowest modulus for trabecular bone (measured at the structural level) was used. By increasing the trabecular bone porosity, stresses and strains were reduced substantially in the lateral tibial cartilage, while they remained relatively constant in the medial tibial plateau. The present results highlight the importance of long bones, in particular, their mechanical properties and porosity, in altering and redistributing forces transmitted through the knee joint.

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.

Prediction of trabecular bone qualitative properties using scanning quantitative ultrasound

Microgravity induced bone loss represents a critical health problem in astronauts, particularly occurred in weight-supporting skeleton, which leads to osteopenia and increase of fracture risk. Lack of suitable evaluation modality makes it difficult for monitoring skeletal status in long term space mission and increases potential risk of complication. Such disuse osteopenia and osteoporosis compromise trabecular bone density, and architectural and mechanical properties. While X-ray based imaging would not be practical in space, quantitative ultrasound may provide advantages to characterize bone density and strength through wave propagation in complex trabecular structure. This study used a scanning confocal acoustic diagnostic and navigation system (SCAN) to evaluate trabecular bone quality in 60 cubic trabecular samples harvested from adult sheep. Ultrasound image based SCAN measurements in structural and strength properties were validated by μCT and compressive mechanical testing. This result indicated a moderately strong negative correlations observed between broadband ultrasonic attenuation (BUA) and μCT-determined bone volume fraction (BV/TV, R²=0.53). Strong correlations were observed between ultrasound velocity (UV) and bone's mechanical strength and structural parameters, i.e., bulk Young's modulus (R²=0.67) and BV/TV (R²=0.85). The predictions for bone density and mechanical strength were significantly improved by using a linear combination of both BUA and UV, yielding R²=0.92 for BV/TV and R²=0.71 for bulk Young's modulus. These results imply that quantitative ultrasound can characterize trabecular structural and mechanical properties through measurements of particular ultrasound parameters, and potentially provide an excellent estimation for bone's structural integrity.

Estimation of multipath transmission parameters for quantitative ultrasound measurements of bone

When applying quantitative ultrasound (QUS) measurements to bone for predicting osteoporotic fracture risk, the multipath transmission of sound waves frequently occurs. In the last 10 years, the interest in separating multipath QUS signals for their analysis awoke, and led to the introduction of several approaches. Here, we compare the performances of the two fastest algorithms proposed for QUS measurements of bone: the modified least-squares Prony method (MLSP), and the space alternating generalized expectation maximization algorithm (SAGE) applied in the frequency domain. In both approaches, the parameters of the transfer functions of the sound propagation paths are estimated. To provide an objective measure, we also analytically derive the Cramér-Rao lower bound of variances for any estimator and arbitrary transmit signals. In comparison with results of Monte Carlo simulations, this measure is used to evaluate both approaches regarding their accuracy and precision. Additionally, with simulations using typical QUS measurement settings, we illustrate the limitations of separating two superimposed waves for varying parameters with focus on their temporal separation. It is shown that for good SNRs around 100 dB, MLSP yields better results when two waves are very close. Additionally, the parameters of the smaller wave are more reliably estimated. If the SNR decreases, the parameter estimation with MLSP becomes biased and inefficient. Then, the robustness to noise of the SAGE clearly prevails. Because a clear influence of the interrelation between the wavelength of the ultrasound signals and their temporal separation is observable on the results, these findings can be transferred to QUS measurements at other sites. The choice of the suitable algorithm thus depends on the measurement conditions.

Arthroscopic ultrasound technique for simultaneous quantitative assessment of articular cartilage and subchondral bone: an in vitro and in vivo feasibility study

Traditional arthroscopic examination is subjective and poorly reproducible. Recently, we introduced an arthroscopic ultrasound method for quantitative diagnostics of cartilage lesions. Here we describe our investigation of the feasibility of ultrasound arthroscopy for simultaneous measurements of articular cartilage and subchondral bone. Human osteochondral samples (n = 13) were imaged using a clinical 9-MHz ultrasound system. Ultrasound reflection coefficients (R, IRC), the ultrasound roughness index (URI) and the apparent integrated backscattering coefficient (AIB) were determined for both tissues. Mechanical testing, histologic analyses and micro-scale computed tomography imaging were the reference methods. Ultrasound arthroscopies were conducted on two patients. The ultrasound reflection coefficient correlated with the Mankin score and Young's modulus of cartilage (|r| > 0.56, p < 0.05). Ultrasound parameters (R, IRC, AIB) for subchondral bone correlated with the bone surface/volume ratio (|r| > 0.70, p < 0.05) and trabecular thickness (|r| > 0.59, p < 0.05). Furthermore, R and subchondral bone mineral density were significantly correlated (|r| > 0.65, p < 0.05). Arthroscopic ultrasound examination provided diagnostically valuable information on cartilage and subchondral bone in vivo.

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.

A backscatter difference technique for ultrasonic bone assessment

Ultrasonic backscatter techniques may offer a useful approach for detecting changes in cancellous bone caused by osteoporosis and other diseases. The goal of this study was to investigate the utility of a backscatter difference technique for ultrasonic bone assessment. Measurements were performed on 22 cube-shaped specimens of human cancellous bone using four broadband transducers with center frequencies 2.25, 5, 7.5, and 10 MHz. The backscatter difference spectrum D(f) was obtained by subtracting power spectra (in dB) from two different portions of the same backscatter signal. D(f) was found to be a monotonically increasing, quasi-linear function of frequency when averaged over multiple measurement sites on multiple specimens. The frequency slope of D(f) demonstrated weak to moderate correlations with specimen density (R = 0.21-0.80). The frequency averaged mean of D(f) demonstrated moderate to good correlations with density (R = 0.70-0.95). These results suggest that parameters based on the frequency averaged mean of the backscatter difference spectrum may be useful for bone assessment purposes.

Computed radiographic and ultrasonic evaluation of bone regeneration during tibial distraction osteogenesis in rabbits

Computed radiography (CR) and a combined ultrasound (US) approach involving two-dimensional (2-D) and three-dimensional (3-D) ultrasonography with ultrasonometry were employed to evaluate their respective efficacies in monitoring bone regeneration during rabbit tibial distraction osteogenesis (DO). Results demonstrated that 2-D and 3-D ultrasonography depicted bone callus growth changes during distraction while CR could not. Evaluation of callus speed of sound, acoustic reflection and attenuation showed significant linear changes over time during early DO stage (p < 0.05). However, surrogate measure of callus density by CR only showed such significant linear changes during consolidation (p < 0.05). Also, callus speed of sound and acoustic reflection during early DO stage showed strong predictions to the bone mineral density and microstructural properties (adjusted-R(2) = 0.43-0.67) of consolidated bone callus measured at the treatment end-point by microcomputed tomography. Findings of the present study indicated a preferred use of the combined US approach over CR in the early monitoring of bone regeneration during DO treatment.

Relationships of quantitative ultrasound parameters with cancellous bone microstructure in human calcaneus in vitro

Ultrasound parameters (attenuation, phase velocity, and backscatter), bone mineral density (BMD), and microarchitectural features were measured on 29 human cancellous calcaneus samples in vitro. Regression analysis was performed to predict ultrasound parameters from BMD and microarchitectural features. The best univariate predictors of the ultrasound parameters were the indexes of bone quantity: BMD and bone volume fraction (BV/TV). The most predictive univariate models for attenuation, phase velocity, and backscatter coefficient yielded adjusted squared correlation coefficients of 0.69-0.73. Multiple regression models yielded adjusted correlation coefficients of 0.74-0.83. Therefore attenuation, phase velocity, and backscatter are primarily determined by bone quantity, but multiple regression models based on bone quantity plus microarchitectural features achieve slightly better predictive performance than models based on bone quantity alone.

Effects of non-optimal focusing on dual-frequency ultrasound measurements of bone

In pulse-echo (PE) ultrasound measurements, the use of focused transducers is desirable for quantitative assessment of bone characteristics because of the attenuation in the overlying soft tissues. However, the variable thickness and composition of the soft tissue overlying bone affect the focal depth of the ultrasound beam and induce errors into the measurements. To compensate for the attenuation-related effects caused by the interfering soft tissue (i.e., fat and lean tissue), a dual-frequency ultrasound (DFUS) technique was recently introduced. The aim of this study was to investigate the effect of non-optimal focal depth of the ultrasound beam on the determination of the integrated reflection coefficient (IRC) of bone when overlaid by an interfering layer composed of oil and water. The feasibility of the DFUS-based correction of the IRC was evaluated through numerical simulations and experimental measurements. Even when the interfering layer-bone interface was out of focus, the total thickness of the interfering layer could be accurately determined with the technique. However, based on the simulations, the errors in the determination of the composition of the interfering layer increased (0.004 to 113.8%) with the increase in distance between the interfering layer-bone interface and the focus of the ultrasound beam. Attenuation compensation, based on the true composition of the interfering layer, resulted in an average relative error of 22.3% in the IRC values calculated from the simulations. Interestingly, the attenuation compensation with the interfering layer composition estimated using the DFUS technique resulted in a smaller average relative error of 14.9% in the IRC values. The simulations suggest that DFUS can reduce the errors induced by soft tissue in bone PE ultrasound measurements. The experimental measurements indicate that the accuracy of the IRC measurements is rather similar when using DFUS correction or correction based on the true composition of the interfering layer. However, the results suggest that accurate determination of soft tissue composition may be difficult without optimal focusing of the ultrasound beam on the soft tissue-bone interface.

Structural parameters of normal and osteoporotic human trabecular bone are affected differently by microCT image resolution

This study employed microCT to investigate whether image resolution affects bone structural parameters differently in healthy normal and osteoporotic trabecular bone. With increasing image voxel size, the originally detected differences between sample groups diminished. The results suggest that structural differences may not be reliably detected with clinical scanners.

INTRODUCTION

Structural parameters of bone reflect its health status, but are highly dependent on the image resolution. We hypothesized that image resolution affects bone structural parameters differently in normal and osteoporotic trabecular bone.

METHODS

Human trabecular bone samples from the iliac crest and the knee were analyzed (normal n = 11, osteoporotic n = 15) using a high-resolution microCT (14 or 18 µm voxel sizes). Images were re-sampled to voxel sizes 1-16 times larger than the original image and thresholded with global or local adaptive algorithms. Absolute and normalized values of each structural parameter were calculated, and the effect of decreasing image resolution was compared between the normal and osteoporotic samples.

RESULTS

Normal and osteoporotic samples had different (p < 0.05) absolute bone volume fractions. However, the normalized values showed that the osteoporotic samples were more prone to errors (p < 0.05) with increased voxel size. The absolute values of trabecular number, trabecular separation, degree of anisotropy, and structure model index were different between the groups at the original voxel size (p < 0.05), but at voxel sizes between 60 and 110 µm, those differences were no longer significant.

CONCLUSIONS

The results suggest that structural differences between osteoporotic and normal trabecular bone may not be reliably detected with clinical CT scanners providing image voxel sizes above 100 µm.

Phase velocity analysis of acoustic propagation in trabecular bone

The trabecular bones are highly dispersive acoustic media composed by randomly shaped trabeculae (considered as scatterers) and surrounded by bone marrow. An acoustic wave generated by an incident ultrasonic pulse with respect to the media under study, shows that its velocity and amplitude are a function of the density size and shaped of the scatterers. Two different methods were used to theoretically study this scattering phenomena: i) a self-consistent model proposed by Sabina and Willis, and ii) the multiple-scattering theory of Waterman-Truell. These methods were able to compute the phase velocity and amplitude as a function of operating frequency, density and size of scatterers. The theoretical results were compared with experimental data already published and the phase veloctiy shows a good agreement for low concentration of scatterers.

Ultrasound backscatter imaging provides frequency-dependent information on structure, composition and mechanical properties of human trabecular bone

The strength as well as the acoustic properties of trabecular bone are determined by its structure and composition. Consequently, tissue structure and compositional properties also affect the ultrasound propagation in bone. The diagnostic potential of ultrasound has not been fully exploited in clinical quantitative ultrasound devices. The aim of this study was to investigate the ability of quantitative ultrasound pulse-echo imaging, conducted over a broad range of frequencies (1 to 5 MHz), to predict the mechanics, composition and microstructure of trabecular bone. Ultrasound reflection and backscatter parameters correlated significantly with the ultimate strength of the trabecular bone and the bone volume fraction (r=0.76-0.90, n=20, p<0.01). Ultrasound backscatter associated significantly (independently of bone structure or mineral content) with the collagen content of the bone matrix (r=0.75, r(adjusted)=0.66, p<0.01). Interestingly, the applied ultrasound frequency seemed to relate the sensitivity of ultrasound backscatter to different properties of trabecular bone. At frequencies ranging from 1 to 3.5 MHz, the ultrasound backscatter associated significantly with the tissue mechanical and structural parameters. At 5MHz, the composition of the bone matrix was a more significant determinant of the measured backscatter. This study provides useful information for optimizing the use of pulse-echo measurements, and thereby further emphasizes the diagnostic potential of the ultrasound backscatter measurements of trabecular bone.

A comparison of principal component analysis and factor analysis strategies for uncovering pleiotropic factors

Principal component analysis (PCA) and factor analysis (FA) are often used to uncover genetic factors that contribute to complex disease phenotypes. The purpose of such an analysis is to distill a genetic signal from a large number of correlated phenotype measurements. That signal can then be used in genetic analyses (e.g. linkage analysis), presumably leading to greater success at finding genes than one would achieve with any one raw trait. Although both PCA and FA have been used this way, there has been no comparison of their performance in the literature. We compared the ability of these two procedures to extract unobserved underlying genetic components from complex simulated data on nuclear families. We first simulated seven underlying genetic and environmentally determined traits. Then we derived two sets of 50 complex (observed) traits using algebraic combinations of the underlying components. We next performed PCA and FA on the complex traits. We assessed two aspects of the performance of the methods: (1) ability to detect the underlying genetic components; (2) whether the methods worked better when applied to raw traits or to residuals (after regressing out significant environmental covariates). Our results indicate that both the methods behave similarly in most cases, although FA generally produced factors that had stronger correlations with the underlying traits. We also found that using residuals in PCA or FA analyses greatly increased the probability that the PCs or factors detected common genetic components instead of common environmental factors, except if there was statistical interaction between genetic and environmental factors.

Effect of bone marrow on acoustic properties of trabecular bone--3D finite difference modeling study

The composition of bone marrow is influenced by many factors, such as age and diseases. The present numerical study investigates the contribution of marrow on the acoustic measurements of trabecular bone. Cylindrical bone samples (n = 11), extracted from three anatomical sites of human cadaver knees, were imaged with a high-resolution microtomography (microCT). Three-dimensional finite difference time domain (FDTD) models (Wave 3000 Pro 2.2, Cyberlogic Inc., NY, USA) were created using the segmented microCT images of each sample. First, we evaluated the effect of voxel size on the computer resource requirements, morphological parameters and acoustic simulations. Second, the effect of bone marrow on ultrasonic measurements was assessed. The simulations were repeated with two voxel sizes before and after substitution of bone marrow (i.e., fat) with water. The voxel size of the FDTD mesh controlled the fine structure of the modeled calcified matrix and significantly affected the simulation results. However, present simulations showed that the effect of bone marrow on ultrasound parameters can be reliably simulated with the applied voxel sizes of 72 and 90 microm. Ultrasound attenuation and speed were found (p < 0.01) to decrease and increase, respectively, when bone marrow was substituted with water. Moreover, reflection from the surface of the sample increased (p < 0.01) and backscatter from internal structures decreased (p < 0.01) after removal of marrow. The effect of bone marrow on the acoustic properties was stronger in samples with low bone volume fraction. The present results indicate that the amount and quality of bone marrow significantly influence the acoustic properties of trabecular bone. Possible interindividual differences in the composition of bone marrow may increase uncertainty in clinical ultrasound diagnostics of osteoporosis. Importantly, the effect is most significant in osteoporotic low-density bone.

Velocity dispersion in trabecular bone: influence of multiple scattering and of absorption

Speed of sound measurements are widely used clinically to assess bone strength. Trabecular bone is an attenuating composite material in which negative values of velocity dispersion have been measured, this behavior remaining poorly explained physically. The aim of this work is to describe the ultrasonic propagation in trabecular bone modeled by infinite cylinders immersed in a saturating matrix, and to derive the physical determinants of velocity dispersion. A homogenization model accounting for the coupling of multiple scattering and absorption phenomena allows the computation of phase velocity and of dispersion while varying bone properties. The present model is adapted from the generalized self-consistent method developed in the work of Yang and Mal [(1994). "Multiple-scattering of elastic waves in a fiber-reinforced composite," J. Mech. Phys. Solids 42, 1945-1968]. It predicts negative values of velocity dispersion, in agreement with experimental results obtained in phantoms mimicking trabecular bone. In trabecular bone, mostly negative and also positive values of velocity dispersion are predicted by the model, which span within the range of values measured experimentally. Scattering effects are responsible for the negative values of dispersion, whereas the frequency dependence of the absorption coefficient in bone marrow and/or in the trabeculae results in an increase in dispersion, which may then become positive.

Ultrasonic characterization of human cancellous bone in vitro using three different apparent backscatter parameters in the frequency range 0.6-15.0 mhz

Ultrasonic techniques based on measurements of apparent backscatter may provide a useful means for diagnosing bone diseases such as osteoporosis. The term "apparent" means that the backscattered signals are not compensated for the frequency-dependent effects of attenuation and diffraction. We performed in vitro apparent backscatter measurements on 23 specimens of human cancellous bone prepared from the left and right femoral heads of seven donors. A mechanical scanning system was used to obtain backscattered signals from each specimen at several sites. Scans were performed using five different ultrasonic transducers with center frequencies of 1, 2.25, 5, 7.5, and 10 MHz. The -6 dB bandwidths of these transducers covered a frequency range of 0.6-15.0 MHz. The backscattered signals were analyzed to determine three ultrasonic parameters: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), and time slope of apparent backscatter (TSAB). Linear regression analysis was used to examine the correlation of these ultrasonic parameters with five measured physical characteristics of the specimens: mass density, X-ray bone mineral density, Young's modulus, yield strength, and ultimate strength. A total of 75 such correlations were examined (3 ultrasonic parameters x 5 specimen characteristics x 5 transducers). Good correlations were observed for AIB using the 5 MHz (r = 0.70 - 0.89) and 7.5 MHz (r = 0.75-0.93) transducers; for FSAB using the 2.25 MHz (r = 0.70 - 0.88), 5 MHz (r = 0.79 - 0.94), and 7.5 MHz (r = 0.80 - 0.92) transducers; and for TSAB using the 5 MHz (r = 0.68 - 0.89), 7.5 MHz (r = 0.75 - 0.89), and 10 MHz (r = 0.75 - 0.92) transducers.

Relationships of trabecular bone structure with quantitative ultrasound parameters: in vitro study on human proximal femur using transmission and backscatter measurements

The present study was designed to assess the relationships between QUS parameters and bone density or microarchitecture on samples of human femoral trabecular bone. The normalized slope of the frequency-dependent attenuation (nBUA), the speed of sound (SOS) and the broadband ultrasound backscatter coefficient (BUB) were measured on 37 specimens of pure trabecular bones removed from upper parts of fresh human femurs. Bone mineral density (BMD) was assessed using a clinical scanner. Finally, 8 mm diameter cylindrical cores were extracted from the specimens and their microarchitecture was reconstructed after synchrotron radiation microtomography experiments (isotropic resolution of 10 microm). A large number of microarchitectural parameters were computed, describing morphology, connectivity and geometry of the specimens. BMD correlated with all the microarchitectural parameters and the number of significant correlations was found among the architectural parameters themselves. All QUS parameters were significantly correlated to BMD (R=0.83 for nBUA, R=0.81 for SOS and R=0.69 for BUB) and to microarchitectural parameters (R=-0.79 between nBUA and Tb.Sp, R=-0.81 between SOS and Tb.Sp, R=-0.65 between BUB and BS/BV). Using multivariate model, it was found that microstructural parameters adds 10%, 19%, and 4% to the respective BMD alone contribution for the three variables BUA, SOS and BUB. Moreover, the RMSE was reduced by up to 50% for SOS, by up to 21% for nBUA and up to 11% when adding structural variables to BMD in explaining QUS results. Given the sample, which is not osteoporosis-enriched, the added contribution is quite substantial. The variability of SOS was indeed completely explained by a multivariate model including BMD and independent structural parameters (R(2)=0.94). The inverse problem on the data presented here has been addressed using simple and multiple linear regressions. It was shown that the predictions (in terms of R(2) or RMSE) of microarchitectural parameters was not enhanced when combining 2 or 3 QUS in multiple linear regressions compared to the prediction obtained with one QUS parameter alone. The best model was found for the prediction of Tb.Th() from BUA (R(2)=0.58, RMSE=17 microm). Given the high values of RMSE, these linear models appear of limited clinical value, suggesting that appropriate models have to be derived in order to solve the inverse problem. In this regard, a very interesting multivariate model was found for nBUA and BUB with Tb.Th and Tb.N, in agreement with single scattering theories by random medium. However, the source of residual variability of nBUA and BUB (15% and 45% respectively) remained unexplained.

Ultrasound transmission technique as a potential tool for physical evaluation of monolithic matrix tablets

The aim of this study was to investigate the effects of tablet porosity and particle size fraction of compacted Starch acetate powders, with and without model drug caffeine, on acoustic properties of tablets. The ultrasound velocity was determined from the transmission measurements. Tablets of starch acetate (SA DS 2.7) powder with two particle size fractions of 0-53 and 0-710 microm were compressed with a compaction simulator. Porosities of tablets varied in the range from 12% to 43% for both particle size fractions. Strong associations were found between the ultrasound velocity and physical properties of the tablets such as porosity and particle size fraction. Interestingly, ultrasound velocity was practically insensitive to inclusion of the model drug caffeine with the concentrations used. Based on this study ultrasound transmission method is a potential non-destructive tool for studying structural changes of tablets and other solid dosage forms.

Development of a method for ultrasound-guided placement of pedicle screws

Abstract-Many forms of spinal fusion involve the placement of long screws through the pedicles into the vertebral body. During the procedure, there is substantial risk of damage to vital neural and vascular structures due to the limited visibility of anatomic landmarks and high anatomic variability. As an alternative to current guidance systems, we have investigated the feasibility of performing ultrasound imaging through cancellous bone for the purpose of pedicle screw guidance. Quantitative ultrasonic characterization and A-mode imaging of seven defatted vertebral cancellous bone specimens was performed along the craniocaudal axis in water with unfocused, 1-MHz and 3.5- MHz broadband transducers. The center frequency attenuation increased considerably from 10.5 +/- 4.6 dB/cm at 1 MHz to 24.1 +/- 7.2 dB/cm at 3.5 MHz, while the speed of sound exhibited moderate positive dispersion, increasing from 1489 +/- 4.7 m/s at 1 MHz to 1494 +/- 4.2 m/s at 3.5 MHz. Despite the high attenuation and large specimen thickness (1.0-1.9 cm), A-mode imaging through cancellous bone to detect an aluminum reflector was possible in 83.2% and 70.1% of the cases at 1 MHz and 3.5 MHz, respectively. Specimen boundaries were identifiable with clinically sufficient average accuracy of 1.1 mm and 0.9 mm in the 1 MHz and 3.5 MHz A-mode images, respectively.

Spatial variation of acoustic properties is related with mechanical properties of trabecular bone

In clinical applications, ultrasound parameters are measured as an average value over a region of interest (ROI) or as a value at a single measurement point. Due to natural adaptation to loading conditions, trabecular bone is structurally, compositionally and mechanically heterogeneous and anisotropic. Thus, spatial variation of ultrasound parameters within ROI may contain valuable information on the mechanical integrity of trabecular bone. However, this issue has not been thoroughly investigated. In the present study, we aimed at investigating the significance of the spatial variation of ultrasound parameters for the prediction of mechanical properties of human trabecular bone. For this aim, parametric maps of apparent integrated backscattering (AIB), integrated reflection coefficient (IRC), speed of sound (SOS), average attenuation (AA) and normalized broadband ultrasound attenuation (nBUA) were calculated for femoral and tibial bone cylinders (n = 19-20). Further, the effect of time window length on the AIB, variation of AIB within ROI and association between AIB and bone mechanical properties were characterized. Based on linear correlation analysis, spatial variation of AIB, assessed as standard deviation of measurements within ROI, was a strong predictor of bone ultimate strength (r = -0.82, n = 19, p < 0.01). Further, the time window length affected absolute values of AIB and strength of correlation between AIB and bone ultimate strength. Interestingly, linear combination of mean IRC and spatial variation of AIB within ROI was the strongest predictor of bone ultimate strength (r = 0.92, n = 19, p < 0.01). In conclusion, our findings suggest that the measurement of two-dimensional parametric maps of ultrasound parameters could yield information on bone status not extractable from single point measurements. This highlights the potential of parametric imaging in osteoporosis diagnostics.

Spectral ratio method to estimate broadband ultrasound attenuation of cortical bonesin vitrousing multiple reflections

Broadband ultrasound attenuation (BUA) is commonly measured by the spectral ratio method. Conventionally BUA is measured in transverse transmission mode where ultrasound signal is recorded with and without the sample. The spectral ratio method was extended to estimate nBUA (BUA normalized by thickness) in axial transmission mode using spectral amplitudes of the primary reflection and multiple reflection, which echoes more than once between the material interfaces within a layer. We performed three experiments. First, reflections were numerically simulated to verify the accuracy of the method. We then applied the method to estimate attenuation of silicon rubber and the cortex of a bovine femur. The center frequency of the transducers is 2.25 MHz. We obtained 93% accuracy for a simulated data set with 10% random noise after bandpass filtering. For the silicon rubber, 15 measurements were collected and the mean attenuation was 6.33 +/- 0.19 dB MHz(-1) cm(-1). For the bovine bone, eight measurements were performed in the middle portion of the femur. The mean attenuation was 4.91 +/- 0.65 dB MHz(-1) cm(-1) and compared well with those reported in the literature. The results demonstrate that the proposed method has the potential to provide a quick, reliable and robust cortical attenuation assessment in vivo.