The mechanical consequences of mineralization in embryonic bone

Mechanoregulatory role of TRPV4 in prenatal skeletal development

Biophysical cues are essential for guiding skeletal development, but the mechanisms underlying the mechanical regulation of cartilage and bone formation are unknown. TRPV4 is a mechanically sensitive ion channel involved in cartilage and bone cell mechanosensing, mutations of which lead to skeletal developmental pathologies. We tested the hypothesis that loading-driven prenatal skeletal development is dependent on TRPV4 activity. We first establish that mechanically stimulating mouse embryo hindlimbs cultured ex vivo stimulates knee cartilage growth, morphogenesis, and expression of TRPV4, which localizes to areas of high biophysical stimuli. We then demonstrate that loading-driven joint cartilage growth and shape are dependent on TRPV4 activity, mediated via control of cell proliferation and matrix biosynthesis, indicating a mechanism by which mechanical loading could direct growth and morphogenesis during joint formation. We conclude that mechanoregulatory pathways initiated by TRPV4 guide skeletal development; therefore, TRPV4 is a valuable target for the development of skeletal regenerative and repair strategies.

Mechanical feedback defines organizing centers to drive digit emergence

During embryonic development, digits gradually emerge in a periodic pattern. Although genetic evidence indicates that digit formation results from a self-organizing process, the underlying mechanisms are still unclear. Here, we find that convergent-extension tissue flows driven by active stresses underlie digit formation. These active stresses simultaneously shape cartilage condensations and lead to the emergence of a compressive stress region that promotes high activin/p-SMAD/SOX9 expression, thereby defining digit-organizing centers via a mechanical feedback. In Wnt5a mutants, such mechanical feedback is disrupted due to the loss of active stresses, organizing centers do not emerge, and digit formation is precluded. Thus, digit emergence does not result solely from molecular interactions, as was previously thought, but requires a mechanical feedback that ensures continuous coupling between phalanx specification and elongation. Our work, which links mechanical and molecular signals, provides a mechanistic context for the emergence of organizing centers that may underlie various developmental processes.

A Direct Comparison of Node and Element-Based Finite Element Modeling Approaches to Study Tissue Growth

Finite element analysis is a useful tool to model growth of biological tissues and predict how growth can be impacted by stimuli. Previous work has simulated growth using node-based or element-based approaches, and this implementation choice may influence predicted growth, irrespective of the applied growth model. This study directly compared node-based and element-based approaches to understand the isolated impact of implementation method on growth predictions by simulating growth of a bone rudiment geometry, and determined what conditions produce similar results between the approaches. We used a previously reported node-based approach implemented via thermal expansion and an element-based approach implemented via osmotic swelling, and we derived a mathematical relationship to relate the growth resulting from these approaches. We found that material properties (modulus) affected growth in the element-based approach, with growth completely restricted for high modulus values relative to the growth stimulus, and no restriction for low modulus values. The node-based approach was unaffected by modulus. Node- and element-based approaches matched marginally better when the conversion coefficient to relate the approaches was optimized based on the results of initial simulations, rather than using the theoretically predicted conversion coefficient (median difference in node position 0.042 cm versus 0.052 cm, respectively). In summary, we illustrate here the importance of the choice of implementation approach for modeling growth, provide a framework for converting models between implementation approaches, and highlight important considerations for comparing results in prior work and developing new models of tissue growth.

Craniofacial skeletal response to encephalization: How do we know what we think we know?

Dramatic changes in cranial capacity have characterized human evolution. Important evolutionary hypotheses, such as the spatial packing hypothesis, assert that increases in relative brain size (encephalization) have caused alterations to the modern human skull, resulting in a suite of traits unique among extant primates, including a domed cranial vault, highly flexed cranial base, and retracted facial skeleton. Most prior studies have used fossil or comparative primate data to establish correlations between brain size and cranial form, but the mechanistic basis for how changes in brain size impact the overall shape of the skull resulting in these cranial traits remains obscure and has only rarely been investigated critically. We argue that understanding how changes in human skull morphology could have resulted from increased encephalization requires the direct testing of hypotheses relating to interaction of embryonic development of the bones of the skull and the brain. Fossil and comparative primate data have thoroughly described the patterns of association between brain size and skull morphology. Here we suggest complementing such existing datasets with experiments focused on mechanisms responsible for producing the observed patterns to more thoroughly understand the role of encephalization in shaping the modern human skull.

