The development of inter-strain variation in cortical and trabecular traits during growth of the mouse lumbar vertebral body

How cortical and trabecular bone co-develop to establish a mechanically functional structure is not well understood. Comparing early postnatal differences in morphology of lumbar vertebral bodies for three inbred mouse strains identified coordinated changes within and between cortical and trabecular traits. These early coordinate changes defined the phenotypic differences among the inbred mouse strains.

INTRODUCTION

Age-related changes in cortical and trabecular traits have been well studied; however, very little is known about how these bone tissues co-develop from day 1 of postnatal growth to establish functional structures by adulthood. In this study, we aimed to establish how cortical and trabecular tissues within the lumbar vertebral body change during growth for three inbred mouse strains that express wide variation in adult bone structure and function.

METHODS

Bone traits were quantified for lumbar vertebral bodies of female A/J, C57BL/6J (B6), and C3H/HeJ (C3H) inbred mouse strains from 1 to 105 days of age (n = 6-10 mice/age/strain).

RESULTS

Inter-strain differences in external bone size were observed as early as 1 day of age. Reciprocal and rapid changes in the trabecular bone volume fraction and alignment in the direction of axial compression were observed by 7 days of age. Importantly, the inter-strain difference in adult trabecular bone volume fraction was established by 7 days of age. Early variation in external bone size and trabecular architecture was followed by progressive increases in cortical area between 28 and 105 days of age, with the greatest increases in cortical area seen in the mouse strain with the lowest trabecular mass.

CONCLUSION

Establishing the temporal changes in bone morphology for three inbred mouse strains revealed that genetic variation in adult trabecular traits were established early in postnatal development. Early variation in trabecular architecture preceded strain-specific increases in cortical area and changes in cortical thickness. This study established the sequence of how cortical and trabecular traits co-develop during growth, which is important for identifying critical early ages to further focus on intervention studies that optimize adult bone strength.

Zoledronate treatment has different effects in mouse strains with contrasting baseline bone mechanical phenotypes

Two strains of mice with distinct bone morphologies and mechanical properties were treated with zoledronate. Our results show a different response to drug treatment in the two strains providing evidence that baseline properties of structure/material may influence response to zoledronate.

INTRODUCTION

Bisphosphonates are highly effective in reducing fracture risk, yet some individuals treated with these agents still experience fracture. The goal of this study was to test the hypothesis that genotype influences the effect of zoledronate on bone mechanical properties.

METHODS

Skeletally mature male mice from genetic backgrounds known to have distinct baseline post-yield properties (C57/B6, high post-yield displacement; A/J, low post-yield displacement) were treated for 8 weeks with saline (VEH) or zoledronate (ZOL, 0.06 mg/kg subcutaneously once every 4 weeks) in a 2 × 2 study design. Ex vivo μCT and mechanical testing (4-pt bending) were conducted on the femur to assess morphological and mechanical differences.

RESULTS

Significant drug and/or genotype effects were found for several mechanical properties and significant drug × genotype interactions were found for measures of strength (ultimate force) and brittleness (total displacement, strain to failure). Treatment with ZOL affected bone biomechanical measures of brittleness (total displacement (-25 %) and strain to failure (-23 %)) in B6 mice significantly differently than in A/J mice. This was driven by unique drug × genotype effects on bone geometry in B6 animals yet likely also reflected changes to the tissue properties.

CONCLUSION

These data may support the concept that properties of the bone geometry and/or tissue at the time of treatment initiation play a role in determining the bone's mechanical response to zoledronate treatment.

Smaller, weaker, and less stiff bones evolve from changes in subsistence strategy

We propose a computational model with which to examine the evolution of bone. Our results indicate that changes in subsistence strategy have influenced the evolution of bone growth and mechanoregulation, and predict that bone size, stiffness, and structural strength may decrease in future generations, bringing increased risk of fracture and prevalence of osteoporosis.

INTRODUCTION

Archeological data suggest that bone size and strength have decreased over evolution. We hypothesize that changing evolutionary pressures and levels of physical activity, both arising from changes in subsistence strategy, have affected the evolution of bone. We propose a computational model with which to examine the evolution of bone growth and mechanoregulation due to the transitions from hunter-gatherer to agricultural to modern lifestyles.

METHODS

The evolution of genes governing growth and mechanoregulation in a population of bones is simulated, where each individual is represented by a 2-D bone cross-section. Genetic variability is assumed to modulate growth through mechanoregulatory factors that direct periosteal expansion, endosteal expansion/infilling, and ash content accretion in response to strains incurred during walking.

