Bone morphology allows estimation of loading history in a murine model of bone adaptation

Elucidating the mechano-molecular dynamics of TRAP activity using CRISPR/Cas9 mediated fluorescent reporter mice

Osteoclasts are essential for bone remodeling by adapting their resorptive activity in response to their mechanical in vivo environment. However, the molecular mechanisms underlying this process remain unclear. Here, we demonstrated the role of tartrate-resistant acid phosphatase (TRAP, Acp5), a key enzyme secreted by osteoclasts, in bone remodeling and mechanosensitivity. Using CRISPR/Cas9 reporter mice, we demonstrated bone cell reporter (BCRIbsp/Acp5) mice feature fluorescent TRAP-deficient osteoclasts and examined their activity during mechanically driven trabecular bone remodeling. Although BCRIbsp/Acp5 mice exhibited trabecular bone impairments and reduced resorption capacity in vitro, RNA sequencing revealed unchanged levels of key osteoclast-associated genes such as Ctsk, Mmp9, and Calcr. These findings, in conjunction with serum carboxy-terminal collagen crosslinks (CTX) and in vivo mechanical loading outcomes collectively indicated an unaltered bone resorption capacity of osteoclasts in vivo. Furthermore, we demonstrated similar mechanoregulation during trabecular bone remodeling in BCRIbsp/Acp5 and wild-type (WT) mice. Hence, this study provides valuable insights into the dynamics of TRAP activity in the context of bone remodeling and mechanosensation.

Cortical and Trabecular Bone Modeling and Implications for Bone Functional Adaptation in the Mammalian Tibia

Bone modeling involves the addition of bone material through osteoblast-mediated deposition or the removal of bone material via osteoclast-mediated resorption in response to perceived changes in loads by osteocytes. This process is characterized by the independent occurrence of deposition and resorption, which can take place simultaneously at different locations within the bone due to variations in stress levels across its different regions. The principle of bone functional adaptation states that cortical and trabecular bone tissues will respond to mechanical stimuli by adjusting (i.e., bone modeling) their morphology and architecture to mechanically improve their mechanical function in line with the habitual in vivo loading direction. This principle is relevant to various research areas, such as the development of improved orthopedic implants, preventative medicine for osteopenic elderly patients, and the investigation of locomotion behavior in extinct species. In the present review, the mammalian tibia is used as an example to explore cortical and trabecular bone modeling and to examine its implications for the functional adaptation of bones. Following a short introduction and an exposition on characteristics of mechanical stimuli that influence bone modeling, a detailed critical appraisal of the literature on cortical and trabecular bone modeling and bone functional adaptation is given. By synthesizing key findings from studies involving small mammals (rodents), large mammals, and humans, it is shown that examining both cortical and trabecular bone structures is essential for understanding bone functional adaptation. A combined approach can provide a more comprehensive understanding of this significant physiological phenomenon, as each structure contributes uniquely to the phenomenon.

The trabecular architecture of the popliteal sesamoid bone (cyamella) from a New Zealand white rabbit (Oryctolagus cuniculus)

Sesamoid bones are ossified structures that are embedded in tendons near articulation. They consist of an inner trabecular bone architecture surrounded by a thin cortical shell. While the formation of sesamoid bones is probably mainly controlled by genetic factors, the proper development and mineralization of a sesamoid bone depends also on mechanical stimulation. While most sesamoid bones are not loaded directly by other bones during locomotion, they still experience forces directed from the tendon in which they are embedded. In cases when the sesamoid bone is experiencing forces only from a single tendon, such as the cyamella in the rabbit, this may give us a tool to study bone functional adaptation in a relatively simple loading setting. This study investigates the internal trabecular architecture of the popliteal sesamoid bone (cyamellae) in New Zealand white (NZW) rabbits (Oryctolagus cuniculus). Five hind limbs of NZW rabbits were micro-computed tomography scanned and the cortical and trabecular architectures of the cyamellae were evaluated. The results revealed that similar to the patella, the cyamella has a thin cortex and a high trabecular bone volume fraction (BV/TV), which is derived mostly from the high trabecular thickness (Tb.Th). Trabecular BV/TV and Tb.Th were not distributed homogeneously, but they were lower at the periphery and higher closer to the proximal and middle of the cyamella, near the musculotendinous junction. The results also demonstrated that trabeculae tend to align along two recognizable orientations, one with the direction of tensile stresses, in line with the popliteal tendon, and the second bridging the narrow space between the cranial and caudal cortical faces of the bone.

Equivalent loads from the life-cycle of acetabular cages in relation to bone-graft transformation

BACKGROUND AND OBJECTIVES

Bone grafts placed behind acetabular cages change their structure in response to mechanical stimuli. The full consideration of lifestyle loads is extremely resource-intensive, so a method using substitutive loads was proposed to reduce the calculation cost. The aim of the study is to present and prove this method.

METHODS

By means of mechanical equations and using the force vectors from the literature which have the same initial point and their relative frequency, while applying a linear model, the average strain energy density distribution for all load cases can be calculated, compiling a matrix from the external loads. From the elements of this matrix, three substitutive load vectors can be calculated, which can be proven to produce the same strain energy density distribution by averaging their effects. The feasibility of using this to model the transformation of bone grafts placed behind acetabular cages is demonstrated with a finite element model, along with a reference calculation.

RESULTS

The substitutive load vectors could be calculated in closed form and the simulations showed that they produced a similar density distribution to the reference model with a numerical calculation error range. Accordingly, the density distribution calculated from bone graft transformation is almost the same.

CONCLUSIONS

In addition to the aforementioned linearity and the same initial point limitations, the applied method is able to produce the substitutive load vectors with which the calculation of the strain energy density distribution and the bone graft's new density distributions can be carried out faster.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Comparison of bone density patterns of the subaxial spine between chimpanzees and gorillas - A case study

Case study on the bone density pattern of subaxial vertebral column in African apes.

INTRODUCTION

African apes have been noted to experience fewer back ailments than humans and to have higher vertebral bone density. Yet, research on the subject is quite limited and has usually included only one or few vertebrae. However, to understand vertebral column as whole and how posture and locomotion might have affected it, we need to know how bone density varies between adjacent vertebrae.

MATERIALS AND METHODS

Bone density in the vertebral body was measured for all subaxial vertebrae of five specimens including two Pan troglodytes (1 male and 1 female) and three Gorilla gorilla (2 males and 1 female) using peripheral quantitative computed tomography (pQCT).

RESULTS

The results tentatively indicated differences between species, especially in the trabecular density of the cervical segment and support the need for further studies on this subject.

Personalized loading conditions for homogenized finite element analysis of the distal sections of the radius

