Non Destructive Characterization of Cortical Bone Micro-Damage by Nonlinear Resonant Ultrasound Spectroscopy

The objective of the study was to evaluate the ability of a nonlinear ultrasound technique, the so-called nonlinear resonant ultrasound spectroscopy (NRUS) technique, for detecting early microdamage accumulation in cortical bone induced by four-point bending fatigue. Small parallelepiped beam-shaped human cortical bone specimens were subjected to cyclic four-point bending fatigue in several steps. The specimens were prepared to control damage localization during four-point bending fatigue cycling and to unambiguously identify resonant modes for NRUS measurements. NRUS measurements were achieved to follow the evolution of the nonlinear hysteretic elastic behavior during fatigue-induced damage. After each fatigue step, a small number of specimens was removed from the protocol and set apart to quantitatively assess the microcrack number density and length using synchrotron radiation micro-computed tomography (SR-µCT). The results showed a significant effect of damage steps on the nonlinear hysteretic elastic behavior. No significant change in the overall length of microcracks was observed in damaged regions compared to the load-free control regions. Only an increased number of shortest microcracks, those in the lowest quartile, was noticed. This was suggestive of newly formed microcracks during the early phases of damage accumulation. The variation of nonlinear hysteretic elastic behavior was significantly correlated to the variation of the density of short microcracks. Our results suggest that the nonlinear hysteretic elastic behavior is sensitive to early bone microdamage. Therefore NRUS technique can be used to monitor fatigue microdamage progression in in vitro experiments.     

Microdamage Caused by Fatigue Loading in Human Cancellous Bone: Relationship to Reductions in Bone Biomechanical Performance

Vertebral fractures associated with osteoporosis are often the result of tissue damage accumulated over time. Microscopic tissue damage (microdamage) generated in vivo is believed to be a mechanically relevant aspect of bone quality that may contribute to fracture risk. Although the presence of microdamage in bone tissue has been documented, the relationship between loading, microdamage accumulation and mechanical failure is not well understood. The aim of the current study was to determine how microdamage accumulates in human vertebral cancellous bone subjected to cyclic fatigue loading. Cancellous bone cores (n = 32) from the third lumbar vertebra of 16 donors (10 male, 6 female, age 76±8.8, mean ± SD) were subjected to compressive cyclic loading at σ/E₀ = 0.0035 (where σ is stress and E₀ is the initial Young’s modulus). Cyclic loading was suspended before failure at one of seven different amounts of loading and specimens were stained for microdamage using lead uranyl acetate. Damage volume fraction (DV/BV) varied from 0.8±0.5% (no loading) to 3.4±2.1% (fatigue-loaded to complete failure) and was linearly related to the reductions in Young’s modulus caused by fatigue loading (r² = 0.60, p<0.01). The relationship between reductions in Young’s modulus and proportion of fatigue life was nonlinear and suggests that most microdamage generation occurs late in fatigue loading, during the tertiary phase. Our results indicate that human vertebral cancellous bone tissue with a DV/BV of 1.5% is expected to have, on average, a Young’s modulus 31% lower than the same tissue without microdamage and is able to withstand 92% fewer cycles before failure than the same tissue without microdamage. Hence, even small amounts of microscopic tissue damage in human vertebral cancellous bone may have large effects on subsequent biomechanical performance.     

Fatigue microcracks that initiate fracture are located near elevated intracortical porosity but not elevated mineralization

