Transcriptional profiling of cortical versus cancellous bone from mechanically-loaded murine tibiae reveals differential gene expression

Bioenergetic programs of cancellous and cortical bone are distinct and differ with age and mechanical loading

Mechanical loading induces bone formation in young rodents, but mechanoresponsiveness is reduced with age. Glycolytic activity and mitochondrial dysfunction increase with age and may change bone mechanotransduction. To evaluate load-induced changes to bioenergetic activity in young and adult animals, we loaded the tibia of 10-wk and 26-wk female C57BL/6J mice and examined transcriptomic responses at the mid-diaphysis, and metaphyseal cortical shell and cancellous core. Across all biological processes, oxidative phosphorylation and mitochondrial pathways were most often enriched with loading and had contrasting enrichment in young and adult animals. Following loading, young animals had temporally-coordinated differential expression of mitochondrial-associated genes, with greatest expression at the mid-diaphysis. In adults, bioenergetic gene expression was lower compared to young animals. To assess individual contributions of glycolysis and pyruvate-mediated oxidative phosphorylation to load-induced bone formation in vivo, we inhibited each pathway therapeutically and loaded the tibia of young and adult female mice for 2 weeks. In both young and adult mice, loading increased cortical bone mass, but inhibition of oxidative phosphorylation reduced cortical area and moment of inertia in both loaded and control limbs. Conversely, load-induced improvements of adult cancellous bone depended on glycolysis. In summary, mechanical loading transcriptionally activated mitochondrial pathways in an age-specific manner and bioenergetic inhibition revealed unique metabolic programs for cortical and cancellous bone.

3D spatial distribution of Sost mRNA and sclerostin protein expression in response to in vivo tibia loading in female mice

Bones adapt to external mechanical loads through a process known as mechanoadaptation. Osteocytes are the bone cells that sense the mechanical environment and initiate a biological response. Investigating the changes in osteocyte molecular expression following mechanical loading has been instrumental in characterizing the regulatory pathways involved in bone adaptation. However, current methods for examining osteocyte molecular expression do not preserve the three-dimensional structure of the bone, which plays a critical role in the mechanical stimuli sensed by the osteocytes and their spatially controlled biological responses. In this study, we used WISH-BONE (Whole-mount In Situ Histology of Bone) to investigate the spatial distribution of Sost-mRNA transcripts and its encoded protein, sclerostin, in 3D mouse tibia midshaft following in vivo tibia loading. Our findings showed a decrease in the percentage of Sost-positive osteocytes, after loading, along the posterior-lateral side of the tibia. The number of sclerostin-positive osteocytes were found to significantly decrease at a very specific 2D location of the tibia after loading. However, in 3D, the total number of sclerostin-positive osteocytes was similar between loaded and control legs. This work is the first to provide a 3D analysis of Sost and sclerostin distribution in loaded versus contralateral mouse tibia midshafts. It also highlights the importance of the bone region analyzed and the method utilized when interpreting mechanoadaptation results. WISH-BONE represents a powerful tool for further characterization of mechanosensitive genes regulation in bone and holds the potential for advancing the development of new treatments targeting mechanosensitivity-related bone disorders.

Spatial transcriptomics in bone mechanomics: Exploring the mechanoregulation of fracture healing in the era of spatial omics

In recent decades, the field of bone mechanobiology has sought experimental techniques to unravel the molecular mechanisms governing the phenomenon of mechanically regulated fracture healing. Each cell within a fracture site resides within different local microenvironments characterized by different levels of mechanical strain; thus, preserving the spatial location of each cell is critical in relating cellular responses to mechanical stimuli. Our spatial transcriptomics-based "mechanomics" platform facilitates spatially resolved analysis of the molecular profiles of cells with respect to their local in vivo mechanical environment by integrating time-lapsed in vivo micro-computed tomography, spatial transcriptomics, and micro-finite element analysis. We investigate the transcriptomic responses of cells as a function of the local strain magnitude by identifying the differential expression of genes in regions of high and low strain within a fracture site. Our platform thus has the potential to address fundamental open questions within the field and to discover mechano-responsive targets to enhance fracture healing.

A proteomics approach to study mouse long bones: examining baseline differences and mechanical loading-induced bone formation in young-adult and old mice

With aging, bone mass declines and the anabolic effects of skeletal loading diminish. While much research has focused on gene transcription, how bone ages and loses its mechanoresponsiveness at the protein level remains unclear. We developed a novel proteomics approach and performed a paired mass spectrometry and RNA-seq analysis on tibias from young-adult (5-month) and old (22-month) mice. We report the first correlation estimate between the bone proteome and transcriptome (Spearman ρ = 0.40), which is in line with other tissues but indicates that a relatively low amount of variation in protein levels is explained by the variation in transcript levels. Of 71 shared targets that differed with age, eight were associated with bone mineral density in previous GWAS, including understudied targets Asrgl1 and Timp2. We used complementary RNA in situ hybridization to confirm that Asrgl1 and Timp2 had reduced expression in osteoblasts/osteocytes in old bones. We also found evidence for reduced TGF-beta signaling with aging, in particular Tgfb2. Next, we defined proteomic changes following mechanical loading. At the protein level, bone differed more with age than with loading, and aged bone had fewer loading-induced changes. Overall, our findings underscore the need for complementary protein-level assays in skeletal biology research.

Piezo1 expression in mature osteocytes is dispensable for the skeletal response to mechanical loading

Mechanical loading is essential for bone growth and homeostasis, with osteocytes considered the primary mechanosensors. Deleting the mechanosensitive ion channel Piezo1 from Dmp1-Cre-targeted cells, which include both osteoblasts and osteocytes, resulted in reduced bone mass and impaired skeletal responses to mechanical stimuli. To specifically isolate Piezo1's role in osteocytes, we used Sost-Cre mice to generate mice with Piezo1 deletion exclusively in mature osteocytes. These mice exhibited lower bone mineral density, decreased cancellous bone mass, and reduced cortical thickness with decrease periosteal expansion. However, unlike mice lacking Piezo1 in both osteoblasts and osteocytes, those with Piezo1 deletion in mature osteocytes responded normally to mechanical loading. These findings suggest that Piezo1 expression in mature osteocytes contributes to bone maintenance under normal physiological condition, but is dispensable for the skeletal response to mechanical loading.

Cortical and Trabecular Bone Modeling and Implications for Bone Functional Adaptation in the Mammalian Tibia

Bone modeling involves the addition of bone material through osteoblast-mediated deposition or the removal of bone material via osteoclast-mediated resorption in response to perceived changes in loads by osteocytes. This process is characterized by the independent occurrence of deposition and resorption, which can take place simultaneously at different locations within the bone due to variations in stress levels across its different regions. The principle of bone functional adaptation states that cortical and trabecular bone tissues will respond to mechanical stimuli by adjusting (i.e., bone modeling) their morphology and architecture to mechanically improve their mechanical function in line with the habitual in vivo loading direction. This principle is relevant to various research areas, such as the development of improved orthopedic implants, preventative medicine for osteopenic elderly patients, and the investigation of locomotion behavior in extinct species. In the present review, the mammalian tibia is used as an example to explore cortical and trabecular bone modeling and to examine its implications for the functional adaptation of bones. Following a short introduction and an exposition on characteristics of mechanical stimuli that influence bone modeling, a detailed critical appraisal of the literature on cortical and trabecular bone modeling and bone functional adaptation is given. By synthesizing key findings from studies involving small mammals (rodents), large mammals, and humans, it is shown that examining both cortical and trabecular bone structures is essential for understanding bone functional adaptation. A combined approach can provide a more comprehensive understanding of this significant physiological phenomenon, as each structure contributes uniquely to the phenomenon.

PTH treatment before cyclic joint loading improves cartilage health and attenuates load-induced osteoarthritis development in mice

Osteoarthritis (OA) treatment is limited by the lack of effective nonsurgical interventions to slow disease progression. Here, we examined the contributions of the subchondral bone properties to OA development. We used parathyroid hormone (PTH) to modulate bone mass before OA initiation and alendronate (ALN) to inhibit bone remodeling during OA progression. We examined the spatiotemporal progression of joint damage by combining histopathological and transcriptomic analyses across joint tissues. The additive effect of PTH pretreatment before OA initiation and ALN treatment during OA progression most effectively attenuated load-induced OA pathology. Individually, PTH directly improved cartilage health and slowed the development of cartilage damage, whereas ALN primarily attenuated subchondral bone changes associated with OA progression. Joint damage reflected early transcriptomic changes. With both treatments, the structural changes were associated with early modulation of immunoregulation and immunoresponse pathways that may contribute to disease mechanisms. Overall, our results demonstrate the potential of subchondral bone-modifying therapies to slow the progression of OA.

Controlled mechanical loading affects the osteocyte transcriptome in porcine trabecular bone in situ

INTRODUCTION