The microstructure of trabecular bone is known to adapt its morphology in response to mechanical loads for achieving a biomechanical homeostasis. Based on this form-function relationship, previous investigators either simulated the remodeling of bone to predict the resulting density and architecture for a specific loading or retraced physiological loading conditions from local density and architecture. The latter inverse approach includes quantifying bone morphology using computed tomography and calculating the relative importance of selected load cases by minimizing the fluctuation of a tissue loading level metric. Along this concept, the present study aims at identifying an optimal, personalized, multiaxial load case at the distal section of the human radius using in vivo HR-pQCT-based isotropic, homogenized finite element (hFE) analysis. The dataset consisted of HR-pQCT reconstructions of the 20 mm most distal section of 21 human fresh-frozen radii. We simulated six different unit canonical load cases (FX palmar-dorsal force, FY ulnar-radial force, FZ distal-proximal force, MX moment about palmar-dorsal, MY moment about ulnar-radial, MZ moment about distal-proximal) using a simplified and efficient hFE method based on a single isotropic bone phase. Once we used a homogeneous mean density (shape model) and once the original heterogeneous density distribution (shape + density model). Using an analytical formulation, we minimized the deviation of the resulting strain tensors ε(x) to a hydrostatic compressive reference strain ε0, once for the 6 degrees of freedom (DOF) optimal (OPT) load case and for all individual 1 DOF load cases (FX, FY, FZ, MX, MY, MZ). All seven load cases were then extended in the nonlinear regime using the scaled displacements of the linear load cases as loading boundary conditions (MAX). We then compared the load cases and models for their objective function (OF) values, the stored energies and their ultimate strength using a specific torsor norm. Both shape and shape + density linear-optimized OPT models were dominated by a positive force in the z-direction (FZ). Transversal force DOFs were close to zero and mean moment DOFs were different depending on the model type. The inclusion of density distribution increased the influence and changed direction of MX and MY, while MZ was small in both models. The OPT load case had 12-15% lower objective function (OF) values than the FZ load case, depending on the model. Stored energies at the optimum were consistently 142-178% higher for the OPT load case than for the FZ load case. Differences in the nonlinear response maximum torsor norm ‖t‖ were heterogeneous, but consistently higher for OPT_MAX than FZ_MAX. We presented the proof of concept of an optimization procedure to estimate patient-specific loading conditions for hFE methods. In contrast to similar models, we included canonical load cases in all six DOFs and used a strain metric that favors hydrostatic compression. Based on a biomechanical analysis of the distal joint surfaces at the radius, the estimated load directions are plausible. For our dataset, the resulting OPT load case is close to the standard axial compression boundary conditions, usually used in HR-pQCT-based FE analysis today. But even using the present simplified hFE model, the optimized linear six DOF load case achieves a more homogeneous tissue loading and can absorb more than twice the energy than the standard uniaxial load case. The ultimate strength calculated with a torsor norm was consistently higher for the 6-DOF nonlinear model (OPT_MAX) than for the 1-DOF nonlinear uniaxial model (FZ_MAX). Defining patient-specific boundary conditions may decrease angulation errors during CT measurements and improve repeatability as well as reproducibility of bone stiffness and strength estimated by HR-pQCT-based hFE analysis. These results encourage the extension of the present method to anisotropic hFE models and their application to repeatability data sets to test the hypothesis of reduced angulation errors during measurement.

Trabecular bone remodeling in the aging mouse: A micro-multiphysics agent-based in silico model using single-cell mechanomics

Bone remodeling is regulated by the interaction between different cells and tissues across many spatial and temporal scales. Notably, in silico models are regarded as powerful tools to further understand the signaling pathways that regulate this intricate spatial cellular interplay. To this end, we have established a 3D multiscale micro-multiphysics agent-based (micro-MPA) in silico model of trabecular bone remodeling using longitudinal in vivo data from the sixth caudal vertebra (CV6) of PolgA(D257A/D257A) mice, a mouse model of premature aging. Our in silico model includes a variety of cells as single agents and receptor-ligand kinetics, mechanomics, diffusion and decay of cytokines which regulate the cells' behavior. We highlighted its capabilities by simulating trabecular bone remodeling in the CV6 of five mice over 4 weeks and we evaluated the static and dynamic morphometry of the trabecular bone microarchitecture. Based on the progression of the average trabecular bone volume fraction (BV/TV), we identified a configuration of the model parameters to simulate homeostatic trabecular bone remodeling, here named basal. Crucially, we also produced anabolic, anti-anabolic, catabolic and anti-catabolic responses with an increase or decrease by one standard deviation in the levels of osteoprotegerin (OPG), receptor activator of nuclear factor kB ligand (RANKL), and sclerostin (Scl) produced by the osteocytes. Our results showed that changes in the levels of OPG and RANKL were positively and negatively correlated with the BV/TV values after 4 weeks in comparison to basal levels, respectively. Conversely, changes in Scl levels produced small fluctuations in BV/TV in comparison to the basal state. From these results, Scl was deemed to be the main driver of equilibrium while RANKL and OPG were shown to be involved in changes in bone volume fraction with potential relevance for age-related bone features. Ultimately, this micro-MPA model provides valuable insights into how cells respond to their local mechanical environment and can help to identify critical pathways affected by degenerative conditions and ageing.

Quantitative Load Dependency Analysis of Local Trabecular Bone Microstructure to Understand the Spatial Characteristics in the Synthetic Proximal Femur

Analysis of the dependency of the trabecular structure on loading conditions is essential for understanding and predicting bone structure formation. Although previous studies have investigated the relationship between loads and structural adaptations, there is a need for an in-depth analysis of this relationship based on the bone region and load specifics. In this study, the load dependency of the trabecular bone microstructure for twelve regions of interest (ROIs) in the synthetic proximal femur was quantitatively analyzed to understand the spatial characteristics under seven different loading conditions. To investigate the load dependency, a quantitative measure, called the load dependency score (LDS), was established based on the statistics of the strain energy density (SED) distribution. The results showed that for the global model and epiphysis ROIs, bone microstructures relied on the multiple-loading condition, whereas the structures in the metaphysis depended on single or double loads. These results demonstrate that a given ROI is predominantly dependent on a particular loading condition. The results confirm that the dependency analysis of the load effects for ROIs should be performed both qualitatively and quantitatively.

Mechanostat parameters estimated from time-lapsed in vivo micro-computed tomography data of mechanically driven bone adaptation are logarithmically dependent on loading frequency

Mechanical loading is a key factor governing bone adaptation. Both preclinical and clinical studies have demonstrated its effects on bone tissue, which were also notably predicted in the mechanostat theory. Indeed, existing methods to quantify bone mechanoregulation have successfully associated the frequency of (re)modeling events with local mechanical signals, combining time-lapsed in vivo micro-computed tomography (micro-CT) imaging and micro-finite element (micro-FE) analysis. However, a correlation between the local surface velocity of (re)modeling events and mechanical signals has not been shown. As many degenerative bone diseases have also been linked to impaired bone (re)modeling, this relationship could provide an advantage in detecting the effects of such conditions and advance our understanding of the underlying mechanisms. Therefore, in this study, we introduce a novel method to estimate (re)modeling velocity curves from time-lapsed in vivo mouse caudal vertebrae data under static and cyclic mechanical loading. These curves can be fitted with piecewise linear functions as proposed in the mechanostat theory. Accordingly, new (re)modeling parameters can be derived from such data, including formation saturation levels, resorption velocity moduli, and (re)modeling thresholds. Our results revealed that the norm of the gradient of strain energy density yielded the highest accuracy in quantifying mechanoregulation data using micro-finite element analysis with homogeneous material properties, while effective strain was the best predictor for micro-finite element analysis with heterogeneous material properties. Furthermore, (re)modeling velocity curves could be accurately described with piecewise linear and hyperbola functions (root mean square error below 0.2 µm/day for weekly analysis), and several (re)modeling parameters determined from these curves followed a logarithmic relationship with loading frequency. Crucially, (re)modeling velocity curves and derived parameters could detect differences in mechanically driven bone adaptation, which complemented previous results showing a logarithmic relationship between loading frequency and net change in bone volume fraction over 4 weeks. Together, we expect this data to support the calibration of in silico models of bone adaptation and the characterization of the effects of mechanical loading and pharmaceutical treatment interventions in vivo.

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models

Inverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from − 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.

How Artificial Intelligence and Machine Learning Is Assisting Us to Extract Meaning from Data on Bone Mechanics?

Dramatic advancements in interdisciplinary research with the fourth paradigm of science, especially the implementation of computer science, nourish the potential for artificial intelligence (AI), machine learning (ML), and artificial neural network (ANN) algorithms to be applied to studies concerning mechanics of bones. Despite recent enormous advancement in techniques, gaining deep knowledge to find correlations between bone shape, material, mechanical, and physical responses as well as properties is a daunting task. This is due to both complexity of the material itself and the convoluted shapes that this complex material forms. Moreover, many uncertainties and ambiguities exist concerning the use of traditional computational techniques that hinders gaining a full comprehension of this advanced biological material. This book chapter offers a review of literature on the use of AI, ML, and ANN in the study of bone mechanics research. A main question as to why to implement AI and ML in the mechanics of bones is fully addressed and explained. This chapter also introduces AI and ML and elaborates on the main features of ML algorithms such as learning paradigms, subtypes, main ideas with examples, performance metrics, training algorithms, and training datasets. As a frequently employed ML algorithm in bone mechanics, feedforward ANNs are discussed to make their taxonomy and working principles more readily comprehensible to researchers. A summary as well as detailed review of papers that employed ANNs to learn from collected data on bone mechanics are presented. Reviewing literature on the use of these data-driven tools is essential since their wider application has the potential to: improve clinical assessments enabling real-time simulations; avoid and/or minimize injuries; and, encourage early detection of such injuries in the first place.

Bone remodelling patterns around orthodontic mini-implants migrating in bone: an experimental study in rat vertebrae

Background

Orthodontic implant migration has been clinically observed in presence of continuous loading forces. Recent studies indicate that osteocytes play a crucial role in this phenomenon.

Objectives

Aim of this study was to investigate local osteocytic gene expression, protein expression, and bone micro-structure in peri-implant regions of pressure and tension.

Material and methods

The present work reports a complementary analysis to a previous micro-computed tomography study. Two customized mini-implants were placed in one caudal rat vertebra and connected by a nickel–titanium contraction spring generating different forces (i.e. 0, 0.5, 1.0, and 1.5 N). Either at 2 or 8 weeks, the vertebrae were harvested and utilized for 1. osteocytic gene expression using laser capture micro-dissection on frozen sections coupled with qPCR, 2. haematoxylin–eosin staining for qualitative and quantitative analyses, 3. immunofluorescence staining and analysis, and 4. bone-to-implant contact on undecalcified samples.

Results

At the two time points for all the performed analyses no significant differences were observed with respect to the applied force magnitudes and cell harvesting localization. However, descriptive histological analysis revealed remarkable bone remodelling at 2 weeks of loading. At 8 weeks the implants were osseointegrated and, especially in 1.0 and 1.5 N groups, newly formed bone presented a characteristic load bearing architecture with trabecula oriented in the direction of the loading.

Conclusions

The present study confirmed that stress-induced bone remodelling is the biological mechanism of orthodontic implant migration. Bone apposition was found at ‘tension’ and ‘pressure’ sites thus limiting implant migration over time.

Bone Mechanoregulation Allows Subject-Specific Load Estimation Based on Time-Lapsed Micro-CT and HR-pQCT in Vivo

Patients at high risk of fracture due to metabolic diseases frequently undergo long-term antiresorptive therapy. However, in some patients, treatment is unsuccessful in preventing fractures or causes severe adverse health outcomes. Understanding load-driven bone remodelling, i.e., mechanoregulation, is critical to understand which patients are at risk for progressive bone degeneration and may enable better patient selection or adaptive therapeutic intervention strategies. Bone microarchitecture assessment using high-resolution peripheral quantitative computed tomography (HR-pQCT) combined with computed mechanical loads has successfully been used to investigate bone mechanoregulation at the trabecular level. To obtain the required mechanical loads that induce local variances in mechanical strain and cause bone remodelling, estimation of physiological loading is essential. Current models homogenise strain patterns throughout the bone to estimate load distribution in vivo, assuming that the bone structure is in biomechanical homoeostasis. Yet, this assumption may be flawed for investigating alterations in bone mechanoregulation. By further utilising available spatiotemporal information of time-lapsed bone imaging studies, we developed a mechanoregulation-based load estimation (MR) algorithm. MR calculates organ-scale loads by scaling and superimposing a set of predefined independent unit loads to optimise measured bone formation in high-, quiescence in medium-, and resorption in low-strain regions. We benchmarked our algorithm against a previously published load history (LH) algorithm using synthetic data, micro-CT images of murine vertebrae under defined experimental in vivo loadings, and HR-pQCT images from seven patients. Our algorithm consistently outperformed LH in all three datasets. In silico-generated time evolutions of distal radius geometries (n = 5) indicated significantly higher sensitivity, specificity, and accuracy for MR than LH (p < 0.01). This increased performance led to substantially better discrimination between physiological and extra-physiological loading in mice (n = 8). Moreover, a significantly (p < 0.01) higher association between remodelling events and computed local mechanical signals was found using MR [correct classification rate (CCR) = 0.42] than LH (CCR = 0.38) to estimate human distal radius loading. Future applications of MR may enable clinicians to link subtle changes in bone strength to changes in day-to-day loading, identifying weak spots in the bone microstructure for local intervention and personalised treatment approaches.

Determination of the representative static loads for cyclically repeated dynamic loads: a case study of bone remodeling simulation with gait loads

BACKGROUND AND OBJECTIVE

Bone has the self-optimizing capability to adjust its structure in order to efficiently support external loads. Bone remodeling simulations have been developed to reflect the above characteristics in a more effective way. In most studies, however, only a set of static loads have been empirically determined although both static and dynamic loads affect bone remodeling phenomenon. The goal of this study is to determine the representative static loads (RSLs) to efficiently consider the statically equivalent effect of cyclically repeated dynamic loads on bone remodeling simulation.

METHODS

Based on the concept of two-scale approach, the RSLs for the gait cycles are determined from five subjects. First, the gait profiles at the hip joint are selected from the public database and then are preprocessed. The finite element model of the proximal femur is constructed from the clinical CT scan data to determine the strain energy distribution during the gait cycles. An optimization problem is formulated to determine the candidate static loads that minimize the errors of the spatial strain energy distribution for five gait profiles. Then, all candidate static loads from five gait profiles are partitioned into multiple clusters. The RSLs and the corresponding coefficients can be determined at the center of the densest cluster. For verification, topology optimization is separately conducted with the whole gait cycle (reference), empirically determined loads (conventional), and the RSLs (proposed). The strain energy density-based bone remodeling simulation is also conducted for another comparison.

