The parathyroid hormone-regulated transcriptome in osteocytes: parallel actions with 1,25-dihydroxyvitamin D3 to oppose gene expression changes during differentiation and to promote mature cell function

Differential Effects of PTH (1-34), PTHrP (1-36), and Abaloparatide on the Murine Osteoblast Transcriptome

Teriparatide (PTH (1-34)), PTHrP (1-36), and abaloparatide (ABL) have been used for the treatment of osteoporosis, but their efficacy long term is significantly limited. The 3 peptides exert time- and dose-dependent differential responses in osteoblasts, leading us to hypothesize they may also differentially modulate the osteoblast transcriptome. Treatment of mouse calvarial osteoblasts with 1 nM of the peptides for 4 hours results in RNA sequencing data with PTH (1-34) regulating 367 genes, including 194 unique genes; PTHrP (1-36) regulating 117 genes, including 15 unique genes; and ABL regulating 179 genes, including 20 unique genes. There were 83 genes shared among all 3 peptides. Gene ontology analyses showed similarities in Wnt signaling, cAMP-mediated signaling, ossification, but differences in morphogenesis of a branching structure in biological processes; receptor ligand activity, transcription factor activity, and cytokine receptor/binding activity in molecular functions. The peptides increased Vdr, Cited1, and Pde10a messenger RNAs (mRNAs) in a pattern similar to Rankl, that is, PTH (1-34) greater than ABL greater than PTHrP (1-36). mRNA abundance of other genes, including Wnt4, Wnt7, Wnt11, Sfrp4, Dkk1, Kcnk10, Hdac4, Epn3, Tcf7, Crem, Fzd5, Ppp2r2a, and Dvl3, showed that some genes were regulated similarly by all 3 peptides; others were not. Finally, small interfering RNA knockdowns of SIK1/2/3 and CRTC1/2/3 in PTH (1-34)-treated cells revealed that Vdr and Wnt4 genes are regulated by salt-inducible kinases (SIKs) and CREB-regulated transcriptional coactivators (CRTCs), while others are not. Although many studies have examined PTH signaling in the osteoblast/osteocyte, ours is the first to compare the global effects of these peptides on the osteoblast transcriptome or to analyze the roles of the SIKs and CRTCs.

Differential effects of PTH (1-34), PTHrP (1-36) and abaloparatide on the murine osteoblast transcriptome

Teriparatide (PTH(1-34)) and its analogs, PTHrP(1-36) and abaloparatide (ABL) have been used for the treatment of osteoporosis, but their efficacy over long-term use is significantly limited. The 3 peptides exert time- and dose-dependent differential responses in osteoblasts, leading us to hypothesize that they may also differentially modulate the osteoblast transcriptome. We show that treatment of mouse calvarial osteoblasts with 1 nM of the 3 peptides for 4 h results in RNA-Seq data with PTH(1-34) regulating 367 genes, including 194 unique genes; PTHrP(1-36) regulating 117 genes, including 15 unique genes; and ABL regulating 179 genes, including 20 unique genes. There were 83 genes shared among all 3 peptides. Gene ontology analyses showed differences in Wnt signaling, cAMP-mediated signaling, bone mineralization, morphogenesis of a branching structure in biological processes; receptor ligand activity, transcription factor activity, cytokine receptor/binding activity and many other actions in molecular functions. The 3 peptides increased Vdr, Cited1 and Pde10a mRNAs in a pattern similar to Rankl , i.e., PTH(1-34) > ABL > PTHrP(1-36). mRNA abundance of other genes based on gene/pathway analyses, including Wnt4, Wnt7, Wnt11, Sfrp4, Dkk1, Kcnk10, Hdac4, Epha3, Tcf7, Crem, Fzd5, Pp2r2a , and Dvl3 showed that some genes were regulated similarly by all 3 peptides; others were not. Finally, siRNA knockdowns of SIK1/2/3 and CRTC1/2/3 in PTH(1-34)-treated cells revealed that Vdr and Wnt4 genes are regulated by SIKs and CRTCs, while others are not. Although many studies have examined PTH signaling in the osteoblast/osteocyte, ours is the first to examine the global effects of these peptides on the osteoblast transcriptome. Further delineation of which signaling events are attributable to PTH(1-34), PTHrP(1-36) or ABL exclusively and which are shared among all 3 will help improve our understanding of the effects these peptides have on the osteoblast and lead to the refinement of PTH-derived treatments for osteoporosis.

Genomic mechanisms controlling renal vitamin D metabolism

