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

Composition and micromechanical properties of the femoral neck compact bone in relation to patient age, sex and hip fracture occurrence

Current clinical methods of bone health assessment depend to a great extent on bone mineral density (BMD) measurements. However, these methods only act as a proxy for bone strength and are often only carried out after the fracture occurs. Besides BMD, composition and tissue-level mechanical properties are expected to affect the whole bone's strength and toughness. While the elastic properties of the bone extracellular matrix (ECM) have been extensively investigated over the past two decades, there is still limited knowledge of the yield properties and their relationship to composition and architecture. In the present study, morphological, compositional and micropillar compression bone data was collected from patients who underwent hip arthroplasty. Femoral neck samples from 42 patients were collected together with anonymous clinical information about age, sex and primary diagnosis (coxarthrosis or hip fracture). The femoral neck cortex from the inferomedial region was analyzed in a site-matched manner using a combination of micromechanical testing (nanoindentation, micropillar compression) together with micro-CT and quantitative polarized Raman spectroscopy for both morphological and compositional characterization. Mechanical properties, as well as the sample-level mineral density, were constant over age. Only compositional properties demonstrate weak dependence on patient age: decreasing mineral to matrix ratio (p = 0.02, R2 = 0.13, 2.6 % per decade) and increasing amide I sub-peak ratio I∼1660/I∼1683 (p = 0.04, R2 = 0.11, 1.5 % per decade). The patient's sex and diagnosis did not seem to influence investigated bone properties. A clear zonal dependence between interstitial and osteonal cortical zones was observed for compositional and elastic bone properties (p < 0.0001). Site-matched microscale analysis confirmed that all investigated mechanical properties except yield strain demonstrate a positive correlation with the mineral fraction of bone. The output database is the first to integrate the experimentally assessed microscale yield properties, local tissue composition and morphology with the available patient clinical information. The final dataset was used for bone fracture risk prediction in-silico through the principal component analysis and the Naïve Bayes classification algorithm. The analysis showed that the mineral to matrix ratio, indentation hardness and micropillar yield stress are the most relevant parameters for bone fracture risk prediction at 70 % model accuracy (0.71 AUC). Due to the low number of samples, further studies to build a universal fracture prediction algorithm are anticipated with the higher number of patients (N > 200). The proposed classification algorithm together with the output dataset of bone tissue properties can be used for the future comparison of existing methods to evaluate bone quality as well as to form a better understanding of the mechanisms through which bone tissue is affected by aging or disease.

Elasticity and material anisotropy of lamellar cortical bone in adult bovine tibia characterized via AFM nanoindentation

The research focuses on the evaluation of the mechanical properties of osteonal cortical bone at the lamellar level. Elastic properties of the mid-diaphysis region of the bovine tibia are investigated via cantilever-based nanoindentation at the submicron length scale utilizing Atomic Force Microscopy, where the force-displacement curves are used for the elastic assessment using the Derjaguin-Muller-Toropov model to calculate indentation modulus. Variations of the modulus and the directional mechanical response of the osteonal bone at different distances from the Haversian canal are investigated. Additionally, the effects of demineralization on the indentation modulus are discussed. It was found that in the axial direction, the first and last untreated thick lamella layers show a significant indentation modulus difference compared to all other layers (4.26 ± 0.4 and 4.6 ± 0.3 GPa vs ∼3.5 GPa). On the other hand, the indentation modulus of transverse thick lamella layers shows a periodic variation between ∼3 ± 0.7 GPa and ∼4 ± 0.3 GPa from near the Haversian canal to near the interstitial bone. A periodic variation in the anisotropy ratio was found. Mineral content was quantified via energy-dispersive X-ray microanalysis at different levels of mineralization and shows a positive correlation with the indentation modulus.

Structural and functional heterogeneity of mineralized fibrocartilage at the Achilles tendon-bone insertion

A demanding task of the musculoskeletal system is the attachment of tendon to bone at entheses. This region often presents a thin layer of fibrocartilage (FC), mineralized close to the bone and unmineralized close to the tendon. Mineralized FC deserves increased attention, owing to its crucial anchoring task and involvement in enthesis pathologies. Here, we analyzed mineralized FC and subchondral bone at the Achilles tendon-bone insertion of rats. This location features enthesis FC anchoring tendon to bone and sustaining tensile loads, and periosteal FC facilitating bone-tendon sliding with accompanying compressive and shear forces. Using a correlative multimodal investigation, we evaluated potential specificities in mineral content, fiber organization and mechanical properties of enthesis and periosteal FC. Both tissues had a lower degree of mineralization than subchondral bone, yet used the available mineral very efficiently: for the same local mineral content, they had higher stiffness and hardness than bone. We found that enthesis FC was characterized by highly aligned mineralized collagen fibers even far away from the attachment region, whereas periosteal FC had a rich variety of fiber arrangements. Except for an initial steep spatial gradient between unmineralized and mineralized FC, local mechanical properties were surprisingly uniform inside enthesis FC while a modulation in stiffness, independent from mineral content, was observed in periosteal FC. We interpreted these different structure-property relationships as a demonstration of the high versatility of FC, providing high strength at the insertion (to resist tensile loading) and a gradual compliance at the periosteal surface (to resist contact stresses). STATEMENT OF SIGNIFICANCE: Mineralized fibrocartilage (FC) at entheses facilitates the integration of tendon in bone, two strongly dissimilar tissues. We focus on the structure-function relationships of two types of mineralized FC, enthesis and periosteal, which have clearly distinct mechanical demands. By investigating them with multiple high-resolution methods in a correlative manner, we demonstrate differences in fiber architecture and mechanical properties between the two tissues, indicative of their mechanical roles. Our results are relevant both from a medical viewpoint, targeting a clinically relevant location, as well as from a material science perspective, identifying FC as high-performance versatile composite.