RESULTS

For the gait loads, the use of the RSLs enables a 99% reduction of the function calls with negligible errors in the bone spatial distribution (6.75% for two representative static loads and 6.24% for three representative static loads) and apparent stiffness (4.84% for two representative static loads and 4.47% for three representative static loads), compared with the use of a whole gait cycle as reference.

CONCLUSION

This study shows the feasibility of the RSLs and provides a theoretical foundation for investigating the relationship between static and dynamic loads in the aspect of bone remodeling simulation.

Application of subject-specific adaptive mechanical loading for bone healing in a mouse tail vertebral defect

Methods to repair bone defects arising from trauma, resection, or disease, continue to be sought after. Cyclic mechanical loading is well established to influence bone (re)modelling activity, in which bone formation and resorption are correlated to micro-scale strain. Based on this, the application of mechanical stimulation across a bone defect could improve healing. However, if ignoring the mechanical integrity of defected bone, loading regimes have a high potential to either cause damage or be ineffective. This study explores real-time finite element (rtFE) methods that use three-dimensional structural analyses from micro-computed tomography images to estimate effective peak cyclic loads in a subject-specific and time-dependent manner. It demonstrates the concept in a cyclically loaded mouse caudal vertebral bone defect model. Using rtFE analysis combined with adaptive mechanical loading, mouse bone healing was significantly improved over non-loaded controls, with no incidence of vertebral fractures. Such rtFE-driven adaptive loading regimes demonstrated here could be relevant to clinical bone defect healing scenarios, where mechanical loading can become patient-specific and more efficacious. This is achieved by accounting for initial bone defect conditions and spatio-temporal healing, both being factors that are always unique to the patient.

Mechano-Regulation of Trabecular Bone Adaptation Is Controlled by the Local in vivo Environment and Logarithmically Dependent on Loading Frequency

It is well-established that cyclic, but not static, mechanical loading has anabolic effects on bone. However, the function describing the relationship between the loading frequency and the amount of bone adaptation remains unclear. Using a combined experimental and computational approach, this study aimed to investigate whether trabecular bone mechano-regulation is controlled by mechanical signals in the local in vivo environment and dependent on loading frequency. Specifically, by combining in vivo micro-computed tomography (micro-CT) imaging with micro-finite element (micro-FE) analysis, we monitored the changes in microstructural as well as the mechanical in vivo environment [strain energy density (SED) and SED gradient] of mouse caudal vertebrae over 4 weeks of either cyclic loading at varying frequencies of 2, 5, or 10 Hz, respectively, or static loading. Higher values of SED and SED gradient on the local tissue level led to an increased probability of trabecular bone formation and a decreased probability of trabecular bone resorption. In all loading groups, the SED gradient was superior in the determination of local bone formation and resorption events as compared to SED. Cyclic loading induced positive net (re)modeling rates when compared to sham and static loading, mainly due to an increase in mineralizing surface and a decrease in eroded surface. Consequently, bone volume fraction increased over time in 2, 5, and 10 Hz (+15%, +21% and +24%, p ≤ 0.0001), while static loading led to a decrease in bone volume fraction (-9%, p ≤ 0.001). Furthermore, regression analysis revealed a logarithmic relationship between loading frequency and the net change in bone volume fraction over the 4 week observation period (R 2 = 0.74). In conclusion, these results suggest that trabecular bone adaptation is regulated by mechanical signals in the local in vivo environment and furthermore, that mechano-regulation is logarithmically dependent on loading frequency with frequencies below a certain threshold having catabolic effects, and those above anabolic effects. This study thereby provides valuable insights toward a better understanding of the mechanical signals influencing trabecular bone formation and resorption in the local in vivo environment.

Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone

Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.

Baby steps towards linking calcaneal trabecular bone ontogeny and the development of bipedal human gait

Trabecular bone structure in adulthood is a product of a process of modelling during ontogeny and remodelling throughout life. Insight into ontogeny is essential to understand the functional significance of trabecular bone structural variation observed in adults. The complex shape and loading of the human calcaneus provides a natural experiment to test the relationship between trabecular morphology and locomotor development. We investigated the relationship between calcaneal trabecular bone structure and predicted changes in loading related to development of gait and body size in growing children. We sampled three main trabecular regions of the calcanei using micro-computed tomography scans of 35 individuals aged between neonate to adult from the Norris Farms #36 site (1300 AD, USA) and from Cambridge (1200-1500 AD, UK). Trabecular properties were calculated in volumes of interest placed beneath the calcaneocuboid joint, plantar ligaments, and posterior talar facet. At birth, thin trabecular struts are arranged in a dense and relatively isotropic structure. Bone volume fraction strongly decreases in the first year of life, whereas anisotropy and mean trabecular thickness increase. Dorsal compressive trabecular bands appear around the onset of bipedal walking, although plantar tensile bands develop prior to predicted propulsive toe-off. Bone volume fraction and anisotropy increase until the age of 8, when gait has largely matured. Connectivity density gradually reduces, whereas trabeculae gradually thicken from birth until adulthood. This study demonstrates that three different regions of the calcaneus develop into distinct adult morphologies through varying developmental trajectories. These results are similar to previous reports of ontogeny in human long bones and are suggestive of a relationship between the mechanical environment and trabecular bone architecture in the human calcaneus during growth. However, controlled experiments combined with more detailed biomechanical models of gait maturation are necessary to establish skeletal markers linking growth to loading. This has the potential to be a novel source of information for understanding loading levels, activity patterns, and perhaps life history in the fossil record.

A novel algorithm to predict bone changes in the mouse tibia properties under physiological conditions

Understanding how bone adapts to mechanical stimuli is fundamental for optimising treatments against musculoskeletal diseases in preclinical studies, but the contribution of physiological loading to bone adaptation in mouse tibia has not been quantified so far. In this study, a novel mechanistic model to predict bone adaptation based on physiological loading was developed and its outputs were compared with longitudinal scans of the mouse tibia. Bone remodelling was driven by the mechanical stimuli estimated from micro-FEA models constructed from micro-CT scans of C57BL/6 female mice (N = 5) from weeks 14 and 20 of age, to predict bone changes in week 16 or 22. Parametric analysis was conducted to evaluate the sensitivity of the models to subject-specific or averaged parameters, parameters from week 14 or week 20, and to strain energy density (SED) or maximum principal strain (ε_maxprinc). The results at week 20 showed no significant difference in bone densitometric properties between experimental and predicted images across the tibia for both stimuli, and 59% and 47% of the predicted voxels matched with the experimental sites in apposition and resorption, respectively. The model was able to reproduce regions of bone apposition in both periosteal and endosteal surfaces (70% and 40% for SED and ε_maxprinc, respectively), but it under-predicted the experimental sites of resorption by over 85%. This study shows for the first time the potential of a subject-specific mechanoregulation algorithm to predict bone changes in a mouse model under physiological loading. Nevertheless, the weak predictions of resorption suggest that a combined stimulus or biological stimuli should be accounted for in the model.

Prediction of load in a long bone using an artificial neural network prediction algorithm

