Differential gene expression from microarray analysis distinguishes woven and lamellar bone formation in the rat ulna following mechanical loading

Formation of woven and lamellar bone in the adult skeleton can be induced through mechanical loading. Although much is known about the morphological appearance and structural properties of the newly formed bone, the molecular responses to loading are still not well understood. The objective of our study was to use a microarray to distinguish the molecular responses between woven and lamellar bone formation induced through mechanical loading. Rat forelimb loading was completed in a single bout to induce the formation of woven bone (WBF loading) or lamellar bone (LBF loading). A set of normal (non-loaded) rats were used as controls. Microarrays were performed at three timepoints after loading: 1 hr, 1 day and 3 days. Confirmation of microarray results was done for a select group of genes using quantitative real-time PCR (qRT-PCR). The micorarray identified numerous genes and pathways that were differentially regulated for woven, but not lamellar bone formation. Few changes in gene expression were evident comparing lamellar bone formation to normal controls. A total of 395 genes were differentially expressed between formation of woven and lamellar bone 1 hr after loading, while 5883 and 5974 genes were differentially expressed on days 1 and 3, respectively. Results suggest that not only are the levels of expression different for each type of bone formation, but that distinct pathways are activated only for woven bone formation. A strong early inflammatory response preceded an increase in angiogenic and osteogenic gene expression for woven bone formation. Furthermore, at later timepoints there was evidence of bone resorption after WBF loading. In summary, the vast coverage of the microarray offers a comprehensive characterization of the early differences in expression between woven and lamellar bone formation.     

The benefit of combining non-mechanical agents with mechanical loading: a perspective based on the Utah Paradigm of Skeletal Physiology

The Utah Paradigm of Skeletal Physiology with its key component, the mechanostat hypothesis, suggest plausible explanations of some of the tissue-level changes occurring from combining selected non-mechanical agents (anabolic and anti-resorptive/( re)modeling agents) with mechanical loading (osteogenic exercise) to increase bone mass and strength. The evidence for combining selected anabolic agents like parathyroid hormone, prostaglandin E(2), growth hormone, etc. with mechanical loading can increase bone mass is strong. Anabolic agents influence loading-related bone formation changes in a permissive manner and modulate (increase) the responsiveness of bone tissue to mechanical loading by changing thresholds for bone formation and resorption. However, any beneficial effect of combining selected anti-resorptive/(re)modeling agents like estrogen with loading is marginal, especially in adult skeletons. Postulated changes in modeling and remodeling thresholds (set points) and known direct effects on bone cells by non-mechanical agents may explain the observed tissue-level changes associated with large and minor increases in bone mass. Although the pharmaceutical industry has avoided considering osteogenic loading in the treatment of osteoporosis, a methodical dose-response study of anabolic agents combined with loading should: (1) provide opportunities for therapeutic intervention to imitate or enhance the osteogenic response to loading in order to correct osteopenias; (2) provide the potential to diminish the dosage of drugs required to induce bone formation in ways that enhanced efficacy and reduced any side effects; and (3) improve the quality of life and reduce the risk of falls by improving balance, gait speed and muscle strength with a non-mechanical agent like GH that could improve both muscle and bone mass and strength. Lastly, more studies are needed which determine bone strength instead of only "mass" in aged skeletons so one can assess how effective such treatments would reduce the risk of fracture in the clinic.     

Long-term potentiation in bone – a role for glutamate in strain-induced cellular memory?

Background

The adaptive response of bone cells to mechanical strain is a primary determinant of skeletal architecture and bone mass. In vivo mechanical loading induces new bone formation and increases bone mineral density whereas disuse, immobilisation and weightlessness induce bone loss. The potency of mechanical strain is such that a single brief period of loading at physiological strain magnitude is able to induce a long-lasting osteogenic response that lasts for days. Although the process of mechanotransduction in bone is incompletely understood, observations that responses to mechanical strain outlast the duration of stimulation necessitate the existence of a form of cellular memory through which transient strain episodes are recorded, interpreted and remembered by bone cells. Recent evidence supports the existence of a complex multicellular glutamate-signalling network in bone that shares functional similarities to glutamatergic neurotransmission in the central nervous system. In neurones, these signalling molecules coordinate synaptic communication required to support learning and memory formation, through a complex process of long-term potentiation.