Stresses and strains on the human fetal skeleton during development

Mechanical forces generated by fetal kicks and movements result in stimulation of the fetal skeleton in the form of stress and strain. This stimulation is known to be critical for prenatal musculoskeletal development; indeed, abnormal or absent movements have been implicated in multiple congenital disorders. However, the mechanical stress and strain experienced by the developing human skeleton in utero have never before been characterized. Here, we quantify the biomechanics of fetal movements during the second half of gestation by modelling fetal movements captured using novel cine-magnetic resonance imaging technology. By tracking these movements, quantifying fetal kick and muscle forces, and applying them to three-dimensional geometries of the fetal skeleton, we test the hypothesis that stress and strain change over ontogeny. We find that fetal kick force increases significantly from 20 to 30 weeks' gestation, before decreasing towards term. However, stress and strain in the fetal skeleton rises significantly over the latter half of gestation. This increasing trend with gestational age is important because changes in fetal movement patterns in late pregnancy have been linked to poor fetal outcomes and musculoskeletal malformations. This research represents the first quantification of kick force and mechanical stress and strain due to fetal movements in the human skeleton in utero, thus advancing our understanding of the biomechanical environment of the uterus. Further, by revealing a potential link between fetal biomechanics and skeletal malformations, our work will stimulate future research in tissue engineering and mechanobiology.

Normal trabecular vertebral bone is formed via rapid transformation of mineralized spicules: A high-resolution 3D ex-vivo murine study

At birth, mouse vertebrae have a reticular fine spongy morphology, yet in the adult animal they exhibit elaborate trabecular architectures. Here, we characterize the physiological microstructural transformations in growing young female mice of the widely used C57BL/6 strain. Extensive architectural changes lead to the establishment of mature cancellous bone in the spine. Vertebrae were mapped in 3D by high resolution microcomputed tomography (µCT), backed by conventional histology. Three different phases are observed in the natural bony biomaterial: In a prenatal templating phase, early vertebrae are composed of foamy, loosely-packed mineralized spicules. During a consolidation phase in the first 7 days after birth, bone material condenses into struts and forms primitive trabeculae accompanied by a significant (>50%) reduction in bone volume/tissue volume ratio (BV/TV). After day 7, the trabeculae expand, reorient and increase in mineral density. Swift growth ensues such that by day 14 the young lumbar spine exhibits all morphological features observed in the mature animal. The greatly varied micro-morphologies of normal trabecular bone observed in 3D within a short timespan are typical for rodent and presumably for other mammalian forming spines. This suggests that fully structured cancellous bone emerges through rapid post-natal restructuring of a foamy mineralized scaffold. STATEMENT OF SIGNIFICANCE: Cancellous bone develops in stages that are not well documented. Using a mouse model, we provide an observer-independent quantification of normal bone formation in the spine. We find that within 14 days, the cancellous bone transforms in 3 phases from a scaffold of spicules into well organized, fully mineralized trabeculae in a functional spine. Detailed knowledge of the physiological restructuring of mineralized material may help to better understand bone formation and may serve as a blueprint for studies of pharmaceuticals effects, tissue healing and regeneration.

Building Finite Element Models to Investigate Zebrafish Jaw Biomechanics