Bone fragility: conceptual framework, therapeutic implications, and COVID-19-related issues

Bone fragility is the susceptibility to fracture even for common loads because of structural, architectural, or material alterations of bone tissue that result in poor bone strength. In osteoporosis, quantitative and qualitative changes in density, geometry, and micro-architecture modify the internal stress state predisposing to fragility fractures. Bone fragility substantially depends on the structural behavior related to the size and shape of the bone characterized by different responses in the load–deformation curve and on the material behavior that reflects the intrinsic material properties of the bone itself, such as yield and fatigue. From a clinical perspective, the measurement of bone density by DXA remains the gold standard for defining the risk of fragility fracture in all population groups. However, non-quantitative parameters, such as macro-architecture, geometry, tissue material properties, and microcracks accumulation can modify the bone’s mechanical strength. This review provides an overview of the role of different contributors to bone fragility and how these factors might be influenced by the use of anti-osteoporotic drugs and by the COVID-19 pandemic.

Closing cones create conical lamellae in secondary osteonal bone

Lamellae are sheets of mineralized collagen 1-20 µm thick, extending over hundreds of µm in bone tissue, occupying bone's structural hierarchy at a level above collagen fibres and osteocytes, and below osteons and trabeculae. Osteons are tubular arrangements of lamellae surrounding central neurovascular canals. Lamellae in osteons are usually described as concentric cylinders based on their annular appearance in transverse section. In this review, I provide a perspective on current understanding of the relationship between geometry of the bone formation front and the shape of lamellae produced at it, reaching the conclusion that the 'closing cone' bone formation front in secondary osteonal remodelling must necessarily result in cone-shaped lamellae in the mature secondary osteon. Secondary osteons replace primary osteons through a tunnelling process of bone turnover, meaning that conical lamellae may become more common in older and damaged bone which is at greatest risk of fracture. Visualization and measurement of three-dimensional lamellar shape over hundreds of microns is needed to provide data for accurate micromechanical simulations. Treating secondary osteonal lamellae as a 'stack of cones' rather than 'nested cylinders' may have important implications for our appreciation of bone's function as a load-bearing tissue and of its behaviour in fracture.

Microscale compressive behavior of hydrated lamellar bone at high strain rates

The increased risk of fracture in the elderly associated with metabolic conditions like osteoporosis poses a significant strain on health care systems worldwide. Due to bone's hierarchical nature, it is necessary to study its mechanical properties and failure mechanisms at several length scales. We conducted micropillar compression experiments on ovine cortical bone to assess the anisotropic mechanical response at the lamellar scale over a wide range of strain rates (10^-4 to 8·10² s^-1). At the microscale, lamellar bone exhibits a strain rate sensitivity similar to what is reported at the macroscale suggesting that it is an intrinsic property of the extracellular matrix. Significant shear band thickening was observed at high strain rates by HRSEM and STEM imaging. This is likely caused by the material's inability to accommodate the imposed deformation by propagation of thin kink bands and shear cracks at high strain rates, leading to shear band thickening and nucleation. The post-yield behavior is strain rate and direction dependent: hardening was observed for transverse oriented micropillars and hardening modulus increases with strain rate by a factor of almost 2, while axially oriented micropillars showed strain softening and an increase of the softening peak width and work to ultimate stress as a function of strain rate. This suggests that for compression at the micrometer scale, energy absorption in bone increases with strain rate. This study highlights the importance of investigating bone strength and post-yield behavior at lower length scales, under hydrated conditions and at clinically relevant strain rates. STATEMENT OF SIGNIFICANCE: We performed micropillar compression experiments of ovine cortical bone at two different orientations and over seven orders of magnitude of strain rate. Experiments were performed under humid condition to mimic the natural conditions of bone in a human body using a newly developed micro-indenter setup. The strain rate sensitivity was found to be of a similar magnitude to what has been reported for higher length scales, suggesting that the strain rate sensitivity is an intrinsic property of the bone extracellular matrix. In addition, localized shear deformation in thick bands was observed for the first time at high strain rates, highlighting the importance of investigating bone under conditions representative of an accident or fall at several length scales.

Cortical and trabecular morphometric properties of the human calvarium