Presentation of the hypothesis

We hypothesise that osteoblasts use a cellular mechanism similar or identical to neuronal long-term potentiation in the central nervous system to mediate long-lasting changes in osteogenesis following brief periods of mechanical strain.

Testing the hypothesis

N-methyl-D-aspartate (NMDA) receptor antagonism should inhibit the saturating response of mechanical strain and reduce the enhanced osteogenicity of segregated loading to that of an equivalent period of uninterrupted loading. Changes in α-amino-3-hydroxy-5-methyl-isoxazole propionate (AMPA) receptor expression, localisation and electrophysiological responses should be induced by mechanical strain and inhibited by modulators of neuronal long-term potentiation.

Implications of the hypothesis

If true, this hypothesis would provide a mechanism through which the skeleton could be pharmacologically primed to enhance or retrieve the normal osteogenic response to exercise. This would form a basis through which novel therapies could be developed to target osteoporosis and other prevalent bone disorders associated with low bone mass.

     

Effects of short-term recovery periods on fluid-induced signaling in osteoblastic cells

It is well known that cyclic mechanical loading can produce an anabolic response in bone. In vivo studies have shown that the insertion of short-term recovery periods (10-15 s) into mechanical loading profiles led to an increased osteogenic response compared to continuous cyclic loading of bone. Although this is suggestive of temporal processing at the bone cell level, there is little evidence to support such a hypothesis. Therefore, the current study investigated the cellular mechanism of bone's response to rest inserted vs. continuous mechanical loading. Cell responses to rest inserted mechanical loading were quantified by applying oscillatory fluid flow (OFF) to osteoblastic cells and quantifying real-time intracellular calcium [Ca2+]i, prostaglandin E2 (PGE2) release, and osteopontin (OPN) mRNA levels. Cells were exposed to OFF (1 Hz) at shear stresses of 1 and 2 Pa with rest periods of 5, 10, and 15s inserted every 10 loading cycles. The insertion of 10 and 15s rest periods into the flow profile resulted in multiple [Ca2+]i responses by individual cells, increased [Ca2+]i response magnitudes, and increased overall percent of cells responding compared to continuously loaded control groups. We determined the source of the multiple calcium responses to be from intracellular stores. In addition, rest inserted OFF led to similar levels of PGE2 release and increased levels of relative OPN mRNA compared to cells exposed to continuous OFF. Our study suggests that the cellular mechanism of bone adaptation to rest inserted mechanical loading may involve modulation of intracellular levels of calcium (frequency, magnitude, percent of cells responding).     

In silico study of bone tissue regeneration in an idealised porous hydrogel scaffold using a mechano-regulation algorithm

Mechanical stimulation, in the form of fluid perfusion or mechanical strain, enhances osteogenic differentiation and overall bone tissue formation by mesenchymal stems cells cultured in biomaterial scaffolds for tissue engineering applications. In silico techniques can be used to predict the mechanical environment within biomaterial scaffolds, and also the relationship between bone tissue regeneration and mechanical stimulation, and thereby inform conditions for bone tissue engineering experiments. In this study, we investigated bone tissue regeneration in an idealised hydrogel scaffold using a mechano-regulation model capable of predicting tissue differentiation, and specifically compared five loading cases, based on known experimental bioreactor regimes. These models predicted that low levels of mechanical loading, i.e. compression (0.5% strain), pore pressure of 10 kPa and a combination of compression (0.5%) and pore pressure (10 kPa), could induce more osteogenic differentiation and lead to the formation of a higher bone tissue fraction. In contrast greater volumes of cartilage and fibrous tissue fractions were predicted under higher levels of mechanical loading (i.e. compression strain of 5.0% and pore pressure of 100 kPa). The findings in this study may provide important information regarding the appropriate mechanical stimulation for in vitro bone tissue engineering experiments.     

Fluid flow in the osteocyte mechanical environment: a fluid–structure interaction approach