Skeletal morphogenesis occurs through tightly regulated cell behaviors during development; many cell types alter their behavior in response to mechanical strain. Skeletal joints are subjected to dynamic mechanical loading. Finite element analysis (FEA) is a computational method, frequently used in engineering that can predict how a material or structure will respond to mechanical input. By dividing a whole system (in this case the zebrafish jaw skeleton) into a mesh of smaller 'finite elements', FEA can be used to calculate the mechanical response of the structure to external loads. The results can be visualized in many ways including as a 'heat map' showing the position of maximum and minimum principal strains (a positive principal strain indicates tension while a negative indicates compression. The maximum and minimum refer the largest and smallest strain). These can be used to identify which regions of the jaw and therefore which cells are likely to be under particularly high tensional or compressional loads during jaw movement and can therefore be used to identify relationships between mechanical strain and cell behavior. This protocol describes the steps to generate Finite Element models from confocal image data on the musculoskeletal system, using the zebrafish lower jaw as a practical example. The protocol leads the reader through a series of steps: 1) staining of the musculoskeletal components, 2) imaging the musculoskeletal components, 3) building a 3 dimensional (3D) surface, 4) generating a mesh of Finite Elements, 5) solving the FEA and finally 6) validating the results by comparison to real displacements seen in movements of the fish jaw.

Modeling the biomechanics of fetal movements

Fetal movements in the uterus are a natural part of development and are known to play an important role in normal musculoskeletal development. However, very little is known about the biomechanical stimuli that arise during movements in utero, despite these stimuli being crucial to normal bone and joint formation. Therefore, the objective of this study was to create a series of computational steps by which the forces generated during a kick in utero could be predicted from clinically observed fetal movements using novel cine-MRI data of three fetuses, aged 20–22 weeks. A custom tracking software was designed to characterize the movements of joints in utero, and average uterus deflection of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6.95 \pm 0.41$$\end{document}6.95±0.41 mm due to kicking was calculated. These observed displacements provided boundary conditions for a finite element model of the uterine environment, predicting an average reaction force of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.52 \pm 0.15$$\end{document}0.52±0.15 N generated by a kick against the uterine wall. Finally, these data were applied as inputs for a musculoskeletal model of a fetal kick, resulting in predicted maximum forces in the muscles surrounding the hip joint of approximately 8 N, while higher maximum forces of approximately 21 N were predicted for the muscles surrounding the knee joint. This study provides a novel insight into the closed mechanical environment of the uterus, with an innovative method allowing elucidation of the biomechanical interaction of the developing fetus with its surroundings.

Finite element modelling predicts changes in joint shape and cell behaviour due to loss of muscle strain in jaw development

Abnormal joint morphogenesis is linked to clinical conditions such as Developmental Dysplasia of the Hip (DDH) and to osteoarthritis (OA). Muscle activity is known to be important during the developmental process of joint morphogenesis. However, less is known about how this mechanical stimulus affects the behaviour of joint cells to generate altered morphology. Using zebrafish, in which we can image all joint musculoskeletal tissues at high resolution, we show that removal of muscle activity through anaesthetisation or genetic manipulation causes a change to the shape of the joint between the Meckel's cartilage and Palatoquadrate (the jaw joint), such that the joint develops asymmetrically leading to an overlap of the cartilage elements on the medial side which inhibits normal joint function. We identify the time during which muscle activity is critical to produce a normal joint. Using Finite Element Analysis (FEA), to model the strains exerted by muscle on the skeletal elements, we identify that minimum principal strains are located at the medial region of the joint and interzone during mouth opening. Then, by studying the cells immediately proximal to the joint, we demonstrate that biomechanical strain regulates cell orientation within the developing joint, such that when muscle-induced strain is removed, cells on the medial side of the joint notably change their orientation. Together, these data show that biomechanical forces are required to establish symmetry in the joint during development.

Effects of normal and abnormal loading conditions on morphogenesis of the prenatal hip joint: application to hip dysplasia

