High-throughput quantification of the mechanical competence of murine femora--a highly automated approach for large-scale genetic studies.

Animal models are widely used to gain insight into the role of genetics on bone structure and function. One of the main strategies to map the genes regulating specific traits is called quantitative trait loci (QTL) analysis, which generally requires a very large number of animals (often more than 1000) to reach statistical significance. QTL analysis for mechanical traits has been mainly based on experimental mechanical testing, which, in view of the large number of animals, is time consuming. Hence, the goal of the present work was to introduce an automated method for large-scale high-throughput quantification of the mechanical properties of murine femora. Specifically, our aims were, first, to develop and validate an automated method to quantify murine femoral bone stiffness. Second, to test its high-throughput capabilities on murine femora from a large genetic study, more specifically, femora from two growth hormone (GH) deficient inbred strains of mice (B6-lit/lit and C3.B6-lit/lit) and their first (F1) and second (F2) filial offsprings. Automated routines were developed to convert micro-computed tomography (micro-CT) images of femora into micro-finite element (micro-FE) models. The method was experimentally validated on femora from C57BL/6J and C3H/HeJ mice: for both inbred strains the micro-FE models closely matched the experimentally measured bone stiffness when using a single tissue modulus of 13.06 GPa. The mechanical analysis of the entire dataset (n=1990) took approximately 44 CPU hours on a supercomputer. In conclusion, our approach, in combination with QTL analysis could help to locate genes directly involved in controlling bone mechanical competence.

Mechanisms of reduced implant stability in osteoporotic bone.

The determining factors for the fixation of uncemented screws in bone are the bone-implant interface and the peri-implant bone. The goal of this work was to explore the role of the peri-implant bone architecture on the mechanics of the bone-implant system. In particular, the specific aims of the study were to investigate: (i) the impact of the different architectural parameters, (ii) the effects of disorder, and (iii) the deformations in the peri-implant region. A three-dimensional beam lattice model to describe trabecular bone was developed. Various microstructural features of the lattice were varied in a systematic way. Implant pull-out tests were simulated, and the stiffness and strength of the bone-implant system were computed. The results indicated that the strongest decrease in pull-out strength was obtained by trabecular thinning, whereas pull-out stiffness was mostly affected by trabecular removal. These findings could be explained by investigating the peri-implant deformation field. For small implant displacements, a large amount of trabeculae in the peri-implant region were involved in the load transfer from implant to bone. Therefore, trabecular removal in this region had a strong negative effect on pull-out stiffness. Conversely, at higher displacements, deformations mainly localized in the trabeculae in contact with the implant; hence, thinning those trabeculae produced the strongest decrease in the strength of the system. Although idealized, the current approach is helpful for a mechanical understanding of the role played by peri-implant bone.

The bone mineralization density distribution as a fingerprint of the mineralization process.

The inhomogeneous mineral content and its topographical distribution on a microscopic scale are major determinants of the mechanical quality of trabecular bone. The kinetics of bone tissue deposition and resorption together with the kinetics of the mineralization process determine the distribution of mineral in the tissue. The heterogeneity of the mineral content is described by the well-established bone mineralization density distribution (BMDD), which is experimentally accessible, e.g., using quantitative electron backscattering imaging (qBEI). In the present work, we demonstrate that the shape of the BMDD histogram of trabecular bone reflects directly the mineralization kinetics. Based on the experimental BMDD data of trabecular bone from healthy human adults and using a mathematical model for the remodeling and the mineralization process, the following main results were obtained. The peaked BMDD reflects necessarily a two-phase mineralization process with a fast primary phase and a slow secondary phase where the corresponding time constants differ three orders of magnitude. The obtained mineralization law, which describes the increase in the mineral content in a bone packet as a function of time, provides information not only about the initial mineralization surge, but also about the slow increase afterwards on the time scale of years. In addition to the mineralization kinetics the turnover rate of the remodeling process has a strong influence on the peak position and the shape of the BMDD. The described theoretical framework opens new possibilities for an analysis of experimentally measured BMDDs with respect to changes caused by diseases or treatments. It allows addressing whether changes in the BMDD have to be attributed to a variation in the turnover rate which consequently affects the density distribution or to a primary disorder in the mineralization process most likely reflecting alterations of the organic matrix. This is of important clinical interest because it helps to find therapeutic approaches directly targeting the primary etiological defects to correct the patients' BMDD towards normal BMDD.