Quantitative trait loci that determine mouse tibial nanoindentation properties in an F2 population derived from C57BL/6J x C3H/HeJ

Regional Nanoindentation Properties in Different Locations on the Mouse Tibia From C57BL/6 and Balb/C Female Mice

The local spatial heterogeneity of the material properties of the cortical and trabecular bone extracted from the mouse tibia is not well-known. Nevertheless, its characterization is fundamental to be able to study comprehensively the effect of interventions and to generate computational models to predict the bone strength preclinically. The goal of this study was to evaluate the nanoindentation properties of bone tissue extracted from two different mouse strains across the tibia length and in different sectors. Left tibiae were collected from four female mice, two C57BL/6, and two Balb/C mice. Nanoindentations with maximum 6 mN load were performed on different microstructures, regions along the axis of the tibiae, and sectors (379 in total). Reduced modulus (Er) and hardness (H) were computed for each indentation. Trabecular bone of Balb/C mice was 21% stiffer than that of C57BL/6 mice (20.8 ± 4.1 GPa vs. 16.5 ± 7.1 GPa). Moreover, the proximal regions of the bones were 13-36% less stiff than the mid-shaft and distal regions of the same bones. No significant differences were found for the different sectors for E _r and H for Balb/C mice. The bone in the medial sector was found to be 8-14% harder and stiffer than the bone in the anterior or posterior sectors for C57BL/6 mice. In conclusion, this study showed that the nanoindentation properties of the mouse tibia are heterogeneous across the tibia length and the trabecular bone properties are different between Balb/C and C57BL/6 mice. These results will help the research community to identify regions where to characterize the mechanical properties of the bone during preclinical optimisation of treatments for skeletal diseases.

The structure, composition and mechanical properties of the skeleton of the naked mole-rat (Heterocephalus glaber)

The naked mole-rat (NMR) is a small rodent with a remarkable array of properties, such as unique physiology, extremely long life-span and unusual social life. However, very little is known regarding its skeleton. The aim of this study was to describe the structure, composition and mechanical properties in an ontogenetic series of naked mole-rat bones. Since common small rodents like mice and rats have an unusual structure of cortical bone, which includes a central region of non-lamellar (disordered) bone, mineralized cartilaginous islands and total lack of remodeling, this study could also determine if these are features of all small rodents. Sixty-one NMRs were included in the study and were divided into the following four age groups: 0-0.5 years old (n = 17), 0.5-3 years old (n = 25), 3-10 years old (n = 13), and >10 years (n = 6). Femora, vertebrae and mandibulae were examined using micro-CT, light microscopy, polarized light microscopy and scanning electron microscopy, thermogravimetric analysis was used to determine their dry ash content and their derived elastic modulus and hardness were determined using micro-indentation. Our findings show that NMR bones are similar in composition and mechanical properties to those of other small rodents. However, in contrast to other small rodents, the cortical bone of NMRs is entirely circumferential-lamellar and lacks mineralized cartilage islands. Furthermore, despite their long life-span, their bones did not show evidence of remodeling at any of the age groups, thus proving that lack of cortical remodeling in small rodents is not caused by their short life-span, but characteristic of this order of mammals.

Lactation-Induced Changes in the Volume of Osteocyte Lacunar-Canalicular Space Alter Mechanical Properties in Cortical Bone Tissue