Osteocytes are believed to be the primary sensor of mechanical stimuli in bone, which orchestrate osteoblasts and osteoclasts to adapt bone structure and composition to meet physiological loading demands. Experimental studies to quantify the mechanical environment surrounding bone cells are challenging, and as such, computational and theoretical approaches have modelled either the solid or fluid environment of osteocytes to predict how these cells are stimulated in vivo. Osteocytes are an elastic cellular structure that deforms in response to the external fluid flow imposed by mechanical loading. This represents a most challenging multi-physics problem in which fluid and solid domains interact, and as such, no previous study has accounted for this complex behaviour. The objective of this study is to employ fluid–structure interaction (FSI) modelling to investigate the complex mechanical environment of osteocytes in vivo. Fluorescent staining of osteocytes was performed in order to visualise their native environment and develop geometrically accurate models of the osteocyte in vivo. By simulating loading levels representative of vigorous physiological activity ( $3,000\,\upmu \upvarepsilon $ compression and 300 Pa pressure gradient), we predict average interstitial fluid velocities $(\sim 60.5\,\upmu \text{ m/s })$ and average maximum shear stresses $(\sim 11\, \text{ Pa })$ surrounding osteocytes in vivo. Interestingly, these values occur in the canaliculi around the osteocyte cell processes and are within the range of stimuli known to stimulate osteogenic responses by osteoblastic cells in vitro. Significantly our results suggest that the greatest mechanical stimulation of the osteocyte occurs in the cell processes, which, cell culture studies have indicated, is the most mechanosensitive area of the cell. These are the first computational FSI models to simulate the complex multi-physics mechanical environment of osteocyte in vivo and provide a deeper understanding of bone mechanobiology.     

Computer modelling of bone’s adaptation: the role of normal strain,shear strain and fluid flow

Bone loss is a serious health problem. In vivo studies have found that mechanical stimulation may inhibit bone loss as elevated strain in bone induces osteogenesis, i.e. new bone formation. However, the exact relationship between mechanical environment and osteogenesis is less clear. Normal strain is considered as a prime stimulus of osteogenic activity; however, there are some instances in the literature where osteogenesis is observed in the vicinity of minimal normal strain, specifically near the neutral axis of bending in long bones. It suggests that osteogenesis may also be induced by other or secondary components of mechanical environment such as shear strain or canalicular fluid flow. As it is evident from the literature, shear strain and fluid flow can be potent stimuli of osteogenesis. This study presents a computational model to investigate the roles of these stimuli in bone adaptation. The model assumes that bone formation rate is roughly proportional to the normal, shear and fluid shear strain energy density above their osteogenic thresholds. In vivo osteogenesis due to cyclic cantilever bending of a murine tibia has been simulated. The model predicts results close to experimental findings when normal strain, and shear strain or fluid shear were combined. This study also gives a new perspective on the relation between osteogenic potential of micro-level fluid shear and that of macro-level bending shear. Attempts to establish such relations among the components of mechanical environment and corresponding osteogenesis may ultimately aid in the development of effective approaches to mitigating bone loss.     

Simulation of bone remodelling in orthodontic treatment

Orthodontic treatments not only displace irregular teeth but also induce responses in surrounding bone tissues. Bone remodelling is regarded as the regulatory mechanism triggered by mechanical loading. This study was aimed at investigating the effect of orthodontic loading on both tooth movement and neighbouring bone density distribution. A set of computational algorithms incorporating both external and internal remodelling mechanisms was implemented into a patient-specific 3D finite element (FE) model to investigate and analyse orthodontic treatment under four typical modes of orthodontic loading. The consequence of orthodontic treatment was reproduced numerically by using this FE-based technique. The results indicated that the diverse modes of orthodontic loading would result in different magnitudes of tooth movement and particular morphology of bone density distribution. It is illuminated that the newly developed algorithms may replicate the clinical situation more closely compared with the previous proposed method.     

A role for the primary cilium in paracrine signaling between mechanically stimulated osteocytes and mesenchymal stem cells