Osteocytes modulate bone adaptation in response to mechanical stimuli imparted by the deforming bone tissue in which they are encased by communicating with osteoclasts and osteoblasts as well as other osteocytes in the lacuna-canalicular network through secreted cytokines and chemokines. Understanding the transcriptional response of osteocytes to mechanical stimulation in situ could identify new targets to inhibit bone loss or enhance bone formation in the presence of diseases like osteoporosis or metastatic cancer. We compared the mechanically regulated transcriptional response of osteocytes in trabecular bone following one or three days of controlled mechanical loading.

METHODS

Porcine trabecular bone explants were cultured in a bioreactor for 48 h and subsequently loaded twice a day for one day or 3 days. RNA was isolated and sequenced, and the Tuxedo suite was used to identify differentially expressed genes and pathway analysis was conducted using Ingenuity Pathway Analysis (IPA).

RESULTS

There were about 4000 differentially expressed genes following in situ culture relative to fresh bone. One hundred six genes were differentially expressed between the loaded and non-loaded groups following one day of loading compared to 913 genes after 3 d of loading. Only 45 of these were coincident between the two time points, indicating an evolving transcriptome. Clustering and principal component analysis indicated differences between the loaded and non-loaded groups after 3 d of loading.

DISCUSSION

With sustained loading, there was a nine-fold increase in the number of differentially expressed genes, suggesting that osteocytes respond to loading through sequential activation of downstream genes in the same pathways. The differentially expressed genes were related to osteoarthritis, osteocyte, and chondrocyte signaling pathways. We noted that NFkB and TNF signaling are affected by early loading and this may drive downstream effects on the mechanobiological response. Moreover, these genes may regulate catabolic effects of mechanical disuse through their actions on pre-osteoclasts in the bone marrow niche.

Transcriptome profile reveals novel candidate genes associated with bone strength in end-of-lay hens

Bone weakness causes many problems such as osteoporosis, bone fractures, and economic loss, especially at the late stage of lay, in laying hen production. However, the genetic factors and molecular mechanism affecting the bone strength is still largely unknown. To elucidate the molecular mechanism and genetic factors affecting bone strength, a total of six cDNA libraries were constructed and used to compare genetic differences between tibia with higher（Group HBS）and lower（Group LBS）breaking strength in Hyline grey layers. A comparison between Groups HBS and LBS revealed nine differentially expressed genes, of which five were upregulated and four were downregulated in the LBS relative to the HBS in tibia. Our results showed novel candidate genes concerned with bone strength in the late laying period. These include transcription factor paired box protein Pax-5 (Pax5), tissue inhibitor of Metallopoteinase-4 (TIMP4), Kelch-like protein 14 (KLHL14), predicted MAGUK p55 subfamily member 7 isoform X4 (MPP7) and Osteoclast-associated Ig-like receptor (OSCAR). Our data provide a vital resource for discovering important candidate genes associated with bone strength and will help further study the molecular mechanisms for bone remodeling.

The transcriptional state and chromatin landscape of cichlid jaw shape variation across species and environments

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, but how they interact remains poorly understood. Here, we utilize the cichlid oral jaw apparatus to better understand these gene-by-environment effects. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. We used two Lake Malawi species adapted to different foraging habitats along the pelagic-benthic ecomorphological axis. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when they were challenged with a pelagic versus benthic foraging mode. Expression data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible. This reduced dataset implicated notable candidate genes including the Hedgehog effector, KIAA0586 and the ETS transcription factor, etv4, which connects environmental stress and craniofacial morphogenesis. Taken together, these data provide novel insights into the epigenetic, genetic and cellular bases of species- and environment-specific bone shapes.

Influence of Scaffold Microarchitecture on Angiogenesis and Regulation of Cell Differentiation during the Early Phase of Bone Healing: A Transcriptomics and Histological Analysis

The early phase of bone healing is a complex and poorly understood process. With additive manufacturing, we can generate a specific and customizable library of bone substitutes to explore this phase. In this study, we produced tricalcium phosphate-based scaffolds with microarchitectures composed of filaments of 0.50 mm in diameter, named Fil050G, and 1.25 mm named Fil125G, respectively. The implants were removed after only 10 days in vivo followed by RNA sequencing (RNAseq) and histological analysis. RNAseq results revealed upregulation of adaptive immune response, regulation of cell adhesion, and cell migration-related genes in both of our two constructs. However, significant overexpression of genes linked to angiogenesis, regulation of cell differentiation, ossification, and bone development was observed solely in Fil050G scaffolds. Moreover, quantitative immunohistochemistry of structures positive for laminin revealed a significantly higher number of blood vessels in Fil050G samples. Furthermore, µCT detected a higher amount of mineralized tissue in Fil050G samples suggesting a superior osteoconductive potential. Hence, different filament diameters and distances in bone substitutes significantly influence angiogenesis and regulation of cell differentiation involved in the early phase of bone regeneration, which precedes osteoconductivity and bony bridging seen in later phases and as consequence, impacts the overall clinical outcome.

The Role of Bone Cell Energetics in Altering Bone Quality and Strength in Health and Disease

PURPOSE OF REVIEW

Bone quality and strength are diminished with age and disease but can be improved by clinical intervention. Energetic pathways are essential for cellular function and drive osteogenic signaling within bone cells. Altered bone quality is associated with changes in the energetic activity of bone cells following diet-based or therapeutic interventions. Energetic pathways may directly or indirectly contribute to changes in bone quality. The goal of this review is to highlight tissue-level and bioenergetic changes in bone health and disease.

RECENT FINDINGS

Bone cell energetics are an expanding field of research. Early literature primarily focused on defining energetic activation throughout the lifespan of bone cells. Recent studies have begun to connect bone energetic activity to health and disease. In this review, we highlight bone cell energetic demands, the effect of substrate availability on bone quality, altered bioenergetics associated with disease treatment and development, and additional biological factors influencing bone cell energetics. Bone cells use several energetic pathways during differentiation and maturity. The orchestration of bioenergetic pathways is critical for healthy cell function. Systemic changes in substrate availability alter bone quality, potentially due to the direct effects of altered bone cell bioenergetic activity. Bone cell bioenergetics may also contribute directly to the development and treatment of skeletal diseases. Understanding the role of energetic pathways in the cellular response to disease will improve patient treatment.

Osteocyte Estrogen Receptor β (Ot-ERβ) Regulates Bone Turnover and Skeletal Adaptive Response to Mechanical Loading Differently in Male and Female Growing and Adult Mice

Age-related bone loss is a failure of balanced bone turnover and diminished skeletal mechanoadaptation. Estrogen receptors, ERα and ERβ, play critical roles in osteoprotective regulation activated by estrogen and mechanical signals. Previous studies mainly focused on ERα and showed that osteocyte-ERα (Ot-ERα) regulated trabecular, but not cortical bone, and played a minor role in load-induced cortical adaptation. However, the role of Ot-ERβ in bone mass regulation remains unrevealed. To address this issue, we characterized bone (re)modeling and gene expression in male and female mice with Ot-ERβ deletion (ERβ-dOT) and littermate control (LC) at 10 weeks (young) or 28 weeks (adult) of age, as well as their responses to in vivo tibial compressive loading. Increased cancellous bone mass appeared in the L4 vertebral body of young male ERβ-dOT mice. At the same time, femoral cortical bone gene expression showed signs consistent with elevated osteoblast and osteoclast activities (type-I collagen, Cat K, RANKL). Upregulated androgen receptor (AR) expression was observed in young male ERβ-dOT mice relative to LC, suggesting a compensatory effect of testosterone on male bone protection. In contrast, bone mass in L4 decreased in adult male ERβ-dOT mice, attributed to potentially increased bone resorption activity (Cat K) with no change in bone formation. There was no effect of ERβ-dOT on bone mass or gene expression in female mice. Sex-dependent regulation of Ot-ERβ also appeared in load-induced cortical responsiveness. Young female ERβ-dOT mice showed an enhanced tibial cortical anabolic adaptation compared with LC. In contrast, an attenuated cortical anabolic response presented at the proximal tibia in male ERβ-dOT mice at both ages. For the first time, our findings suggest that Ot-ERβ regulates bone (re)modeling and the response to mechanical signals through different mechanisms in males and females. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Protein Kinase G2 Is Essential for Skeletal Homeostasis and Adaptation to Mechanical Loading in Male but Not Female Mice