RESULTS

The model predicts decreases in cortical area and section modulus (a measure of structural strength) and increases in maximum compressive strain over the course of the simulation, meaning evolution of smaller, less strong, and less stiff bones is predicted for the population average. The model predicts small but continued decreases in size, strength, and stiffness in modern populations, despite the absence of a strong evolutionary advantage to efficient bones during this phase.

CONCLUSION

In conclusion, our results show that changing loading regimes and evolutionary pressures may have influenced the evolution of bone growth and mechanoregulation, and predict that bone size and strength may continue to decrease in future generations, bringing increased risk of fracture and prevalence of osteoporosis.

Activation of bone remodeling after fatigue: differential response to linear microcracks and diffuse damage

Recent experiments point to two predominant forms of fatigue microdamage in bone: linear microcracks (tens to a few hundred microns in length) and "diffuse damage" (patches of diffuse stain uptake in fatigued bone comprised of clusters of sublamellar-sized cracks). The physiological relevance of diffuse damage in activating bone remodeling is not known. In this study microdamage amount and type were varied to assess whether linear or diffuse microdamage has similar effects on the activation of intracortical resorption. Activation of resorption was correlated to the number of linear microcracks (Cr.Dn) in the bone (R(2)=0.60, p<0.01). In contrast, there was no activation of resorption in response to diffuse microdamage alone. Furthermore, there was no significant change in osteocyte viability in response to diffuse microdamage, suggesting that osteocyte apoptosis, which is known to activate remodeling at typical linear microcracks in bone, does not result from sublamellar damage. These findings indicate that inability of diffuse microdamage to activate resorption may be due to lack of a focal injury response. Finally, we found that duration of loading does not affect the remodeling response. In conclusion, our data indicate that osteocytes activate resorption in response to linear microcracks but not diffuse microdamage, perhaps due to lack of a focal injury-induced apoptotic response.

Genetic variation in the structural pattern of osteoclast activity during post-natal growth of mouse femora

While the spatial activity of osteoblasts has been associated with modeling of bones during development, few studies have examined if variation in the spatial activity of osteoclasts also contributes to the morphogenesis of skeletal tissues. We examined this question by histomorphometric analysis and reconstructing the three-dimensional spatial distribution of osteoclasts in the femora of three inbred strains of male mice (A/J, C57BL/6J [B6], and C3H/HeJ [C3H]) that have differing skeletal, structural, and material properties. Our data show that total osteoclast surface area and osteoclast numbers are related to the overall bone density, but not related to the development of bone diameter or overall cortical area. The analysis of the spatial distribution of the osteoclasts showed that the asymmetrical mid-diaphyseal distribution of osteoclasts in A/J and B6 compared to the more uniform distribution of these cells around the circumference in the C3H mice was consistent with the more ellipsoid shape of A/J and B6 femora compared to the more circular mid-diaphyseal shape of the femora in the C3H mice. The statistically 2- to 3-fold fewer cells on the periosteal surface in the C3H compared to either the B6 or A/J mice is also consistent with the greater cortical thickness that is seen for the C3H mice compared to either B6 or A/J strains. In vitro studies of osteoclastogenesis and the expression of numerous phenotypic properties of osteoclasts prepared from the three strains of mice showed that A/J and B6 mice developed statistically greater numbers of tartrate resistant acid phosphatase (TRAP) positive cells and expressed statistically higher levels of multiple mRNAs that are unique to differentiated osteoclasts than those isolated from the C3H strain. In summary, the 3D reconstructions and histomorphometric analysis suggest that genetic differences lead to spatial variation in the distribution of osteoclasts. These variations in spatial distribution of osteoclasts in turn contribute in part to the development of the structural variations of the femora that are seen in the three strains of mice. In vitro studies suggest that intrinsic genetic variation in osteoclastogenesis and their phenotypic expression may contribute to the differences in their functional activities that give rise to the unique spatial distributions of these cells in bones.

Osteocyte apoptosis and control of bone resorption following ovariectomy in mice

INTRODUCTION

Osteocyte apoptosis has been linked to bone resorption resulting from estrogen depletion and other resorptive stimuli; however, precise spatial and temporal relationships between the two events have not been clearly established. The purpose of this study was to characterize the patterns of osteocyte apoptosis in relation to bone resorption following ovariectomy to test whether osteocyte apoptosis occurs preferentially in areas known to activate resorption. Moreover, we report that osteocyte apoptosis is necessary to initiate endocortical remodeling in response to estrogen withdrawal.