Bone turnover is a mechanically regulated process, coordinated in part by the network of mechanosensitive osteocytes residing within the tissue. The recruitment and bone forming activity of the mesenchymal derived osteoblast is determined by numerous factors including mechanical loading. It is therefore somewhat surprising that although mechanically regulated signaling between the coordinating osteocytes and mesenchymal stem cells (MSCs) should exist, to date it has not been directly demonstrated. In this study, conditioned media from mechanically stimulated osteocytes (MLO-Y4 cell line) was collected and added to MSCs (C3H10T1/2 cell line). The addition of mechanically stimulated osteocyte conditioned media resulted in a significant upregulation of the osteogenic genes OPN and COX-2 in MSCs compared to statically cultured conditioned media, demonstrating a novel paracrine signaling mechanism between the two cell types. The same mechanically conditioned media did not alter gene expression in osteoblasts (MC3T3 cell line), and mechanically stimulated osteoblast conditioned media did not alter gene expression in MSCs demonstrating that this signaling is unique to osteocytes and MSCs. Finally, the upregulation in osteogenic genes in MSCs was not observed if primary cilia formation was inhibited prior to mechanical stimulation of the osteocyte. In summary, the results of this study indicate that soluble factors secreted by osteocytes in response to mechanical stimulation can enhance osteogenic gene expression in MSCs demonstrating a novel, unique signaling mechanism and introduces a role for the primary cilium in flow mediated paracrine signaling in bone thereby highlighting the cilium as a potential target for therapeutics aimed at enhancing bone formation.     

The effects of misalignment during in vivo loading of bone: Techniques to detect the proximity of objects in three-dimensional models

Theories of mechanical adaptation of bone suggest that mechanical loading causes bone formation at discrete locations within bone microstructure experiencing the greatest mechanical stress/strain. Experimental testing of such theories requires in vivo loading experiments and high-resolution finite element models to determine the distribution of mechanical stresses. Finite element models of in vivo loading experiments typically assume idealized boundary conditions with applied load perfectly oriented on the bone, however small misalignments in load orientation during an in vivo experiment are unavoidable, and potentially confound the ability of finite element models to predict locations of bone formation at the scale of micrometers. Here we demonstrate two different three-dimensional spatial correlation methods to determine the effects of misalignment in load orientation on the locations of high mechanical stress/strain in the rodent tail loading model. We find that, in cancellous bone, the locations of tissue with high stress are maintained under reasonable misalignments in load orientation (p<0.01). In cortical bone, however, angular misalignments in the dorsal direction can alter the locations of high mechanical stress, but the locations of tissue with high stress are maintained under other misalignments (p<0.01). We conclude that, when using finite element models of the rodent tail loading model, small misalignments in loading orientation do not affect the predicted locations of high mechanical stress within cancellous bone.     

Novel oxysterols have pro-osteogenic and anti-adipogenic effects in vitro and induce spinal fusion in vivo

Stimulation of bone formation by osteoinductive materials is of great clinical importance in spinal fusion surgery, repair of bone fractures, and in the treatment of osteoporosis. We previously reported that specific naturally occurring oxysterols including 20(S)-hydroxycholesterol (20S) induce the osteogenic differentiation of pluripotent mesenchymal cells, while inhibiting their adipogenic differentiation. Here we report the characterization of two structural analogues of 20S, Oxy34 and Oxy49, which induce the osteogenic and inhibit the adipogenic differentiation of bone marrow stromal cells (MSC) through activation of Hedgehog (Hh) signaling. Treatment of M2-10B4 MSC with Oxy34 or Oxy49 induced the expression of osteogenic differentiation markers Runx2, Osterix (Osx), alkaline phosphatase (ALP), bone sialoprotein (BSP), and osteocalcin (OCN), as well as ALP enzymatic activity and robust mineralization. Treatment with oxysterols together with PPARγ activator, troglitazone (Tro), inhibited mRNA expression for adipogenic genes PPARγ, LPL, and aP2, and inhibited the formation of adipocytes. Efficacy of Oxy34 and Oxy49 in stimulating bone formation in vivo was assessed using the posterolateral intertransverse process rat spinal fusion model. Rats receiving collagen implants with Oxy 34 or Oxy49 showed comparable osteogenic efficacy to BMP2/collagen implants as measured by radiography, MicroCT, and manual inspection. Histological analysis showed trabecular and cortical bone formation by oxysterols and rhBMP2 within the fusion mass, with robust adipogenesis in BMP2-induced bone and significantly less adipocytes in oxysterol-induced bone. These data suggest that Oxy34 and Oxy49 are effective novel osteoinductive molecules and may be suitable candidates for further development and use in orthopedic indications requiring local bone formation.     