Vitamin D metabolism centers on regulation in the kidney of CYP27B1 induction by PTH, suppression by FGF23 and 1,25(OH)2D3, and reciprocal CYP24A1 suppression by PTH, and induction by FGF23 and 1,25(OH)2D3. This coordinated genomic regulation through enhancer modules results in the production and dynamic maintenance of circulating endocrine 1,25(OH)2D3 which, together with PTH and FGF23, controls mineral homeostasis. We discovered enhancers near Cyp27b1 in the mouse kidney located within intronic regions of Mettl1 and Mettl21b genes. These kidney-specific enhancers ("M1", "M21") control Cyp27b1. Through CRISPR/Cas deletion, we found that PTH activation of Cyp27b1 is lost with deletion of M1, whereas FGF23 suppression is lost with deletion of M21. The combination of both deletions (M1/M21-DIKO) eliminated the suppression by 1,25(OH)2D3. Cyp24a1 activation by 1,25(OH)2D3 is controlled by a promoter proximal pair of VDREs as well as a distal region - 35 to - 37 kb (DS2). We also found that FGF23 activation and PTH suppression of Cyp24a1 was located in a region - 21 to - 37 kb downstream (DS1). More recently, using in vivo ChIP-seq in mouse kidney, we demonstrate that PTH activation rapidly induces increased recruitment of pCREB and its coactivators, CBP and CRTC2, to the M1 and M21 enhancers near the Cyp27b1 gene. At distal enhancers of the Cyp24a1 gene, PTH suppression promotes dismisses CBP with only minor changes in pCREB and CRTC2 occupancy, all of which correlate with a suppression of basal histone acetylation across this locus and reduced transcripts. Surprisingly, we find that 1,25(OH)2D3 suppression increases the occupancy of CRTC2 in the M1 enhancer, a novel observation for CRTC2 and/or 1,25(OH)2D3 action. The suppressive actions of 1,25(OH)2D3 and FGF23 at the Cyp27b1 gene are associated with a reduction in CBP recruitment at these enhancers. Although FGF23-regulated transcription factors remain unknown, we hypothesize that VDR occupancy induced at the M1 and M21 enhancers by 1,25(OH)2D3 likely disrupts or competes with the active conformation of these CREB modules thereby preventing full induction by PTH. Our findings show coactivators such as CRTC2 and CBP contribute to Cyp27b1 and Cyp24a1 transcription and provide molecular insight into the coordinated mechanistic actions of PTH, FGF23, and 1,25(OH)2D3 in the kidney that regulate mineral homeostasis.

Increased PHOSPHO1 and alkaline phosphatase expression during the anabolic bone response to intermittent parathyroid hormone delivery

The administration of intermittent parathyroid hormone (iPTH) is anabolic to the skeleton. Recent studies with cultured osteoblasts have revealed that the expression of PHOSPHO1, a bone-specific phosphatase essential for the initiation of mineralisation, is regulated by PTH. Therefore, this study sought to determine whether the bone anabolic response to iPTH involves modulation of expression of Phospho1 and of other enzymes critical for bone matrix mineralisation. To mimic iPTH treatment, primary murine osteoblasts were challenged with 50 nM PTH for 6 h in every 48 h period for 8 days (4 cycles), 14 days (7 cycles) and 20 days (10 cycles) in total. The expression of both Phospho1 and Smpd3 was almost completely inhibited after 4 cycles, whereas 10 cycles were required to stimulate a similar response in Alpl expression. To explore the in vivo role of PHOSPHO1 in PTH-mediated osteogenesis, the effects of 14- and 28-day iPTH (80 µg/kg/day) administration was assessed in male wild-type (WT) and Phospho1-/- mice. The expression of Phospho1, Alpl, Smpd3, Enpp1, Runx2 and Trps1 expression was enhanced in the femora of WT mice following iPTH administration but remained unchanged in the femora of Phospho1-/- mice. After 28 days of iPTH administration, the anabolic response in the femora of WT was greater than that noted in Phospho1-/- mice. Specifically, cortical and trabecular bone volume/total volume, as well as cortical thickness, were increased in femora of iPTH-treated WT but not in iPTH-treated Phospho1-/- mice. Trabecular bone osteoblast number was also increased in iPTH-treated WT mice but not in iPTH-treated Phospho1-/-  mice. The increased levels of Phospho1, Alpl, Enpp1 and Smpd3 in WT mice in response to iPTH administration is consistent with their contribution to the potent anabolic properties of iPTH in bone. Furthermore, as the anabolic response to iPTH was attenuated in mice deficient in PHOSPHO1, this suggests that the osteoanabolic effects of iPTH are at least partly mediated via bone mineralisation processes.

MDS cells impair osteolineage differentiation of MSCs via extracellular vesicles to suppress normal hematopoiesis

Myelodysplastic syndrome (MDS) is a clonal disorder of hematopoietic stem cells (HSCs), characterized by ineffective hematopoiesis and frequent progression to leukemia. It has long remained unresolved how MDS cells, which are less proliferative, inhibit normal hematopoiesis and eventually dominate the bone marrow space. Despite several studies implicating mesenchymal stromal or stem cells (MSCs), a principal component of the HSC niche, in the inhibition of normal hematopoiesis, the molecular mechanisms underlying this process remain unclear. Here, we demonstrate that both human and mouse MDS cells perturb bone metabolism by suppressing the osteolineage differentiation of MSCs, which impairs the ability of MSCs to support normal HSCs. Enforced MSC differentiation rescues the suppressed normal hematopoiesis in both in vivo and in vitro MDS models. Intriguingly, the suppression effect is reversible and mediated by extracellular vesicles (EVs) derived from MDS cells. These findings shed light on the novel MDS EV-MSC axis in ineffective hematopoiesis.

Vitamins D2 and D3 Have Overlapping But Different Effects on the Human Immune System Revealed Through Analysis of the Blood Transcriptome