Osteocytes can remove and remodel small amounts of their surrounding bone matrix through osteocytic osteolysis, which results in increased volume occupied by lacunar and canalicular space (LCS). It is well established that cortical bone stiffness and strength are strongly and inversely correlated with vascular porosity, but whether changes in LCS volume caused by osteocytic osteolysis are large enough to affect bone mechanical properties is not known. In the current studies we tested the hypotheses that (1) lactation and postlactation recovery in mice alter the elastic modulus of bone tissue, and (2) such local changes in mechanical properties are related predominantly to alterations in lacunar and canalicular volume rather than bone matrix composition. Mechanical testing was performed using microindentation to measure modulus in regions containing solely osteocytes and no vascular porosity. Lactation caused a significant (∼13%) reduction in bone tissue-level elastic modulus (p < 0.001). After 1 week postweaning (recovery), bone modulus levels returned to control levels and did not change further after 4 weeks of recovery. LCS porosity tracked inversely with changes in cortical bone modulus. Lacunar and canalicular void space increased 7% and 15% with lactation, respectively (p < 0.05), then returned to control levels at 1 week after weaning. Neither bone mineralization (assessed by high-resolution backscattered scanning electron microscopy) nor mineral/matrix ratio or crystallinity (assessed by Raman microspectroscopy) changed with lactation. Thus, changes in bone mechanical properties induced by lactation and recovery appear to depend predominantly on changes in osteocyte LCS dimensions. Moreover, this study demonstrates that tissue-level cortical bone mechanical properties are rapidly and reversibly modulated by osteocytes in response to physiological challenge. These data point to a hitherto unappreciated role for osteocytes in modulating and maintaining local bone mechanical properties. © 2016 American Society for Bone and Mineral Research.

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.

Nanomechanics insights into the performance of healthy and osteoporotic bones

In situ nanoscopic observations of healthy and osteoporotic bone nanopillars under compression were performed. The structural-mechanical property relationship at the atomic scale suggests that cortical bone performance is correlated to the feature, arrangement, movement, distortion, and fracture of hydroxyapatite nanocrystals. Healthy bone comprising tightly bound mineral nanocrystals shows high structural stability with nanoscopic lattice distortions and dislocation activities. On the other hand, osteoporotic bone exhibits brittleness owing to the movements of dispersed minerals in and intergranular fracture along a weak organic matrix.

Genomic dissection and prioritizing of candidate genes of QTL for regulating spontaneous arthritis on chromosome 1 in mice deficient for interleukin-1 receptor antagonist

Rheumatoid arthritis is a heterogeneous disease with clinical and biological polymorphisms. IL-1RN is a protein that binds to interleukin-1 (IL-1) receptors and inhibits the binding of IL-1-alpha and IL-1-beta. IL-1RN levels are elevated in the blood of patients with a variety of infectious, immune, and traumatic conditions. Balb/c mice deficient in IL-1ra (mouse gene of IL-1RN) develop spontaneous autoimmune arthritis while DBA/1 mice deficient in IL-1ra do not. Previously, we identified a major QTL that regulates the susceptibility to arthritis in Balb/c mice with IL-1ra deficiency. In this study, we found that the QTL may contain two peaks that are regulated by two sets of candidate genes. By haplotype analysis, the total genomic regions of candidate genes were reduced from about 19 Mbp to approximately 9 Mbp. The total number of candidate genes was reduced from 208 to 21.

Assessment of lamellar level properties in mouse bone utilizing a novel spherical nanoindentation data analysis method

In this work, we demonstrate the viability of using our recently developed data analysis procedures for spherical nanoindentation in conjunction with Raman spectroscopy for studying lamellar-level correlations between the local composition and local mechanical properties in mouse bone. Our methodologies allow us to convert the raw load-displacement datasets to much more meaningful indentation stress-strain curves that accurately capture the loading and unloading elastic moduli, the indentation yield points, as well as the post-yield characteristics in the tested samples. Using samples of two different inbred mouse strains, A/J and C57BL/6J (B6), we successfully demonstrate the correlations between the mechanical information obtained from spherical nanoindentation measurements to the local composition measured using Raman spectroscopy. In particular, we observe that a higher mineral-to-matrix ratio correlated well with a higher local modulus and yield strength in all samples. Thus, new bone regions exhibited lower moduli and yield strengths compared to more mature bone. The B6 mice were also found to exhibit lower modulus and yield strength values compared to the more mineralized A/J strain.

FINDING OF MICRODAMAGE MORPHOLOGY DIFFERENCES IN MOUSE FEMORAL BONES WITH DISTINCT MINERALIZATION LEVELS