The hierarchical nature of bone makes it a difficult material to fully comprehend. The equine third metacarpal (MC3) bone experiences nonuniform surface strains, which are a measure of displacement induced by loads. This paper investigates the use of an artificial neural network expert system to quantify MC3 bone loading. Previous studies focused on determining the response of bone using load, bone geometry, mechanical properties, and constraints as input parameters. This is referred to as a forward problem and is generally solved using numerical techniques such as finite element analysis (FEA). Conversely, an inverse problem has to be solved to quantify load from the measurements of strain and displacement. Commercially available FEA packages, without manipulating their underlying algebraic formulae, are incapable of completing a solution to the inverse problem. In this study, an artificial neural network (ANN) was employed to quantify the load required to produce the MC3 displacement and surface strains determined experimentally. Nine hydrated MC3 bones from thoroughbred horses were loaded in compression in an MTS machine. Ex-vivo experiments measured strain readings from one three-gauge rosette and three distinct single-element gauges at different locations on the MC3 midshaft, associated displacement, and load exposure time. Horse age and bone side (left or right limb) were also recorded for each MC3 bone. This information was used to construct input variables for the ANN model. The ability of this expert system to predict the MC3 loading was investigated. The ANN prediction offered excellent reliability for the prediction of load in the MC3 bones investigated, i.e. R² ≥ 0.98.

Are we crying Wolff? 3D printed replicas of trabecular bone structure demonstrate higher stiffness and strength during off-axis loading

Roux's principle of bone functional adaptation postulates that bone tissue, and particularly trabecular bone tissue, responds to mechanical stimuli by adjusting (modeling) its architecture accordingly. Hence, it predicts that the new modeled trabecular structure is mechanically improved (stiffer and stronger) in line with the habitual in vivo loading direction. While previous studies found indirect evidence to support this theory, direct support was so far unattainable. This is attributed to the fact that each trabecular bone is unique, and that trabecular bone tissue tends to be damaged during mechanical testing. Consequently, a unique modeled trabecular structure can be mechanically tested only along one direction and a comparison to other directions for that specific structure is impossible. To address this issue, we have 3D printed 10 replicas of a trabecular structure from a sheep talus cropped along the 3 principal axes of the bone and in line with the principal direction of loading (denoted on-axis model). Next, we have rotated the same cropped trabecular structure in increments of 10° up to 90° to the bone principal direction of loading (denoted off-axis models) and printed 10 replicas of each off-axis model. Finally, all on-axis and off-axis 3D printed replicas were loaded in compression until failure and trabecular structure stiffness and strength were calculated. Contrary to our prediction, and conflicting with Roux's principle of bone functional adaptation, we found that a trabecular structure loaded off-axis tended to have higher stiffness and strength values when compared to the same trabecular structure loaded on-axis. These unexpected results may not disprove Roux's principle of bone functional adaptation, but they do imply that trabecular bone adaptation may serve additional purposes than simply optimizing bone structure to one principal loading scenario and this suggests that we still don't fully understand bone modeling in its entirety.

Inverse remodelling algorithm identifies habitual manual activities of primates based on metacarpal bone architecture

Previously, a micro-finite element (micro-FE)-based inverse remodelling method was presented in the literature that reconstructs the loading history of a bone based on its architecture alone. Despite promising preliminary results, it remains unclear whether this method is sensitive enough to detect differences of bone loading related to pathologies or habitual activities. The goal of this study was to test the sensitivity of the inverse remodelling method by predicting joint loading histories of metacarpal bones of species with similar anatomy but clearly distinct habitual hand use. Three groups of habitual hand use were defined using the most representative primate species: manipulation (human), suspensory locomotion (orangutan), and knuckle-walking locomotion (bonobo, chimpanzee, gorilla). Nine to ten micro-computed tomography scans of each species ([Formula: see text] in total) were used to create micro-FE models of the metacarpal head region. The most probable joint loading history was predicted by optimally scaling six load cases representing joint postures ranging from [Formula: see text] (extension) to [Formula: see text] (flexion). Predicted mean joint load directions were significantly different between knuckle-walking and non-knuckle-walking groups ([Formula: see text]) and in line with expected primary hand postures. Mean joint load magnitudes tended to be larger in species using their hands for locomotion compared to species using them for manipulation. In conclusion, this study shows that the micro-FE-based inverse remodelling method is sensitive enough to detect differences of joint loading related to habitual manual activities of primates and might, therefore, be useful for palaeoanthropologists to reconstruct the behaviour of extinct species and for biomedical applications such as detecting pathological joint loading.

Virtual supersampling as post-processing step preserves the trabecular bone morphometry in human peripheral quantitative computed tomography scans

In the clinical field of diagnosis and monitoring of bone diseases, high-resolution peripheral quantitative computed tomography (HR-pQCT) is an important imaging modality. It provides a resolution where quantitative bone morphometry can be extracted in vivo on patients. It is known that HR-pQCT provides slight differences in morphometric indices compared to the current standard approach micro-computed tomography (micro-CT). The most obvious reason for this is the restriction of the radiation dose and with this a lower image resolution. With advances in micro-CT evaluation techniques such as patient-specific remodeling simulations or dynamic bone morphometry, a higher image resolution would potentially also allow the application of such novel evaluation techniques to clinical HR-pQCT measurements. Virtual supersampling as post-processing step was considered to increase the image resolution of HR-pQCT scans. The hypothesis was that this technique preserves the structural bone morphometry. Supersampling from 82 μm to virtual 41 μm by trilinear interpolation of the grayscale values of 42 human cadaveric forearms resulted in strong correlations of structural parameters (R²: 0.96–1.00). BV/TV was slightly overestimated (4.3%, R²: 1.00) compared to the HR-pQCT resolution. Tb.N was overestimated (7.47%; R²: 0.99) and Tb.Th was slightly underestimated (-4.20%; R²: 0.98). The technique was reproducible with PE_%CV between 1.96% (SMI) and 7.88% (Conn.D). In a clinical setting with 205 human forearms with or without fracture measured at 82 μm resolution HR-pQCT, the technique was sensitive to changes between groups in all parameters (p < 0.05) except trabecular thickness. In conclusion, we demonstrated that supersampling preserves the bone morphometry from HR-pQCT scans and is reproducible and sensitive to changes between groups. Supersampling can be used to investigate on the resolution dependency of HR-pQCT images and gain more insight into this imaging modality.

Cancellous bone and theropod dinosaur locomotion. Part II-a new approach to inferring posture and locomotor biomechanics in extinct tetrapod vertebrates

This paper is the second of a three-part series that investigates the architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and therefore has the potential to provide insight into locomotor biomechanics in extinct tetrapod vertebrates such as dinosaurs. Here in Part II, a new biomechanical modelling approach is outlined, one which mechanistically links cancellous bone architectural patterns with three-dimensional musculoskeletal and finite element modelling of the hindlimb. In particular, the architecture of cancellous bone is used to derive a single 'characteristic posture' for a given species-one in which bone continuum-level principal stresses best align with cancellous bone fabric-and thereby clarify hindlimb locomotor biomechanics. The quasi-static approach was validated for an extant theropod, the chicken, and is shown to provide a good estimate of limb posture at around mid-stance. It also provides reasonable predictions of bone loading mechanics, especially for the proximal hindlimb, and also provides a broadly accurate assessment of muscle recruitment insofar as limb stabilization is concerned. In addition to being useful for better understanding locomotor biomechanics in extant species, the approach hence provides a new avenue by which to analyse, test and refine palaeobiomechanical hypotheses, not just for extinct theropods, but potentially many other extinct tetrapod groups as well.