MATERIALS AND METHODS

Adult female C57BL/6J mice (17 weeks old) underwent either bilateral ovariectomy (OVX), or sham surgery (SHAM) and were euthanized on days 3, 7, 14, or 21 days after OVX. Diaphyseal cross-sections were stained by immunohistochemistry for activated caspase-3 as a marker of apoptosis. The percentages of caspase-positive stained osteocytes (Casp+Ot.) were measured along major and minor anatomical axes around the femoral diaphysis to evaluate the distribution of osteocyte apoptosis after estrogen loss; resorption surface was measured at the adjacent endocortical regions. In a second study to test whether osteocyte apoptosis plays a regulatory role in the initiation of bone resorption, a group of OVX mice received the pan-caspase inhibitor, QVDOPh, to inhibit osteocyte apoptosis. Remaining experimental and sham groups received either QVD or Vehicle.

RESULTS

OVX increased osteocyte apoptosis in a non-uniform distribution throughout the femoral diaphyses. Increases in Casp+osteocytes were predominantly located in the posterior diaphyseal cortex. Here, the number of apoptotic osteocytes 4- to 7-fold higher than sham controls (p<0.005) by day 3 post-OVX and remained elevated. Increases in resorption post-OVX also occurred along the posterior endocortical surface overlying the region of osteocyte apoptosis, but these increases occurred only at 14 and 21 days post-OVX (p<0.002) well after the increases in osteocyte apoptosis. Treatment with QVD in OVX animals suppressed osteocyte apoptosis, with levels in QVD-treated samples equivalent to baseline. Moreover, the increases in osteoclastic resorption normally observed after estrogen loss did not occur in OVX mice treated with QVD.

CONCLUSIONS

The results of this study demonstrate that osteocyte apoptosis following estrogen loss occur regionally, rather than uniformly throughout the cortex. We also showed that estrogen loss increased osteocyte apoptosis. Apoptotic osteocytes were overwhelmingly localized within the posterior cortical region, the location where endocortical resorption was subsequently activated in ovariectomized mice. Finally, the increases in osteoclastic resorption normally observed after estrogen withdrawal did not occur in the absence of osteocyte apoptosis indicating that this apoptosis is necessary to activate endocortical remodeling following estrogen loss.

Lack of sexual dimorphism in femora of the eusocial and hypogonadic naked mole-rat: a novel animal model for the study of delayed puberty on the skeletal system

Sex steroid hormones are major determinants of bone morphology and quality and are responsible for sexually dimorphic skeletal traits. Hypogonadism results in suboptimal skeletal development and may lead to an increased risk of bone fracture later in life. The etiology of delayed puberty and/or hypothalamic amenorrhea is poorly understood, and experimental animal models addressing this issue are predominantly based upon short-term experimental induction of hormonal suppression via gonadotropin releasing hormone antagonists (GnRH-a). This acute change in hormone profile does not necessarily emulate the natural progression of hypogonadic bone disorders. We propose a novel animal model with which to explore the effects of chronic hypogonadism on bone quality, the naked mole-rat (NMR; Heterocephalus glaber). This mouse-size rodent may remain reproductively suppressed throughout its life, if it remains as a subordinate within the eusocial mole-rat colony. NMRs live in large colonies with a single dominant breeding female. She, primarily by using aggressive social contact, naturally suppresses the hypothalamic gonadotropic axis of subordinate NMRs and thereby their reproductive expression. However, should an NMR be separated from the dominant breeder, within less than a week reproductive hormones may become elevated and the animal attains breeding status. We questioned if sexual suppression of subordinates impact upon the development and maintenance of the femora and lead to a sexually indistinct monomorphic skeleton. Femora were obtained from male and female NMRs that were either non-breeders (subordinate) or breeders at the time of sacrifice. Diaphyseal cross-sectional morphology, metaphyseal trabecular micro-architecture and tissue mineral density of the femur were measured using microcomputed tomography and diaphyseal mechanical properties were assessed by four-point bending tests to failure. Subordinates were sexually monomorphic and showed no significant differences in body weight or femoral bone structure and quality between males and females. Femora of subordinate females differed significantly from that of breeding animals, whereas in males, the divergent trend among breeders and non-breeders did not reach statistical significance. Subordinate NMRs, naturally suppressed from entering puberty, may prove to be a useful model to tease apart the relationship between bone morphology and hypogonadism and evaluate skeletal development during pubertal maturation.

The mechanical consequences of mineralization in embryonic bone