We previously showed that the NO/cGMP/protein kinase G (PKG) signaling pathway positively regulates osteoblast proliferation, differentiation, and survival in vitro, and that cGMP-elevating agents have bone-anabolic effects in mice. Here, we generated mice with an osteoblast-specific (OB) knockout (KO) of type 2 PKG (gene name Prkg2) using a Col1a1(2.3 kb)-Cre driver. Compared to wild type (WT) littermates, 8-week-old male OB Prkg2-KO mice had fewer osteoblasts, reduced bone formation rates, and lower trabecular and cortical bone volumes. Female OB Prkg2-KO littermates showed no bone abnormalities, despite the same degree of PKG2 deficiency in bone. Expression of osteoblast differentiation- and Wnt/β-catenin-related genes was lower in primary osteoblasts and bones of male KO but not female KO mice compared to WT littermates. Osteoclast parameters were unaffected in both sexes. Since PKG2 is part of a mechano-sensitive complex in osteoblast membranes, we examined its role during mechanical loading. Cyclical compression of the tibia increased cortical thickness and induced mechanosensitive and Wnt/β-catenin-related genes to a similar extent in male and female WT mice and female OB Prkg2-KO mice, but loading had a minimal effect in male KO mice. We conclude that PKG2 drives bone acquisition and adaptation to mechanical loading via the Wnt/β-catenin pathway in male mice. The striking sexual dimorphism of OB Prkg2-KO mice suggests that current U.S. Food and Drug Administration-approved cGMP-elevating agents may represent novel effective treatment options for male osteoporosis. © 2022 American Society for Bone and Mineral Research (ASBMR).

From mesenchymal niches to engineered in vitro model systems: Exploring and exploiting biomechanical regulation of vertebrate hedgehog signalling

Tissue patterning is the result of complex interactions between transcriptional programs and various mechanical cues that modulate cell behaviour and drive morphogenesis. Vertebrate Hedgehog signalling plays key roles in embryogenesis and adult tissue homeostasis, and is central to skeletal development and the osteogenic differentiation of mesenchymal stem cells. The expression of several components of the Hedgehog signalling pathway have been reported to be mechanically regulated in mesodermal tissue patterning and osteogenic differentiation in response to external stimulation. Since a number of bone developmental defects and skeletal diseases, such as osteoporosis, are directly linked to aberrant Hedgehog signalling, a better knowledge of the regulation of Hedgehog signalling in the mechanosensitive bone marrow-residing mesenchymal stromal cells will present novel avenues for modelling these diseases and uncover novel opportunities for extracellular matrix-targeted therapies. In this review, we present a brief overview of the key molecular players involved in Hedgehog signalling and the basic concepts of mechanobiology, with a focus on bone development and regeneration. We also highlight the correlation between the activation of the Hedgehog signalling pathway in response to mechanical cues and osteogenesis in bone marrow-derived mesenchymal stromal cells. Finally, we propose different tissue engineering strategies to apply the expanding knowledge of 3D material-cell interactions in the modulation of Hedgehog signalling in vitro for fundamental and translational research applications.

Molecular Identification of Spatially Distinct Anabolic Responses to Mechanical Loading in Murine Cortical Bone

Osteoporosis affects over 200 million women worldwide, one-third of whom are predicted to suffer from an osteoporotic fracture in their lifetime. The most promising anabolic drugs involve administration of expensive antibodies. Because mechanical loading stimulates bone formation, our current data, using a mouse model, replicates the anabolic effects of loading in humans and may identify novel pathways amenable to oral treatment. Murine tibial compression produces axially varying deformations along the cortical bone, inducing highest strains at the mid-diaphysis and lowest at the metaphyseal shell. To test the hypothesis that load-induced transcriptomic responses at different axial locations of cortical bone would vary as a function of strain magnitude, we loaded the left tibias of 10-week-old female C57Bl/6 mice in vivo in compression, with contralateral limbs as controls. Animals were euthanized at 1, 3, or 24 hours post-loading or loaded for 1 week (n = 4-5/group). Bone marrow and cancellous bone were removed, cortical bone was segmented into the metaphyseal shell, proximal diaphysis, and mid-diaphysis, and load-induced differential gene expression and enriched biological processes were examined for the three segments. At each time point, the mid-diaphysis (highest strain) had the greatest transcriptomic response. Similarly, biological processes regulating bone formation and turnover increased earlier and to the greatest extent at the mid-diaphysis. Higher strain induced greater levels of osteoblast and osteocyte genes, whereas expression was lower in osteoclasts. Among the top differentially expressed genes at 24-hours post-loading, 17 had known functions in bone biology, of which 12 were present only in osteoblasts, 3 exclusively in osteoclasts, and 2 were present in both cell types. Based on these results, we conclude that murine tibial loading induces spatially unique transcriptomic responses correlating with strain magnitude in cortical bone. © 2022 American Society for Bone and Mineral Research (ASBMR).

Loading-induced bone formation is mediated by Wnt1 induction in osteoblast-lineage cells

Mechanical loading on the skeleton stimulates bone formation. Although the exact mechanism underlying this process remains unknown, a growing body of evidence indicates that the Wnt signaling pathway is necessary for the skeletal response to loading. Recently, we showed that Wnts produced by osteoblast lineage cells mediate the osteo-anabolic response to tibial loading in adult mice. Here, we report that Wnt1 specifically plays a crucial role in mediating the mechano-adaptive response to loading. Independent of loading, short-term loss of Wnt1 in the Osx-lineage resulted in a decreased cortical bone area in the tibias of 5-month-old mice. In females, strain-matched loading enhanced periosteal bone formation in Wnt1F/F controls, but not in Wnt1F/F; OsxCreERT2 knockouts. In males, strain-matched loading increased periosteal bone formation in both control and knockout mice; however, the periosteal relative bone formation rate was 65% lower in Wnt1 knockouts versus controls. Together, these findings show that Wnt1 supports adult bone homeostasis and mediates the bone anabolic response to mechanical loading.

Proteomic profiling of human bone from different anatomical sites - A pilot study

PURPOSE

The study aim is a comparative proteome-based analysis of different autologous bone entities (alveolar bone [AB], iliac cortical [IC] bone, and iliac spongiosa [IS]) used for alveolar onlay grafting.

EXPERIMENTAL DESIGN

Site-matched bone samples of AB, IC, and IS were harvested during alveolar onlay grafting. Proteins were extracted using a detergent-based (sodium dodecyl sulfate) strategy and trypsinized. Proteome analysis was performed using liquid chromatography-tandem mass spectrometry (LC-MS/MS). MaxQuant was used for peptide-to-spectrum matching, peak detection, and quantitation. Linear models for microarray analysis (LIMMA) were used to detect differentially abundant peptides and proteins.

RESULTS

A total of 1730 different proteins were identified across the 15 samples at a false discovery rate of 1%. Partial least-squares discriminant analysis approved segregation of AB, IC, and IS protein profiles. LIMMA statistics highlighted 66 proteins that were more abundant in AB then in IC (vs. 92 proteins were enriched in IC over AB). Gene Ontology enrichment analysis revealed a matrisomal versus an immune-related proteome fingerprint in AB versus IC.

CONCLUSION AND CLINICAL RELEVANCE

This pilot study demonstrates an ECM protein-related proteome fingerprint in AB and an immune-related proteome fingerprint in IS and IC.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.

Osteoblast-Specific Wnt Secretion Is Required for Skeletal Homeostasis and Loading-Induced Bone Formation in Adult Mice

Wnt signaling is critical to many aspects of skeletal regulation, but the importance of Wnt ligands in the bone anabolic response to mechanical loading is not well established. Recent transcriptome profiling studies by our laboratory and others show that mechanical loading potently induces genes encoding Wnt ligands, including Wnt1 and Wnt7b. Based on these findings, we hypothesized that mechanical loading stimulates adult bone formation by inducing Wnt ligand expression. To test this hypothesis, we inhibited Wnt ligand secretion in adult (5 months old) mice using a systemic (drug) and a bone-targeted (conditional gene knockout) approach, and subjected them to axial tibial loading to induce lamellar bone formation. Mice treated with the Wnt secretion inhibitor WNT974 exhibited a decrease in bone formation in non-loaded bones as well as a 54% decline in the periosteal bone formation response to tibial loading. Next, osteoblast-specific Wnt secretion was inhibited by dosing 5-month-old Osx-CreERT2; WlsF/F mice with tamoxifen. Within 1 to 2 weeks of Wls deletion, skeletal homeostasis was altered with decreased bone formation and increased resorption, and the anabolic response to loading was reduced 65% compared to control (WlsF/F ). Together, these findings show that Wnt ligand secretion is required for adult bone homeostasis, and furthermore establish a role for osteoblast-derived Wnts in mediating the bone anabolic response to tibial loading. © 2021 American Society for Bone and Mineral Research (ASBMR).

Bone adaptation to mechanical loading in mice is affected by circadian rhythms