The relationship between oral tori and bite force

OBJECTIVE

The relationship between bite force and torus palatinus or mandibularis remains to be explained. The major aim of this study was to determine the correlation between bite force and oral tori.

METHODS

The bite force of 345 patients was measured with a bite force recorder; impressions of the shape and size of the oral tori were taken on plaster models prior to orthodontic treatments. Subsequently, the relationship between oral tori and bite force was analyzed.

RESULTS

The size, shape, and incidence of torus palatinus was not significantly correlated with bite force. However, the size of torus mandibularis increased significantly in proportion to the bite force (p = 0.020). The occurrence of different types of oral tori was not correlated with the bite force.

DISCUSSION

The size of torus mandibularis provides information about bite force and can thus be used to clinically assess occlusal stress.

Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur

A micro-finite element-based method to estimate the bone loading history based on bone architecture was recently presented in the literature. However, a thorough investigation of the parameter sensitivity and plausibility of this method to predict joint loads is still missing. The goals of this study were (1) to analyse the parameter sensitivity of the joint load predictions at one proximal femur and (2) to assess the plausibility of the results by comparing load predictions of ten proximal femora to in vivo hip joint forces measured with instrumented prostheses (available from www.orthoload.com ). Joint loads were predicted by optimally scaling the magnitude of four unit loads (inclined [Formula: see text] to [Formula: see text] with respect to the vertical axis) applied to micro-finite element models created from high-resolution computed tomography scans ([Formula: see text]m voxel size). Parameter sensitivity analysis was performed by varying a total of nine parameters and showed that predictions of the peak load directions (range 10[Formula: see text]-[Formula: see text]) are more robust than the predicted peak load magnitudes (range 2344.8-4689.5 N). Comparing the results of all ten femora with the in vivo loading data of ten subjects showed that peak loads are plausible both in terms of the load direction (in vivo: [Formula: see text], predicted: [Formula: see text]) and magnitude (in vivo: [Formula: see text], predicted: [Formula: see text]). Overall, this study suggests that micro-finite element-based joint load predictions are both plausible and robust in terms of the predicted peak load direction, but predicted load magnitudes should be interpreted with caution.

Computational modeling of long-term effects of prophylactic vertebroplasty on bone adaptation

Cement augmentation in vertebrae (vertebroplasty) is usually used to restore mechanical strength after spinal fracture but could also be used as a prophylactic treatment. So far, the mechanical competence has been determined immediately post-treatment, without considering long-term effects of bone adaptation. In this work, we investigated such long-term effects of vertebroplasty on the stiffness of the augmented bone by means of computational simulation of bone adaptation. Using micro-finite element analysis, we determined sites of increased mechanical stress (stress raisers) and stress shielding and, based on the simulations, regions with increased or decreased bone loss due to augmentation. Cement volumes connecting the end plates led to increased stress shielding and bone loss. The increased stiffness due to the augmentation, however, remained constant over the simulation time of 30 years. If the intervention was performed at an earlier time point, it did lead to more bone loss, but again, it did not affect long-term stability as this loss was compensated by bone gains in other areas. In particular, around the augmentation cement, bone structures were preserved, suggesting a long-term integration of the cement in the augmented bone. We conclude that, from a biomechanical perspective, the impact of vertebroplasty on the bone at the microstructural level is less detrimental than previously thought.

In vivo Visualisation and Quantification of Bone Resorption and Bone Formation from Time-Lapse Imaging

PURPOSE OF REVIEW

Mechanoregulation of bone cells was proposed over a century ago, but only now can we visualise and quantify bone resorption and bone formation and its mechanoregulation. In this review, we show how the newest advances in imaging and computational methods paved the way for this breakthrough.

RECENT FINDINGS

Non-invasive in vivo assessment of bone resorption and bone formation was demonstrated by time-lapse micro-computed tomography in animals, and by high-resolution peripheral quantitative computed tomography in humans. Coupled with micro-finite element analysis, the relationships between sites of bone resorption and bone formation and low and high tissue loading, respectively, were shown. Time-lapse in vivo imaging and computational methods enabled visualising and quantifying bone resorption and bone formation as well as its mechanoregulation. Future research includes visualising and quantifying mechanoregulation of bone resorption and bone formation from molecular to organ scales, and translating the findings into medicine using personalised bone health prognosis.

Biomaterials and Tissue Biomechanics: A Match Made in Heaven?

Biomaterials and tissue biomechanics have been traditionally separate areas of research with relatively little overlap in terms of methodological approaches. Recent advances in both fields on the one hand and developments in fabrication techniques and design approaches on the other have prepared the ground for joint research efforts by both communities. Additive manufacturing and rational design are examples of the revolutionary fabrication techniques and design methodologies that could facilitate more intimate collaboration between biomaterial scientists and biomechanists. This editorial article highlights the various ways in which the research on tissue biomechanics and biomaterials are related to each other and could benefit from each other’s results and methodologies.

A modelling approach demonstrating micromechanical changes in the tibial cemented interface due to in vivo service

Post-operative changes in trabecular bone morphology at the cement-bone interface can vary depending on time in service. This study aims to investigate how micromotion and bone strains change at the tibial bone-cement interface before and after cementation. This work discusses whether the morphology of the post-mortem interface can be explained by studying changes in these mechanical quantities. Three post-mortem cement-bone interface specimens showing varying levels of bone resorption (minimal, extensive and intermediate) were selected for this study Using image segmentation techniques, masks of the post-mortem bone were dilated to fill up the mould spaces in the cement to obtain the immediately post-operative situation. Finite element (FE) models of the post-mortem and post-operative situation were created from these segmentation masks. Subsequent removal of the cement layer resulted in the pre-operative situation. FE micromotion and bone strains were analyzed for the interdigitated trabecular bone. For all specimens micromotion increased from the post-operative to the post-mortem models (distally, in specimen 1: 0.1 to 0.5µm; specimen 2: 0.2 to 0.8µm; specimen 3: 0.27 to 1.62µm). Similarly bone strains were shown to increase from post-operative to post-mortem (distally, in specimen 1: -185 to -389µε; specimen 2: -170 to -824µε; specimen 3: -216 to -1024µε). Post-mortem interdigitated bone was found to be strain shielded in comparison with supporting bone indicating that failure of bone would occur distal to the interface. These results indicate that stress shielding of interdigitated trabeculae is a plausible explanation for resorption patterns observed in post-mortem specimens.

Load-adaptive bone remodeling simulations reveal osteoporotic microstructural and mechanical changes in whole human vertebrae