Vitamin D is best known for its role in maintaining bone health and calcium homeostasis. However, it also exerts a broad range of extra-skeletal effects on cellular physiology and on the immune system. Vitamins D₂ and D₃ share a high degree of structural similarity. Functional equivalence in their vitamin D-dependent effects on human physiology is usually assumed but has in fact not been well defined experimentally. In this study we seek to redress the gap in knowledge by undertaking an in-depth examination of changes in the human blood transcriptome following supplementation with physiological doses of vitamin D₂ and D₃. Our work extends a previously published randomized placebo-controlled trial that recruited healthy white European and South Asian women who were given 15 µg of vitamin D₂ or D₃ daily over 12 weeks in wintertime in the UK (Nov-Mar) by additionally determining changes in the blood transcriptome over the intervention period using microarrays. An integrated comparison of the results defines both the effect of vitamin D₃ or D₂ on gene expression, and any influence of ethnic background. An important aspect of this analysis was the focus on the changes in expression from baseline to the 12-week endpoint of treatment within each individual, harnessing the longitudinal design of the study. Whilst overlap in the repertoire of differentially expressed genes was present in the D₂ or D₃-dependent effects identified, most changes were specific to either one vitamin or the other. The data also pointed to the possibility of ethnic differences in the responses. Notably, following vitamin D₃ supplementation, the majority of changes in gene expression reflected a down-regulation in the activity of genes, many encoding pathways of the innate and adaptive immune systems, potentially shifting the immune system to a more tolerogenic status. Surprisingly, gene expression associated with type I and type II interferon activity, critical to the innate response to bacterial and viral infections, differed following supplementation with either vitamin D₂ or vitamin D₃, with only vitamin D₃ having a stimulatory effect. This study suggests that further investigation of the respective physiological roles of vitamin D₂ and vitamin D₃ is warranted.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

New Approaches to Assess Mechanisms of Action of Selective Vitamin D Analogues

Recent studies of transcription have revealed an advanced set of overarching principles that govern vitamin D action on a genome-wide scale. These tenets of vitamin D transcription have emerged as a result of the application of now well-established techniques of chromatin immunoprecipitation coupled to next-generation DNA sequencing that have now been linked directly to CRISPR-Cas9 genomic editing in culture cells and in mouse tissues in vivo. Accordingly, these techniques have established that the vitamin D hormone modulates sets of cell-type specific genes via an initial action that involves rapid binding of the VDR-ligand complex to multiple enhancer elements at open chromatin sites that drive the expression of individual genes. Importantly, a sequential set of downstream events follows this initial binding that results in rapid histone acetylation at these sites, the recruitment of additional histone modifiers across the gene locus, and in many cases, the appearance of H3K36me3 and RNA polymerase II across gene bodies. The measured recruitment of these factors and/or activities and their presence at specific regions in the gene locus correlate with the emerging presence of cognate transcripts, thereby highlighting sequential molecular events that occur during activation of most genes both in vitro and in vivo. These features provide a novel approach to the study of vitamin D analogs and their actions in vivo and suggest that they can be used for synthetic compound evaluation and to select for novel tissue- and gene-specific features. This may be particularly useful for ligand activation of nuclear receptors given the targeting of these factors directly to genetic sites in the nucleus.

Inhibition of miR-29-3p isoforms via tough decoy suppresses osteoblast function in homeostasis but promotes intermittent parathyroid hormone-induced bone anabolism

miRNAs play a vital role in post-transcriptional regulation of gene expression in osteoblasts and osteoclasts, and the miR-29 family is expressed in both lineages. Using mice globally expressing a miR-29-3p tough decoy, we demonstrated a modest 30-60% decrease all three miR-29-3p isoforms: miR-29a, miR-29b, and miR-29c. While the miR-29-3p decoy did not impact osteoclast number or function, the tough decoy decreased bone formation in growing mice, which led to decreased trabecular bone volume in mature animals. These data support previous in vitro studies suggesting that miR-29-3p is a positive regulator of osteoblast differentiation. In contrast, when mice were treated with intermittent parathyroid hormone (PTH1-34), inhibition of miR-29-3p augmented the effect of PTH on cortical bone anabolism, increased bone formation rate and osteoblast surface, and increased levels of Ctnnb1/βcatenin mRNA, which is a miR-29 target. These findings highlight differences in the mechanisms controlling basal level bone formation and bone formation induced by intermittent PTH. Overall, the global miR-29-3p tough decoy model represents a modest loss-of-function, which could be a relevant tool for assessing the possible impact of systemically administered miR-29-3p inhibitors. Our studies provide a potential rationale for co-administration of PTH1-34 and miR-29-3p inhibitors, to boost bone formation in severely affected osteoporosis patients, particularly in the cortical compartment.

Nfil3, a target of the NACA transcriptional coregulator, affects osteoblast and osteocyte gene expression differentially