Physical forces are critical for successful function of many organs including bone. Interestingly, the timing of exercise during the day alters physiology and gene expression in many organs due to circadian rhythms. Circadian clocks in tissues, such as bone, express circadian clock genes that target tissue-specific genes, resulting in tissue-specific rhythmic gene expression (clock-controlled genes). We hypothesized that the adaptive response of bone to mechanical loading is regulated by circadian rhythms. First, mice were sham loaded and sacrificed 8 h later, which amounted to tissues being collected at zeitgeber time (ZT)2, 6, 10, 14, 18, and 22. Cortical bone of the tibiae collected from these mice displayed diurnal expression of core clock genes and key osteocyte and osteoblast-related genes, such as the Wnt-signaling inhibitors Sost and Dkk1, indicating these are clock-controlled genes. Serum bone turnover markers did not display rhythmicity. Second, mice underwent a single bout of in vivo loading at either ZT2 or ZT14 and were sacrificed 1, 8, or 24 h after loading. Loading at ZT2 resulted in Sost upregulation, while loading at ZT14 led to Sost and Dkk1 downregulation. Third, mice underwent daily in vivo tibial loading over 2 weeks administered either in the morning, (ZT2, resting phase) or evening (ZT14, active phase). In vivo microCT was performed at days 0, 5, 10, and 15 and conventional histomorphometry was performed at day 15. All outcome measures indicated a robust response to loading, but only microCT-based time-lapse morphometry showed that loading at ZT14 resulted in a greater endocortical bone formation response compared to mice loaded at ZT2. The decreased Sost and Dkk1 expression coincident with the modest, but significant time-of-day specific increase in adaptive bone formation, suggests that circadian clocks influence bone mechanoresponse.

Osteocytic Pericellular Matrix (PCM): Accelerated Degradation under In Vivo Loading and Unloading Conditions Using a Novel Imaging Approach

The proteoglycan-containing pericellular matrix (PCM) controls both the biophysical and biochemical microenvironment of osteocytes, which are the most abundant cells embedded and dispersed in bones. As a molecular sieve, osteocytic PCMs not only regulate mass transport to and from osteocytes but also act as sensors of external mechanical environments. The turnover of osteocytic PCM remains largely unknown due to technical challenges. Here, we report a novel imaging technique based on metabolic labeling and “click-chemistry,” which labels de novo PCM as “halos” surrounding osteocytes in vitro and in vivo. We then tested the method and showed different labeling patterns in young vs. old bones. Further “pulse-chase” experiments revealed dramatic difference in the “half-life” of PCM of cultured osteocytes (~70 h) and that of osteocytes in vivo (~75 d). When mice were subjected to either 3-week hindlimb unloading or 7-week tibial loading (5.1 N, 4 Hz, 3 d/week), PCM half-life was shortened (~20 d) and degradation accelerated. Matrix metallopeptidase MMP-14 was elevated in mechanically loaded osteocytes, which may contribute to PCM degradation. This study provides a detailed procedure that enables semi-quantitative study of the osteocytic PCM remodeling in vivo and in vitro.

The Osteocyte Transcriptome: Discovering Messages Buried Within Bone

PURPOSE OF THE REVIEW

Osteocytes are cells embedded within the bone matrix, but their function and specific patterns of gene expression remain only partially defined; this is beginning to change with recent studies using transcriptomics. This unbiased approach can generate large amounts of data and is now being used to identify novel genes and signalling pathways within osteocytes both at baseline conditions and in response to stimuli. This review outlines the methods used to isolate cell populations containing osteocytes, and key recent transcriptomic studies that used osteocyte-containing preparations from bone tissue.

RECENT FINDINGS

Three common methods are used to prepare samples to examine osteocyte gene expression: digestion followed by sorting, laser capture microscopy, and the isolation of cortical bone shafts. All these methods present challenges in interpreting the data generated. Genes previously not known to be expressed by osteocytes have been identified and variations in osteocyte gene expression have been reported with age, sex, anatomical location, mechanical loading, and defects in bone strength. A substantial proportion of newly identified transcripts in osteocytes remain functionally undefined but several have been cross-referenced with functional data. Future work and improved methods (e.g. scRNAseq) likely provide useful resources for the study of osteocytes and important new information on the identity and functions of this unique cell type within the skeleton.

The WNT1G177C mutation specifically affects skeletal integrity in a mouse model of osteogenesis imperfecta type XV

The recent identification of homozygous WNT1 mutations in individuals with osteogenesis imperfecta type XV (OI-XV) has suggested that WNT1 is a key ligand promoting the differentiation and function of bone-forming osteoblasts. Although such an influence was supported by subsequent studies, a mouse model of OI-XV remained to be established. Therefore, we introduced a previously identified disease-causing mutation (G177C) into the murine Wnt1 gene. Homozygous Wnt1^G177C/G177C mice were viable and did not display defects in brain development, but the majority of 24-week-old Wnt1^G177C/G177C mice had skeletal fractures. This increased bone fragility was not fully explained by reduced bone mass but also by impaired bone matrix quality. Importantly, the homozygous presence of the G177C mutation did not interfere with the osteoanabolic influence of either parathyroid hormone injection or activating mutation of LRP5, the latter mimicking the effect of sclerostin neutralization. Finally, transcriptomic analyses revealed that short-term administration of WNT1 to osteogenic cells induced not only the expression of canonical WNT signaling targets but also the expression of genes encoding extracellular matrix modifiers. Taken together, our data demonstrate that regulating bone matrix quality is a primary function of WNT1. They further suggest that individuals with WNT1 mutations should profit from existing osteoanabolic therapies.

RNAseq and RNA molecular barcoding reveal differential gene expression in cortical bone following hindlimb unloading in female mice

Disuse-induced bone loss is seen following spinal cord injury, prolonged bed rest, and exposure to microgravity. We performed whole transcriptomic profiling of cortical bone using RNA sequencing (RNAseq) and RNA molecular barcoding (NanoString) on a hindlimb unloading (HLU) mouse model to identify genes whose mRNA transcript abundances change in response to disuse. Eleven-week old female C57BL/6 mice were exposed to ambulatory loading or HLU for 7 days (n = 8/group). Total RNA from marrow-flushed femoral cortical bone was analyzed on HiSeq and NanoString platforms. The expression of several previously reported genes associated with Wnt signaling and metabolism was altered by HLU. Furthermore, the increased abundance of transcripts, such as Pfkfb3 and Mss51, after HLU imply these genes also have roles in the cortical bone’s response to altered mechanical loading. Our study demonstrates that an unbiased approach to assess the whole transcriptomic profile of cortical bone can reveal previously unidentified mechanosensitive genes and may eventually lead to novel targets to prevent disuse-induced osteoporosis.

Mechanosensitive miR-100 coordinates TGFβ and Wnt signaling in osteocytes during fluid shear stress

Organism scale mechanical forces elicit cellular scale changes through coordinated regulation of multiple signaling pathways. The mechanisms by which cells integrate signaling to generate a unified biological response remains a major question in mechanobiology. For example, the mechanosensitive response of bone and other tissues requires coordinated signaling by the transforming growth factor beta (TGFβ) and Wnt pathways through mechanisms that are not well‐defined. Here we report a new microRNA‐dependent mechanism that mediates mechanosensitive crosstalk between TGFβ and Wnt signaling in osteocytes exposed to fluid shear stress (FSS). From 60 mechanosensitive microRNA (miRs) identified by small‐RNAseq, miR100 expression is suppressed by in vivo hindlimb loading in the murine tibia and by cellular scale FSS in OCY454 cells. Though FSS activates both TGFβ and Wnt signaling in osteocytes, only TGFβ represses miR‐100 expression. miR‐100, in turn, antagonizes Wnt signaling by targeting and inhibiting expression of Frizzled receptors (FZD5/FZD8). Accordingly, miR‐100 inhibition blunts FSS‐ and TGFβ‐inducible Wnt signaling. Therefore, our results identify FSS‐responsive miRNAs in osteocytes, including one that integrates the mechanosensitive function of two essential signaling pathways in the osteoanabolic response of bone to mechanical load.

Gene expression of intracortical bone demonstrates loading-induced increases in Wnt1 and Ngf and inhibition of bone remodeling processes

Osteocytes are the primary mechanosensitive cells in bone. However, their location in mineralized matrix has limited the in vivo study of osteocytic genes induced by mechanical loading. Laser Capture Microdissection (LCM) allows isolation of intracortical bone (Intra-CB), enriched for osteocytes, from bone tissue for gene expression analysis. We used microarray to analyze gene expression from mouse tibial Intra-CB dissected using LCM 4 h after a single loading bout or after 5 days of loading. Osteocyte enrichment was supported by greater expression of Sost, Dmp1, Dkk1, and Mepe in Intra-CB regions vs. Mixed regions containing periosteum and muscle (fold-change (FC) = 3.4, 2.2, 5.1, 3.0, respectively). Over 150 differentially expressed genes (DEGs) due to loading (loaded vs. contralateral control) in Intra-CB were found on Day 1 and Day 5, but only 10 genes were differentially expressed on both days, including Ngf (Day 1 FC = 13.5, Day 5 FC = 11.1) and Wnt1 (Day 1 FC = 1.5, Day 5 FC = 5.1). The expression of Ngf and Wnt1 within Intra-CB was confirmed by in situ hybridization, and a significant increase in number of Wnt1 mRNA molecules occurred on day 1. We also found changes in extracellular matrix remodeling with Timp1 (FC = 3.1) increased on day 1 and MMP13 (FC = 0.3) decreased on day 5. Supporting this result, IHC for osteocytic MMP13 demonstrated a marginal decrease due to loading on day 5. Gene Ontology (GO) biological processes for loading DEGs indicated regulation of vasculature, neuronal and immune processes while cell-type specific gene lists suggested regulation of osteoclast, osteoblast, and endothelial related genes. In summary, microarray analysis of microdissected Intra-CB revealed differential regulation of Ngf, Wnt1, and MMP13 due to loading in osteocytes.