Microdamage progression in bone is dependent on the ultrastructure of the tissue. Thus, any pathological changes in bone ultrastructure may be reflected in the pattern and capacity of microdamage accumulation. A previous numerical study of microdamage progression in bone using a probabilistic failure analysis approach predicts that the microdamage morphology (either linear microcracks or diffuse damage) is very sensitive to the level of mineralization in bone, which is also implicated in some experimental observations. To examine the prediction, femurs from two strains of mice (C57BL/6J, N = 10 and C3H/HeJ, N = 11) that have distinct mineralization levels were fatigued under four-point bending to create damage in the bone tissues. After testing, the microdamage morphology of the bone samples was examined using bulk-staining technique with basic fuchsin. The results demonstrate that more linear microcracks are observed in femurs of C3H/HeJ (higher mineralization), whereas more diffuse-like damage is found in C57BL/6J femurs (less mineralized). Compared with linear microcracks, the formation of diffuse damage tends to dissipate more energy and help bone to avoid catastrophic failures. Therefore, results from this study may help explain why highly mineralized bone tends to be more brittle. Observations from this study are consistent with the numerical prediction from the previous study, suggesting that mineralization has a significant effect on the microdamage morphology of bone.

Identifying a major locus that regulates spontaneous arthritis in IL-1ra-deficient mice and analysis of potential candidates

To identify genetic loci that regulate spontaneous arthritis in interleukin-1 receptor antagonist (IL-1ra)-deficient mice, an F2 population was created from a cross between Balb/c IL-1ra-deficient mice and DBA/1 IL-1ra-deficient mice. Spontaneous arthritis in the F2 population was examined and recorded. Genotypes of those F2 mice were determined using microsatellite markers. Quantitative trail locus (QTL) analysis was conducted with R/qtlbim. Functions of genes within QTL chromosomal regions were evaluated using a bioinformatics tool, PGMapper, and microarray analysis. Potential candidate genes were further evaluated using GeneNetwork. A total of 137 microsatellite markers with an average of 12 cM spacing along the whole genome were used for determining the correlation of arthritis phenotypes with genotypes of 191 F2 progenies. By whole-genome mapping, we obtained QTLs on chromosomes 1 and 6 that were above the significance threshold for strong Bayesian evidence. The QTL on chromosome 1 had a peak near D1Mit55 and D1Mit425 at 82·6 cM. It may account for as much as 12% of the phenotypic variation in susceptibility to spontaneous arthritis. The QTL region contained 208 known transcripts. According to their functions, Mr1, Pla2g4a and Fasl are outstanding candidate genes. From microarray analysis, 11 genes were selected as favourable candidates based on their function and expression profiles. Three of those 11 genes, Prg4, Ptgs2 and Mr1, correlated with the IL-1ra pathway. Those genes were considered to be the best candidates.

Genetically determined phenotype covariation networks control bone strength

To identify genes affecting bone strength, we studied how genetic variants regulate components of a phenotypic covariation network that was previously shown to accurately characterize the compensatory trait interactions involved in functional adaptation during growth. Quantitative trait loci (QTLs) regulating femoral robustness, morphologic compensation, and mineralization (tissue quality) were mapped at three ages during growth using AXB/BXA Recombinant Inbred (RI) mouse strains and adult B6-i(A) Chromosome Substitution Strains (CSS). QTLs for robustness were identified on chromosomes 8, 12, 18, and 19 and confirmed at all three ages, indicating that genetic variants established robustness postnatally without further modification. A QTL for morphologic compensation, which was measured as the relationship between cortical area and body weight, was identified on chromosome 8. This QTL limited the amount of bone formed during growth and thus acted as a setpoint for diaphyseal bone mass. Additional QTLs were identified from the CSS analysis. QTLs for robustness and morphologic compensation regulated bone structure independently (ie, in a nonpleiotropic manner), indicating that each trait may be targeted separately to individualize treatments aiming to improve strength. Multiple regression analyses showed that variation in morphologic compensation and tissue quality, not bone size, determined femoral strength relative to body weight. Thus an individual inheriting slender bones will not necessarily inherit weak bones unless the individual also inherits a gene that impairs compensation. This systems genetic analysis showed that genetically determined phenotype covariation networks control bone strength, suggesting that incorporating functional adaptation into genetic analyses will advance our understanding of the genetic basis of bone strength.

Quantitative trait loci for tibial bone strength in C57BL/6J and C3H/HeJ inbred strains of mice