In vivo microcracks in cortical bone are typically observed within more highly mineralized interstitial tissue, but postmortem investigations are inherently limited to cracks that did not lead to fracture which may be misleading with respect to understanding fracture mechanisms. We hypothesized that the one fatigue microcrack which initiates fracture is located spatially adjacent to elevated intracortical porosity but not elevated mineralization. Therefore, the spatial correlation between intracortical porosity, elevated mineralization, and fatigue microdamage was investigated by combining, for the first time, sequential, nondestructive, three-dimensional micro-computed tomography (micro-CT) measurements of each in cortical bone specimens subjected to compressive fatigue loading followed by a tensile overload to fracture. Fatigue loading resulted in significant microdamage accumulation and compromised mechanical properties upon tensile overload compared to control specimens. The microdamage that initiated fracture upon tensile overload was able to be identified in all fatigue-loaded specimens using contrast-enhanced micro-CT and registered images. Two-point (or pair) correlation functions revealed a spatial correlation between microdamage at the fracture initiation site and intracortical porosity, but not highly mineralized tissue, confirming the hypothesis. This difference was unique to the fracture initiation site. Intracortical porosity and highly mineralized tissue exhibited a significantly lower and higher probability, respectively, of being located spatially adjacent to all sites of microdamage compared to the fracture initiation site. Therefore, the results of this study suggest that human cortical bone is tolerant of most microcracks, which are generally compartmentalized within the more highly mineralized interstitial tissue, but a single microcrack of sufficient size located in spatial proximity to intracortical porosity can compromise fracture resistance.     

Low bone mass may not be the only cause of skeletal fragility in osteoporosis

Skeletal fragility in postmenopausal osteoporosis is not due solely to reduction in bone mass. This fact explains some of the overlap in bone mineral measurements observed between patients who are fracturing and age- and sex-matched normals who are not. Changes in skeletal architecture and bone remodeling occur with age which can account for some of the fragility. These changes are exaggerated in patients with postmenopausal osteoporosis who are suffering spine fractures. Three abnormalities have been described by histomorphometric methods which can account for skeletal fragility out of proportion to the degree of bone loss. They are: (i) loss of trabecular connectivity such that vertical weight-bearing bars lose their cross-attachments with each other, thus becoming susceptible to buckling; (ii) inefficient and prolonged microdamage repair due to periods of pause in the formation phase of remodeling; and (iii) accumulation of unrepaired microdamage in unremodeled bone tissue in the central part of trabeculae due to reduced osteon wall thickness coupled with maintenance of trabecular thickness. Recognition of these abnormalities should broaden our approach to the study of skeletal fragility in the syndrome of postmenopausal osteoporosis.     

Simulation of vertebral trabecular bone loss using voxel finite element analysis

Trabecular bone loss in human vertebral bone is characterised by thinning and eventual perforation of the horizontal trabeculae. Concurrently, vertical trabeculae are completely lost with no histological evidence of significant thinning. Such bone loss results in deterioration in apparent modulus and strength of the trabecular core. In this study, a voxel-based finite element program was used to model bone loss in three specimens of human vertebral trabecular bone. Three sets of analyses were completed. In Set 1, strain adaptive resorption was modelled, whereby elements which were subject to the lowest mechanical stimulus (principal strain) were removed. In Set 2, both strain adaptive and microdamage mechanisms of bone resorption were included. Perforation of vertical trabeculae occurred due to microdamage resorption of elements with strains that exceeded a damage threshold. This resulted in collapse of the trabecular network under compression loading for two of the specimens tested. In Set 3, the damage threshold strain was gradually increased as bone loss progressed, resulting in reduced levels of microdamage resorption. This mechanism resulted in trabecular architectures in which vertical trabeculae had been perforated and which exhibited similar apparent modulus properties compared to experimental values reported in the literature. Our results indicate that strain adaptive remodelling alone does not explain the deterioration in mechanical properties that have been observed experimentally. Our results also support the hypothesis that horizontal trabeculae are lost principally by strain adaptive resorption, while vertical trabeculae may be lost due to perforation from microdamage resorption followed by rapid strain adaptive resorption of the remaining unloaded trabeculae.     

Theoretical analysis of alendronate and risedronate effects on canine vertebral remodeling and microdamage