The purpose of this study was to examine the effect of mineralization on the mechanical properties of embryonic bone rudiments. For this purpose, four-point bending experiments were performed on unmineralized and mineralized embryonic mouse ribs at 16 and 17 days of gestational age. Young's modulus was calculated using force-displacement data from the experiment in combination with finite element analysis (FEA). For the unmineralized specimens, a calculated average for the Young's modulus of 1.11 (+/- 0.62) MPa was established after corrections for sticking to the four-point bending device and aspect ratio, which is the ratio between the length of the bone and its diameter. For the mineralized specimens, the value was 117 (+/- 62) MPa after corrections. Hence, Young's moduli of embryonic bone rudiments increase by two orders of magnitude within 1 day, during endochondral ossification. As an effect, the hypertrophic chondrocytes in the calcifying cartilage experience a significant change in their mechanical environment. The chondrocytes are effectively stress shielded, which means that they do not carry stresses since stresses are supported by the stiffest parts of the tissue, which are in this case the diaphyseal cortex and the calcified matrix. The deformability of the hypertrophic chondrocytes is, therefore, severely reduced. Since the transition is so sudden and enormous, it can be seen as a process of 'catastrophic' proportion for the hypertrophic chondrocytes. The subsequent resorption of calcified cartilage and the expansion of the marrow cavity could be consequential to stress shielding.

Cyclic hydrostatic pressure enhances the chondrogenic phenotype of human mesenchymal progenitor cells differentiated in vitro

Much attention has been given to the influences of bioactive factors on mesenchymal progenitor cell differentiation and proliferation, but few studies have examined the effect of mechanical factors on these cells. This study examined the effects of cyclic hydrostatic pressure on human bone marrow-derived mesenchymal progenitor cells undergoing chondrogenic differentiation. Aggregates of bone marrow-derived mesenchymal progenitor cells were cultured in a defined chondrogenic medium and were subjected to cyclic hydrostatic pressure. Aggregates were loaded at various time points: single (day 1 or 3) or multiple (days 1-7). At 14 and 28 days, aggregates were harvested for histology, immunohistochemistry, and quantitative DNA and matrix macromolecule analysis. The aggregates loaded for a single day did not demonstrate significant changes in proteoglycan and collagen contents compared with the non-loaded controls. In contrast, for the multi-day loaded aggregates, statistically significant increases in proteoglycan and collagen contents were found on both day 14 and day 28. Aggregates loaded for seven days were larger and histological staining indicated a greater matrix/cell ratio. This study indicates that hydrostatic pressure enhances the cartilaginous matrix formation of mesenchymal progenitor cells differentiated in vitro, and suggests that mechanical forces may play an important role in cartilage repair and regeneration in vivo.

Bone brittleness varies with genetic background in A/J and C57BL/6J inbred mice

The contribution of genetic and environmental factors to variations in bone quality are understood poorly. We tested whether bone brittleness varies with genetic background using the A/J and C57BL/6J inbred mouse strains. Whole bone four-point bending tests revealed a 70% decrease in postyield deflection of A/J femurs compared with C57BL/6J, indicating that A/J femurs failed in a significantly more brittle manner. Cyclic loading studies indicated that A/J femurs accumulated damage differently than C57BL/6J femurs, consistent with their increased brittleness. Differences in matrix composition also were observed between the two mouse strains. A/J femurs had a 4.5% increase in ash content and an 11.8% decrease in collagen content. Interestingly, a reciprocal relationship was observed between femoral geometry and material stiffness; this relationship may have contributed to the brittle phenotype of A/J femurs. A/J femurs are more slender than those of C57BL/6J femurs; however, their 47% smaller moment of inertia appeared to be compensated by an increased tissue stiffness at the expense of altered tissue damageability. Importantly, these differences in whole bone mechanical properties between A/J and C57BL/6J femurs could not have been predicted from bone mass or density measures alone. The results indicated that bone brittleness is a genetically influenced trait and that it is associated with genetically determined differences in whole bone architecture, bone matrix composition, and mechanisms of cyclical damage accumulation.

The role of the lamellar interface during torsional yielding of human cortical bone