Three-point bending technology has been widely used in the measurement of bone strength. Quantitative trait loci (QTLs) for bone strength have been identified using mouse femurs. In this study, we investigate the use of mouse tibiae in identification of QTLs that regulate bone strength. Mouse tibiae were from a F(2) population derived from C57BL/6J (B6) and C3H/HeJ (C3H). Three-point bending was measured using ISO 4049, with the support width adjustable to accommodate specimen sizes outside the scope of ISO 4049. The strain rate is selectable from 0.05 to 500 mm per min. All stress strain diagrams are recorded and retrieved in digital electronic form. Genome scan was performed in The Jackson Laboratory (TJL). QTL mapping was conducted using Map Manager QTX software. Data show that (i) both elastic modulus (stiffness) and maximum loading (strength) value appear as normal distributions, suggesting that multiple genetic factors control the bone strength; (ii) 11 QTLs, accounting for 90% of variation for strength, have been detected. More than half QTLs of three-point bending are located on the same locations of bone density earlier identified from mouse femurs; (iii) a major QTL of femoral and vertebral bone mineral density (BMD) was not detected for bone strength of tibiae; (iv) the QTL on chromosome 4 has extremely high LOD score of 31.8 and represents 60% of the variation of bone strength; and (v) four QTLs of stiffness (chromosomes 2, 11, 15 and 19) have been identified.

Linkage mapping of femoral material properties in a reciprocal intercross of HcB-8 and HcB-23 recombinant mouse strains

Skeletal fragility is an important health problem with a large genetic component. We performed a 603 animal F2 reciprocal intercross of the recombinant congenic strains HcB-8 and HcB-23 to genetically map quantitative trait loci (QTLs) for tissue-level femoral biomechanical performance. These included elastic and post-yield strain, Young's modulus, elastic and maximum stress, and toughness and were calculated from 3-point bend testing of femora by the application of standard beam equations. We mapped these with R/qtl and QTL Cartographer and established significance levels empirically by permutation testing. Significant QTLs for at least one trait are present on chromosomes 1, 6, and 10 in the full F2 population, with additional QTLs evident in subpopulations defined by sex and cross direction. On chromosome 10, we find a QTL for post-yield strain and toughness, phenotypes that have not been mapped previously. Notably, the HcB-8 allele at this QTL increases post-yield strain and toughness, but decreases bone mineral density (BMD), while the material property QTLs on chromosomes 1, 6, and at a second chromosome 10 QTL are independent of BMD. We find significant sex x QTL and cross direction x QTL interactions. A robust, pleiotropic chromosome 4 QTL that we previously reported at the whole-bone level showed no evidence of linkage at the tissue-level, supporting our interpretation that modeling capacity is its primary phenotype. Our data demonstrate an inverse relationship between femoral perimeter and Young's modulus, with R(2)=0.27, supporting the view that geometric and material bone properties are subject to an integrated set of regulatory mechanisms. Mapping QTLs for tissue-level biomechanical performance advances understanding of the genetic basis of bone quality.

Genotype-phenotype associations in understanding the role of corticosteroid-binding globulin in health and disease animal models

Corticosteroid-binding globulin (CBG) is a plasma glycoprotein discovered more than 60 years ago for its high-affinity for glucocorticoids. Although its molecular structure and its biochemical properties have been described, its various biological roles and its importance are not yet fully understood. This review focuses first on studies that have used no-hypothesis-driven genetic approaches in animal models to reveal the higher than expected importance of CBG in particular in glucocorticoid stress responses. Then the dissection of some CBG physiological roles in an animal model of genetic CBG deficiency is reported. Finally, studies on the role of CBG genetic variability in human obesity traits are reviewed and discussed.

Variation in within-bone stiffness measured by nanoindentation in mice bred for high levels of voluntary wheel running