There is currently a gap in the literature that quantitatively describes the complex bone microarchitecture within the diploë (trabecular bone) and cortical layers of the human calvarium. The purpose of this study was to determine the morphometric properties of the diploë and cortical tables of the human calvarium in which key interacting factors of sex, location on the calvarium, and layers of the sandwich structure were considered. Micro-computed tomography (micro-CT) was utilized to capture images at 18 μm resolution of male (n = 26) and female (n = 24) embalmed calvarium specimens in the frontal and parietal regions (N = 50). All images were post-processed and analyzed using vendor bundled CT-Analyzer software to determine the morphometric properties of the diploë and cortical layers. A two-way mixed (repeated measures) analysis of variance (ANOVA) was used to determine diploë morphometric properties accounting for factors of sex and location. A three-way mixed ANOVA was performed to determine cortical morphometric properties accounting for factors of cortical layer (inner and outer table), sex, and location. The study revealed no two-way interaction effects between sex and location on the diploë morphometry except for fractal dimension. Trabecular thickness and separation in the diploë were significantly greater in the male specimens; however, females showed a greater number of trabeculae and fractal dimension on average. Parietal specimens revealed a greater porosity, trabecular separation, and deviation from an ideal plate structure, but a lesser number of trabeculae and connectivity compared to the frontal location. Additionally, the study observed a lower density and greater porosity in the inner cortical layer than the outer which may be due to clear distinctions between each layer's physiological environment. The study provides valuable insight into the quantitative morphometry of the calvarium in which finite element modelers of the skull can refer to when designing detailed heterogenous or subject-specific skull models to effectively predict injury. Furthermore, this study contributes towards the recent developments on physical surrogate models of the skull which require approximate measures of calvarium bone architecture in order to effectively fabricate a model and then accurately simulate a traumatic head impact event.

Microtensile properties and failure mechanisms of cortical bone at the lamellar level

Bone features a remarkable combination of toughness and strength which originates from its complex hierarchical structure and motivates its investigation on multiple length scales. Here, in situ microtensile experiments were performed on dry ovine osteonal bone for the first time at the length scale of a single lamella. The micromechanical response was brittle and revealed larger ultimate tensile strength compared to the macroscale (factor of 2.3). Ultimate tensile strength for axial and transverse specimens was 0.35 ± 0.05 GPa and 0.13 ± 0.02 GPa, respectively. A significantly greater strength anisotropy relative to compression was observed (axial to transverse strength ratio of 2.7:1 for tension, 1.3:1 for compression). Fracture surface and transmission electron microscopic analysis suggested that this may be rationalized by a change in failure mode from fibril-matrix interfacial shearing for axial specimens to fibril-matrix debonding in the transverse direction. An improved version of the classic Hashin's composite failure model was applied to describe lamellar bone strength as a function of fibril orientation. Together with our experimental observations, the model suggests that cortical bone strength at the lamellar level is remarkably tolerant to variations of fibrils orientation of about ±30°. This study highlights the importance of investigating bone's hierarchical organization at several length scales for gaining a deeper understanding of its macroscopic fracture behavior. STATEMENT OF SIGNIFICANCE: Understanding bone deformation and failure behavior at different length scales of its hierarchical structure is fundamental for the improvement of bone fracture prevention, as well as for the development of multifunctional bio-inspired materials combining toughness and strength. The experiments reported in this study shed light on the microtensile properties of dry primary osteonal bone and establish a baseline from which to start further investigations in more physiological conditions. Microtensile specimens were stronger than their macroscopic counterparts by a factor of 2.3. Lamellar bone strength seems remarkably tolerant to variations of the sub-lamellar fibril orientation with respect to the loading direction (±30°). This study underlines the importance of studying bone on all length scales for improving our understanding of bone's macroscopic mechanical response.

Combining polarized Raman spectroscopy and micropillar compression to study microscale structure-property relationships in mineralized tissues

Bone is a natural composite possessing outstanding mechanical properties combined with a lightweight design. The key feature contributing to this unusual combination of properties is the bone hierarchical organization ranging from the nano- to the macro-scale. Bone anisotropic mechanical properties from two orthogonal planes (along and perpendicular to the main bone axis) have already been widely studied. In this work, we demonstrate the dependence of the microscale compressive mechanical properties on the angle between loading direction and the mineralized collagen fibril orientation in the range between 0° and 82°. For this, we calibrated polarized Raman spectroscopy for quantitative collagen fibril orientation determination and validated the method using widely used techniques (small angle X-ray scattering, micro-computed tomography). We then performed compression tests on bovine cortical bone micropillars with known mineralized collagen fibril angles. A strong dependence of the compressive micromechanical properties of bone on the fibril orientation was found with a high degree of anisotropy for both the elastic modulus (E_a/E_t=3.80) and the yield stress (σ_a^y/σ_t^y=2.54). Moreover, the post-yield behavior was found to depend on the MCF orientation with a transition between softening to hardening behavior at approximately 50°. The combination of methods described in this work allows to reliably determine structure-property relationships of bone at the microscale, which may be used as a measure of bone quality.

Compositional and mechanical properties of growing cortical bone tissue: a study of the human fibula

Human cortical bone contains two types of tissue: osteonal and interstitial tissue. Growing bone is not well-known in terms of its intrinsic material properties. To date, distinctions between the mechanical properties of osteonal and interstitial regions have not been investigated in juvenile bone and compared to adult bone in a combined dataset. In this work, cortical bone samples obtained from fibulae of 13 juveniles patients (4 to 18 years old) during corrective surgery and from 17 adult donors (50 to 95 years old) were analyzed. Microindentation was used to assess the mechanical properties of the extracellular matrix, quantitative microradiography was used to measure the degree of bone mineralization (DMB), and Fourier transform infrared microspectroscopy was used to evaluate the physicochemical modifications of bone composition (organic versus mineral matrix). Juvenile and adult osteonal and interstitial regions were analyzed for DMB, crystallinity, mineral to organic matrix ratio, mineral maturity, collagen maturity, carbonation, indentation modulus, indicators of yield strain and tissue ductility using a mixed model. We found that the intrinsic properties of the juvenile bone were not all inferior to those of the adult bone. Mechanical properties were also differently explained in juvenile and adult groups. The study shows that different intrinsic properties should be used in case of juvenile bone investigation.