Intermittent administration of PTH(1-34) has a profound osteoanabolic effect on the skeleton. At the cellular level, osteoblasts and osteocytes are two crucial cell types that respond to PTH stimulation in bone. The transcriptional cofactor Nascent polypeptide Associated Complex and coregulator alpha (NACA) is a downstream target of the PTH-Gαs-PKA axis in osteoblasts. NACA functions as a transcriptional cofactor affecting bZIP factor-mediated transcription of target promoters in osteoblasts, such as Osteocalcin (Bglap2). Here, we used RNA-Seq and ChIP-Seq against NACA in PTH-treated MC3T3-E1 osteoblastic cells to identify novel targets of the PTH-activated NACA. Our approach identified Nuclear factor interleukin-3-regulated (Nfil3) as a target promoter of this pathway. Knockdown of Naca reduced the response of Nfil3 to PTH(1-34) stimulation. In silico analysis of the Nfil3 promoter revealed potential binding sites for NACA (located within the ChIP fragment) and CREB. We show that following PTH stimulation, phosphorylated-CREB binds the proximal promoter of Nfil3 in osteoblasts. The activity of the Nfil3 promoter (-818/+182 bp) is regulated by CREB and this activation relies on the presence of NACA. In addition, we show that knockdown of Nfil3 enhances the expression of osteoblastic differentiation markers in MC3T3-E1 cells while it represses osteocytic marker gene expression in IDG-SW3 cells. These results show that the PTH-induced NACA axis regulates Nfil3 expression and suggest that NFIL3 acts as a transcriptional repressor in osteoblasts while it exhibits differential activity as an activator in osteocytes.

Parathyroid hormone ameliorates temporomandibular joint osteoarthritic-like changes related to age

OBJECTIVES

Ageing could be a contributing factor to the progression of temporomandibular joint osteoarthritis (TMJ OA), whereas its pathogenesis and potential therapeutic strategy have not been comprehensively investigated.

MATERIALS AND METHODS

We generated ageing mouse models (45-week and 60-week; 12-week mice as control) and intermittently injected 45-week mice with parathyroid hormone (PTH(1-34)) or vehicle for 4 weeks. Cartilage and subchondral bone of TMJ were analysed by microCT, histological and immunostaining. Western blot, qRT-PCR, ChIP, ELISA and immunohistochemical analysis were utilized to examination the mechanism of PTH(1-34)'s function.

RESULTS

We showed apparent OA-like phenotypes in ageing mice. PTH treatment could ameliorate the degenerative changes and improve bone microarchitecture in the subchondral bone by activating bone remodelling. Moreover, PTH inhibited phosphorylation level of Smad3, which can combine with p16^ink4a gene promoter region, resulting in reduced senescent cells accumulation and increased cellular proliferation of marrow mesenchymal stem cells (MSCs). ELISA also showed relieved levels of specific senescent-associated secretory phenotype (SASP) in ageing mice after PTH treatment.

CONCLUSIONS

In summary, PTH may reduce the accumulation of senescent cells in subchondral bone by inhibiting p16^ink4a and improve bone marrow microenvironment to active bone remodelling process, indicating PTH administration could be a potential preventative and therapeutic treatment for age-related TMJ OA.

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.

How To Build a Bone: PHOSPHO1, Biomineralization, and Beyond

Since its characterization two decades ago, the phosphatase PHOSPHO1 has been the subject of an increasing focus of research. This work has elucidated PHOSPHO1's central role in the biomineralization of bone and other hard tissues, but has also implicated the enzyme in other biological processes in health and disease. During mineralization PHOSPHO1 liberates inorganic phosphate (P_i) to be incorporated into the mineral phase through hydrolysis of its substrates phosphocholine (PCho) and phosphoethanolamine (PEA). Localization of PHOSPHO1 within matrix vesicles allows accumulation of P_i within a protected environment where mineral crystals may nucleate and subsequently invade the organic collagenous scaffold. Here, we examine the evidence for this process, first discussing the discovery and characterization of PHOSPHO1, before considering experimental evidence for its canonical role in matrix vesicle–mediated biomineralization. We also contemplate roles for PHOSPHO1 in disorders of dysregulated mineralization such as vascular calcification, along with emerging evidence of its activity in other systems including choline synthesis and homeostasis, and energy metabolism. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Parathyroid hormone(1-34) and its analogs differentially modulate osteoblastic Rankl expression via PKA/SIK2/SIK3 and PP1/PP2A-CRTC3 signaling

Osteoporosis can result from the loss of sex hormones and/or aging. Abaloparatide (ABL), an analog of parathyroid hormone-related protein (PTHrP(1-36)), is the second osteoanabolic therapy approved by the United States Food and Drug Administration after teriparatide (PTH(1-34)). All three peptides bind PTH/PTHrP receptor type 1 (PTHR1), but the effects of PTHrP(1-36) or ABL in the osteoblast remain unclear. We show that, in primary calvarial osteoblasts, PTH(1-34) promotes a more robust cAMP response than PTHrP(1-36) and ABL and causes a greater activation of protein kinase A (PKA) and cAMP response element-binding protein (CREB). All three peptides similarly inhibited sclerostin (Sost). Interestingly, the three peptides differentially modulated two other PKA target genes, c-Fos and receptor activator of NF-κB ligand (Rankl), and the latter both in vitro and in vivo Knockdown of salt-inducible kinases (SIKs) 2 and 3 and CREB-regulated transcription coactivator 3 (CRTC3), indicated that all three are part of the pathway that regulates osteoblastic Rankl expression. We also show that the peptides differentially regulate the nuclear localization of CRTC2 and CRTC3, and that this correlates with PKA activation. Moreover, inhibition of protein phosphatases 1 and 2A (PP1/PP2A) activity revealed that they play a major role in both PTH-induced Rankl expression and the effects of PTH(1-34) on CRTC3 localization. In summary, in the osteoblast, the effects of PTH(1-34), PTHrP(1-36), and ABL on Rankl are mediated by differential stimulation of cAMP/PKA signaling and by their downstream effects on SIK2 and -3, PP1/PP2A, and CRTC3.