Fragility fractures are a result of alterations in bone quantity, tissue properties, applied loads, or a combination of these factors. The current study addresses the contribution of cortical bone tissue properties to skeletal fragility by characterizing the shear damage accumulation processes which occur during torsional yielding in normal bone. Samples of human femoral cortical bone were loaded in torsion and damaged at a post-yield twist level. The number of microcracks within osteons, interstitial tissue, and along cement lines were assessed using basic fuchsin staining. Damage density measures (number of cracks/mm2) were correlated with stiffness degradation and changes in relaxation. Damaged samples exhibited a wide variation in total microcrack density, ranging from 1.1 to 43.3 cracks/mm2 with a mean density of 19.7 +/- 9.8 cracks/mm2. Lamellar interface cracks comprised more than 75% of the total damage, indicating that the lamellar interface is weak in shear and is a principal site of shear damage accumulation. Damage density was positively correlated with secant stiffness degradation, but only explained 22% of the variability in degradation. In contrast, damage density was uncorrelated with the changes in relaxation, indicating that a simple crack counting measure such as microcrack density was not an appropriate measure of relaxation degradation. Finally, a nonuniform microcrack density distribution was observed, suggesting that internal shear stresses were redistributed within the torsion samples during post-yield loading. The results suggested that the lamellar interface in human cortical bone plays an important role in torsional yielding by keeping cracks physically isolated from each other and delaying microcrack coalescence in order to postpone the inevitable formation of the fatal crack.

Lumican regulates collagen fibril assembly: skin fragility and corneal opacity in the absence of lumican

Lumican, a prototypic leucine-rich proteoglycan with keratan sulfate side chains, is a major component of the cornea, dermal, and muscle connective tissues. Mice homozygous for a null mutation in lumican display skin laxity and fragility resembling certain types of Ehlers-Danlos syndrome. In addition, the mutant mice develop bilateral corneal opacification. The underlying connective tissue defect in the homozygous mutants is deregulated growth of collagen fibrils with a significant proportion of abnormally thick collagen fibrils in the skin and cornea as indicated by transmission electron microscopy. A highly organized and regularly spaced collagen fibril matrix typical of the normal cornea is also missing in these mutant mice. This study establishes a crucial role for lumican in the regulation of collagen assembly into fibrils in various connective tissues. Most importantly, these results provide a definitive link between a necessity for lumican in the development of a highly organized collagenous matrix and corneal transparency.

Type I collagen mutation alters the strength and fatigue behavior of Mov13 cortical tissue

Despite advances in understanding the molecular basis of Osteogenesis Imperfecta, the mechanisms by which type I collagen mutations compromise whole bone function are not well understood. Previously, we have shown that a heterozygous type I collagen mutation is associated with increased brittleness of long bones from Mov13 transgenic mice, a model of the mild form of Osteogenesis Imperfecta. In the current study, we investigated tissue-level damage processes by testing the hypothesis that the fatigue properties of Mov13 tissue were significantly compromised relative to littermate controls. We also quantified tissue structure and mineral content to explain variations in the fatigue behavior. Micro-beam specimens were machined from the anterior and posterior quadrants of Mov13 and control femurs and subjected to cyclic bending at one of four stress levels. Mov13 tissue exhibited a 22-25% reduction in tissue bending strength and a similar reductions in fatigue life and the stress level at which damage was apparent. These results provided tissue-level evidence that damage accumulation mechanisms were significantly compromised in Mov13 cortical tissue. Given that significant alterations in tissue structure were observed in Mov13 femurs, the results of this study support the idea that Mov13 femurs were brittle because alterations in tissue structure associated with the mutation interfered with normal damage processes. These results provide new insight into the pathogenesis of Osteogenesis Imperfecta and are consistent with bone behaving as a damaging composite material, where damage accumulation is central to bone fracture.

Comparison of damage accumulation measures in human cortical bone

Elastic modulus degradation, strength reduction, and energy dissipation have traditionally been the properties of choice to monitor the damage process in cortical bone. However, these properties only provide limited insight into the damage process given the complex mechanical nature of bone. In the current study, alternative measures of the damage process were investigated for machined human cortical bone specimens loaded under torsion. Seventy-two bone specimens from 6 human femurs were subjected to a series of torsional relaxation cycles in which damage was induced during a single relaxation cycle and the effects of damage on the elastic, yield, viscous, and failure properties were determined from pre- and post-damage relaxation cycles. The results revealed that degradation of all torsion properties exhibited a significant twist magnitude effect. However, the yield stress and strain, the relaxation rate, and the total relaxation exhibited 5-10 fold greater degradation than both strength and modulus, when residual strength tests were conducted at high shear strain rates. For the loading conditions examined in this study, the results indicated that the relaxation and yield properties of cortical bone are more sensitive to shear damage accumulation and better measures of the damage process than either strength or modulus. Further, the results reveal an important interaction between damage and the viscous behavior of bone which provides new insight into the effects of damage on bone mechanical properties.