The Cellular Choreography of Osteoblast Angiotropism in Bone Development and Homeostasis

Interaction between endothelial cells and osteoblasts is essential for bone development and homeostasis. This process is mediated in large part by osteoblast angiotropism, the migration of osteoblasts alongside blood vessels, which is crucial for the homing of osteoblasts to sites of bone formation during embryogenesis and in mature bones during remodeling and repair. Specialized bone endothelial cells that form "type H" capillaries have emerged as key interaction partners of osteoblasts, regulating osteoblast differentiation and maturation and ensuring their migration towards newly forming trabecular bone areas. Recent revolutions in high-resolution imaging methodologies for bone as well as single cell and RNA sequencing technologies have enabled the identification of some of the signaling pathways and molecular interactions that underpin this regulatory relationship. Similarly, the intercellular cross talk between endothelial cells and entombed osteocytes that is essential for bone formation, repair, and maintenance are beginning to be uncovered. This is a relatively new area of research that has, until recently, been hampered by a lack of appropriate analysis tools. Now that these tools are available, greater understanding of the molecular relationships between these key cell types is expected to facilitate identification of new drug targets for diseases of bone formation and remodeling.

The TrkA agonist gambogic amide augments skeletal adaptation to mechanical loading

The periosteal and endosteal surfaces of mature bone are densely innervated by sensory nerves expressing TrkA, the high-affinity receptor for nerve growth factor (NGF). In previous work, we demonstrated that administration of exogenous NGF significantly increased load-induced bone formation through the activation of Wnt signaling. However, the translational potential of NGF is limited by the induction of substantial mechanical and thermal hyperalgesia in mice and humans. Here, we tested the effect of gambogic amide (GA), a recently identified robust small molecule agonist for TrkA, on hyperalgesia and load-induced bone formation. Behavioral analysis was used to assess pain up to one week after axial forelimb compression. Contrary to our expectations, GA treatment was not associated with diminished use of the loaded forelimb or sensitivity to thermal stimulus. Furthermore, dynamic histomorphometry revealed a significant increase in relative periosteal bone formation rate as compared to vehicle treatment. Additionally, we found that GA treatment was associated with an increase in the number of osteoblasts per bone surface in loaded limbs as well as a significant increase in the fold change of Ngf, Wnt7b, and Axin2 mRNA expression as compared to vehicle (control). To test the effect of GA on osteoblasts directly, we cultured MC3T3-E1 cells for up to 21 days in osteogenic differentiation media containing NGF, GA, or vehicle (control). Media containing GA induced the significant upregulation of the osteoblastic differentiation markers Runx2, Bglap2, and Sp7 in a dose-dependent manner, whereas treatment with NGF was not associated with any significant increases in these markers. Furthermore, consistent with our in vivo findings, we observed that administration of 50 nM of GA upregulated expression of Ngf at both Day 3 and Day 7. However, cells treated with the highest dose of GA (500 nM) had significantly increased apoptosis and impaired cell proliferation. In conclusion, our study indicates GA may be useful for augmenting skeletal adaptation to mechanical forces without inducing hyperalgesia.

A collagen extraction and deuterium oxide stable isotope tracer method for the quantification of bone collagen synthesis rates in vivo

The development of safe and practical strategies to prevent weakening of bone tissue is vital, yet attempts to achieve this have been hindered by a lack of understanding of the short-term (days-weeks) physiology of bone collagen turnover. To address this, we have developed a method to quantify bone collagen synthesis in vivo, using deuterium oxide (D₂ O) tracer incorporation techniques combined with gas chromatography pyrolysis isotope-ratio mass spectrometry (GC-pyrolysis-IRMS). Forty-six male and female rats from a selectively bred model ingested D₂ O for 3 weeks. Femur diaphyses (FEM), tibia proximal (T-PRO), and distal (T-DIS) epiphyses-metaphyses and tibia mid-shaft diaphyses (T-MID) were obtained from all rats after necropsy. After demineralisation, collagen proteins were isolated and hydrolysed and collagen fractional synthetic rates (FSRs) determined by incorporation of deuterium into protein-bound alanine via GC-pyrolysis-IRMS. The collagen FSR for the FEM (0.131 ± 0.078%/day; 95% CI [0.106-0.156]) was greater than the FSR at T-MID (0.055 ± 0.049%/day; 95% CI [0.040-0.070]; p < 0.001). The T-PRO site had the highest FSR (0.203 ± 0.123%/day; 95% CI [0.166-0.241]) and T-DIS the lowest (0.027 ± 0.015%/day; 95% CI [0.022-0.031]). The three tibial sites exhibited different FSRs (p < 0.001). Herein, we have developed a sensitive method to quantify in vivo bone collagen synthesis and identified site-specific rates of synthesis, which could be applicable to studies of human bone collagen turnover.

Cancellous Bone May Have a Greater Adaptive Strain Threshold Than Cortical Bone

Strain magnitude has a controlling influence on bone adaptive response. However, questions remain as to how and if cancellous and cortical bone tissues respond differently to varied strain magnitudes, particularly at a molecular level. The goal of this study was to characterize the time‐dependent gene expression, bone formation, and structural response of the cancellous and cortical bone of female C57Bl/6 mice to mechanical loading by applying varying load levels (low: −3.5 N; medium: −5.2 N; high: −7 N) to the skeleton using a mouse tibia loading model. The loading experiment showed that cortical bone mass at the tibial midshaft was significantly enhanced following all load levels examined and bone formation activities were particularly elevated at the medium and high loads applied. In contrast, for the proximal metaphyseal cancellous bone, only the high load led to significant increases in bone mass and bone formation indices. Similarly, expression of genes associated with inhibition of bone formation (e.g., Sost) was altered in the diaphyseal cortical bone at all load levels, but in the metaphyseal cortico‐cancellous bone only by the high load. Finite element analysis determined that the peak tensile or compressive strains that were osteogenic for the proximal cancellous bone under the high load were significantly greater than those that were osteogenic for the midshaft cortical tissues under the low load. These results suggest that the magnitude of the strain stimulus regulating structural, cellular, and molecular responses of bone to loading may be greater for the cancellous tissues than for the cortical tissues. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

New Advances in Osteocyte Mechanotransduction

PURPOSE OF REVIEW

Skeletal adaptation to mechanical loading plays a critical role in bone growth and the maintenance of bone homeostasis. Osteocytes are postulated to serve as a hub orchestrating bone remodeling. The recent findings on the molecular mechanisms by which osteocytes sense mechanical loads and the downstream bone-forming factors are reviewed.

RECENT FINDINGS

Calcium channels have been implicated in mechanotransduction in bone cells for a long time. Efforts have been made to identify a specific calcium channel mediating the skeletal response to mechanical loads. Recent studies have revealed that Piezo1, a mechanosensitive ion channel, is critical for normal bone growth and is essential for the skeletal response to mechanical loading. Identification of mechanosensors and their downstream effectors in mechanosensing bone cells is essential for new strategies to modulate regenerative responses and develop therapies to treat the bone loss related to disuse or advanced age.

Hypergravity as a gravitational therapy mitigates the effects of knee osteoarthritis on the musculoskeletal system in a murine model

Insights into the effects of osteoarthritis (OA) and physical interventions on the musculoskeletal system are limited. Our goal was to analyze musculoskeletal changes in OA mice and test the efficacy of 8-week exposure to hypergravity, as a replacement of physical activity. 16-week-old male (C57BL/6J) mice allocated to sham control and OA groups not centrifuged (Ctrl 1g and OA 1g, respectively) or centrifuged at 2g acceleration (Ctrl 2g and OA 2g). OA 1g displayed decreased trabecular bone in the proximal tibia metaphysis and increased osteoclastic activity and local TNFα gene expression, all entirely prevented by 2g gravitational therapy. However, while cortical bone of tibia midshaft was preserved in OA 1g (vs. ctrl), it is thinner in OA 2g (vs. OA 1g). In the hind limb, OA at 1g increased fibers with lipid droplets by 48% in the tibialis anterior, a fact fully prevented by 2g. In Ctrl, 2g increased soleus, tibialis anterior and gastrocnemius masses. In the soleus of both Ctrl and OA, 2g induced larger fibers and a switch from type-II to type-I fiber. Catabolic (myostatin and its receptor activin RIIb and visfatine) and anabolic (FNDC5) genes dramatically increased in Ctrl 2g and OA 2g (p<0.01 vs 1g). Nevertheless, the overexpression of FNDC5 (and follistatine) was smaller in OA 2g than in Ctrl 2g. Thus, hypergravity in OA mice produced positive effects for trabecular bone and muscle typology, similar to resistance exercises, but negative effects for cortical bone.