Joint morphogenesis is an important phase of prenatal joint development during which the opposing cartilaginous rudiments acquire their reciprocal and interlocking shapes. At an early stage of development, the prenatal hip joint is formed of a deep acetabular cavity that almost totally encloses the head. By the time of birth, the acetabulum has become shallower and the femoral head has lost substantial sphericity, reducing joint coverage and stability. In this study, we use a dynamic mechanobiological simulation to explore the effects of normal (symmetric), reduced and abnormal (asymmetric) prenatal movements on hip joint shape, to understand their importance for postnatal skeletal malformations such as developmental dysplasia of the hip (DDH). We successfully predict the physiological trends of decreasing sphericity and acetabular coverage of the femoral head during fetal development. We show that a full range of symmetric movements helps to maintain some of the acetabular depth and femoral head sphericity, while reduced or absent movements can lead to decreased sphericity and acetabular coverage of the femoral head. When an abnormal movement pattern was applied, a deformed joint shape was predicted, with an opened asymmetric acetabulum and the onset of a malformed femoral head. This study provides evidence for the importance of fetal movements in the prevention and manifestation of congenital musculoskeletal disorders such as DDH.

Mechanical factors explain development of cam-type deformity

OBJECTIVE

A cam-type deformity drastically increases the risk of hip osteoarthritis (OA). Since this type of skeletal anomaly is more prevalent among young active adults, it is hypothesized that the loading conditions experienced during certain types of vigorous physical activities stimulates formation of cam-type deformity. We further hypothesize that the growth plate shape modulates the influence of mechanical factors on the development of cam-type deformity.

DESIGN

We used finite element (FE) models of the proximal femur with an open growth plate to study whether mechanical factors could explain the development of cam-type deformity in adolescents. Four different loading conditions (representing different types of physical activities) and three different levels of growth plate extension towards the femoral neck were considered. Mechanical stimuli at the tissue level were calculated by means of the osteogenic index (OI) for all loading conditions and growth plate shape variations.

RESULTS

Loading conditions and growth plate shape influence the distribution of OI in hips with an open growth plate, thereby driving the development of cam-type deformity. In particular, specific types of loads experienced during physical activities and a larger growth plate extension towards the femoral neck increase the chance of cam-type deformity.

CONCLUSIONS

Specific loading patterns seem to stimulate the development of cam-type deformity by modifying the distribution of the mechanical stimulus. This is in line with recent clinical studies and reveals mechanobiological mechanisms that trigger the development of cam-type deformity. Avoiding these loading patterns during skeletal growth might be a potential preventative strategy for future hip OA.

Skeletal biology over the life span: a view from the surfaces

The biocultural interpretation of skeletal remains is based upon the foundation of skeletal biology. In this review we examine the current state of skeletal biology research outside of the mainstream anthropology literature. The focus is on the structural changes of bone development and growth, and modeling and repair in the four bone surfaces: periosteal, Haversian, endosteal, and trabecular. The pattern of skeletal changes is placed within the framework of the human life span. New perspectives and direction of research on the environmental, biological, and genetic influences on modeling and remodeling processes are discussed chronologically at each bone surface. Implications for biological anthropologists are considered. This approach emphasizes variation in skeletal biology as a dynamic record of development, maturity, and aging.

Fetal and postnatal mouse bone tissue contains more calcium than is present in hydroxyapatite

It has been shown for developing enamel and zebrafish fin that hydroxyapatite (HA) is preceded by an amorphous precursor, motivating us to examine the mineral development in mammalian bone, particularly femur and tibia of fetal and young mice. Mineral particle thickness and arrangement were characterized by (synchrotron) small-angle X-ray scattering (SAXS) combined with wide-angle X-ray diffraction (WAXD) and X-ray fluorescence (XRF) analysis. Simultaneous measurements of the local calcium content and the HA content via XRF and WAXD, respectively, revealed the total calcium contained in HA crystals. Interestingly, bones of fetal as well as newborn mice contained a certain fraction of calcium which is not part of the HA crystals. Mineral deposition could be first detected in fetal tibia at day 16.5 by environmental scanning electron microscopy (ESEM). SAXS revealed a complete lack of orientation in the mineral particles at this stage, whereas 1day after birth particles were predominantly aligned parallel to the longitudinal bone axis, with the highest degree of alignment in the midshaft. Moreover, we found that mineral particle length increased with age as well as the thickness, while fetal particles were thicker but much shorter. In summary, this study revealed strong differences in size and orientation of the mineral particles between fetal and postnatal bone, with bulkier, randomly oriented particles at the fetal stage, and highly aligned, much longer particles after birth. Moreover, a part of the calcium seems to be present in other form than HA at all stages of development.