Functional disuse initiates medullary endosteal micro-architectural impairment in cortical bone characterized by nanoindentation

In this study, we evaluated the effect of functional disuse-induced bone remodeling on its mechanical properties, individually at periosteum and medullary endosteum regions of the cortical bone. Left middle tibiae were obtained from 5-month-old female Sprague-Dawley rats for the baseline control as well as hindlimb suspended (disuse) groups. Micro-nano-mechanical elastic moduli (at lateral region) was evaluated along axial (Z), circumferential (C) and radial (R) orientations using nanoindentation. Results indicated an anisotropic microstructure with axial orientation having the highest and radial orientation with the lowest moduli at periosteum and medullary endosteum for both baseline control as well as disuse groups. Between the groups: at periosteum, an insignificant difference was evaluated for each of the orientations (p > 0.05) and at endosteum, a significant decrease of elastic moduli in the radial (p < 0.0001), circumferential (p < 0.001) and statistically insignificant difference in axial (p > 0.05) orientation. For the moduli ratios between groups: at periosteum, only significant difference in the Z/R (p < 0.05) anisotropy ratio, whereas at endosteum, a statistically significant difference in Z/C (p < 0.001), and Z/R (p < 0.001), as well as C/R (p < 0.05) anisotropy ratios, was evaluated. The results suggested initial bone remodeling impaired bone micro-architecture predominantly at the medullary endosteum with possible alterations in the geometric orientations of collagen and mineral phases inside the bone. The findings could be significant for studying the mechanotransduction pathways involved in maintaining the bone micro-architecture and possibly have high clinical significance for drug use against impairment from functional disuse.

Application of the Johnson-Cook plasticity model in the finite element simulations of the nanoindentation of the cortical bone

The mechanical behavior of the cortical bone in nanoindentation is a complicated mechanical problem. The finite element analysis has commonly been assumed to be the most appropriate approach to this issue. One significant problem in nanoindentation modeling of the elastic-plastic materials is pile-up deformation, which is not observed in cortical bone nanoindentation testing. This phenomenon depends on the work-hardening of materials; it doesn't occur for work-hardening materials, which suggests that the cortical bone could be considered as a work-hardening material. Furthermore, in a recent study [59], a plastic hardening until failure was observed on the micro-scale of a dry ovine osteonal bone samples subjected to micropillar compression. The purpose of the current study was to apply an isotropic hardening model in the finite element simulations of the nanoindentation of the cortical bone to predict its mechanical behavior. The Johnson-Cook (JC) model was chosen as the constitutive model. The finite element modeling in combination with numerical optimization was used to identify the unknown material constants and then the finite element solutions were compared to the experimental results. A good agreement of the numerical curves with the target loading curves was found and no pile-up was predicted. A Design Of Experiments (DOE) approach was performed to evaluate the linear effects of the material constants on the mechanical response of the material. The strain hardening modulus and the strain hardening exponent were the most influential parameters. While a positive effect was noticed with the Young's modulus, the initial yield stress and the strain hardening modulus, an opposite effect was found with the Poisson's ratio and the strain hardening exponent. Finally, the JC model showed a good capability to describe the elastoplastic behavior of the cortical bone.

Mechanics of linear microcracking in trabecular bone

Microcracking in trabecular bone is responsible both for the mechanical degradation and remodeling of the trabecular bone tissue. Recent results on trabecular bone mechanics have demonstrated that bone tissue microarchitecture, tissue elastic heterogeneity and tissue-level mechanical anisotropy all should be considered to obtain detailed information on the mechanical stress state. The present study investigated the influence of tissue microarchitecture, tissue heterogeneity in elasticity and material separation properties and tissue-level anisotropy on the microcrack formation process. Microscale bone models were executed with the extended finite element method. It was demonstrated that anisotropy and heterogeneity of the bone tissue contribute significantly to bone tissue toughness and the resistance of trabecular bone to microcrack formation. The compressive strain to microcrack initiation was computed to increase by a factor of four from an assumed homogeneous isotropic tissue to an assumed anisotropic heterogenous tissue.

Modal analysis of nanoindentation data, confirming that reduced bone turnover may cause increased tissue mineralization/elasticity

It is widely believed that the activities of bone cells at the tissue scale not only govern the size of the vascular pore spaces (and hence, the amount of bone tissue available for actually carrying the loads), but also the characteristics of the extracellular bone matrix itself. In this context, increased mechanical stimulation (in mediolateral regions of human femora, as compared to anteroposterior regions) may lead to increased bone turnover, lower bone matrix mineralization, and therefore lower tissue modulus. On the other hand, resorption-only processes (in endosteal versus periosteal regions) may have the opposite effect. A modal analysis of nanoindentation data obtained on femurs from the Melbourne Femur Research Collection (MFRC) indeed confirms that bone is stiffer in endosteal regions compared to periosteal regions (E̅_endost = 29.34 ± 0.75 GPa >E̅_periost = 24.67 ± 1.63 GPa), most likely due to the aging-related increase in resorption modeling on endosteal surfaces resulting in trabecularization of cortical bone. The results also show that bone is stiffer along the anteroposterior direction compared the mediolateral direction (E̅_anteropost = 28.89 ± 1.08 GPa >E̅_mediolat = 26.03 ± 2.31 GPa), the former being aligned with the neutral bending axis of the femur and, thus, undergoing more resorption modeling and consequently being more mineralized.

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.