Osteoporosis is a major medical burden and its impact is expected to increase in our aging society. It is associated with low bone density and microstructural deterioration. Treatments are available, but the critical factor is to define individuals at risk from osteoporotic fractures. Computational simulations investigating not only changes in net bone tissue volume, but also changes in its microstructure where osteoporotic deterioration occur might help to better predict the risk of fractures. In this study, bone remodeling simulations with a mechanical feedback loop were used to predict microstructural changes due to osteoporosis and their impact on bone fragility from 50 to 80 years of age. Starting from homeostatic bone remodeling of a group of seven, mixed sex whole vertebrae, five mechanostat models mimicking different biological alterations associated with osteoporosis were developed, leading to imbalanced bone formation and resorption with a total net loss of bone tissue. A model with reduced bone formation rate and cell sensitivity led to the best match of morphometric indices compared to literature data and was chosen to predict postmenopausal osteoporotic bone loss in the whole group. Thirty years of osteoporotic bone loss were predicted with changes in morphometric indices in agreement with experimental measurements, and only showing major deviations in trabecular number and trabecular separation. In particular, although being optimized to match to the morphometric indices alone, the predicted bone loss revealed realistic changes on the organ level and on biomechanical competence. While the osteoporotic bone was able to maintain the mechanical stability to a great extent, higher fragility towards error loads was found for the osteoporotic bones.

Voxel size dependency, reproducibility and sensitivity of an in vivo bone loading estimation algorithm

A bone loading estimation algorithm was previously developed that provides in vivo loading conditions required for in vivo bone remodelling simulations. The algorithm derives a bone's loading history from its microstructure as assessed by high-resolution (HR) computed tomography (CT). This reverse engineering approach showed accurate and realistic results based on micro-CT and HR-peripheral quantitative CT images. However, its voxel size dependency, reproducibility and sensitivity still need to be investigated, which is the purpose of this study. Voxel size dependency was tested on cadaveric distal radii with micro-CT images scanned at 25 µm and downscaled to 50, 61, 75, 82, 100, 125 and 150 µm. Reproducibility was calculated with repeated in vitro as well as in vivo HR-pQCT measurements at 82 µm. Sensitivity was defined using HR-pQCT images from women with fracture versus non-fracture, and low versus high bone volume fraction, expecting similar and different loading histories, respectively. Our results indicate that the algorithm is voxel size independent within an average (maximum) error of 8.2% (32.9%) at 61 µm, but that the dependency increases considerably at voxel sizes bigger than 82 µm. In vitro and in vivo reproducibility are up to 4.5% and 10.2%, respectively, which is comparable to other in vitro studies and slightly higher than in other in vivo studies. Subjects with different bone volume fraction were clearly distinguished but not subjects with and without fracture. This is in agreement with bone adapting to customary loading but not to fall loads. We conclude that the in vivo bone loading estimation algorithm provides reproducible, sensitive and fairly voxel size independent results at up to 82 µm, but that smaller voxel sizes would be advantageous.

3D Local in vivo Environment (LivE) imaging for single cell protein analysis of bone tissue

The molecular processes behind pathological bone remodelling seen in diseases such as osteoporosis are unclear. However, a recently developed methodological platform known as Local in vivo Environment (LivE) imaging has been used to link cellular expression data to the local remodelling and mechanical environment in 2D sections of bone tissue. The method therefore can be used to give insight into which proteins are important for pathological bone remodelling. However, the cells within bone tissue exist as a 3D network. Therefore extension of LivE to accommodate 3D data may provide additional physiologically relevant information that is not possible to determine using 2D analysis alone. This will have implications for the further understanding of the cellular basis that underlies bone diseases such as osteoporosis. Here the LivE imaging technique is expanded to incorporate data from cells in a three dimensional manner via a serial sectioning technique. The methodological steps involved in the LivE imaging approach are defined and the optimisation steps performed are explained in detail.

Subject-specific musculoskeletal loading of the tibia: Computational load estimation

The systematic development of subject-specific computer models for the analysis of personalized treatments is currently a reality. In fact, many advances have recently been developed for creating virtual finite element-based models. These models accurately recreate subject-specific geometries and material properties from recent techniques based on quantitative image analysis. However, to determine the subject-specific forces, we need a full gait analysis, typically in combination with an inverse dynamics simulation study. In this work, we aim to determine the subject-specific forces from the computer tomography images used to evaluate bone density. In fact, we propose a methodology that combines these images with bone remodelling simulations and artificial neural networks. To test the capability of this novel technique, we quantify the personalized forces for five subject-specific tibias using our technique and a gait analysis. We compare both results, finding that similar vertical loads are estimated by both methods and that the dominant part of the load can be reliably computed. Therefore, we can conclude that the numerical-based technique proposed in this work has great potential for estimating the main forces that define the mechanical behaviour of subject-specific bone.

Estimation of Local Bone Loads for the Volume of Interest

Computational bone remodeling simulations have recently received significant attention with the aid of state-of-the-art high-resolution imaging modalities. They have been performed using localized finite element (FE) models rather than full FE models due to the excessive computational costs of full FE models. However, these localized bone remodeling simulations remain to be investigated in more depth. In particular, applying simplified loading conditions (e.g., uniform and unidirectional loads) to localized FE models have a severe limitation in a reliable subject-specific assessment. In order to effectively determine the physiological local bone loads for the volume of interest (VOI), this paper proposes a novel method of estimating the local loads when the global musculoskeletal loads are given. The proposed method is verified for the three VOI in a proximal femur in terms of force equilibrium, displacement field, and strain energy density (SED) distribution. The effect of the global load deviation on the local load estimation is also investigated by perturbing a hip joint contact force (HCF) in the femoral head. Deviation in force magnitude exhibits the greatest absolute changes in a SED distribution due to its own greatest deviation, whereas angular deviation perpendicular to a HCF provides the greatest relative change. With further in vivo force measurements and high-resolution clinical imaging modalities, the proposed method will contribute to the development of reliable patient-specific localized FE models, which can provide enhanced computational efficiency for iterative computing processes such as bone remodeling simulations.

Trabecular bone loss at a distant skeletal site following noninvasive knee injury in mice

Traumatic injuries can have systemic consequences, as the early inflammatory response after trauma can lead to tissue destruction at sites not affected by the initial injury. This systemic catabolism may occur in the skeleton following traumatic injuries such as anterior cruciate ligament (ACL) rupture. However, bone loss following injury at distant,unrelated skeletal sites has not yet been established. In the current study, we utilized a mouse knee injury model to determine whether acute knee injury causes a mechanically significant trabecular bone loss at a distant, unrelated skeletal site (L5 vertebral body).Knee injury was noninvasively induced using either high-speed (HS; 500 mm/s) or lowspeed(LS; 1 mm/s) tibial compression overload. HS injury creates an ACL rupture by midsubstance tear, while LS injury creates an ACL rupture with an associated avulsion bone fracture. At 10 days post-injury, vertebral trabecular bone structure was quantified using high-resolution microcomputed tomography (lCT), and differences in mechanical properties were determined using finite element modeling (FEM) and compressive mechanical testing. We hypothesized that knee injury would initiate a loss of trabecular bone structure and strength at the L5 vertebral body. Consistent with our hypothesis, we found significant decreases in trabecular bone volume fraction (BV/TV) and trabecular number at the L5 vertebral body in LS injured mice compared to sham (8.8% and 5.0%, respectively), while HS injured mice exhibited a similar, but lower magnitude response (5.1% and 2.5%, respectively). Contrary to our hypothesis, this decrease intrabecular bone structure did not translate to a significant deficit in compressive stiffness or ultimate load of the full trabecular body assessed by mechanical testing or FEM. However,we were able to detect significant decreases in compressive stiffness in both HS and LS injured specimens when FE models were loaded directly through the trabecular bone region (9.9% and 8.1%, and 3, respectively). This finding may be particularly important for osteoporotic fracture risk, as damage within vertebral bodies has been shown to initiate within the trabecular bone compartment. Altogether, these data point to a systemic trabecular bone loss as a consequence of fracture or traumatic musculoskeletal injury, which may be an underlying mechanism contributing to increased risk of refracture following an initial injury. This finding may have consequences for treatment of acute musculoskeletal injuries and the prevention of future bone fragility.