Type-I collagen mutation compromises the post-yield behavior of Mov13 long bone

Despite recent advances in our understanding of the molecular basis of skeletal fragility, little is known about how these molecular alterations lead to whole bone brittleness. In the current study, we investigated the relationship between a type-I collagen mutation and post-yield behavior of whole bone in Mov13 transgenic mice by considering tissue-level organizational issues known to be important for normal bone fracture. Mechanical assays revealed that the post-yield deflection of Mov13 femurs was reduced by 61% relative to littermate controls. Fractographic images revealed that lamellar interfaces which were important for dissipating energy during the failure process of control femurs, were not effective in Mov13 mice. Further investigation revealed that a 22% reduction in bone collagen content, a 2-fold increase in tissue porosity, and significant alterations in collagen organization interfered with normal energy dissipation mechanisms of Mov13 microstructure. Collectively, the results provided the first evidence that the reduced ductility associated with a type-I collagen mutation was mediated by alterations in intermediate structures that normally contribute to the post-yield behavior of cortical bone. The results suggest that, to better understand the pathogenesis of skeletal fragility, it is important to consider the effects of molecular alterations on higher-level structures, particularly those structures that contribute to the failure mechanisms in normal bone.

Osteoporosis induced in mice by overproduction of interleukin 4

Osteoporosis is a common disease in which loss of bone mass results in skeletal fragility. The development of therapies for this disorder has been hampered by the lack of a convenient animal model. Here we describe a disorder in bone homeostasis in transgenic mice that inappropriately express the cytokine interleukin 4 (IL-4) under the direction of the lymphocyte-specific proximal promoter for the lck gene. Bone disease in lck-IL-4 mice appeared to result from markedly decreased bone formation by osteoblasts, features strikingly similar to those observed in cases of severe low-turnover human involutional osteoporosis. By 2 months of age, female and male lck-IL-4 mice invariably developed severe osteoporosis of both cortical and trabecular bone. Osteoporosis was observed in two independently derived founder animals, indicating that this phenotype was directly mediated by the IL-4 transgene.

A murine skeletal adaptation that significantly increases cortical bone mechanical properties. Implications for human skeletal fragility

Mov13 mice carry a provirus that prevents transcription initiation of the alpha 1(I) collagen gene. Mutant mice homozygous for the null mutation produce no type I collagen and die at mid-gestation, whereas heterozygotes survive to adulthood. Dermal fibroblasts from heterozygous mice produce approximately 50% less type I collagen than normal littermates, and the partial deficiency in collagen production results in a phenotype similar to osteogenesis imperfecta type I (an inherited form of skeletal fragility). In this study, we have identified an adaptation of Mov13 skeletal tissue that significantly improves the bending strength of long bone. The adaptive response occurred over a 2-mo period, during which time a small number of newly proliferated osteogenic cells produced a significant amount of matrix components and thus generated new bone along periosteal surfaces. New bone deposition resulted in a measurable increase in cross-sectional geometry which, in turn, led to a dramatic increase in long bone bending strength.

Application of homogenization theory to the study of trabecular bone mechanics

It is generally accepted that the strength and stiffness of trabecular bone is strongly affected by trabecular microstructure. It has also been hypothesized that stress induced adaptation of trabecular bone is affected by trabecular tissue level stress and/or strain. At this time, however, there is no generally accepted (or easily accomplished) technique for predicting the effect of microstructure on trabecular bone apparent stiffness and strength or estimating tissue level stress or strain. In this paper, a recently developed mechanics theory specifically designed to analyze microstructured materials, called the homogenization theory, is presented and applied to analyze trabecular bone mechanics. Using the homogenization theory it is possible to perform microstructural and continuum analyses separately and then combine them in a systematic manner. Stiffness predictions from two different microstructural models of trabecular bone show reasonable agreement with experimental results, depending on metaphyseal region, (R2 greater than 0.5 for proximal humerus specimens, R2 less than 0.5 for distal femur and proximal tibia specimens). Estimates of both microstructural strain energy density (SED) and apparent SED show that there are large differences (up to 30 times) between apparent SED (as calculated by standard continuum finite element analyses) and the maximum microstructural or tissue SED. Furthermore, a strut and spherical void microstructure gave very different estimates of maximum tissue SED for the same bone volume fraction (BV/TV). The estimates from the spherical void microstructure are between 2 and 20 times greater than the strut microstructure at 10-20% BV/TV.