Autoimmune arthritis deteriorates bone quantity and quality of periarticular bone in a mouse model of rheumatoid arthritis

This study showed that autoimmune arthritis induces especially severe osteoporosis in the periarticular region adjacent to inflamed joints, suggesting that arthritis increases the fragility fracture risk near inflamed joints, which is frequently observed in patients with RA.

INTRODUCTION

Periarticular osteoporosis near inflamed joints is a hallmark of early rheumatoid arthritis (RA). Here we show that rheumatic inflammation deteriorates the bone quality and bone quantity of periarticular bone, thereby decreasing bone strength and toughness in a mouse model of RA.

METHODS

Female BALB/c mice and SKG mice, a mutant mouse model of autoimmune arthritis on the BALB/c background, were used. At 12 weeks of age, BALB/c mice underwent either Sham surgery or bilateral ovariectomy (OVX), and SKG mice underwent intraperitoneal injection of mannan to induce arthritis. Eight weeks later, the mice were killed and the femurs and tibias were subjected to micro-computed tomography, Fourier transform infrared (FTIR) spectroscopic imaging, X-ray diffraction, histology, and mechanical testing.

RESULTS

SKG mice developed significant trabecular bone loss in both the distal metaphysis of the femur and the lumbar vertebral body, but the extent of the bone loss was more severe in the distal metaphysis. Neither SKG nor OVX mice exhibited changes in the geometry and matrix properties of the diaphysis of the femur, whereas SKG mice, but not OVX mice, did exhibit changes in these properties in the distal metaphysis of the femur. Bone strength and fracture toughness of the distal metaphysis of the tibia adjacent to the inflamed ankle joint were significantly decreased in SKG mice.

CONCLUSIONS

Autoimmune arthritis induces periarticular osteoporosis, characterized by deterioration of cortical bone geometry and quality as well as by trabecular bone loss, leading to severe bone fragility in periarticular bone adjacent to inflamed joints.

Mechanical properties of cortical bone and their relationships with age, gender, composition and microindentation properties in the elderly

The growing incidence of skeletal fractures poses a significant challenge to ageing societies. Since a major part of physiological loading in the lower limbs is carried by cortical bone, it would be desirable to better understand the structure-mechanical property relationships and scale effects in this tissue. This study aimed at assessing whether microindentation properties combined with chemical and morphological information are usable to predict macroscopic elastic and strength properties in a donor- and site-matched manner. Specimens for quasi-static macroscopic tests in tension, compression, and torsion and microindentation were prepared from a cohort of 19 male and 20 female donors (46 to 99 years). All tests were performed under fully hydrated conditions. The chemical composition of the extra-cellular matrix was investigated with Raman spectroscopy. The results of the micro-mechanical tests were combined with morphological and compositional properties using a power law relationship to predict the macro-mechanical results. Microindentation properties were not gender dependent, remarkably constant over age, and showed an overall small variation with standard deviations of approximately 10 %. Similar results were obtained for chemical tissue composition. Macro-mechanical stiffness and strength were significantly related to porosity for all load cases (p<0.05). In case of macroscopic yield strain and work-to-failure this was only true in torsion and compression, respectively. The correlations of macro-mechanical with micro-mechanical, morphological, and chemical properties showed no significance for cement line density, mineralisation, or variations in the microindentation results and were dominated by porosity with a moderate explanatory power of predominately less than 50 %. The results confirm that age, with minor exceptions gender, and small variations in average mineralisation have negligible effect on the tissue microindentation properties of human lamellar bone in the elderly. Furthermore, our findings suggest that microindentation experiments are suitable to predict macroscopic mechanical properties in the elderly only on average and not on a one to one basis. The presented data may help to form a better understanding of the mechanisms of ageing in bone tissue and of the length scale at which they are active. This may be used for future prediction of fracture risk in the elderly.

TECHNOLOGIES FOR STRAIN ASSESSMENT FROM WHOLE BONE TO MINERALIZED OSTEOID LEVEL: A CRITICAL REVIEW

Bone has distinctive structures and mechanical properties at the whole bone, perilacunar and mineralized osteoid levels. A systematic understanding of bone strain magnitudes at different anatomical levels and their internal interactions is the prerequisite to advances in bone mechanobiology. However, due to the intrinsic shortcomings of the strain-measuring technologies, the systematic assessment of bone strain at different anatomical levels under physiological conditions and a deep understanding of their internal interactions are still restricted. To promote technological advances and provide systematic and valuable information for mechanical engineers and bone biomechanical researchers, the most useful methods for measuring bone strain at different anatomical levels are demonstrated in this review, and suggestions for the future development of the technologies and their potential integrated applications are proposed.

NANOINDENTATION AND ASH CONTENT STUDY OF AGE DEPENDENT CHANGES IN PORCINE CORTICAL BONE

Composition-structure-property relations of bone provide fundamental understanding of bone quality. The objective of this paper was to investigate age dependent changes in the composition, structure and mechanical properties of porcine femoral cortical bone at mid-diaphysis region from six age groups (1, 3.5, 6, 12, 30, 48 months). This study was motivated by the fact that limited data is available in the literature on young porcine cortical bone. Nanoindentation technique with Berkovich fluid cell tip was employed to measure the elastic modulus and hardness. Individual lamellae were indented in the longitudinal direction of bone in different microstructural components (osteonal, interstitial and plexiform bone). A grid of indentations was also made on one bone sample to obtain spatial variations in the elastic modulus and hardness. Ash and water content tests were performed to measure water, organic and mineral contents of bone as a function of age. Finally, high resolution micro-computed tomography was used to measure porosity and visualize three-dimensional void structures. We found that the elastic modulus and hardness of bone increased with age but at different rates in each microstructural component. The mineral content increased correspondingly with age while the porosity decreased. The obtained structure, composition, and mechanical properties data give new insights on the age related changes in young cortical bone and can serve as inputs for and validation of multiscale models of bone.