Regulation of Bone Remodeling by Parathyroid Hormone

Parathyroid hormone (PTH) exerts profound effects on skeletal homeostasis through multiple cellular and molecular mechanisms. Continuous hyperparathyroidism causes net loss of bone mass, despite accelerating bone formation by osteoblasts. Intermittent treatment with PTH analogs represents the only Food and Drug Administration (FDA)-approved bone anabolic osteoporosis treatment strategy. Functional PTH receptors are present on cells of the osteoblast lineage, ranging from early skeletal stem cells to matrix-embedded osteocytes. In addition, bone remodeling by osteoclasts liberates latent growth factors present within bone matrix. Here, we will provide an overview of the multiple cellular and molecular mechanisms through which PTH influences bone homeostasis. Notably, net skeletal effects of continuous versus intermittent can differ significantly. Where possible, we will highlight mechanisms through which continuous hyperparathyroidism leads to bone loss, and through which intermittent hyperparathyroidism boosts bone mass. Given the therapeutic usage of intermittent PTH (iPTH) treatment for osteoporosis, particular attention will be paid toward mechanisms underlying the bone anabolic effects of once daily PTH administration.

A Novel Distal Enhancer Mediates Inflammation-, PTH-, and Early Onset Murine Kidney Disease-Induced Expression of the Mouse Fgf23 Gene

Fibroblast growth factor 23 (FGF23) production is regulated by both calciotropic hormones and inflammation. Consistent with this, elevated FGF23 levels are associated with inflammatory markers as well as parathyroid hormone (PTH) in various disease states, including chronic kidney disease (CKD). However, the molecular mechanisms underpinning Fgf23 transcription in response to these regulators are largely unknown. We therefore utilized chromatin immunoprecipitation followed by DNA sequencing (ChIP‐seq) data from an osteocyte cell line to identify potential regulatory regions of the Fgf23 gene. Based on ChIP‐seq analysis of enhancer‐associated histone modifications, including H3K4 methylation and H3K9 acetylation, we discovered several potential enhancers for Fgf23, one of which was located 16kb upstream of the gene's transcriptional start site. Deletion of this putative enhancer from the mouse genome using CRISPR‐Cas9 technology led to lower bone, thymus, and spleen expression of Fgf23 mRNA without altering circulating levels of the intact hormone, although as previously reported, only bone displayed significant basal expression. Nevertheless, lack of the −16kb enhancer blunted FGF23 upregulation in a tissue‐specific manner by the acute inflammatory inducers lipopolysaccharide (LPS), interleukin‐1‐beta (IL‐1β), and tumor necrosis factor‐alpha (TNFα) in bone, non‐osseous tissues, and in circulation. Lack of the −16kb enhancer also inhibited PTH‐induced bone Fgf23 mRNA. Moreover, the absence of this Fgf23 enhancer in an oxalate diet‐induced murine CKD model prevented the early onset induction of osseous, renal, and thymic Fgf23 mRNA levels and led to a significant blunting of elevated circulating intact FGF23 levels. These results suggest that −16kb enhancer mediates the induction of Fgf23 by inflammation and PTH and facilitates the increase in FGF23 expression in a murine model of CKD. As exemplified herein, these Fgf23 enhancer‐deleted mice will provide a unique model in which to study the role of FGF23 expression in inflammatory diseases. © 2017 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Parathyroid Hormone Signaling in Osteocytes

Osteocytes are the most abundant cell type in bone and play a central role in orchestrating skeletal remodeling, in part by producing paracrine‐acting factors that in turn influence osteoblast and osteoclast activity. Recent evidence has indicated that osteocytes are crucial cellular targets of parathyroid hormone (PTH). Here, we will review the cellular and molecular mechanisms through which PTH influences osteocyte function. Two well‐studied PTH target genes in osteocytes are SOST and receptor activator of NF‐κB ligand (RANKL). The molecular mechanisms through which PTH regulates expression of these two crucial target genes will be discussed. Beyond SOST and RANKL, PTH/PTH‐related peptide (PTHrP) signaling in osteocytes may directly influence the way osteocytes remodel their perilacunar environment to influence bone homeostasis in a cell‐autonomous manner. Here, I will highlight novel, additional mechanisms used by PTH and PTHrP to modulate bone homeostasis through effects in osteocytes. © 2017 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Spotlight on vitamin D receptor, lipid metabolism and mitochondria: Some preliminary emerging issues

Transcriptional control and modulation of calcium fluxes underpin the differentiating properties of vitamin D (1,25(OH)₂D₃). In the latest years however few studies have pointed out the relevance of the mitochondrial effects of the hormone. It is now time to focus on the metabolic results of vitamin D receptor (VDR) action in mitochondria, which can explain the pleiotropic effects of 1,25(OH)₂D₃ and may elucidate few contrasting aspects of its activity. The perturbation of lipid metabolism described in VDR knockout mice and vitamin D deficient animals can be revisited based on the newly identified mechanism of action of 1,25(OH)₂D₃ in mitochondria. From the same point of view, the controversial role of 1,25(OH)₂D₃ in adipogenesis can be better interpreted.