The hierarchical structure of bone, involving micro-scale organization and interaction of material components, is a critical determinant of macro-scale mechanics. Changes in whole-bone morphology in response to the actions of individual genes, physiological loading during life, or evolutionary processes, may be accompanied by alterations in underlying mineralization or architecture. Here, we used nanoindentation to precisely measure compressive stiffness in the femoral mid-diaphysis of mice that had experienced 37 generations of selective breeding for high levels of voluntary wheel running (HR). Mice (n = 48 total), half from HR lines and half from non-selected control (C) lines, were divided into two experimental groups, one with 13-14 weeks of access to a running wheel and one housed without wheels (n = 12 in each group). At the end of the experiment, gross and micro-computed tomography (microCT)-based morphometric traits were measured, and reduced elastic modulus (E(r)) was estimated separately for four anatomical quadrants of the femoral cortex: anterior, posterior, lateral, and medial. Two-way, mixed-model analysis of covariance (ancova) showed that body mass was a highly significant predictor of all morphometric traits and that structural change is more apparent at the microCT level than in conventional morphometrics of whole bones. Both line type (HR vs. C) and presence of the mini-muscle phenotype (caused by a Mendelian recessive allele and characterized by a approximately 50% reduction in mass of the gastrocnemius muscle complex) were significant predictors of femoral cortical cross-sectional anatomy. Measurement of reduced modulus obtained by nanoindentation was repeatable within a single quadrant and sensitive enough to detect inter-individual differences. Although we found no significant effects of line type (HR vs. C) or physical activity (wheel vs. no wheel) on mean stiffness, anterior and posterior quadrants were significantly stiffer (P < 0.0001) than medial and lateral quadrants (32.67 and 33.09 GPa vs. 29.78 and 30.46 GPa, respectively). Our findings of no significant difference in compressive stiffness in the anterior and posterior quadrants agree with previous results for mice, but differ from those for large mammals. Integrating these results with others from ongoing research on these mice, we hypothesize that the skeletons of female HR mice may be less sensitive to the effects of chronic exercise, due to decreased circulating leptin levels and potentially altered endocannabinoid signaling.

ENU induced single mutation locus on chr 16 leads to high-frequency hearing loss in mice

The hallmark of age-related (presbycusis) and noise-induced hearing loss is high-frequency (> 20 kHz) hearing loss. Through a collaborative study with TMGC (Tennessee Mouse Genome Consortium), seventeen ENU-induced mouse mutation strains with high-frequency hearing loss have been identified, but affected genes are yet identified. As a first step in identifying the gene/s underlying the ENU mutations, we created a F2 population between a representative mutation strain, 118 TNE and a wild type strain, CAST/EJ (CAST). Phenotypic analysis showed that there is a 3:1 ratio of segregation between normal and hearing loss in the F2 population, suggestion a single locus regulation. However, the linkage mapping identified 2 QTLs, each on chromosomes 15 and 16. Further statistical analysis of marker segregation patterns revealed that the locus on Chr 16 was ENU induced while the one on Chr 15 was derived from the parental strain, CAST.

Determining the elastic modulus of mouse cortical bone using electronic speckle pattern interferometry (ESPI) and micro computed tomography: a new approach for characterizing small-bone material properties

Mice phenotypes are invaluable for understanding bone formation and function, as well as bone disease. The elastic modulus is an important property of bones that can provide insights into bone quality. The determination of the elastic modulus of mouse cortical bone is complicated by the small dimensions of the bones. Whole bone bending tests are known to under estimate the elastic modulus compared to nanoindentation tests. The latter however provides information on extremely localized areas that do not necessarily correspond to the bulk elastic modulus under compression. This study presents a novel method for determining the bulk or effective elastic modulus of mouse cortical bone using the femur. We use Electronic Speckle Pattern Interferometry (ESPI), an optical method that enables the measurement of displacements on the bone surface, as it is compressed under water. This data is combined with geometric information obtained from micro-CT to calculate the elastic modulus. Roughly tubular cortical bone segments (2 mm) were cut from the diaphyses of femora of four week old C57BL/6 (B6) female mice and compressed axially using a mechanical tension-compression device. Displacements in the loading direction were mapped on the bone surface after loading the specimen. A linear regression of the displacement vs. axial-position enabled the calculation of the effective strain. Effective stress was calculated using force (N) data from the system's load cell and the mean cross-sectional area of the sample as determined by micro-CT. The effective elastic modulus (E) was calculated from the stress to strain ratio. The method was shown to be accurate and precise using a standard material machined to similar dimensions as those of the mouse femoral segments. Diaphyses of mouse femora were shown to have mean elastic moduli of 10.4+/-0.9 GPa for femora frozen for eight months, 8.6+/-1.4 GPa for femora frozen for two weeks and 8.9+/-1.1 GPa for the fresh femora. These values are much higher than those measured using three-point bending, and lower than values reported in the literature based on nanoindentation tests from mice bones of the same age. We show that this method can be used to accurately and precisely measure the effective elastic modulus of mouse cortical bone.