3D scanning SAXS: a novel method for the assessment of bone ultrastructure orientation

The arrangement and orientation of the ultrastructure plays an important role for the mechanical properties of inhomogeneous and anisotropic materials, such as polymers, wood, or bone. However, there is a lack of techniques to spatially resolve and quantify the material's ultrastructure orientation in a macroscopic context. In this study, a new method is presented, which allows deriving the ultrastructural 3D orientation in a quantitative and spatially resolved manner. The proposed 3D scanning small-angle X-ray scattering (3D sSAXS) method was demonstrated on a thin trabecular bone specimen of a human vertebra. A micro-focus X-ray beam from a synchrotron radiation source was used to raster scan the sample for different rotation angles. Furthermore, a mathematical framework was developed, validated and employed to describe the relation between the SAXS data for the different rotation angles and the local 3D orientation and degree of orientation (DO) of the bone ultrastructure. The resulting local 3D orientation was visualized by a 3D orientation map using vector fields. Finally, by applying the proposed 3D scanning SAXS method on consecutive bone sections, a 3D map of the local orientation of a complete trabecular element could be reconstructed for the first time. The obtained 3D orientation map provided information on the bone ultrastructure organization and revealed links between trabecular bone microarchitecture and local bone ultrastructure. More specifically, we observed that trabecular bone ultrastructure is organized in orientation domains of tens of micrometers in size. In addition, it was observed that domains with a high DO were more likely to be found near the surface of the trabecular structure, and domains with lower DO (or transition zones) were located in-between the domains with high DO. The method reproducibility was validated by comparing the results obtained when scanning the sample under different sample tilt angles. 3D orientation maps such as the ones created using 3D scanning SAXS will help to quantify and understand structure-function relationships between bone ultrastructure and bone mechanics. Beyond that, the proposed method can also be used in other research fields such as material sciences, with the aim to locally determine the 3D orientation of material components.

In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

Application of a micro-Brillouin scattering technique to characterize bone in the GHz range

The evaluation of elastic properties of bone matrix has been investigated using several techniques such as nanoindentation and scanning acoustic microscopy (SAM). These techniques make use of good spatial resolution, which can prevent effects due to microstructures at the level of several hundreds of microns. In this paper, micro-Brillouin scattering (μ-BR) is introduced as another possible technique to characterize the elastic properties of bone. This technique is well known as a non-contact and non-destructive method to evaluate viscoelastic properties of transparent materials in the GHz range. Using thin, translucent bone specimens with thicknesses of around 100 μm, and the reflection induced optical geometry, ultrasonic wave velocities in the GHz range were obtained. Because this technique optically measures thermal phonons in the specimen, we can easily measure in-plane anisotropy of wave velocities by rotating the specimen. In a single trabecula, the site matched data between SAM and μ-BR showed good correlation, revealing the applicability of this technique to characterize material properties of bone. Some recent results on the anisotropy in a trabecula and the elasticity evaluation of newly and matured bones are also introduced.

Thermal influence of radiofrequency ablation for bone: an experimental study in normal rabbit bone

OBJECTIVE

To evaluate the heat effects of radiofrequency ablation (RFA) on normal bone by mechanical testing, MRI, and histology.

MATERIALS AND METHODS

The institutional animal care and use committee approved the animal study. Thirty-two adult Japanese white rabbits were included in our study. Bone biopsy needles were inserted from the distal end of the right (RFA side) and the left (control side) femurs using a fluoroscopic guide. A 17-gauge internally cooled RFA electrode with a 2-cm active tip was inserted through the needle to the right femur, and RFA was performed for 12 min using a 200-W generator. Animals were divided into four groups and 8 animals from each group were euthanized on day 1, day 7, day 30, and day 60 after RFA. MRI was performed prior to euthanasia. Three-point bending test was performed to measure flexural strength. Student's t test was used to evaluate for significant differences between RFA and control side for each group. Femurs underwent histological examination by hematoxylin and eosin staining after the bending test.

RESULTS

MRI showed a high-intensity rim around the bone on T2WI. Three-point bending test showed no statistically significant differences (P < 0.05) between the RFA and the control side in any of the groups. Histologically, osteocytes of cortical bone showed cell death, but the lamellar structure was preserved in all groups and bone remodeling was observed.

CONCLUSION

Heat by RFA did not change normal bone strength within 2 months, despite the heat effects in the cortical bone and cell death.

EXPERIMENTAL VALIDATION OF A MEAN FILED MODEL OF MINERALIZED COLLAGEN FIBER ARRAYS AT TWO LEVELS OF HIERARCHY

In the course of this study, stiffness of a fibril array of mineralized collagen fibrils modeled with a mean field method was validated experimentally at site-matched two levels of tissue hierarchy using mineralized turkey leg tendons (MTLT). The applied modeling approaches allowed to model the properties of this unidirectional tissue from nanoscale (mineralized collagen fibrils) to macroscale (mineralized tendon).