Bisphosphonates suppress bone remodeling activity, increase bone volume, and significantly reduce fracture risk in individuals with osteoporosis and other metabolic bone diseases. The objectives of the current study were to develop a mathematical model that simulates control and 1 year experimental results following bisphosphonate treatment (alendronate or risedronate) in the canine fourth lumbar vertebral body, validate the model by comparing simulation predictions to 3 year experimental results, and then use the model to predict potential long term effects of bisphosphonates on remodeling and microdamage accumulation. To investigate the effects of bisphosphonates on bone volume and microdamage, a mechanistic biological model was modified from previous versions to simulate remodeling in a representative volume of vertebral trabecular bone in dogs treated with various doses of alendronate or risedronate, including doses equivalent to those used for treatment of post-menopausal osteoporosis in humans. Bisphosphonates were assumed to affect remodeling by suppressing basic multicellular unit activation and reducing resorption area. Model simulation results for trabecular bone volume fraction, microdamage, and activation frequency following 1 year of bisphosphonate treatment are consistent with experimental measurements. The model predicts that trabecular bone volume initially increases rapidly with 1 year of bisphosphonate treatment, and continues to slowly rise between 1 and 3 years of treatment. The model also predicts that microdamage initially increases rapidly, 0.5–1.5-fold for alendronate or risedronate during the first year of treatment, and reaches its maximum value by 2.5 years before trending downward for all dosages. The model developed in this study suggests that increasing bone volume fraction with long term bisphosphonate treatment may sufficiently reduce strain and damage formation rate so that microdamage does not accumulate above that which is initiated in the first two years of treatment.     

Detection of trabecular bone microdamage by micro-computed tomography

Microdamage is an important component of bone quality and affects bone remodeling. Improved techniques to assess microdamage without the need for histological sectioning would provide insight into the role of microdamage in trabecular bone strength by allowing the spatial distribution of damage within the trabecular microstructure to be measured. Nineteen cylindrical trabecular bone specimens were prepared and assigned to two groups. The specimens in group I were damaged to 3% compressive strain and labeled with BaSO(4). Group II was not loaded, but was labeled with BaSO(4). Micro-computed tomography (Micro-CT) images of the specimens were obtained at 10 microm resolution. The median intensity of the treated bone tissue was compared between groups. Thresholding was also used to measure the damaged area fraction in the micro-CT scans. The histologically measured damaged area fraction, the median CT intensity, and the micro-CT measured damaged area fraction were all higher in the loaded group than in the unloaded group, indicating that the micro-CT images could differentiate the damaged specimen group from the unloaded specimens. The histologically measured damaged area fraction was positively correlated with the micro-CT measured damaged area fraction and with the median CT intensity of the bone, indicating that the micro-CT images can detect microdamage in trabecular bone with sufficient accuracy to differentiate damage levels between samples. This technique provides a means to non-invasively assess the three-dimensional distribution of microdamage within trabecular bone test specimens and could be used to gain insight into the role of trabecular architecture in microdamage formation.     

Microdamage: a cell transducing mechanism based on ruptured osteocyte processes

As a result of underlying pathological diseases, such as osteoporosis, osteopenia, or due to altered loading after joint replacements, bones become more susceptible to microdamage accumulation than those of normal human beings, as are those of athletes who undertake strenuous exercise [Stromsoe, 2004. Fracture fixation problems in osteoporosis. Injury 35, 107-113]. Experimental evidence has linked bone adaptation to microdamage, and to increased cell activity. In this work, we investigated whether microcrack detection is related to rupturing of the cellular material itself due to crack face displacements. Using specific cell staining techniques, it was confirmed that relative crack displacements are capable of tearing cell processes between neighbouring osteocytes. No ruptured cell processes were found near the crack tip where the displacements are less. Rupturing of cell processes due to crack opening and shear displacement is a feasible new mechanism by which bone can detect and estimate the size of a microcrack. Ruptured cell processes may directly secrete passive and active components in the extracellular matrix, triggering a repair response.     

Synchrotron radiation micro-CT at the micrometer scale for the analysis of the three-dimensional morphology of microcracks in human trabecular bone