Bone structural changes in osteoarthritis as a result of mechanoregulated bone adaptation: a modeling approach

OBJECTIVE

There are strong indications that subchondral bone may play an important role in osteoarthritis (OA), making it an interesting target for medical therapies. The subchondral bone structure changes markedly during OA, and it has long been assumed that this occurs secondary to cartilage degeneration. However, for various conditions that are associated with OA, it is known that they may also induce bone structural changes in the absence of cartilage degeneration. We therefore aimed to investigate if OA bone structural changes can result from mechanoregulated bone adaptation, independent of cartilage degeneration.

METHOD

With a bone adaptation model, we simulated various conditions associated with OA -without altering the articular cartilage- and we evaluated if mechanoregulated bone remodeling by itself could lead to OA-like bone structural changes.

RESULTS

For each of the conditions, the predicted changes in bone structural parameters (bone fraction, trabecular thickness, trabecular number, and trabecular separation) were similar to those observed in OA.

CONCLUSION

This indicates that bone adaptation in OA can be mechanoregulated with structural changes occurring independent of cartilage degeneration.

Biophysical stimuli induced by passive movements compensate for lack of skeletal muscle during embryonic skeletogenesis

In genetically modified mice with abnormal skeletal muscle development, bones and joints are differentially affected by the lack of skeletal muscle. We hypothesise that unequal levels of biophysical stimuli in the developing humerus and femur can explain the differential effects on these rudiments when muscle is absent. We find that the expression patterns of four mechanosensitive genes important for endochondral ossification are differentially affected in muscleless limb mutants, with more extreme changes in the expression in the humerus than in the femur. Using finite element analysis, we show that the biophysical stimuli induced by muscle forces are similar in the humerus and femur, implying that the removal of muscle contractile forces should, in theory, affect the rudiments equally. However, simulations in which a displacement was applied to the end of the limb, such as could be caused in muscleless mice by movements of the mother or normal littermates, predicted higher biophysical stimuli in the femur than in the humerus. Stimuli induced by limb movement were much higher than those induced by the direct application of muscle forces, and we propose that movements of limbs caused by muscle contractions, rather than the direct application of muscle forces, provide the main mechanical stimuli for normal skeletal development. In muscleless mice, passive movement induces unequal biophysical stimuli in the humerus and femur, providing an explanation for the differential effects seen in these mice. The significance of these results is that forces originating external to the embryo may contribute to the initiation and progression of skeletal development when muscle development is abnormal.

Three-dimensional engineered bone from bone marrow stromal cells and their autogenous extracellular matrix

Most bone tissue-engineering research uses porous three-dimensional (3D) scaffolds for cell seeding. In this work, scaffold-less 3D bone-like tissues were engineered from rat bone marrow stromal cells (BMSCs) and their autogenous extracellular matrix (ECM). The BMSCs were cultured on a 2D substrate in medium that induced osteogenic differentiation. After reaching confluence and producing a sufficient amount of their own ECM, the cells contracted their tissue monolayer around two constraint points, forming scaffold-less cylindrical engineered bone-like constructs (EBCs). The EBCs exhibited alizarin red staining for mineralization and alkaline phosphatase activity and contained type I collagen. The EBCs developed a periosteum characterized by fibroblasts and unmineralized collagen on the periphery of the construct. Tensile tests revealed that the EBCs in culture had a tangent modulus of 7.5 +/- 0.5 MPa at 7 days post-3D construct formation and 29 +/- 9 MPa at 6 weeks after construct formation. Implantation of the EBCs into rats 7 days after construct formation resulted in further bone development and vascularization. Tissue explants collected at 4 weeks contained all three cell types found in native bone: osteoblasts, osteocytes, and osteoclasts. The resulting engineered tissues are the first 3D bone tissues developed without the use of exogenous scaffolding.

