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

Local anisotropy in mineralized fibrocartilage and subchondral bone beneath the tendon-bone interface

The enthesis allows the insertion of tendon into bone thanks to several remarkable strategies. This complex and clinically relevant location often features a thin layer of fibrocartilage sandwiched between tendon and bone to cope with a highly heterogeneous mechanical environment. The main purpose of this study was to investigate whether mineralized fibrocartilage and bone close to the enthesis show distinctive three-dimensional microstructural features, possibly to enable load transfer from tendon to bone. As a model, the Achilles tendon-calcaneus bone system of adult rats was investigated with histology, backscattered electron imaging and micro-computed tomography. The microstructural porosity of bone and mineralized fibrocartilage in different locations including enthesis fibrocartilage, periosteal fibrocartilage and bone away from the enthesis was characterized. We showed that calcaneus bone presents a dedicated protrusion of low porosity where the tendon inserts. A spatially resolved analysis of the trabecular network suggests that such protrusion may promote force flow from the tendon to the plantar ligament, while partially relieving the trabecular bone from such a task. Focusing on the tuberosity, highly specific microstructural aspects were highlighted. Firstly, the interface between mineralized and unmineralized fibrocartilage showed the highest roughness at the tuberosity, possibly to increase failure resistance of a region carrying large stresses. Secondly, fibrochondrocyte lacunae inside mineralized fibrocartilage, in analogy with osteocyte lacunae in bone, had a predominant alignment at the enthesis and a rather random organization away from it. Finally, the network of subchondral channels inside the tuberosity was highly anisotropic when compared to contiguous regions. This dual anisotropy of subchondral channels and cell lacunae at the insertion may reflect the alignment of the underlying collagen network. Our findings suggest that the microstructure of fibrocartilage may be linked with the loading environment. Future studies should characterize those microstructural aspects in aged and or diseased conditions to elucidate the poorly understood role of bone and fibrocartilage in enthesis-related pathologies.

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.

Predicting bone strength from CT data: Clinical applications

In this review we summarise over 15 years of research and development around the prediction of whole bones strength from Computed Tomography data, with particular reference to the prediction of the risk of hip fracture in osteoporotic patients. We briefly discuss the theoretical background, and then provide a summary of the laboratory and clinical validation of these modelling technologies. We then discuss the three current clinical applications: in clinical research, in clinical trials, and in clinical practice. On average the strength predicted with finite element models (QCT-FE) based on computed tomography is 7% more accurate that that predicted with areal bone mineral density from Dual X-ray Absorptiometry (DXA-aBMD), the current standard of care, both in term of laboratory validation on cadaver bones and in terms of stratification accuracy on clinical cohorts of fractured and non-fractured women. This improved accuracy makes QCT-FE superior to DXA-aBMD in clinical research and in clinical trials, where the its use can cut in half the number of patients to be enrolled to get the same statistical power. For routine clinical use to decide who to treat with antiresorptive drugs, QCT-FE is more accurate but less cost-effective than DXA-aBMD, at least when the decision is on first line treatment like bisphosphonates. But the ability to predict skeletal strength from medical imaging is now opening a number of other applications, for example in paediatrics and oncology.

Weak Epistasis Generally Stabilizes Phenotypes in a Mouse Intercross

The extent and strength of epistasis is commonly unresolved in genetic studies, and observed epistasis is often difficult to interpret in terms of biological consequences or overall genetic architecture. We investigated the prevalence and consequences of epistasis by analyzing four body composition phenotypes—body weight, body fat percentage, femoral density, and femoral circumference—in a large F2 intercross of B6-lit/lit and C3.B6-lit/lit mice. We used Combined Analysis of Pleiotropy and Epistasis (CAPE) to examine interactions for the four phenotypes simultaneously, which revealed an extensive directed network of genetic loci interacting with each other, circulating IGF1, and sex to influence these phenotypes. The majority of epistatic interactions had small effects relative to additive effects of individual loci, and tended to stabilize phenotypes towards the mean of the population rather than extremes. Interactive effects of two alleles inherited from one parental strain commonly resulted in phenotypes closer to the population mean than the additive effects from the two loci, and often much closer to the mean than either single-locus model. Alternatively, combinations of alleles inherited from different parent strains contribute to more extreme phenotypes not observed in either parental strain. This class of phenotype-stabilizing interactions has effects that are close to additive and are thus difficult to detect except in very large intercrosses. Nevertheless, we found these interactions to be useful in generating hypotheses for functional relationships between genetic loci. Our findings suggest that while epistasis is often weak and unlikely to account for a large proportion of heritable variance, even small-effect genetic interactions can facilitate hypotheses of underlying biology in well-powered studies.

The role of statistical epistasis in the genetic architecture of complex traits has been of great interest to the genetics community since Fisher introduced the concept in 1918. However, assessing epistasis in human and model organism populations has been impeded by limited statistical power. To mitigate this limitation, we analyzed bone and body composition traits in an unusually large mouse intercross population of over 2000 mice, paired with a recently-developed computational approach that leverages information to detect interactions across multiple phenotypes. We discovered a large network of highly significant genetic interactions between variants that influence complex body composition traits. Although epistasis was abundant, the interaction network was dominated by epistasis that stabilizes phenotypes by reducing phenotypic deviation from the parent strains. Nevertheless, the observed network provides an overview of genetic architecture and specific hypotheses of how QTL combine to affect phenotypes. These findings suggest that epistatic effects are generally of lesser magnitude than main QTL effects, and therefore are unlikely to account for major components of variance, but also reinforce genetic interaction analysis as a potent tool for dissecting the biology of complex traits.