Systems analysis of bone

The genetic variants contributing to variability in skeletal traits has been well studied, and several hundred QTLs have been mapped and several genes contributing to trait variation have been identified. However, many questions remain unanswered. In particular, it is unclear whether variation in a single gene leads to alterations in function. Bone is a highly adaptive system and genetic variants affecting one trait are often accompanied by compensatory changes in other traits. The functional interactions among traits, which is known as phenotypic integration, has been observed in many biological systems, including bone. Phenotypic integration is a property of bone that is critically important for establishing a mechanically functional structure that is capable of supporting the forces imparted during daily activities. In this paper, bone is reviewed as a system and primarily in the context of functionality. A better understanding of the system properties of bone will lead to novel targets for future genetic analyses and the identification of genes that are directly responsible for regulating bone strength. This systems analysis has the added benefit of leaving a trail of valuable information about how the skeletal system works. This information will provide novel approaches to assessing skeletal health during growth and aging and for developing novel treatment strategies to reduce the morbidity and mortality associated with fragility fractures.

The use of nanoindentation for characterizing the properties of mineralized hard tissues: state-of-the art review

The use of nanoindentation to determine nanomechanical properties of mineralized tissues has been investigated extensively. A detailed, critical, and comprehensive review of this literature is the subject of the present work. After stating the motivation for the review, a succinct presentation of the challenges, advantages, and disadvantages of the various quasi-static nanoindentation test methods (to obtain elastic modulus, E, and hardness, H) and dynamic test methods (to obtain storage and loss moduli and/or loss/damping factor) is given in the form of a primer. Explicative summaries of literature reports on various intrinsic and extrinsic factors that significantly influence E and H, followed by 15 suggested topics for future research, are included additionally. This review is designed to present a compact guide to the principles of the nanoindentation technique and to emphasize considerations when determining material properties of mineralized tissues.

Tissue modulus calculated from beam theory is biased by bone size and geometry: implications for the use of three-point bending tests to determine bone tissue modulus

Current practice to determine bone tissue modulus of murine cortical bone is to estimate it from three-point bending tests, using Euler-Bernoulli beam theory. However, murine femora are not perfect beams; hence, results can be inaccurate. Our aim was to assess the accuracy of beam theory, which we tested for two commonly used inbred strains of mice, C57BL/6 (B6) and C3H/He (C3H). We measured the three-dimensional structure of male and female B6 and C3H femora (N=20/group) by means of micro-computed tomography. For each femur five micro-finite element (micro-FE) models were created that simulated three-point bending tests with varying distances between the supports. Tissue modulus was calculated from beam theory using micro-FE results. The accuracy of beam theory was assessed by comparing the beam theory-derived moduli with the modulus as used in the micro-FE analyses. An additional set of fresh-frozen femora (10 B6 and 12 C3H) was biomechanically tested and subjected to the same micro-FE analyses. These combined experimental-computational analyses enabled an unbiased assessment of specimen-specific tissue modulus. We found that by using beam theory, tissue modulus was underestimated for all femora. Femoral geometry and size had strong effects on beam theory-derived tissue moduli. Owing to their relatively thin cortex, underestimation was markedly higher for B6 than for C3H. Underestimation was dependent on support width in a strain-specific manner. From our combined experimental-computational approach we calculated tissue moduli of 12.0+/-1.3 GPa and 13.4+/-2.1 GPa for B6 and C3H, respectively. We conclude that tissue moduli in murine femora are strongly underestimated when calculated from beam theory. Using image-based micro-FE analyses we could precisely quantify this underestimation. We showed that previously reported murine inbred strain-specific differences in tissue modulus are largely an effect of geometric differences, not accounted for by beam theory. We suggest a re-evaluation of the tissue properties obtained from three-point bending tests, especially in mouse genetics.