Bone quality is an important concept to explain bone fragility in addition to bone mass. Among bone quality factors, microdamage which appears in daily life is thought to have a marked impact on bone strength and plays a major role in the repair process. The starting point for all studies designed to further our understanding of how bone microdamage initiate or dissipate energy, or to investigate the impact of age, gender or disease, remains reliable observation and measurement of microdamage. In this study, 3D Synchrotron Radiation (SR) micro-CT at the micrometric scale was coupled to image analysis for the three-dimensional characterization of bone microdamage in human trabecular bone specimens taken from femoral heads. Specimens were imaged by 3D SR micro-CT with a voxel size of 1.4 μm. A new tailored 3D image analysis technique was developed to segment and quantify microcracks. Microcracks from human trabecular bone were observed in different tomographic sections as well as from 3D renderings. New 3D quantitative measurements on the microcrack density and morphology are reported on five specimens. The 3D microcrack density was found between 3.1 and 9.4/mm3 corresponding to a 2D density between 0.55 and 0.76 /mm2. The microcrack length and width measured in 3D on five selected microcrack ranged respectively from 164 μm to 209 μm and 100 μm to 120 μm. This is the first time that various microcracks in unloaded human trabecular bone--from the simplest linear crack to more complex cross-hatch cracks--have been examined and quantified by 3D imaging at this scale. The suspected complex morphology of microcracks is here considerably more evident than in the 2D observations. In conclusion, this technique opens new perspective for the 3D investigation of microcracks and the impact of age, disease or treatment.     

Determination of ultrasound phase velocity in trabecular bone using time dependent phase tracking technique

Ultrasound velocity is one of the key acoustic parameters for noninvasive diagnosis of osteoporosis. Ultrasound phase velocity can be uniquely measured from the phase of the ultrasound signal at a specified frequency. Many previous studies used fast Fourier transform (FFT) to determine the phase velocity, which may cause errors due to the limitations of FFT. The new phase tracking technique applied an adaptive tracking algorithm to detect the time dependent phase and amplitude of the ultrasound signal at a specified frequency. This overcame the disadvantages of FFT to ensure the accuracy of the ultrasound phase velocity. As a result, the new method exhibited high accuracy in the measurement of ultrasound phase velocity of two phantom blocks with the error less than 0.4%. 41 cubic trabecular samples from sheep femoral condyles were used in the study. The phase velocity of the samples using the new method had significantly high correlation to the bulk stiffness of the samples (r = 0.84) compared to the phase velocity measured using fast Fourier transform FFT (r = 0.14). In conclusion, the new method provided an accurate measurement of the ultrasound phase velocity in bone.     

Probabilistic failure analysis of bone using a finite element model of mineral–collagen composites

Microdamage accumulation is a major pathway for energy dissipation during the post-yield deformation of bone. In this study, a two-dimensional probabilistic finite element model of a mineral–collagen composite was developed to investigate the influence of the tissue and ultrastructural properties of bone on the evolution of microdamage from an initial defect in tension. The probabilistic failure analyses indicated that the microdamage progression would be along the plane of the initial defect when the debonding at mineral–collagen interfaces was either absent or limited in the vicinity of the defect. In this case, the formation of a linear microcrack would be facilitated. However, the microdamage progression would be scattered away from the initial defect plane if interfacial debonding takes place at a large scale. This would suggest the possible formation of diffuse damage. In addition to interfacial debonding, the sensitivity analyses indicated that the microdamage progression was also dependent on the other material and ultrastructural properties of bone. The intensity of stress concentration accompanied with microdamage progression was more sensitive to the elastic modulus of the mineral phase and the nonlinearity of the collagen phase, whereas the scattering of failure location was largely dependent on the mineral to collagen ratio and the nonlinearity of the collagen phase. The findings of this study may help understanding the post-yield behavior of bone at the ultrastructural level and shed light on the underlying mechanism of bone fractures.     

Perforation of cancellous bone trabeculae by damage-stimulated remodelling at resorption pits: a computational analysis