At the microlevel, the indentation moduli obtained with a mean field homogenization scheme were compared to the experimental ones obtained with microindentation. At the macrolevel, the macroscopic stiffness predicted with micro finite element (μFE) models was compared to the experimental stiffness measured with uniaxial tensile tests. Elastic properties of the elements in μFE models were injected from the mean field model or two-directional microindentations.

Quantitatively, the indentation moduli can be properly predicted with the mean-field models. Local stiffness trends within specific tissue morphologies are very weak, suggesting additional factors responsible for the stiffness variations. At macrolevel, the μFE models underestimate the macroscopic stiffness, as compared to tensile tests, but the correlations are strong.

Local ultrasonic wave velocities in trabeculae measured by micro-Brillouin scattering

Ultrasonic wave velocities in trabeculae of distal end of bovine femurs were measured using micro-Brillouin scattering (μ-BR). μ-BR allows the measurement of wave velocities in a small area (diameter, 10 μm). Trabecular structure and alignment were evaluated with x-ray micro-computed tomography techniques before μ-BR measurements. Wave velocities in rod-type trabeculae [4.90 × 10(3) m/s with standard deviation (SD) of 0.05 × 10(3) m/s] were higher than those in plate-type trabeculae (4.79 × 10(3) m/s with SD of 0.05 × 10(3) m/s). The elastic properties of trabeculae appeared to change with trabecular type and direction of trabecular alignment.

Investigation of the three-dimensional orientation of mineralized collagen fibrils in human lamellar bone using synchrotron X-ray phase nano-tomography

We investigate the three-dimensional (3-D) organization of mineralized collagen fibrils in human cortical bone based on synchrotron X-ray phase nano-tomography images. In lamellar bone the collagen fibrils are assumed to have a plywood-like arrangement, but due to experimental limitations the 3-D fibril structure has only been deduced from section surfaces so far and the findings have been controversial. Breakthroughs in synchrotron tomographic imaging have given access to direct 3-D information on the bone structure at the nanoscale level. Using an autocorrelation-based orientation measure we confirm that the fibrils are unidirectional in quasi-planes of sub-lamellae and find two specific dominant patterns, oscillating and twisted plywoods coexisting in a single osteon. Both patterns exhibit smooth orientation changes between adjacent quasi-planes. Moreover, we find that the periodic changes in collagen fibril orientation are independent of fluctuations in local mass density. These data improve our understanding of the lamellar arrangement in bone and allow more detailed investigations of structure-function relationships at this scale, providing templates for bio-inspired materials. The presented methodology can be applied to non-destructive 3-D characterization of the sub-micron scale structure of other natural and artificial mineralized biomaterials.

Computational and experimental methodology for site-matched investigations of the influence of mineral mass fraction and collagen orientation on the axial indentation modulus of lamellar bone

Relationships between mineralization, collagen orientation and indentation modulus were investigated in bone structural units from the mid-shaft of human femora using a site-matched design. Mineral mass fraction, collagen fibril angle and indentation moduli were measured in registered anatomical sites using backscattered electron imaging, polarized light microscopy and nano-indentation, respectively. Theoretical indentation moduli were calculated with a homogenization model from the quantified mineral densities and mean collagen fibril orientations. The average indentation moduli predicted based on local mineralization and collagen fibers arrangement were not significantly different from the average measured experimentally with nanoindentation (p=0.9). Surprisingly, no substantial correlation of the measured indentation moduli with tissue mineralization and/or collagen fiber arrangement was found. Nano-porosity, micro-damage, collagen cross-links, non-collagenous proteins or other parameters affect the indentation measurements. Additional testing/simulation methods need to be considered to properly understand the variability of indentation moduli, beyond the mineralization and collagen arrangement in bone structural units.

•

Site-matched assessment of nanoindentation modulus, mineral mass fraction and collagen fibers orientation in human cortical bone sections.

•

Comparison of experimental nanoindentation modulus with its computed equivalent based on the site-matched morphological data.

•

While mean experimental and computed nanoindentation moduli match well, their variations exhibit very weak correlations.

•

Considering factors like nano-porosity and damage may be necessary to understand variability of lamellar stiffness of bone structural units.

•

This is not in conflict with the well known anisotropy associated with the rotated plywood model at the sublamellar scale.

Indentation of poroviscoelastic vocal fold tissue using an atomic force microscope

The elastic properties of the vocal folds (VFs) vary as a function of depth relative to the epithelial surface. The poroelastic anisotropic properties of porcine VFs, at various depths, were measured using atomic force microscopy (AFM)-based indentation. The minimum tip diameter to effectively capture the local properties was found to be 25µm, based on nonlinear laser scanning microscopy data and image analysis. The effects of AFM tip dimensions and AFM cantilever stiffness were systematically investigated. The indentation tests were performed along the sagittal and coronal planes for an evaluation of the VF anisotropy. Hertzian contact theory was used along with the governing equations of linear poroelasticity to calculate the diffusivity coefficient of the tissue from AFM indentation creep testing. The permeability coefficient of the porcine VF was found to be 1.80±0.32×10(-15)m(4)/Ns.

Osteonal lamellae elementary units: lamellar microstructure, curvature and mechanical properties