Old Mice Have Less Transcriptional Activation But Similar Periosteal Cell Proliferation Compared to Young-Adult Mice in Response to in vivo Mechanical Loading

Mechanical loading is a potent strategy to induce bone formation, but with aging, the bone formation response to the same mechanical stimulus diminishes. Our main objectives were to (i) discover the potential transcriptional differences and (ii) compare the periosteal cell proliferation between tibias of young-adult and old mice in response to strain-matched mechanical loading. First, to discover potential age-related transcriptional differences, we performed RNA sequencing (RNA-seq) to compare the loading responses between tibias of young-adult (5-month) and old (22-month) C57BL/6N female mice following 1, 3, or 5 days of axial loading (loaded versus non-loaded). Compared to young-adult mice, old mice had less transcriptional activation following loading at each time point, as measured by the number of differentially expressed genes (DEGs) and the fold-changes of the DEGs. Old mice engaged fewer pathways and gene ontology (GO) processes, showing less activation of processes related to proliferation and differentiation. In tibias of young-adult mice, we observed prominent Wnt signaling, extracellular matrix (ECM), and neuronal responses, which were diminished with aging. Additionally, we identified several targets that may be effective in restoring the mechanoresponsiveness of aged bone, including nerve growth factor (NGF), Notum, prostaglandin signaling, Nell-1, and the AP-1 family. Second, to directly test the extent to which periosteal cell proliferation was diminished in old mice, we used bromodeoxyuridine (BrdU) in a separate cohort of mice to label cells that divided during the 5-day loading interval. Young-adult and old mice had an average of 15.5 and 16.7 BrdU+ surface cells/mm, respectively, suggesting that impaired proliferation in the first 5 days of loading does not explain the diminished bone formation response with aging. We conclude that old mice have diminished transcriptional activation following mechanical loading, but periosteal proliferation in the first 5 days of loading does not differ between tibias of young-adult and old mice. © 2020 American Society for Bone and Mineral Research.

Murine Axial Compression Tibial Loading Model to Study Bone Mechanobiology: Implementing the Model and Reporting Results

In vivo, tibial loading in mice is increasingly used to study bone adaptation and mechanotransduction. To achieve standardized and defined experimental conditions, loading parameters and animal-related factors must be considered when performing in vivo loading studies. In this review, we discuss these loading and animal-related experimental conditions, present methods to assess bone adaptation, and suggest reporting guidelines. This review originated from presentations by each of the authors at the workshop "Developing Best Practices for Mouse Models of In Vivo Loading" during the Preclinical Models Section at the Orthopaedic Research Society Annual Meeting, San Diego, CA, March 2017. Following the meeting, the authors engaged in detailed discussions with consideration of relevant literature. The guidelines and recommendations in this review are provided to help researchers perform in vivo loading experiments in mice, and thus further our knowledge of bone adaptation and the mechanisms involved in mechanotransduction. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:233-252, 2020.

Perlecan/Hspg2 deficiency impairs bone's calcium signaling and associated transcriptome in response to mechanical loading

Perlecan, a heparan sulfate proteoglycan, acts as a mechanical sensor for bone to detect external loading. Deficiency of perlecan increases the risk of osteoporosis in patients with Schwartz-Jampel Syndrome (SJS) and attenuates loading-induced bone formation in perlecan deficient mice (Hypo). Considering that intracellular calcium [Ca²⁺]_i is an ubiquitous messenger controlling numerous cellular processes including mechanotransduction, we hypothesized that perlecan deficiency impairs bone's calcium signaling in response to loading. To test this, we performed real-time [Ca²⁺]_i imaging on in situ osteocytes of adult murine tibiae under cyclic loading (8N). Relative to wild type (WT), Hypo osteocytes showed decreases in the overall [Ca²⁺]_i response rate (-58%), calcium peaks (-33%), cells with multiple peaks (-53%), peak magnitude (-6.8%), and recovery speed to baseline (-23%). RNA sequencing and pathway analysis of tibiae from mice subjected to one or seven days of unilateral loading demonstrated that perlecan deficiency significantly suppressed the calcium signaling, ECM-receptor interaction, and focal adhesion pathways following repetitive loading. Defects in the endoplasmic reticulum (ER) calcium cycling regulators such as Ryr1/ryanodine receptors and Atp2a1/Serca1 calcium pumps were identified in Hypo bones. Taken together, impaired calcium signaling may contribute to bone's reduced anabolic response to loading, underlying the osteoporosis risk for the SJS patients.

Conditional Activation of NF-κB Inducing Kinase (NIK) in the Osteolineage Enhances Both Basal and Loading-Induced Bone Formation

Studies from global loss-of-function mutants suggest that alternative NF-κB downstream of NF-κB inducing kinase (NIK) is a cell-intrinsic negative regulator of osteogenesis. However, the interpretation of the osteoblast and/or osteocyte contribution to the bone phenotype is complicated by simultaneous osteoclast defects in these models. Therefore, we turned to a transgenic mouse model to investigate the direct role of NIK in the osteolineage. Osx-Cre;NT3 animals (NT3-Cre +), which bear a constitutively active NIK allele (NT3) driven by Osx-Cre, were compared with their Cre-negative, Control (Ctrl) littermates. NT3-Cre + mice had elevated serum P1NP and CTX levels. Despite this high turnover state, µCT showed that constitutive activation of NIK resulted in a net increase in basal bone mass in both cortical and cancellous compartments. Furthermore, NT3-Cre + mice exhibited a greater anabolic response following mechanical loading compared with controls. We next performed RNA-Seq on nonloaded and loaded tibias to elucidate possible mechanisms underlying the increased bone anabolism seen in NT3-Cre + mice. Hierarchical clustering revealed two main transcriptional programs: one loading-responsive and the other NT3 transgene-driven. Gene ontology (GO) analysis indicated a distinct upregulation of receptor, kinase, and growth factor activities including Wnts, as well as a calcium-response signature in NT3-Cre + limbs. The promoters of these GO-term associated genes, including many known to be bone-anabolic, were highly enriched for multiple κB recognition elements (κB-RE) relative to the background frequency in the genome. The loading response in NT3-Cre + mice substantially overlapped (>90%) with Ctrl. Surprisingly, control animals had 10-fold more DEGs in response to loading. However, most top DEGs shared between genotypes had a high incidence of multiple κB-RE in their promoters. Therefore, both transcriptional programs (loading-responsive and NT3 transgene-driven) are modulated by NF-κB. Our studies uncover a previously unrecognized role for NF-κB in the promotion of both basal and mechanically stimulated bone formation. © 2019 American Society for Bone and Mineral Research.

Stimulation of Piezo1 by mechanical signals promotes bone anabolism

Mechanical loading, such as caused by exercise, stimulates bone formation by osteoblasts and increases bone strength, but the mechanisms are poorly understood. Osteocytes reside in bone matrix, sense changes in mechanical load, and produce signals that alter bone formation by osteoblasts. We report that the ion channel Piezo1 is required for changes in gene expression induced by fluid shear stress in cultured osteocytes and stimulation of Piezo1 by a small molecule agonist is sufficient to replicate the effects of fluid flow on osteocytes. Conditional deletion of Piezo1 in osteoblasts and osteocytes notably reduced bone mass and strength in mice. Conversely, administration of a Piezo1 agonist to adult mice increased bone mass, mimicking the effects of mechanical loading. These results demonstrate that Piezo1 is a mechanosensitive ion channel by which osteoblast lineage cells sense and respond to changes in mechanical load and identify a novel target for anabolic bone therapy.

Transcriptional profiling of intramembranous and endochondral ossification after fracture in mice

Bone fracture repair represents an important clinical challenge with nearly 1 million non-union fractures occurring annually in the U.S. Gene expression differs between non-union and healthy repair, suggesting there is a pattern of gene expression that is indicative of optimal repair. Despite this, the gene expression profile of fracture repair remains incompletely understood. In this work, we used RNA-seq of two well-established murine fracture models to describe gene expression of intramembranous and endochondral bone formation. We used top differentially expressed genes, enriched gene ontology terms and pathways, callus cellular phenotyping, and histology to describe and contrast these bone formation processes across time. Intramembranous repair, as modeled by ulnar stress fracture, and endochondral repair, as modeled by femur full fracture, exhibited vastly different transcriptional profiles throughout repair. Stress fracture healing had enriched differentially expressed genes associated with bone repair and osteoblasts, highlighting the strong osteogenic repair process of this model. Interestingly, the PI3K-Akt signaling pathway was one of only a few pathways uniquely enriched in stress fracture repair. Full fracture repair involved a higher level of inflammatory and immune cell related genes than did stress fracture repair. Full fracture repair also differed from stress fracture in a robust downregulation of ion channel genes following injury, the role of which in fracture repair is unclear. This study offers a broad description of gene expression in intramembranous and endochondral ossification across several time points throughout repair and suggests several potentially intriguing genes, pathways, and cells whose role in fracture repair requires further study.