Loss of trabeculae in cancellous bone is often attributed to a general decline in the bone mass leading to fracture of the thin trabeculae. It has never been investigated whether trabecular perforation may have any other biomechanical mechanism. In this paper, an alternative hypothesis is proposed and tested using a computational model. Taking it as given that osteoclastic resorption is targeted to microdamage, it is hypothesised that the creation of a resorption cavity during normal bone remodelling could cause a stress-concentration in the bone tissue. If the resorption cavities were excessively deep, as is seen during osteoporosis, then this stress concentration may be sufficient to generate more microdamage so that osteoclasts "chase" newly formed damage leading to perforation. If this were true then we should find that, for a given trabecular thickness, there is a critical depth of resorption cavity such that smaller cavities refill whereas deeper cavities cause microdamage accumulation, continued osteoclast activity, and eventual trabecular perforation. Computer simulation is used to test this hypothesis. Using a remodelling stimulus calculated from both strain and damage and a simplified finite element model of a trabeculum with cavities of different sizes, it is predicted that such a critical depth of resorption cavity does indeed exist. Therefore we suggest that an increase in resorption depth relative to the thickness of trabeculae may be responsible for trabecular perforation during osteoporosis, rather than simply trabecular fracture due to insufficient strength.     

Heterogeneous glycation of cancellous bone and its association with bone quality and fragility

Non-enzymatic glycation (NEG) and enzymatic biochemical processes create crosslinks that modify the extracellular matrix (ECM) and affect the turnover of bone tissue. Because NEG affects turnover and turnover at the local level affects microarchitecture and formation and removal of microdamage, we hypothesized that NEG in cancellous bone is heterogeneous and accounts partly for the contribution of microarchitecture and microdamage on bone fragility. Human trabecular bone cores from 23 donors were subjected to compression tests. Mechanically tested cores as well as an additional 19 cores were stained with lead-uranyl acetate and imaged to determine microarchitecture and measure microdamage. Post-yield mechanical properties were measured and damaged trabeculae were extracted from a subset of specimens and characterized for the morphology of induced microdamage. Tested specimens and extracted trabeculae were quantified for enzymatic and non-enzymatic crosslink content using a colorimetric assay and Ultra-high Performance Liquid Chromatography (UPLC). Results show that an increase in enzymatic crosslinks was beneficial for bone where they were associated with increased toughness and decreased microdamage. Conversely, bone with increased NEG required less strain to reach failure and were less tough. NEG heterogeneously modified trabecular microarchitecture where high amounts of NEG crosslinks were found in trabecular rods and with the mechanically deleterious form of microdamage (linear microcracks). The extent of NEG in tibial cancellous bone was the dominant predictor of bone fragility and was associated with changes in microarchitecture and microdamage.     

Detection of fatigue microdamage in whole rat femora using contrast-enhanced micro-computed tomography

Microdamage in bone tissue is typically studied using destructive, two-dimensional histological techniques. Contrast-enhanced micro-computed tomography (micro-CT) was recently demonstrated to enable non-destructive, three-dimensional (3-D) detection of microdamage in machined cortical and trabecular bone specimens in vitro. However, the accumulation of microdamage in whole bones is influenced by variations in the magnitude and mode of loading due to the complex whole bone morphology. Therefore, the objective of this study was to detect the presence, spatial location, and accumulation of fatigue microdamage in whole rat femora in vitro using micro-CT with a BaSO₄ contrast agent. Microdamage was detected and observed to accumulate at specific spatial locations within the cortex of femora loaded in cyclic three-point bending to a 5% or 10% reduction in secant modulus. The ratio of the segmented BaSO₄ stain volume (SV) to the total volume (TV) of cortical bone was adopted as a measure of damage. The amount of microdamage measured by micro-CT (SV/TV) was significantly greater for both loaded groups compared to the control group (p<0.05), but the difference between loaded groups was not statistically significant. At least one distinct region of microdamage, as indicated by the segmented SV, was observed in 85% of loaded specimens. A specimen-specific finite element model confirmed elevated tensile principal strains localized in regions of tissue corresponding to the accumulated microdamage. These regions were not always located where one might expect a priori based upon Euler–Bernoulli beam theory, demonstrating the utility of contrast-enhanced micro-CT for non-destructive, 3-D detection of fatigue microdamage in whole bones in vitro.     