Sensitivity of tissue differentiation and bone healing predictions to tissue properties

Computational models are employed as tools to investigate possible mechano-regulation pathways for tissue differentiation and bone healing. However, current models do not account for the uncertainty in input parameters, and often include assumptions about parameter values that are not yet established. The aim was to clarify the importance of the assumed tissue material properties in a computational model of tissue differentiation during bone healing. An established mechano-biological model was employed together with a statistical approach. The model included an adaptive 2D finite element model of a fractured long bone. Four outcome criteria were quantified: (1) ability to predict sequential healing events, (2) amount of bone formation at specific time points, (3) total time until healing, and (4) mechanical stability at specific time points. Statistical analysis based on fractional factorial designs first involved a screening experiment to identify the most significant tissue material properties. These seven properties were studied further with response surface methodology in a three-level Box-Behnken design. Generally, the sequential events were not significantly influenced by any properties, whereas rate-dependent outcome criteria and mechanical stability were significantly influenced by Young's modulus and permeability. Poisson's ratio and porosity had minor effects. The amount of bone formation at early, mid and late phases of healing, the time until complete healing and the mechanical stability were all mostly dependent on three material properties; permeability of granulation tissue, Young's modulus of cartilage and permeability of immature bone. The consistency between effects of the most influential parameters was high. To increase accuracy and predictive capacity of computational models of bone healing, the most influential tissue mechanical properties should be accurately quantified.

Mechanical modulation of osteochondroprogenitor cell fate

Mesenchymal cells are natural tissue builders. They exhibit an extraordinary capacity to metamorphize into differentiated cells, using extrinsic spatial and temporal inputs and intrinsic algorithms, as well as to build and adapt their own habitat. In addition to providing a habitat for osteoprogenitor cells, tissues of the skeletal system provide mechanical support and protection for the multiple organs of vertebrate organisms. This review examines the role of mechanics on determination of cell fate during pre-, peri- and postnatal development of the skeleton as well as during tissue genesis and repair in postnatal life. The role of cell mechanics is examined and brought into context of intrinsic cues during mesenchymal condensation. Remarkable new insights regarding structure function relationships in mesenchymal stem cells, and their influence on determination of cell fate are integrated in the context of de novo tissue generation and postnatal repair. Key differences in the formation of osteogenic and chondrogenic condensations are discussed in relation to direct intramembranous and indirect endochondral ossification. New approaches are discussed to elucidate and exploit extrinsic cues to generate tissues in the laboratory and in the clinic.

Determination of bone's mechanical matrix properties by nanoindentation

Osteoporosis is a devastating disease that is characterized not only by a reduction in bone quantity but also by deterioration in bone quality. The quality of bone tissue is greatly influenced by its mechanical properties and, therefore, investigations into the etiology and enhanced detection of osteoporosis, or the efficacy of interventions, may require the assessment of bone's mechanical properties at the level of the tissue. Nanoindentation is a relatively new technique that is capable of evaluating bone's quasi-static and dynamic mechanical properties on extremely small volumes of tissue. These data can be used directly to describe the pre-yield properties of the matrix, but can also be combined with imaging techniques and mechanical models to extrapolate the mechanical properties from the level of the tissue to that of the organ.

Accretion of bone quantity and quality in the developing mouse skeleton

In this work, we found that bone mineral formation proceeded very rapidly in mice by 1 day of age, where the degree of mineralization, the tissue mineral density, and the mineral crystallinity reached 36%, 51%, and 87% of the adult values, respectively. However, even though significant mineralization had occurred, the elastic modulus of 1-day-old bone was only 14% of its adult value, indicating that the intrinsic stiffening of the bone lags considerably behind the initial mineral formation.