RNA-seq in Skeletal Biology

PURPOSE OF REVIEW

The goal of this paper is to review state-of-the-art transcriptome profiling methods and their recent applications in the field of skeletal biology.

RECENT FINDINGS

Next-generation sequencing of mRNA (RNA-seq) methods have been established and routinely used in skeletal biology research. RNA-seq has led to the identification of novel genes and transcription factors involved in skeletal development and disease, through its application in small and large animal models, as well as human tissue and cells. With the availability of advanced techniques such as single-cell RNA-seq, novel cell types in skeletal tissues are being identified. As the sequencing technologies are rapidly evolving, the exciting discoveries supported by transcriptomics will continue to emerge and improve our understanding of the biology of the skeleton.

NOTCH Signaling Is Activated through Mechanical Strain in Human Bone Marrow-Derived Mesenchymal Stromal Cells

Skeletal development and remodeling of adult bone are critically controlled by activated NOTCH signaling in genetically modified mice. It is yet unclear whether NOTCH signaling is activated by mechanical strain sensed by bone cells. We found that expression of specific NOTCH target genes is induced after in vivo tibial mechanical loading in wild-type mice. We further applied mechanical strain through cyclic stretching in human bone marrow-derived mesenchymal stromal cells (BMSCs) in vitro by using a bioreactor system and detected upregulation of NOTCH target gene expression. Inhibition of the NOTCH pathway in primary BMSCs as well as telomerase-immortalized human BMSCs (hMSC-TERT) through the gamma-secretase inhibitor GSI XII blocked mechanotransduction and modulated actin cytoskeleton organization. Short-hairpin RNA gene silencing identified NOTCH2 as the key receptor mediating NOTCH effects on hMSC-TERT cells. Our data indicate a functional link between NOTCH activation and mechanotransduction in human BMSCs. We suggest that NOTCH signaling is an important contributor to molecular mechanisms that mediate the bone formation response to mechanical strain.

Mechanically stimulated ATP release from murine bone cells is regulated by a balance of injury and repair

Bone cells sense and actively adapt to physical perturbations to prevent critical damage. ATP release is among the earliest cellular responses to mechanical stimulation. Mechanical stimulation of a single murine osteoblast led to the release of 70 ± 24 amole ATP, which stimulated calcium responses in neighboring cells. Osteoblasts contained ATP-rich vesicles that were released upon mechanical stimulation. Surprisingly, interventions that promoted vesicular release reduced ATP release, while inhibitors of vesicular release potentiated ATP release. Searching for an alternative ATP release route, we found that mechanical stresses induced reversible cell membrane injury in vitro and in vivo. Ca²⁺/PLC/PKC-dependent vesicular exocytosis facilitated membrane repair, thereby minimizing cell injury and reducing ATP release. Priming cellular repair machinery prior to mechanical stimulation reduced subsequent membrane injury and ATP release, linking cellular mechanosensitivity to prior mechanical exposure. Thus, our findings position ATP release as an integrated readout of membrane injury and repair.

Athletes' skeletons get stronger with training, while bones weaken in people who cannot move or in astronauts experiencing weightlessness. This is because bone cells thrive when exposed to forces. When a bone cell is exposed to a physical force, the first thing that happens is the release of the energy-rich molecule called ATP into the space outside the cell. This molecule then binds to the neighboring cell to unleash a cascade of responses. ATP can exit the cell either through special canals in the cell membrane or released in tiny pod-like structures called vesicles. It is known that strong forces can injure the cell membrane and cause ATP to spill out. However, it is less clear how ATP is released when cells are subjected to regular forces.

Mikolajewicz et al. investigated whether ATP exits through injured membranes of cells experiencing regular forces. Bone cells grown in the laboratory were gently poked with a glass needle or placed in a turbulent fluid to simulate forces experienced in the body. Dyes and fluorescent imaging techniques were used to observe the movement of vesicles and calculate the concentration of ATP in these cells.

The experiments showed that regular forces in the body do indeed injure the cell membranes and cause ATP to spill out. But importantly, the cells repaired the injuries quickly by releasing vesicles that patch the wound. As soon as the membrane is sealed, ATP stops coming out. From the first injury, cells adapted and quickly strengthened their membrane and repair system to be more resilient against future forces. This process was also seen in the shin bones of mice.

These results are important because knowing how bone cells sense, respond and convert physical forces can help us develop treatments for astronauts, the injured and aged.

Comparison of knee injury threshold during tibial compression based on limb orientation in mice

Our previous studies used tibial compression overload to induce anterior cruciate ligament (ACL) rupture in mice, while others have applied similar or greater compressive magnitudes without injury. The causes of these differences in injury threshold are not known. In this study, we compared knee injury thresholds using a "prone configuration" and a "supine configuration" that differed with respect to hip, knee, and ankle flexion, and utilized different fixtures to stabilize the knee. Right limbs of female and male C57BL/6 mice were loaded using the prone configuration, while left limbs were loaded using the supine configuration. Mice underwent progressive loading from 2 to 20 N, or cyclic loading at 9 N or 14 N (n = 9-11/sex/loading method). Progressive loading with the prone configuration resulted in ACL rupture at an average of 10.2 ± 0.9 N for females and 11.4 ± 0.7 N for males. In contrast, progressive loading with the supine configuration resulted in ACL rupture in only 36% of female mice and 50% of male mice. Cyclic loading with the prone configuration resulted in ACL rupture after 15 ± 8 cycles for females and 24 ± 27 cycles for males at 9 N, and always during the first cycle for both sexes at 14 N. In contrast, cyclic loading with the supine configuration was able to complete 1,200 cycles at 9 N without injury for both sexes, and an average of 45 ± 41 cycles for females and 49 ± 25 cycles for males at 14 N before ACL rupture. These results show that tibial compression configurations can strongly affect knee injury thresholds during loading.

Sclerostin Antibody Augments the Anabolic Bone Formation Response in a Mouse Model of Mechanical Tibial Loading

Decreased activity or expression of sclerostin, an endogenous inhibitor of Wnt/β-catenin signaling, results in increased bone formation and mass. Antibodies targeting and neutralizing sclerostin (Scl-Ab) have been shown to increase bone mass and reduce fracture risk. Sclerostin is also important in modulating the response of bone to changes in its biomechanical environment. However, the effects of Scl-Ab on mechanotransduction are unclear, and it was speculated that the loading response may be altered for individuals receiving Scl-Ab therapy. To address this, we carried out a 2-week study of tibial cyclic compressive loading on C57Bl/6 mice treated with vehicle or 100 mg/kg/wk Scl-Ab. Increases in bone volume, density, and dynamic bone formation were found with loading, and the anabolic response was further increased by the combination of load and Scl-Ab. To investigate the underlying mechanism, gene profiling by RNA sequencing (RNAseq) was performed on tibias isolated from mice from all four experimental groups. Major alterations in Wnt/β-catenin gene expression were found with tibial loading, however not with Scl-Ab treatment alone. Notably, the combination of load and Scl-Ab elicited a synergistic response from a number of specific Wnt-related and mechanotransduction factors. An unexpected finding was significant upregulation of factors in the Rho GTPase signaling pathway with combination treatment. In summary, combination therapy had a more profound anabolic response than either Scl-Ab or loading treatment alone. The Wnt/β-catenin and Rho GTPase pathways were implicated within bone mechanotransduction and support the concept that bone mechanotransduction is likely to encompass a number of interconnected signaling pathways. © 2017 American Society for Bone and Mineral Research.

Temporal mechanically-induced signaling events in bone and dorsal root ganglion neurons after in vivo bone loading