The vitamin D receptor functions as a transcription regulator in the absence of 1,25-dihydroxyvitamin D3

The vitamin D receptor (VDR) is a critical mediator of the biological actions of 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). As a nuclear receptor, ligand activation of the VDR leads to the protein's binding to specific sites on the genome that results in the modulation of target gene expression. The VDR is also known to play a role in the hair cycle, an action that appears to be 1,25(OH)₂D₃-independent. Indeed, in the absence of the VDR as in hereditary 1,25-dihydroxyvitamin D resistant rickets (HVDRR) both skin defects and alopecia emerge. Recently, we generated a mouse model of HVDRR without alopecia wherein a mutant human VDR lacking 1,25(OH)₂D₃-binding activity was expressed in the absence of endogenous mouse VDR. While 1,25(OH)₂D₃ failed to induce gene expression in these mice, resulting in an extensive skeletal phenotype, the receptor was capable of restoring normal hair cycling. We also noted a level of secondary hyperparathyroidism that was much higher than that seen in the VDR null mouse and was associated with an exaggerated bone phenotype as well. This suggested that the VDR might play a role in parathyroid hormone (PTH) regulation independent of 1,25(OH)₂D₃. To evaluate this hypothesis further, we contrasted PTH levels in the HVDRR mouse model with those seen in Cyp27b1 null mice where the VDR was present but the hormone was absent. The data revealed that PTH was indeed higher in Cyp27b1 null mice compared to VDR null mice. To evaluate the mechanism of action underlying such a hypothesis, we measured the expression levels of a number of VDR target genes in the duodena of wildtype mice and in transgenic mice expressing either normal or hormone-binding deficient mutant VDRs. We also compared expression levels of these genes between VDR null mice and Cyp27b1 null mice. In a subset of cases, the expression of VDR target genes was lower in mice containing the VDR as opposed to mice that did not. We suggest that the VDR may function as a selective suppressor/de-repressor of gene expression in the absence of 1,25(OH)₂D₃.

The Expression of PHOSPHO1, nSMase2 and TNAP is Coordinately Regulated by Continuous PTH Exposure in Mineralising Osteoblast Cultures

Sustained exposure to high levels of parathyroid hormone (PTH), as observed in hyperparathyroidism, is catabolic to bone. The increase in the RANKL/OPG ratio in response to continuous PTH, resulting in increased osteoclastogenesis, is well established. However, the effects of prolonged PTH exposure on key regulators of skeletal mineralisation have yet to be investigated. This study sought to examine the temporal expression of PHOSPHO1, TNAP and nSMase2 in mineralising osteoblast-like cell cultures and to investigate the effects of continuous PTH exposure on the expression of these enzymes in vitro. PHOSPHO1, nSMase2 and TNAP expression in cultured MC3T3-C14 cells significantly increased from day 0 to day 10. PTH induced a rapid downregulation of Phospho1 and Smpd3 gene expression in MC3T3-C14 cells and cultured hemi-calvariae. Alpl was differentially regulated by PTH, displaying upregulation in cultured MC3T3-C14 cells and downregulation in hemi-calvariae. PTH was also able to abolish the stimulatory effects of bone morphogenic protein 2 (BMP-2) on Smpd3 and Phospho1 expression. The effects of PTH on Phospho1 expression were mimicked with the cAMP agonist forskolin and blocked by the PKA inhibitor PKI (5-24), highlighting a role for the cAMP/PKA pathway in this regulation. The potent down-regulation of Phospho1 and Smpd3 in osteoblasts in response to continuous PTH may provide a novel explanation for the catabolic effects on the skeleton of such an exposure. Furthermore, our findings support the hypothesis that PHOSPHO1, nSMase2 and TNAP function cooperatively in the initiation of skeletal mineralisation.

Epigenetic Plasticity Drives Adipogenic and Osteogenic Differentiation of Marrow-derived Mesenchymal Stem Cells

Terminal differentiation of multipotent stem cells is achieved through a coordinated cascade of activated transcription factors and epigenetic modifications that drive gene transcription responsible for unique cell fate. Within the mesenchymal lineage, factors such as RUNX2 and PPARγ are indispensable for osteogenesis and adipogenesis, respectively. We therefore investigated genomic binding of transcription factors and accompanying epigenetic modifications that occur during osteogenic and adipogenic differentiation of mouse bone marrow-derived mesenchymal stem cells (MSCs). As assessed by ChIP-sequencing and RNA-sequencing analyses, we found that genes vital for osteogenic identity were linked to RUNX2, C/EBPβ, retinoid X receptor, and vitamin D receptor binding sites, whereas adipocyte differentiation favored PPARγ, retinoid X receptor, C/EBPα, and C/EBPβ binding sites. Epigenetic marks were clear predictors of active differentiation loci as well as enhancer activities and selective gene expression. These marrow-derived MSCs displayed an epigenetic pattern that suggested a default preference for the osteogenic pathway; however, these patterns were rapidly altered near the Adipoq, Cidec, Fabp4, Lipe, Plin1, Pparg, and Cebpa genes during adipogenic differentiation. Surprisingly, we found that these cells also exhibited an epigenetic plasticity that enabled them to trans-differentiate from adipocytes to osteoblasts (and vice versa) after commitment, as assessed by staining, gene expression, and ChIP-quantitative PCR analysis. The osteogenic default pathway may be subverted during pathological conditions, leading to skeletal fragility and increased marrow adiposity during aging, estrogen deficiency, and skeletal unloading. Taken together, our data provide an increased mechanistic understanding of the epigenetic programs necessary for multipotent differentiation of MSCs that may prove beneficial in the development of therapeutic strategies.