INTRODUCTION

To meet the mechanical challenges during early development, the skeleton requires the rapid accretion of bone quality and bone quantity. Here, we describe early bone development in the mouse skeleton and test the hypothesis that specific compositional properties determine the stiffness of the tissue.

MATERIALS AND METHODS

Tibias of female BALB mice were harvested at eight time-points (n = 4 each) distributed between 1 and 40 days of age and subjected to morphometric (muCT), chemical (Fourier transform infrared microspectroscopy), and mechanical (nanoindentation) analyses. Tibias of 450-day-old mice served as fully mineralized control specimens.

RESULTS

Bone growth proceeded very rapidly; at 1 day of age, the degree of mineralization (phosphate/protein ratio), the density of mineralized bone (TMD), and mineral crystallinity had reached 36%, 51%, and 87% of the adult (450 days) values, respectively. Spatially, the variability in mineralization across the mid-diaphysis was very high for the early time-points and declined over time. In contrast to the notable changes in mineralization, carbonate substitution into the mineral lattice (carbonate/phosphate ratio) and collagen cross-linking did not show any significant changes over this time period. Even though significant mineralization had occurred, the elastic modulus of 1-day-old bone was only 14% of the adult value and increased to 89% (of its adult value) after 40 days. Between samples of different time-points, significant positive correlations were observed between the elastic modulus and TMD (r(2) = 0.84), phosphate/protein ratio (r(2) = 0.59), and crystallinity (r(2) = 0.23), whereas collagen cross-linking showed a small but significant negative correlation (r(2) = 0.15).

CONCLUSIONS

These data indicate that specific chemical and morphometric properties modulate bone's stiffness during early growth. The intrinsic stiffening of the bone, however, lags considerably behind the initial mineral formation, emphasizing the importance of bone mineral quality for optimizing matrix integrity.

Mechanobiology of embryonic limb development

Considerable evidence exists to support the hypothesis that mechanical forces have an essential role in healthy embryonic skeletal development. Clinical observations and experimental data indicate the importance of muscle contractions for limb development. However, the influence of these forces is seldom referred to in biological descriptions of bone development, and perhaps this is due to the fact that the hypothesis that mechanical forces are essential for normal embryonic skeletal development is difficult to test and elaborate experimentally in vivo, particularly in humans. Computational modeling has the potential to address this issue by simulating embryonic growth under a range of loading conditions but the potential of such models has yet to be fully exploited. In this article, we review the literature on mechanobiology of limb development in three main sections: (a) experimental alteration of the mechanical environment, (b) mechanical properties of embryonic tissues, and (c) the use of computational models. Then we analyze the main issues, and suggest how experimental and computational fields could work closer together to enhance our understanding of mechanobiology of the embryonic skeleton.

The PTHrP-Ihh feedback loop in the embryonic growth plate allows PTHrP to control hypertrophy and Ihh to regulate proliferation

Growth plate and long bone development is governed by biochemical signaling pathways of which the PTHrP-Ihh system is the best known. Other factors, such as BMPs, FGFs and mechanical loading, may interact with this system. This study aims at elucidating the relative importance of PTHrP and Ihh for controlling proliferation, and hypertrophy in fetal growth plate cartilage. We assessed the question why reduced Ihh expression leads to more pronounced effects on the number of non-hypertrophic cells and total bone formation, compared to PTHrP down-regulation. Using few basic equations, constituted from literature data, this paper shows how the PTHrP-Ihh feedback system can control different aspects of tissue differentiation at distinct locations. In particular, it is shown that (mechanical or biochemical) perturbations will affect proliferation via Ihh-related parameters, whereas changes in PTHrP-related parameters selectively interact with hypertrophy. This is contra-intuitive, since PTHrP acts to keep cells proliferating. In this context, the critical PTHrP level for keeping cells proliferating has been reconsidered. In addition, an explanation is provided for the aforementioned difference in effect between reduced Ihh and PTHrP expression.