A potential mechanism for allometric trabecular bone scaling in terrestrial mammals

Trabecular bone microstructural parameters, including trabecular thickness, spacing, and number, have been reported to scale with animal size with negative allometry, whereas bone volume fraction is animal size-invariant in terrestrial mammals. As for the majority of scaling patterns described in animals, its underlying mechanism is unknown. However, it has also been found that osteocyte density is inversely related to animal size, possibly adapted to metabolic rate, which shows a negative relationship as well. In addition, the signalling reach of osteocytes is limited by the extent of the lacuno-canalicular network, depending on trabecular dimensions and thus also on animal size. Here we propose animal size-dependent variations in osteocyte density and their signalling influence distance as a potential mechanism for negative allometric trabecular bone scaling in terrestrial mammals. Using an established and tested computational model of bone modelling and remodelling, we run simulations with different osteocyte densities and influence distances mimicking six terrestrial mammals covering a large range of body masses. Simulated trabecular structures revealed negative allometric scaling for trabecular thickness, spacing, and number, constant bone volume fraction, and bone turnover rates inversely related to animal size. These results are in agreement with previous observations supporting our proposal of osteocyte density and influence distance variation as a potential mechanism for negative allometric trabecular bone scaling in terrestrial mammals. The inverse relationship between bone turnover rates and animal size further indicates that trabecular bone scaling may be linked to metabolic rather than mechanical adaptations.

Patient-specific bone modeling and analysis: the role of integration and automation in clinical adoption

Patient-specific analysis of bones is considered an important tool for diagnosis and treatment of skeletal diseases and for clinical research aimed at understanding the etiology of skeletal diseases and the effects of different types of treatment on their progress. In this article, we discuss how integration of several important components enables accurate and cost-effective patient-specific bone analysis, focusing primarily on patient-specific finite element (FE) modeling of bones. First, the different components are briefly reviewed. Then, two important aspects of patient-specific FE modeling, namely integration of modeling components and automation of modeling approaches, are discussed. We conclude with a section on validation of patient-specific modeling results, possible applications of patient-specific modeling procedures, current limitations of the modeling approaches, and possible areas for future research.

Bone adaptation to cyclic loading in murine caudal vertebrae is maintained with age and directly correlated to the local micromechanical environment

The ability of the skeleton to adapt to mechanical stimuli (mechanosensitivity) has most often been investigated at the whole-bone level, but less is known about the local mechanoregulation of bone remodeling at the bone surface, especially in context of the aging skeleton. The aim of this study was to determine the local and global mechanosensitivity of the sixth caudal vertebra during cyclic loading (8 N, three times per week, for six weeks) in mice aged 15, 52, and 82 weeks at the start of loading. Bone adaptation was monitored with in vivo micro-computed tomography. Strain energy density (SED), assumed as the mechanical stimulus for bone adaptation, was determined with micro-finite element models. Mechanical loading had a beneficial effect on the bone microstructure and bone stiffness in all age groups. Mineralizing surface was on average 13% greater (p<0.05) in loaded than control groups in 15- and 82-week-old mice, but not for 52-week-old mice. SED at the start of loading correlated to the change in bone volume fraction in the following 6 weeks for loaded groups (r(2)=0.69-0.85) but not control groups. At the local level, SED was 14-20% greater (p<0.01) at sites of bone formation, and 15-20% lower (p<0.01) at sites of bone resorption compared to quiescent bone surfaces for all age groups, indicating SED was a stimulus for bone adaptation. Taken together, these results support that mechanosensitivity is maintained with age in caudal vertebrae of mice at a local and global level. Since age-related bone loss was not observed in caudal vertebrae, results from the current study might not be translatable to aged humans.

Assessing bone's adaptive capacity around dental implants: a literature review

BACKGROUND

Increased stress (force) on prostheses induces strain (deformation) in the peri-implant bone. Elevated stress and strain could result in the failure of implants that support prostheses. However, the survival rate of implants supporting prostheses under increased stress is high. Either the bone is stronger than expected or it adapts to increased stress. Concepts regarding bone's adaptive capacity continue to evolve and are the focus of this literature review.

TYPES OF STUDIES REVIEWED

The authors searched the literature to find studies that addressed the bone's capacity to adjust to increased stress and strain. They assessed experimental and clinical trials in which investigators monitored healing after placement of dental implants.

RESULTS

The data indicate that forces greater than the bone's adaptive ability can induce loss of osseointegration, as well as osseous resorption. In contrast, it is possible that increased stress on prostheses initiates a reparative process, thereby facilitating retention of implants experiencing increased stress. Numerous lines of evidence support the concept that bone can modify itself to withstand increased mechanical forces.

PRACTICAL IMPLICATIONS

The authors provide an explanation for the high success rate of prostheses and implants in bone that are exposed to increased stress and strain.

Statistical shape and appearance models of bones

When applied to bones, statistical shape models (SSM) and statistical appearance models (SAM) respectively describe the mean shape and mean density distribution of bones within a certain population as well as the main modes of variations of shape and density distribution from their mean values. The availability of this quantitative information regarding the detailed anatomy of bones provides new opportunities for diagnosis, evaluation, and treatment of skeletal diseases. The potential of SSM and SAM has been recently recognized within the bone research community. For example, these models have been applied for studying the effects of bone shape on the etiology of osteoarthritis, improving the accuracy of clinical osteoporotic fracture prediction techniques, design of orthopedic implants, and surgery planning. This paper reviews the main concepts, methods, and applications of SSM and SAM as applied to bone.

Trabecular evidence for a human-like gait in Australopithecus africanus

Although the earliest known hominins were apparently upright bipeds, there has been mixed evidence whether particular species of hominins including those in the genus Australopithecus walked with relatively extended hips, knees and ankles like modern humans, or with more flexed lower limb joints like apes when bipedal. Here we demonstrate in chimpanzees and humans a highly predictable and sensitive relationship between the orientation of the ankle joint during loading and the principal orientation of trabecular bone struts in the distal tibia that function to withstand compressive forces within the joint. Analyses of the orientation of these struts using microCT scans in a sample of fossil tibiae from the site of Sterkfontein, of which two are assigned to Australopithecus africanus, indicate that these hominins primarily loaded their ankles in a relatively extended posture like modern humans and unlike chimpanzees. In other respects, however, trabecular properties in Au africanus are distinctive, with values that mostly fall between those of chimpanzees and humans. These results indicate that Au. africanus, like Homo, walked with an efficient, extended lower limb.