The mechanical and structural properties of the sublayers of osteonal lamellae were studied. Young's modulus (E) of adjacent individual lamellae was measured by nanoindentation of parallel slices every 1-3 μm, in planes parallel and perpendicular to the osteon axis (OA). In planes parallel to the OA, the modulus of a lamella could vary significantly between sequential slices. Significant modulus variations were also sometimes found on opposing sides of the osteonal canal for the same lamella. These results are rationalized by considerations involving the microstructural organization of the collagen fibrils in the lamellae. Scanning electron microscope imaging of freeze fractured surfaces revealed that the substructure of a single lamella can vary significantly on the opposing sides of the osteonal axis. Using a serial surface view method, parallel planes were exposed every 8-10 nm using a dual-beam microscope. Analysis of the orientations of fibrils revealed that the structure is rotated plywood like, consisting of unidirectional sublayers of fibrils of several orientations, with occasional randomly oriented sublayers. The dependence of the measured mechanical properties of the lamellae on the indentation location may be explained by the observed structure, as well as by the curvature of the osteonal lamellae through simple geometrical-structural considerations. Mechanical advantages arising from the curved laminate structure are discussed.

Symmetrically reduced stiffness and increased extensibility in compression and tension at the mineralized fibrillar level in rachitic bone

In metabolic bone diseases, the alterations in fibrillar level bone-material quality affecting macroscopic mechanical competence are not well-understood quantitatively. Here, we quantify the fibrillar level deformation in cantilever bending in a mouse model for hereditary rickets (Hpr). Microfocus in-situ synchrotron small-angle X-ray scattering (SAXS) combined with cantilever bending was used to resolve nanoscale fibril strain in tensile- and compressive tissue regions separately, with quantitative backscattered scanning electron microscopy used to measure microscale mineralization. Tissue-level flexural moduli for Hpr mice were significantly (p<0.01) smaller compared to wild-type (~5 to 10-fold reduction). At the fibrillar level, the fibril moduli within the tensile and compressive zones were significantly (p<0.05) lower by ~3- to 5-fold in Hpr mice compared to wild-type mice. Hpr mice have a lower mineral content (24.2±2.1Cawt.% versus 27.4±3.3Ca wt.%) and its distribution was more heterogeneous compared to wild-type animals. However, the average effective fibril modulus did not differ significantly (p>0.05) over ages (4, 7 and 10weeks) between tensile and compressive zones. Our results indicate that incompletely mineralized fibrils in Hpr mice have greater deformability and lower moduli in both compression and tension, and those compressive and tensile zones have similar moduli at the fibrillar level.

Influence of osteon area fraction and degree of orientation of HAp crystals on mechanical properties in bovine femur

Cortical bone has a hierarchical structure, spanning from the macrostructure at several millimeters or whole bone level, the microstructure at several hundred micrometers level, to the nanostructure at hydroxyapatite (HAp) crystals and collagen fibrils levels. The aim of the study is to understand the relationship between the HAp crystal orientation and the elastic modulus and the relationship between the osteon area fraction and the deformation behavior of HAp crystals in cortical bone. In the experiments, five strip specimens (40×2×1mm(3)) aligned with the bone axis were taken from the cortical bone of a bovine femur. The degree of c-axis orientation of HAp crystals in the specimens was measured with the X-ray diffraction technique with the imaging plate. To measure the deformation behavior of HAp crystals in the specimens, tensile tests under X-ray irradiation were conducted. The specimens were cut at the X-ray measurement positions and osteon area fraction and porosity at the transverse cross-sections were observed. Further, the volume fraction of HAp of the specimens was measured. Results showed the degree of c-axis orientation of HAp crystals was positively correlated with the elastic modulus of the specimens (r=0.94). The volume fraction of HAp and the porosity showed no statistical correlation with the elastic modulus and the tensile strength. The HAp crystal strain ε(H) increased linearly with the bone tissue strain ε. The average value of ε(H)/ε was 0.69±0.13 and there was no correlation between the osteon area fraction and ε(H)/ε (r=-0.27, p=0.33). The results suggest that the degree of c-axis orientation of HAp crystals affects the elastic modulus and the magnitude of HAp crystal strain does not depend on the osteon area fraction.

Comparative investigation of elastic properties in a trabecula using micro-Brillouin scattering and scanning acoustic microscopy

Micro-Brillouin scattering (μ-BR) and a 200 MHz scanning acoustic microscope (SAM) with similar spatial resolutions were applied to evaluate tissue elastic properties in two directions in a trabecula. Acoustic impedance measured by SAM was in the range of 5-9 Mrayl. Wave velocities determined by μ-BR were in the range of (4.75-5.11) × 10(3) m/s. Both exhibited a similar trend of variation across the trabecula and were significantly correlated (R(2) = 0.63-0.67, p < 0.01). μ-BR is useful for the evaluation of tissue stiffness within a trabecula. Combined with SAM or nanoindentation, it can provide additional information to assess elastic anisotropy at the micro-scale.

Correlative light and scanning X-ray scattering microscopy of healthy and pathologic human bone sections

Scanning small and wide angle X-ray scattering (scanning SWAXS) experiments were performed on healthy and pathologic human bone sections. Via crystallographic tools the data were transformed into quantitative images and as such compared with circularly polarized light (CPL) microscopy images. SWAXS and CPL images allowed extracting information of the mineral nanocrystalline phase embedded, with and without preferred orientation, in the collagen fibrils, mapping local changes at sub-osteon resolution. This favorable combination has been applied for the first time to biopsies of dwarfism syndrome and Paget's disease to shed light onto the cortical structure of natural bone in healthy and pathologic sections.