Modulation of endochondral development of the distal femoral condyle by mechanical loading

Although previous theoretical modeling studies have predicted that various mechanical stresses accelerate or inhibit the ossification process of the neonatal chondroepiphysis, there is a paucity of experimental data to verify these models. The present study was designed to provide experimental evidence on whether the ossification of the chondroepiphysis is modulated by mechanical loading on the distal femoral condyle explant of the neonatal (5-day-old) rabbit in organ culture. Upon aseptic dissection, the right condyle explant was immersed in and fixated to an organ culture system, and received cyclic forces at 200 mN and 1 Hz for 12 h (N=8) directly on its slightly convex articular surface, whereas the contralateral, left condyle explant was immersed separately in organ culture (N=8). Subsequently, both loaded and control explants were placed in a bioreactor rotating at 20 rpm for 72 h. In each mechanically loaded specimen, a structure reminiscent of the secondary ossification center (SOC) appeared with an average area of 1.17+/-0.13 mm(2), or 15.2+/-8.2% of the total epiphysis area. In contrast, no SOC was detected in any of the unloaded contralateral control specimens. The SOC in mechanically loaded specimens was stained intensively with fast green, whereas either the rest of the loaded epiphysis or the entire control epiphysis was stained intensely to safranin-O but lacked fast green staining. Immunolocalization revealed that the SOC of the mechanically loaded specimens expressed Run x 2 and osteopontin, both of which were absent in the unloaded control specimens. Type X collagen was expressed surrounding hypertrophic chondrocytes adjacent to the SOC, but was absent in the control specimen. Type II collagen and decorin were absent in the SOC of the loaded specimen, but were expressed throughout the rest of the loaded epiphysis and the unloaded control epiphysis. The intensity of type II collagen and decorin expression was significantly stronger among hypertrophic chondrocytes surrounding the SOC than the control. The numbers of hypertrophic chondrocytes surrounding the SOC and superior to metaphyseal bone were significantly higher in the loaded specimens than the unloaded controls. Taken together, mechanical stresses accelerate the formation of the secondary ossification center, and therefore modulate endochondral ossification.

Rapid establishment of chemical and mechanical properties during lamellar bone formation

The development of prophylaxes and treatments of bone diseases that can effectively increase the strength of bone as a structure necessitates a better understanding of the time course by which chemical properties define the stiffness of the material during primary and secondary mineralization. It was hypothesized that these processes would be relatively slow in the actively growing skeleton. Seven-week-old Sprague-Dawley female rats (n = 8) were injected with multiple fluorochrome labels over a time span of 3 weeks and killed. Chemical and mechanical properties of the tibial mid-diaphysis were spatially characterized between the endocortical and periosteal surface by in situ infrared microspectroscopy and nanoindentation. The phosphate-to-protein ratio of bone 2-6 days old was 20% smaller at the periosteal surface and 22% smaller at the endocortical surface (P < 0.05 each) compared to older intracortical regions. The ratios of carbonate to protein, crystallinity, type A/type B carbonate, collagen cross-linking, and bone elastic modulus did not differ significantly between bone 2-6, 10-14, and 8-22 days old and intracortical regions. Intracortical properties of 10-week-old rats, except for the carbonate-to-protein ratio which was 23% smaller (P < 0.01), were not significantly different from intracortical matrix properties of young adult rats (5 months, n = 4). Spatially, the phosphate-to-protein ratio (R(2) = 0.33) and the phosphate-to-carbonate ratio (R(2) = 0.55) were significantly correlated with bone material stiffness, while the combination of all chemical parameters raised the R(2) value to 0.83. These data indicate that lamellar bone has the ability to quickly establish its mechanical and chemical tissue properties during primary and secondary mineralization even when the skeleton experiences rapid growth.