PTH Signaling and Epigenetic Control of Bone Remodeling

As our understanding of the mechanisms that govern bone development advance, the role of epigenetic modifications in these processes become increasingly evident. Interestingly, in parathyroid hormone (PTH)-induced bone metabolism and remodeling, recent evidence shows that PTH signaling employs a particular facet of the epigenetic machinery to elicit its desired effects. In this review, we briefly discuss the known epigenetic events occurring in cells of the osteoblast lineage. More specifically, we elaborate on current findings that reveal the utilization of histone deacetylating enzymes (HDACs) in PTH-regulated modulation of gene expression in bone.

Deletion of the Distal Tnfsf11 RL-D2 Enhancer That Contributes to PTH-Mediated RANKL Expression in Osteoblast Lineage Cells Results in a High Bone Mass Phenotype in Mice

Receptor activator of nuclear factor-κB ligand (RANKL) is a tumor necrosis factor (TNF)-like cytokine that is necessary for osteoclast formation and survival. Elevated RANKL synthesis is associated with both increased osteoclast number and bone resorption. Earlier studies identified an enhancer 76 kb upstream of the Tnfsf11 transcriptional start site (TSS) termed RL-D5 or the distal control region (DCR) that modulates RANKL expression in response to PTH, 1,25(OH)2D3,, and an array of cytokines. Mice lacking RL-D5 exhibit high bone mass associated with decreased RANKL expression in bone, spleen, and thymus. In addition to RL-D5, genome-wide studies have identified 9 additional Tnfsf11 enhancers residing upstream of the gene's TSS, which provide RANKL cell type-specificity and responsiveness to local and systemic factors. ChIP-chip analyses has revealed inducible vitamin D receptor (VDR) and cAMP response element-binding protein (CREB) binding at an enhancer termed RL-D2 23 kb upstream of the Tnfsf11 TSS in osteoblastic ST2 cells. Herein, we use ChIP-seq analyses to confirm this finding and then delete this enhancer from the mouse genome to determine its physiological role in vivo. RL-D2(-/-) primary stromal cells showed decreased RANKL-induction by both forskolin and 1,25(OH)2D3 ex vivo. Consistent with this, the parathyroid hormone (PTH) induction of RANKL expression was significantly blunted in RL-D2(-/-) mice in vivo. In contrast, lack of RL-D2 had no effect on 1,25(OH)2D3 induction of RANKL in vivo. Similar to the results found in RL-D5(-/-) mice, lack of RL-D2 led to decreased skeletal RANKL expression, resulting in decreased osteoclast numbers and a progressive increase in bone mineral density. Lack of RL-D2 increased cancellous bone mass in femur and spine but did not alter femoral cortical bone thickness. These results highlight the role of distal enhancers in the regulation of RANKL expression by PTH and perhaps 1,25(OH)2D3 and suggest that the RL-D2 and RL-D5 enhancers contribute in either an additive or synergistic manner to regulate bone remodeling.

Epigenetic histone modifications and master regulators as determinants of context dependent nuclear receptor activity in bone cells

Genomic annotation of unique and combinatorial epigenetic modifications along with transcription factor occupancy is having a profound impact on our understanding of the genome. These studies have led to a better appreciation of the dynamic nature of the epigenetic and transcription factor binding components that reveal overarching principles of the genome as well as tissue specificity. In this minireview, we discuss the presence and potential functions of several of these features across the genome in osteoblast lineage cells. We examine how these features are modulated during cellular maturation, affect transcriptional output and phenotype, and how they alter the ability of cells to respond to systemic signals directed by calcemic hormones such as 1,25-dihydroxyvitamin D3 and PTH. In particular, we describe recent experiments which indicate that progressive stages of bone cell differentiation affect RUNX2 binding to the genome, modify and restrict patterns of gene expression, and dramatically alter cellular response to the vitamin D hormone. These studies expand our understanding of mechanisms that govern steroid hormone regulation of gene expression, while highlighting the increasing complexity that is evident relative to these basic cellular processes. The results also have profound implications with respect to the impact of skeletal diseases on transcriptional outcomes as well. This article is part of a Special Issue entitled Epigenetics and Bone.

Genomic Determinants of Vitamin D-Regulated Gene Expression