Mechanical signals play an integral role in the regulation of bone mass and functional adaptation to bone loading. The osteocyte has long been considered the principle mechanosensory cell type in bone, although recent evidence suggests the sensory nervous system may play a role in mechanosensing. The specific signaling pathways responsible for functional adaptation of the skeleton through modeling and remodeling are not clearly defined. In vitro studies suggest involvement of intracellular signaling through mitogen-activated protein kinase (MAPK), phosphatidylinositol 3-kinase (PI3K)/protein kinase B (Akt), and mammalian target of rapamycin (mTOR). However, anabolic signaling responses to bone loading using a whole animal in vivo model have not been studied in detail. Therefore, we examined mechanically-induced signaling events at five time points from 0 to 24 hours after loading using the rat in vivo ulna end-loading model. Western blot analysis of bone for MAPK's, PI3K/Akt, and mTOR signaling, and quantitative reverse transcription polymerase chain reaction (qRT-PCR) to estimate gene expression of calcitonin gene-related protein alpha (CGRP-α), brain-derived neurotrophic factor (BDNF), nerve growth factor (NGF), c-jun, and c-fos in dorsal root ganglion (DRG) of the brachial intumescence were performed. There was a significant increase in signaling through MAPK's including extracellular signal-related kinase (ERK) and c-Jun N-terminal kinase (JNK) in loaded limbs at 15 minutes after mechanical loading. Ulna loading did not significantly influence expression of the genes of interest in DRG neurons. Bone signaling and DRG gene expression from the loaded and contralateral limbs was correlated (SR>0.40, P<0.05). However, bone signaling did not correlate with expression of the genes of interest in DRG neurons. These results suggest that signaling through the MAPK pathway may be involved in load-induced bone formation in vivo. Further characterization of the molecular events involved in regulation of bone adaptation is needed to understand the timing and impact of loading events, and the contribution of the neuronal signaling to functional adaptation of bone.

Collagen XI mutation lowers susceptibility to load-induced cartilage damage in mice

Interactions among risk factors for osteoarthritis (OA) are not well understood. We investigated the combined impact of two prevalent risk factors: mechanical loading and genetically abnormal cartilage tissue properties. We used cyclic tibial compression to simulate mechanical loading in the cho/+ (Col11a1 haploinsufficient) mouse, which has abnormal collagen fibrils in cartilage due to a point mutation in the Col11a1 gene. We hypothesized that the mutant collagen would not alter phenotypic bone properties and that cho/+ mice, which develop early onset OA, would develop enhanced load-induced cartilage damage compared to their littermates. To test our hypotheses, we applied cyclic compression to the left tibiae of 6-month-old cho/+ male mice and wild-type (WT) littermates for 1, 2, and 6 weeks at moderate (4.5 N) and high (9.0 N) peak load magnitudes. We then characterized load-induced cartilage and bone changes by histology, microcomputed tomography, and immunohistochemistry. Prior to loading, cho/+ mice had less dense, thinner cortical bone compared to WT littermates. In addition, in loaded and non-loaded limbs, cho/+ mice had thicker cartilage. With high loads, cho/+ mice experienced less load-induced cartilage damage at all time points and displayed decreased matrix metalloproteinase (MMP)-13 levels compared to WT littermates. The thinner, less dense cortical bone and thicker cartilage were unexpected and may have contributed to the reduced severity of load-induced cartilage damage in cho/+ mice. Furthermore, the spontaneous proteoglycan loss resulting from the mutant collagen XI was not additive to cartilage damage from mechanical loading, suggesting that these risk factors act through independent pathways. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:711-720, 2018.

Exercise initiated after the onset of insulin resistance improves trabecular microarchitecture and cortical bone biomechanics of the tibia in hyperphagic Otsuka Long Evans Tokushima Fatty rats

The present study extends our previous findings that exercise, which prevents the onset of insulin resistance and type 2 diabetes (T2D), also prevents the detrimental effects of T2D on whole-bone and tissue-level strength. Our objective was to determine whether exercise improves bone's structural and material properties if insulin resistance is already present in the Otsuka Long-Evans Tokushima Fatty (OLETF) rat. The OLETF rat is hyperphagic due to a loss-of-function mutation in cholecystokinin-1 receptor (CCK-1 receptor), which leads to progressive obesity, insulin resistance and T2D after the majority of skeletal growth is complete. Because exercise reduces body mass, which is a significant determinant of bone strength, we used a body-mass-matched caloric-restricted control to isolate body-mass-independent effects of exercise on bone. Eight-wk old, male OLETF rats were fed ad libitum until onset of hyperglycemia (20weeks of age), at which time they were randomly assigned to three groups: ad libitum fed, sedentary (O-SED); ad libitum fed, treadmill running (O-EX); or, sedentary, mild caloric restriction to match body mass of O-EX (O-CR). Long-Evans Tokushima Otsuka rats served as the normophagic, normoglycemic controls (L-SED). At 32weeks of age, O-SED rats had T2D as evidenced by hyperglycemia and a significant reduction in fasting insulin compared to OLETFs at 20weeks of age. O-SED rats also had reduced total body bone mineral content (BMC), increased C-terminal telopeptide of type I collagen (CTx)/tartrate resistant acid phosphatase isoform 5b (TRAP5b), decreased N-terminal propeptide of type I procollagen (P1NP), reduced percent cancellous bone volume (BV/TV), trabecular number (Tb.N) and increased trabecular separation (Tb.Sp) and structural model index (SMI) of the proximal tibia compared to L-SED. T2D also adversely affected biomechanical properties of the tibial diaphysis, and serum sclerostin was increased and β-catenin, runt-related transcription factor 2 (Runx2) and insulin-like growth factor-I (IGF-I) protein expression in bone were reduced in O-SED vs. L-SED. O-EX or O-CR had greater total body bone mineral density (BMD) and BMC, and BV/TV, Tb.N, Tb.Sp, and SMI compared to O-SED. O-EX had lower CTx and CR greater P1NP relative to O-SED. O-EX, not O-CR, had greater cortical thickness and area, and improved whole-bone and tissue-level biomechanical properties associated with a 4-fold increase in cortical bone β-catenin protein expression vs. O-SED. In summary, EX or CR initiated after the onset of insulin resistance preserved cancellous bone volume and structure, and EX elicited additional benefits in cortical bone.

Bone mechanobiology in mice: toward single-cell in vivo mechanomics

Mechanically driven bone (re)modeling is a multiscale process mediated through complex interactions between multiple cell types and their microenvironments. However, the underlying mechanisms of how cells respond to mechanical signals are still unclear and are at the focus of the field of bone mechanobiology. Traditionally, this complex process has been addressed by reducing the system to single scales and cell types. It is only recently that more integrative approaches have been established to study bone mechanobiology across multiple scales in which mechanical load at the organ level is related to molecular responses at the cellular level. The availability of mouse loading models and imaging techniques with improved spatial and temporal resolution has made it possible to track dynamic bone (re)modeling at the tissue and cellular level in vivo. Coupled with advanced computational models, the (re)modeling activities at the tissue scale can be associated with the mechanical microenvironment. However, methods are lacking to link the molecular responses of different cell types to their local mechanical microenvironment and bone (re)modeling activities occurring at the tissue scale. With recent improvements in "omics" technologies and single-cell molecular biology, it is now possible to sequence the complete genome and transcriptome of single cells. These technologies offer unique opportunities to comprehensively investigate the cellular transcriptional profiles within their specific microenvironment. By combining single-cell "omics" technologies with well-established tissue-scale models of bone mechanobiology, we propose a mechanomics approach to locally analyze the transcriptome of single cells with respect to their local 3D mechanical in vivo environment.

Sclerostin antibody inhibits skeletal deterioration in mice exposed to partial weight-bearing

Whereas much is known regarding the musculoskeletal responses to full unloading, little is known about the physiological effects and response to pharmacological agents in partial unloading (e.g. Moon and Mars) environments. To address this, we used a previously developed ground-based model of partial weight-bearing (PWB) that allows chronic exposure to reduced weight-bearing in mice to determine the effects of murine sclerostin antibody (SclAbII) on bone microstructure and strength across different levels of mechanical unloading. We hypothesize that treatment with SclAbII would improve bone mass, microarchitecture and strength in all loading conditions, but that there would be a greater skeletal response in the normally loaded mice than in partially unloaded mice suggesting the importance of combined countermeasures for exploration-class long duration spaceflight missions. Eleven-week-old female mice were assigned to one of four loading groups: normal weight-bearing controls (CON) or weight-bearing at 20% (PWB20), 40% (PWB40) or 70% (PWB70) of normal. Mice in each group received either SclAbII (25mg/kg) or vehicle (VEH) via twice weekly subcutaneous injection for 3 weeks. In partially-unloaded VEH-treated groups, leg BMD decreased -5 to -10% in a load-dependent manner. SclAbII treatment completely inhibited bone deterioration due to PWB, with bone properties in SclAbII-treated groups being equal to or greater than those of CON, VEH-treated mice. SclAbII treatment increased leg BMD from +14 to +18% in the PWB groups and 30 ± 3% in CON (p< 0.0001 for all). Trabecular bone volume, assessed by μCT at the distal femur, was lower in all partially unloaded VEH-treated groups vs. CON-VEH (p< 0.05), and was 2-3 fold higher in SclAbII-treated groups (p< 0.001). Midshaft femoral strength was also significantly higher in SclAbII vs. VEH-groups in all-loading conditions. These results suggest that greater weight bearing leads to greater benefits of SclAbII on bone mass, particularly in the trabecular compartment. Altogether, these results demonstrate the efficacy of sclerostin antibody therapy in preventing astronaut bone loss during terrestrial solar system exploration.