Effect of cell seeding and mechanical loading on vascularization and tissue formation inside a scaffold: A mechano-biological model using a lattice approach to simulate cell activity

Achieving successful vascularization remains one of the main problems in bone tissue engineering. After scaffold implantation, the growth of capillaries into the porous construct may be too slow to provide adequate nutrients to the cells in the scaffold interior and this inhibits tissue formation in the scaffold core. Often, prior to implantation, a controlled cell culture environment is used to stimulate cell proliferation and, once in place, the mechanical environment acting on the tissue construct is determined by the loading conditions at the implantation site. To what extent do cell seeding conditions and the construct loading environment have an effect on scaffold vascularization and tissue growth? In this study, a mechano-biological model for tissue differentiation and blood vessel growth was used to determine the influence of cell seeding on vascular network development and tissue growth inside a regular-structured bone scaffold under different loading conditions. It is predicted that increasing the number of cells seeded homogeneously reduces the rate of vascularization and the maximum penetration of the vascular network, which in turn reduces bone tissue formation. The seeding of cells in the periphery of the scaffold was predicted to be beneficial for vascularization and therefore for bone growth; however, tissue formation occurred more slowly during the first weeks after implantation compared to homogeneous seeding. Low levels of mechanical loading stimulated bone formation while high levels of loading inhibited bone formation and capillary growth. This study demonstrates the feasibility of computational design approaches for bone tissue engineering.     

Increased osteogenic response to exercise in metaphyseal versus diaphyseal cortical bone

Recent experimental data suggest that the anabolic response of bone to changes in physical activity and mechanical loading may vary among different skeletal elements, and even within different regions of the same bone. In order to better understand site-specific variation in bone modeling we used an experimental protocol in which locomotor activity was increased in laboratory mice with regular treadmill exercise for only 30 min/day. We predicted that the regular muscle contractions that occur during exercise would significantly increase cortical bone formation in these animals, and that the increase in cortical bone mass would vary between metaphyseal and diaphyseal regions. Cortical bone mass, density, and bone geometry were compared between these two regions using pQCT technology. Results indicate that exercise increases bone mineral content (BMC) in the mid-diaphysis by approximately 20%, whereas bone mass in the metaphyseal region is increased by approximately 35%. Endosteal and periosteal circumference at the midshaft are increased with exercise, whereas increased periosteal circumference is accompanied by marked endosteal contraction at the metaphysis, resulting in an increase in cortical area of more than 50%. These findings suggest that the osteogenic response of cortical bone to exercise varies significantly along the length of a bone, and more distal regions appear most likely to exhibit morphologic changes when loading conditions are altered.     

Evolution of mechanoregulation of bone growth will lead to non-optimal bone phenotypes

Mechanical forces acting on the bones during growth affect their final shape and strength. Mechanoregulation of bone growth may be recognized in embryogenesis, and also in the adaptation of the adult skeleton to changes in mechanical loading. Mechanoregulatory responses for tissues have arisen during evolution, but does evolution give rise to responses that produce optimal skeletal phenotypes? In this paper, we investigate the emergence of an optimal mechanoregulation response in a population. By combining equations describing long bone growth with a genetic algorithm to describe evolutionary change, we created a computational model to simulate the evolution of mechanoregulation in bone growth. A population of individuals is created where each individual is assigned a diploid gene set which controls the growth and remodelling of the bone. At maturity, each bone is assessed and its 'fitness' calculated; fitness is quantified as bone strength relative to bone mass. The simulation continues for many generations, and includes mutations and a varying environment. The genes present in the population are tracked and the evolution of parameters governing mechanoregulation is calculated. The results indicate that a population may converge to one bone growth algorithm but, more usually, a range of mechanoregulation algorithms for different individuals will persist after many generations. Even if the population converges to one mechanoregulation law, convergence to the 'optimum' bone was never found. Although many researchers propose that natural selection has pushed skeletal structure towards an optimum, our computational model suggests that this is unlikely to be the case.