Insight into mechanisms that link the actions of 1,25-dihydroxyvitamin D3 (1,25(OH)2D3) to the regulation of gene expression has evolved extensively since the initial discovery of a nuclear protein known as the vitamin D receptor (VDR). Perhaps most important was the molecular cloning of this receptor which enabled its inclusion within the nuclear receptor gene family and further studies of both its structure and regulatory function. Current studies are now refocused on the vitamin D hormone's action at the genome, where VDR together with other transcription factors coordinates the recruitment of chromatin active coregulatory complexes that participate directly in the modification of gene output. These studies highlight the role of chromatin in the expression of genes and the dynamic impact of the epigenetic landscape that contextualizes individual gene loci thus influencing the VDR's transcriptional actions. In this chapter, we summarize advances made over the past few years in understanding vitamin D action on a genome-wide scale, focusing on overarching principles that have emerged at this level. Of particular significance is the finding that dynamic changes that occur to the genome during cellular differentiation at both genetic and epigenetic levels profoundly alter the ability of 1,25(OH)2D3 and its receptor to regulate gene expression. We address the broad impact of differentiation on specific epigenetic histone modifications that occur across the genome and the ability of the VDR to influence this activity at selected gene loci as well. These studies advance our understanding of not only vitamin D action but also of the complex and dynamic role played by the genome itself as a major determinant of VDR activity.

PTH and Vitamin D Repress DMP1 in Cementoblasts

A complex feedback mechanism between parathyroid hormone (PTH), 1,25(OH)2D3 (1,25D), and fibroblast growth factor 23 (FGF-23) maintains mineral homeostasis, in part by regulating calcium and phosphate absorption/reabsorption. Previously, we showed that 1,25D regulates mineral homeostasis by repressing dentin matrix protein 1 (DMP1) via the vitamin D receptor pathway. Similar to 1,25D, PTH may modulate DMP1, but the underlying mechanism remains unknown. Immortalized murine cementoblasts (OCCM.30), similar to osteoblasts and known to express DMP1, were treated with PTH (1-34). Real-time quantitative polymerase chain reaction (PCR) and Western blot revealed that PTH decreased DMP1 gene transcription (85%) and protein expression (30%), respectively. PTH mediated the downregulation of DMP1 via the cAMP/protein kinase A (PKA) pathway. Immunohistochemistry confirmed the decreased localization of DMP1 in vivo in cellular cementum and alveolar bone of mice treated with a single dose (50 µg/kg) of PTH (1-34). RNA-seq was employed to further identify patterns of gene expression shared by PTH and 1,25D in regulating DMP1, as well as other factors involved in mineral homeostasis. PTH and 1,25D mutually upregulated 36 genes and mutually downregulated 27 genes by ≥2-fold expression (P ≤ 0.05). Many identified genes were linked with the regulation of bone/tooth homeostasis, cell growth and differentiation, calcium signaling, and DMP1 transcription. Validation of RNA-seq results via PCR array confirmed a similar gene expression pattern in response to PTH and 1,25D treatment. Collectively, these results suggest that PTH and 1,25D share complementary effects in maintaining mineral homeostasis by mutual regulation of genes/proteins associated with calcium and phosphate metabolism while also exerting distinct roles on factors modulating mineral metabolism. Furthermore, PTH may modulate phosphate homeostasis by downregulating DMP1 expression via the cAMP/PKA pathway. Targeting genes/proteins mutually governed by PTH and 1,25D may be a viable approach for designing new therapies for preserving mineralized tissue health.

Parathyroid Hormone Induces Bone Cell Motility and Loss of Mature Osteocyte Phenotype through L-Calcium Channel Dependent and Independent Mechanisms

Parathyroid Hormone (PTH) can exert both anabolic and catabolic effects on the skeleton, potentially through expression of the PTH type1 receptor (PTH1R), which is highly expressed in osteocytes. To determine the cellular and molecular mechanisms responsible, we examined the effects of PTH on osteoblast to osteocyte differentiation using primary osteocytes and the IDG-SW3 murine cell line, which differentiate from osteoblast to osteocyte-like cells in vitro and express GFP under control of the dentin matrix 1 (Dmp1) promoter. PTH treatment resulted in an increase in some osteoblast and early osteocyte markers and a decrease in mature osteocyte marker expression. The gene expression profile of PTH-treated Day 28 IDG-SW3 cells was similar to PTH treated primary osteocytes. PTH treatment induced striking changes in the morphology of the Dmp1-GFP positive cells in IDG-SW3 cultures and primary cells from Dmp1-GFP transgenic mice. The cells changed from a more dendritic to an elongated morphology and showed increased cell motility. E11/gp38 has been shown to be important for cell migration, however, deletion of the E11/gp38/podoplanin gene had no effect on PTH-induced motility. The effects of PTH on motility were reproduced using cAMP, but not with protein kinase A (PKA), exchange proteins activated by cAMP (Epac), protein kinase C (PKC) or phosphatidylinositol-4,5-bisphosphonate 3-kinase (Pi3K) agonists nor were they blocked by their antagonists. However, the effects of PTH were mediated through calcium signaling, specifically through L-type channels normally expressed in osteoblasts but decreased in osteocytes. PTH was shown to increase expression of this channel, but decrease the T-type channel that is normally more highly expressed in osteocytes. Inhibition of L-type calcium channel activity attenuated the effects of PTH on cell morphology and motility but did not prevent the downregulation of mature osteocyte marker expression. Taken together, these results show that PTH induces loss of the mature osteocyte phenotype and promotes the motility of these cells. These two effects are mediated through different mechanisms. The loss of phenotype effect is independent and the cell motility effect is dependent on calcium signaling.