High-throughput phenotyping and genetic linkage of cortical bone microstructure in the mouse

Background

Understanding cellular structure and organization, which plays an important role in biological systems ranging from mechanosensation to neural organization, is a complicated multifactorial problem depending on genetics, environmental factors, and stochastic processes. Isolating these factors necessitates the measurement and sensitive quantification of many samples in a reliable, high-throughput, unbiased manner. In this manuscript we present a pipelined approach using a fully automated framework based on Synchrotron-based X-ray Tomographic Microscopy (SRXTM) for performing a full 3D characterization of millions of substructures.

Results

We demonstrate the framework on a genetic study on the femur bones of in-bred mice. We measured 1300 femurs from a F2 cross experiment in mice without the growth hormone (which can confound many of the smaller structural differences between strains) and characterized more than 50 million osteocyte lacunae (cell-sized hollows in the bone). The results were then correlated with genetic markers in a process called quantitative trait localization (QTL). Our findings provide a mapping between regions of the genome (all 19 autosomes) and observable phenotypes which could explain between 8–40 % of the variance using between 2–10 loci for each trait. This map shows 4 areas of overlap with previous studies looking at bone strength and 3 areas not previously associated with bone.

Conclusions

The mapping of microstructural phenotypes provides a starting point for both structure-function and genetic studies on murine bone structure and the specific loci can be investigated in more detail to identify single gene candidates which can then be translated to human investigations. The flexible infrastructure offers a full spectrum of shape, distribution, and connectivity metrics for cellular networks and can be adapted to a wide variety of materials ranging from plant roots to lung tissue in studies requiring high sample counts and sensitive metrics such as the drug-gene interactions and high-throughput screening.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1617-y) contains supplementary material, which is available to authorized users.

Characterization of complex, co-adapted skeletal biomechanics phenotypes: a needed paradigm shift in the genetics of bone structure and function

The genetic architecture of skeletal biomechanical performance has tremendous potential to advance our knowledge of the biological mechanisms that drive variation in skeletal fragility and osteoporosis risk. Research using traditional approaches that focus on specific gene pathways is increasing our understanding of how and to what degree those pathways may affect population-level variation in fracture susceptibility, and shows that known pathways may affect bone fragility through unsuspected mechanisms. Non-traditional approaches that incorporate a new appreciation for the degree to which bone traits co-adapt to functional loading environments, using a wide variety of redundant compensatory mechanisms to meet both physiological and mechanical demands, represent a radical departure from the dominant reductionist paradigm and have the potential to rapidly advance our understanding of bone fragility and identification of new targets for therapeutic intervention.

Modeling microdamage behavior of cortical bone

Bone is a complex material which exhibits several hierarchical levels of structural organization. At the submicron-scale, the local tissue porosity gives rise to discontinuities in the bone matrix which have been shown to influence damage behavior. Computational tools to model the damage behavior of bone at different length scales are mostly based on finite element (FE) analysis, with a range of algorithms developed for this purpose. Although the local mechanical behavior of bone tissue is influenced by microstructural features such as bone canals and osteocyte lacunae, they are often not considered in FE damage models due to the high computational cost required to simulate across several length scales, i.e., from the loads applied at the organ level down to the stresses and strains around bone canals and osteocyte lacunae. Hence, the aim of the current study was twofold: First, a multilevel FE framework was developed to compute, starting from the loads applied at the whole bone scale, the local mechanical forces acting at the micrometer and submicrometer level. Second, three simple microdamage simulation procedures based on element removal were developed and applied to bone samples at the submicrometer-scale, where cortical microporosity is included. The present microdamage algorithm produced a qualitatively analogous behavior to previous experimental tests based on stepwise mechanical compression combined with in situ synchrotron radiation computed tomography. Our results demonstrate the feasibility of simulating microdamage at a physiologically relevant scale using an image-based meshing technique and multilevel FE analysis; this allows relating microdamage behavior to intracortical bone microstructure.

The application of digital volume correlation (DVC) to study the microstructural behaviour of trabecular bone during compression

Digital Volume Correlation (DVC) has been emerged recently as an innovative approach to full volume (i.e. internal) displacement and strain field measurement in materials and structures, particularly in conjunction with high resolution X-ray computed tomography (CT). As a relatively novel technique certain aspects of precision, accuracy and the breadth of application are yet to be fully established. This study has applied DVC to volume images of porcine trabecular bone assessing the effect of noise and sub-volume size on strain measurement. Strain resolutions ranging between 70 and 800με were obtained for the optimum sub-volume size of 64 voxels with a 50% overlap for metrological studies conducted. These values allowed the mechanical behaviour of porcine trabecular bone during compression to be investigated. During compression a crushed layer formed adjacent to the boundary plate which increased in thickness as the specimen was further deformed. The structure of the crushed layer was altered to such an extent that it confounded the correlation method. While investigating this factor, it was found that for reliable strain calculations a correlation coefficient of 0.90 or above was required between the sub-volumes in the reference and the deformed volumes. Good agreements between the results and published bone strain failures were obtained. Using the full field strain measurements, Poisson's ratio was identified for each compression step using a dedicated inverse method called the virtual fields method (VFM). It was found that for a given region outside of the crushed zone the Poisson ratio decreased from 0.32 to 0.21 between the first and the final compression steps, which was hypothesised to be due to the bone geometry and its resulting deformation behaviour. This study demonstrates that volumetric strain measurement can be obtained successfully using DVC, making it a useful tool for quantitatively investigating the micro-mechanical behaviour of macroscale bone specimens.