Observations of microdamage around osteocyte lacunae in bone

Tensile microdamage was examined using laser scanning confocal microscopy in beam specimens of bovine, equine and human long bones loaded in vitro and whole specimens of rat ulnae loaded in vivo. Microcracks were observed to initiate frequently at osteocyte lacunae. The implication is that osteocyte lacunae act as stress concentrating features in bone. This association provides a potential mechanism for the detection of strain and/ or damage by osteocytes in bone.     

Fatigue microdamage in bovine trabecular bone

Microdamage, in the form of small cracks, may accumulate in trabecular bone loaded in fatigue. Specimens of bovine trabecular bone were loaded in compressive fatigue at one of four normalized stresses and loading was stopped after the specimens reached one of six maximum strains. Microdamage was identified using a fluorochrome staining technique, and microdamage parameters, including the number of damaged trabeculae and the damaged area fraction, were measured. No microdamage was observed during loading to strains below the yield strain; at higher strains, all microdamage parameters increased with increasing maximum compressive strain. Few significant differences were observed in the type or amount of microdamage accumulation between specimens loaded to the same maximum strain at different normalized stresses; however, more trabecular fractures were observed at high numbers of cycles, which corresponded to low normalized stresses.     

Postfailure modulus strongly affects microcracking and mechanical property change in human iliac cancellous bone: a study using a 2D nonlinear finite element method

A two-dimensional (2D) finite element (FE) method was used to estimate the ability of bone tissue to sustain damage as a function of postfailure modulus. Briefly, 2D nonlinear compact-tension FE models were created from quantitative back-scattered electron images taken of human iliac crest bone specimens. The effects of different postfailure moduli on predicted microcrack propagation were examined. The 2D FE models were used as surrogates for real bone tissues. The crack number was larger in models with higher postfailure modulus, while mean crack length and area were smaller in these models. The rate of stiffness reduction was greater in the models with lower postfailure modulus. Hence, the current results supported the hypothesis that hard tissue postfailure properties have strong effects on bone microdamage morphology and the rate of change in apparent mechanical properties.     

Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.     

A comparison of the fatigue behavior of human trabecular and cortical bone tissue

The fatigue properties of trabecular bone tissue (single trabeculae) and similarly sized cortical bone specimens from human tibia were experimentally determined on a microstructural level using four-point bending cyclic tests, and they were compared based on modulus, mineral density, and microstructural characteristics. The results showed that trabecular specimens had significantly lower moduli and lower fatigue strength than cortical specimens, despite their higher mineral density values. Fracture surface and microdamage analyses illustrated different fracture and damage patterns between trabecular and cortical bone tissue, depending upon their microstructural characteristics. Based on the results from mechanical tests and qualitative observations, a possible mechanical role of the cement lines in trabecular tissue microfracture was suggested.     

Nonlinear acoustics in pant hoots of common chimpanzees (Pan troglodytes): frequency jumps, subharmonics, biphonation, and deterministic chaos

The pant hoot calls produced by common chimpanzees (Pan troglodytes) are multi-call vocalizations that have figured prominently in investigations of acoustic communication in this species. Although pant hoots are predominantly harmonically structured, they can exhibit an acoustic complexity that has recently been linked to nonlinearity in the vocal-fold dynamics underlying typical mammalian sound production. We examined the occurrence of these sorts of nonlinear phenomena in pant hoot vocalizations, contrasting quieter and lower-pitched "introduction" components with loud and high-pitched "climax" calls in the same bouts. Spectrographic evidence revealed four kinds of nonlinear phenomena, including discrete frequency jumps, subharmonics, biphonation, and deterministic chaos. While these events were virtually never observed during the introduction, they occurred in more than half of the climax calls. Biphonation was by far the most common phenomenon, followed by subharmonics, chaos, and frequency jumps. Individual callers varied in the degree to which their climax calls exhibited nonlinear phenomena, but were consistent in showing more biphonation than other forms. These outcomes show that nonlinear phenomena are routinely present in chimpanzee pant hoots, and help lay the foundation for investigating the function of such events.