Genomic determinants of gene regulation by 1,25-dihydroxyvitamin D3 during osteoblast-lineage cell differentiation

Epigenetic regulators controlling osteogenic lineage commitment and bone formation

Bone formation and homeostasis are controlled by environmental factors and endocrine regulatory cues that initiate intracellular signaling pathways capable of modulating gene expression in the nucleus. Bone-related gene expression is controlled by nucleosome-based chromatin architecture that limits the accessibility of lineage-specific gene regulatory DNA sequences and sequence-specific transcription factors. From a developmental perspective, bone-specific gene expression must be suppressed during the early stages of embryogenesis to prevent the premature mineralization of skeletal elements during fetal growth in utero. Hence, bone formation is initially inhibited by gene suppressive epigenetic regulators, while other epigenetic regulators actively support osteoblast differentiation. Prominent epigenetic regulators that stimulate or attenuate osteogenesis include lysine methyl transferases (e.g., EZH2, SMYD2, SUV420H2), lysine deacetylases (e.g., HDAC1, HDAC3, HDAC4, HDAC7, SIRT1, SIRT3), arginine methyl transferases (e.g., PRMT1, PRMT4/CARM1, PRMT5), dioxygenases (e.g., TET2), bromodomain proteins (e.g., BRD2, BRD4) and chromodomain proteins (e.g., CBX1, CBX2, CBX5). This narrative review provides a broad overview of the covalent modifications of DNA and histone proteins that involve hundreds of enzymes that add, read, or delete these epigenetic modifications that are relevant for self-renewal and differentiation of mesenchymal stem cells, skeletal stem cells and osteoblasts during osteogenesis.

The neutrophil-osteogenic cell axis promotes bone destruction in periodontitis

The immune-stromal cell interactions play a key role in health and diseases. In periodontitis, the most prevalent infectious disease in humans, immune cells accumulate in the oral mucosa and promote bone destruction by inducing receptor activator of nuclear factor-κB ligand (RANKL) expression in osteogenic cells such as osteoblasts and periodontal ligament cells. However, the detailed mechanism underlying immune-bone cell interactions in periodontitis is not fully understood. Here, we performed single-cell RNA-sequencing analysis on mouse periodontal lesions and showed that neutrophil-osteogenic cell crosstalk is involved in periodontitis-induced bone loss. The periodontal lesions displayed marked infiltration of neutrophils, and in silico analyses suggested that the neutrophils interacted with osteogenic cells through cytokine production. Among the cytokines expressed in the periodontal neutrophils, oncostatin M (OSM) potently induced RANKL expression in the primary osteoblasts, and deletion of the OSM receptor in osteogenic cells significantly ameliorated periodontitis-induced bone loss. Epigenomic data analyses identified the OSM-regulated RANKL enhancer region in osteogenic cells, and mice lacking this enhancer showed decreased periodontal bone loss while maintaining physiological bone metabolism. These findings shed light on the role of neutrophils in bone regulation during bacterial infection, highlighting the novel mechanism underlying osteoimmune crosstalk.

Neuropilin 2 in osteoblasts regulates trabecular bone mass in male mice

Introduction

Neuropilin 2 (NRP2) mediates the effects of class 3 semaphorins and vascular endothelial growth factor and is implicated in axonal guidance and angiogenesis. Moreover, NRP2 expression is suggested to be involved in the regulation of bone homeostasis. Indeed, osteoblasts and osteoclasts express NRP2 and male and female global Nrp2 knockout mice have a reduced bone mass accompanied by reduced osteoblast and increased osteoclast counts.

Methods

We first examined the in vitro effect of the calciotropic hormone 1,25-dihydroxyvitamin D3 [1,25(OH)2D3] on Nrp2 transcription in osteoblasts. We next generated mice with a conditional deletion of Nrp2 in the osteoblast cell lineage under control of the paired related homeobox 1 promoter and mice with a conditional Nrp2 knockdown in osteoclasts under control of the Lysozyme promoter. Mice were examined under basal conditions or after treatment with either the bone anabolic vitamin D3 analog WY 1048 or with 1,25(OH)2D3.

Results and discussion

We show that Nrp2 expression is induced by 1,25(OH)2D3 in osteoblasts and is associated with enrichment of the vitamin D receptor in an intronic region of the Nrp2 gene. In male mice, conditional deletion of Nrp2 in osteoblast precursors and mature osteoblasts recapitulated the bone phenotype of global Nrp2 knockout mice, with a reduced cortical cross-sectional tissue area and lower trabecular bone content. However, female mice with reduced osteoblastic Nrp2 expression display a reduced cross-sectional tissue area but have a normal trabecular bone mass. Treatment with the vitamin D3 analog WY 1048 (0.4 μg/kg/d, 14 days, ip) resulted in a similar increase in bone mass in both genotypes and genders. Deleting Nrp2 from the osteoclast lineage did not result in a bone phenotype, even though in vitro osteoclastogenesis of hematopoietic cells derived from mutant mice was significantly increased. Moreover, treatment with a high dose of 1,25(OH)2D3 (0.5 μg/kg/d, 6 days, ip), to induce osteoclast-mediated bone resorption, resulted in a similar reduction in trabecular and cortical bone mass. In conclusion, osteoblastic Nrp2 expression is suggested to regulate bone homeostasis in a sex-specific manner.

Molecular insights into mineralotropic hormone inter-regulation

The regulation of mineral homeostasis involves the three mineralotropic hormones PTH, FGF23 and 1,25-dihydroxyvitamin D3 (1,25(OH)2D3). Early research efforts focused on PTH and 1,25(OH)2D3 and more recently on FGF23 have revealed that each of these hormones regulates the expression of the other two. Despite early suggestions of transcriptional processes, it has been only recently that research effort have begun to delineate the genomic mechanisms underpinning this regulation for 1,25(OH)2D3 and FGF23; the regulation of PTH by 1,25(OH)2D3, however, remains obscure. We review here our molecular understanding of how PTH induces Cyp27b1 expression, the gene encoding the enzyme responsible for the synthesis of 1,25(OH)2D3. FGF23 and 1,25(OH)2D3, on the other hand, function by suppressing production of 1,25(OH)2D3. PTH stimulates the PKA-induced recruitment of CREB and its coactivator CBP at CREB occupied sites within the kidney-specific regulatory regions of Cyp27b1. PKA activation also promotes the nuclear translocation of SIK bound coactivators such as CRTC2, where it similarly interacts with CREB occupied Cyp27b1 sites. The negative actions of both FGF23 and 1,25(OH)2D3 appear to suppress Cyp27b1 expression by opposing the recruitment of CREB coactivators at this gene. Reciprocal gene actions are seen at Cyp24a1, the gene encoding the enzyme that degrades 1,25(OH)2D3, thereby contributing to the overall regulation of blood levels of 1,25(OH)2D3. Relative to PTH regulation, we summarize what is known of how 1,25(OH)2D3 regulates PTH suppression. These studies suggest that it is not 1,25(OH)2D3 that controls PTH levels in healthy subjects, but rather calcium itself. Finally, we describe current progress using an in vivo approach that furthers our understanding of the regulation of Fgf23 expression by PTH and 1,25(OH)2D3 and provide the first evidence that P may act to induce Fgf23 expression via a complex transcriptional mechanism in bone. It is clear, however, that additional advances will need to be made to further our understanding of the inter-regulation of each of these hormonal genes.

Genomically anchored vitamin D receptor mediates an abundance of bioprotective actions elicited by its 1,25-dihydroxyvitamin D hormonal ligand

The nuclear vitamin D receptor (VDR) mediates the actions of its physiologic 1,25-dihydroxyvitamin D3 (1,25D) ligand produced in kidney and at extrarenal sites during times of physiologic and cellular stress. The ligand-receptor complex transcriptionally controls genes encoding factors that regulate calcium and phosphate sensing/transport, bone remodeling, immune function, and nervous system maintenance. With the aid of parathyroid hormone (PTH) and fibroblast growth factor 23 (FGF23), 1,25D/VDR primarily participates in an intricate network of feedback controls that govern extracellular calcium and phosphate concentrations, mainly influencing bone formation and mineralization, ectopic calcification, and indirectly supporting many fundamental roles of calcium. Beyond endocrine and intracrine effects, 1,25D/VDR signaling impacts multiple biochemical phenomena that potentially affect human health and disease, including autophagy, carcinogenesis, cell growth/differentiation, detoxification, metabolic homeostasis, and oxidative stress mitigation. Several health advantages conferred by 1,25D/VDR appear to be promulgated by induction of klotho, an anti-aging renal peptide hormone which functions as a co-receptor for FGF23 and, like 1,25D, regulates nrf2, foxo, mTOR and other cellular protective pathways. Among hundreds of genes for which expression is modulated by 1,25D/VDR either primarily or secondarily in a cell-specific manner, the resulting gene products (in addition to those expressed in the classic skeletal mineral regulatory tissues kidney, intestine, and bone), fall into multiple biochemical categories including apoptosis, cholesterol homeostasis, glycolysis, hypoxia, inflammation, p53 signaling, unfolded protein response and xenobiotic metabolism. Thus, 1,25D/VDR is a bone mineral control instrument that also signals the maintenance of multiple cellular processes in the face of environmental and genetic challenges.

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.

Epigenetic regulation during 1,25-dihydroxyvitamin D3-dependent gene transcription

Multiple evidence accumulated over the years, demonstrates that vitamin D-dependent physiological control in vertebrates occurs primarily through the regulation of target gene transcription. In addition, there has been an increasing appreciation of the role of the chromatin organization of the genome on the ability of the active form of vitamin D, 1,25(OH)₂D₃, and its specific receptor VDR to regulate gene expression. Chromatin structure in eukaryotic cells is principally modulated through epigenetic mechanisms including, but not limited to, a wide number of post-translational modifications of histone proteins and ATP-dependent chromatin remodelers, which are operative in different tissues during response to physiological cues. Hence, there is necessity to understand in depth the epigenetic control mechanisms that operate during 1,25(OH)₂D₃-dependent gene regulation. This chapter provides a general overview about epigenetic mechanisms functioning in mammalian cells and discusses how some of these mechanisms represent important components during transcriptional regulation of the model gene system CYP24A1 in response to 1,25(OH)₂D₃.

p57Kip2 is an essential regulator of vitamin D receptor-dependent mechanisms

A cyclin-dependent kinase (CDK) inhibitor, p57Kip2, is an important molecule involved in bone development; p57Kip2-deficient (p57-/-) mice display neonatal lethality resulting from abnormal bone formation and cleft palate. The modulator 1α,25-dihydroxyvitamin D3 (l,25-(OH)2VD3) has shown the potential to suppress the proliferation and induce the differentiation of normal and tumor cells. The current study assessed the role of p57Kip2 in the 1,25-(OH)2VD3-regulated differentiation of osteoblasts because p57Kip2 is associated with the vitamin D receptor (VDR). Additionally, 1,25-(OH)2VD3 treatment increased p57KIP2 expression and induced the colocalization of p57KIP2 with VDR in the osteoblast nucleus. Primary p57-/- osteoblasts exhibited higher proliferation rates with Cdk activation than p57+/+ cells. A lower level of nodule mineralization was observed in p57-/- osteoblasts than in p57+/+ cells. In p57+/+ osteoblasts, 1,25-(OH)2VD3 upregulated the p57Kip2 and opn mRNA expression levels, while the opn expression levels were significantly decreased in p57-/- cells. The osteoclastogenesis assay performed using bone marrow cocultured with 1,25-(OH)2VD3-treated osteoblasts revealed a decreased efficiency of 1,25-(OH)2VD3-stimulated osteoclastogenesis in p57-/- cells. Based on these results, p57Kip2 might function as a mediator of 1,25-(OH)2VD3 signaling, thereby enabling sufficient VDR activation for osteoblast maturation.

PTH, FGF-23, Klotho and Vitamin D as regulators of calcium and phosphorus: Genetics, epigenetics and beyond

The actions of several bone-mineral ion regulators, namely PTH, FGF23, Klotho and 1,25(OH)2 vitamin D (1,25(OH)2D), control calcium and phosphate metabolism, and each of these molecules has additional biological effects related to cell signaling, metabolism and ultimately survival. Therefore, these factors are tightly regulated at various levels - genetic, epigenetic, protein secretion and cleavage. We review the main determinants of mineral homeostasis including well-established genetic and post-translational regulators and bring attention to the epigenetic mechanisms that affect the function of PTH, FGF23/Klotho and 1,25(OH)2D. Clinically relevant epigenetic mechanisms include methylation of cytosine at CpG-rich islands, histone deacetylation and micro-RNA interference. For example, sporadic pseudohypoparathyroidism type 1B (PHP1B), a disease characterized by resistance to PTH actions due to blunted intracellular cAMP signaling at the PTH/PTHrP receptor, is associated with abnormal methylation at the GNAS locus, thereby leading to reduced expression of the stimulatory G protein α-subunit (Gsα). Post-translational regulation is critical for the function of FGF-23 and such modifications include glycosylation and phosphorylation, which regulate the cleavage of FGF-23 and hence the proportion of available FGF-23 that is biologically active. While there is extensive data on how 1,25(OH)2D and the vitamin D receptor (VDR) regulate other genes, much more needs to be learned about their regulation. Reduced VDR expression or VDR mutations are the cause of rickets and are thought to contribute to different disorders. Epigenetic changes, such as increased methylation of the VDR resulting in decreased expression are associated with several cancers and infections. Genetic and epigenetic determinants play crucial roles in the function of mineral factors and their disorders lead to different diseases related to bone and beyond.

miR-155-5p can be involved in acquisition of osseointegration on titanium surface

Dental implants made of titanium are commonly used. Although titanium implants succeed by osseointegration with bone, the detailed molecular mechanism of osseointegration is unclear. To clarify the involvement of microRNA (miRNA) in the acquisition of osseointegration on titanium, here we compared the miRNA expression profiles of mouse osteoblast-like cells (MC3T3-E1) cultured on titanium-, gold-, and stainless steel-coating glass dishes by microarray analysis. Three kinds of metals, namely titanium, gold, and stainless steel, were coated on the surface of the glass dishes by sputtering with similar roughness and shape of their surface. After MC3T3-E1 cells were cultured on the dishes without coating or coating with titanium, gold, or stainless steel for 6 h, total RNA was extracted, and miRNA expression was analyzed by microarray. To confirm the expression of the selected miRNA during osteogenic differentiation of MC3T3-E1 cells, real-time PCR analysis was performed. Furthermore, the effects of selected miRNA were examined by ectopic overexpression in MC3T3-E1 cells. The microarray analysis revealed that the expressions of miR-155-5p and miR-7023-3p were significantly increased in MC3T3-E1 cells cultured on titanium-coating glass dishes, compared to non-coating, gold-, and stainless steel-coating glass dishes. Interestingly, miR-155-5p was upregulated during osteogenic differentiation of MC3T3-E1 cells. Furthermore, overexpression of miR-155-5p enhanced the expression of Runx2 and Col1a1. In this study, miR-155-5p may be involved in the acquisition of osseointegration on titanium implant via upregulating osteogenic differentiation-related genes.

Vitamin D and Its Target Genes

The vitamin D metabolite 1α,25-dihydroxyvitamin D₃ is the natural, high-affinity ligand of the transcription factor vitamin D receptor (VDR). In many tissues and cell types, VDR binds in a ligand-dependent fashion to thousands of genomic loci and modulates, via local chromatin changes, the expression of hundreds of primary target genes. Thus, the epigenome and transcriptome of VDR-expressing cells is directly affected by vitamin D. Vitamin D target genes encode for proteins with a large variety of physiological functions, ranging from the control of calcium homeostasis, innate and adaptive immunity, to cellular differentiation. This review will discuss VDR’s binding to genomic DNA, as well as its genome-wide locations and interaction with partner proteins, in the context of chromatin. This information will be integrated into a model of vitamin D signaling, explaining the regulation of vitamin D target genes.

SLC38A2 provides proline to fulfil unique synthetic demands arising during osteoblast differentiation and bone formation

Cellular differentiation is associated with the acquisition of a unique protein signature which is essential to attain the ultimate cellular function and activity of the differentiated cell. This is predicted to result in unique biosynthetic demands that arise during differentiation. Using a bioinformatic approach, we discovered osteoblast differentiation is associated with increased demand for the amino acid proline. When compared to other differentiated cells, osteoblast-associated proteins including RUNX2, OSX, OCN and COL1A1 are significantly enriched in proline. Using a genetic and metabolomic approach, we demonstrate that the neutral amino acid transporter SLC38A2 acts cell autonomously to provide proline to facilitate the efficient synthesis of proline-rich osteoblast proteins. Genetic ablation of SLC38A2 in osteoblasts limits both osteoblast differentiation and bone formation in mice. Mechanistically, proline is primarily incorporated into nascent protein with little metabolism observed. Collectively, these data highlight a requirement for proline in fulfilling the unique biosynthetic requirements that arise during osteoblast differentiation and bone formation.

Transcriptional Regulation of 25-Hydroxyvitamin D-24-Hydroxylase (CYP24A1) by Calcemic Factors in Keratinocytes

CYP24A1 regulates serum vitamin D (VD) levels by inactivating 25(OH)₂D₃, which is the precursor of the active form of VD [1α,25(OH)₂D₃], and CYP24A1 expression is controlled by multiple calcemic factors such as 1α,25(OH)₂D₃, calcium, and phosphate. A major phosphaturic factor, FGF23, has also been identified as a regulator of serum VD levels by affecting renal CYP24A1 gene expression; however, its effect on CYP24A1 in extrarenal cells remains largely unstudied. Therefore, the direct effect of FGF23 on CYP24A1 was examined in a human keratinocyte cell line (HaCaT). In this cell line, significant induction of CYP24A1 gene expression by 1α,25(OH)₂D₃ was seen within 4 h by qRT-PCR, and this was mediated by the VD receptor, as shown in a mutant cell line genetically deficient in this receptor. However, FGF23 treatment up to 12 h did not induce CYP24A1 expression, although the expected activation of the downstream MAPK signaling pathway was seen. High calcium and phosphate treatments were also ineffective in inducing CYP24A1 gene expression. Furthermore, a luciferase assay showed no activation of a VD-sensitive proximal CYP24A1 promoter in response to the calcium and phosphate treatments, suggesting that the effect of FGF23 on dermal CYP24A1 gene expression is indirect. From these findings, we speculate that CYP24A1 gene regulation by FGF23 occurs mainly in renal cells, but not in extrarenal cells, at least not in keratinocytes.

Inhibition of Nrf2 degradation alleviates age-related osteoporosis induced by 1,25-Dihydroxyvitamin D deficiency

Previous studies have shown that 1,25(OH)₂D plays an anti-osteoporosis role by an anti-aging mechanism. Oxidative stress is a key mediator of aging and bone loss; however, whether 1,25(OH)₂D can exert its anti-osteoporosis effect by inhibiting oxidative stress is unclear. In this study, osteoporosis and the bone aging phenotype induced by 1,25(OH)₂D deficiency in male mice were significantly rescued in vivo upon the supplementation of oltipraz, an inhibitor of Nrf2 degradation. Increased oxidative stress, cellular senescence and reduced osteogenesis of BM-MSCs from VDR knockout mice were also significantly rescued when the cells were pre-treated with oltipraz. We found that 1,25(OH)₂D₃ promoted Nrf2 accumulation by inhibiting its ubiquitin-proteasome degradation, thus facilitating Nrf2 activation of its transcriptional targets. Mechanistically, 1,25(OH)₂D₃ enhances VDR-mediated recruitment of Ezh2 and facilitation of H3K27me3 action at the promoter region of Keap1, thus transcriptionally repressing Keap1. To further validate that the Nrf2-Keap1 pathway serves as the key mediator in the anabolic effect of 1,25(OH)₂D₃ on bone, Nrf2^-/- mice, or hBM-MSCs with shRNA-mediated Nrf2-knockdown, were treated with 1,25(OH)₂D₃; we found that Nrf2 knockout largely blocked the bone anabolic effect of 1,25(OH)₂D₃ in vivo and ex vivo, and Nrf2 knockdown in hBM-MSCs markedly blocked the role of 1,25(OH)₂D₃ in inhibiting oxidative stress and promoting osteogenic differentiation and bone formation. This study provides insight into the mechanism whereby 1,25(OH)₂D₃ postpones age-related osteoporosis via VDR-mediated activation of Nrf2-antioxidant signaling and inhibition of oxidative stress, and thus provides evidence for oltipraz as a potential reagent for clinical prevention and treatment of age-related osteoporosis.

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.

Systemic Jak1 activation causes extrarenal calcitriol production and skeletal alterations provoking stunted growth

Mineral homeostasis is regulated by a complex network involving endocrine actions by calcitriol, parathyroid hormone (PTH), and FGF23 on several organs including kidney, intestine, and bone. Alterations of mineral homeostasis are found in chronic kidney disease and other systemic disorders. The interplay between the immune system and the skeletal system is not fully understood, but cytokines play a major role in modulating calcitriol production and function. One of the main cellular signaling pathways mediating cytokine function is the Janus kinase (JAK)--signal transducer and activator of transcription (STAT) pathway. Here, we used a mouse model (Jak1^S645P+/- ) that resembles a constitutive activating mutation of the Jak1/Stat3 signaling pathway in humans, and shows altered mineral metabolism, with higher fibroblast growth factor 23 (FGF23) levels, lower PTH levels, and higher calcitriol levels. The higher calcitriol levels are probably due to extrarenal calcitriol production. Furthermore, systemic Jak1/Stat3 activation led to growth impairment and skeletal alterations. The growth plate in long bones showed decreased chondrocyte proliferation rates and reduced height of terminal chondrocytes. Furthermore, we demonstrate that Jak1 is also involved in bone remodeling early in life. Jak1^S645P+/- animals have decreased bone and cortical volume, imbalanced bone remodeling, reduced MAP kinase signaling, and local inflammation. In conclusion, Jak1 plays a major role in bone health probably both, directly and systemically by regulating mineral homeostasis. Understanding the role of this signaling pathway will contribute to a better knowledge in bone growth and in mineral physiology, and to the development of selective Jak inhibitors as osteoprotective agents.

Vitamin K and D Supplementation and Bone Health in Chronic Kidney Disease-Apart or Together?

Vitamin K (VK) and vitamin D (VD) deficiency/insufficiency is a common feature of chronic kidney disease (CKD), leading to impaired bone quality and a higher risk of fractures. CKD patients, with disturbances in VK and VD metabolism, do not have sufficient levels of these vitamins for maintaining normal bone formation and mineralization. So far, there has been no consensus on what serum VK and VD levels can be considered sufficient in this particular population. Moreover, there are no clear guidelines how supplementation of these vitamins should be carried out in the course of CKD. Based on the existing results of preclinical studies and clinical evidence, this review intends to discuss the effect of VK and VD on bone remodeling in CKD. Although the mechanisms of action and the effects of these vitamins on bone are distinct, we try to find evidence for synergy between them in relation to bone metabolism, to answer the question of whether combined supplementation of VK and VD will be more beneficial for bone health in the CKD population than administering each of these vitamins separately.

Genomic Mechanisms Governing Mineral Homeostasis and the Regulation and Maintenance of Vitamin D Metabolism

Our recent genomic studies identified a complex kidney‐specific enhancer module located within the introns of adjacent Mettl1 (M1) and Mettl21b (M21) genes that mediate basal and PTH induction of Cyp27b1, as well as suppression by FGF23 and 1,25‐dihydroxyvitamin D₃ [1,25(OH)₂D₃]. The tissue specificity for this regulatory module appears to be localized exclusively to renal proximal tubules. Gross deletion of these segments in mice has severe consequences on skeletal health, and directly affects Cyp27b1 expression in the kidney. Deletion of both the M1 and M21 submodules together almost completely eliminates basal Cyp27b1 expression in the kidney, creating a renal specific pseudo‐null mouse, resulting in a systemic and skeletal phenotype similar to that of the Cyp27b1‐KO mouse caused by high levels of both 25‐hydroxyvitamin D₃ [25(OH)D₃] and PTH and depletion of 1,25(OH)₂D₃. Cyp24a1 levels in the double KO mouse also decrease because of compensatory downregulation of the gene by elevated PTH and reduced FGF23 that is mediated by an intergenic module located downstream of the Cyp24a1 gene. Outside of the kidney in nonrenal target cells (NRTCs), expression of Cyp27b1 in these mutant mice was unaffected. Dietary normalization of calcium, phosphate, PTH, and FGF23 rescues the aberrant phenotype of this mouse and normalizes the skeleton. In addition, both the high levels of 25(OH)D₃ were reduced and the low levels of 1,25(OH)₂D₃ were fully eliminated in these mutant mice as a result of the rescue‐induced normalization of renal Cyp24a1. Thus, these hormone‐regulated enhancers for both Cyp27b1 and Cyp24a1 in the kidney are responsible for the circulating levels of 1,25(OH)₂D₃ in the blood. The retention of Cyp27b1 and Cyp24a1 expression in NRTCs of these endocrine 1,25(OH)₂D₃‐deficient mice suggests that this Cyp27b1 pseudo‐null mouse will provide a model for the future exploration of the role of NRTC‐produced 1,25(OH)₂D₃ in the hormone's diverse noncalcemic actions in both health and disease. © 2020 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

RNF40 exerts stage-dependent functions in differentiating osteoblasts and is essential for bone cell crosstalk

The role of histone ubiquitination in directing cell lineage specification is only poorly understood. Our previous work indicated a role of the histone 2B ubiquitin ligase RNF40 in controlling osteoblast differentiation in vitro. Here, we demonstrate that RNF40 has a stage-dependent function in controlling osteoblast differentiation in vivo. RNF40 expression is essential for early stages of lineage specification, but is dispensable in mature osteoblasts. Paradoxically, while osteoblast-specific RNF40 deletion led to impaired bone formation, it also resulted in increased bone mass due to impaired bone cell crosstalk. Loss of RNF40 resulted in decreased osteoclast number and function through modulation of RANKL expression in OBs. Mechanistically, we demonstrate that Tnfsf11 (encoding RANKL) is an important target gene of H2B monoubiquitination. These data reveal an important role of RNF40-mediated H2B monoubiquitination in bone formation and remodeling and provide a basis for exploring this pathway for the treatment of conditions such as osteoporosis or cancer-associated osteolysis.

Downregulation of the Ca2+-activated K+ channel KCa3.1 in mouse preosteoblast cells treated with vitamin D receptor agonist

The maturity of osteoblasts by proliferation and differentiation in preosteoblasts is essential for maintaining bone homeostasis. The beneficial effects of vitamin D on bone homeostasis in mammals have been demonstrated experimentally and clinically. However, the direct actions of vitamin D on preosteoblasts remain to be fully elucidated. In this study, we found that the functional activity of intermediate-conductance Ca²⁺-activated K⁺ channels (K_Ca3.1) positively regulated cell proliferation in MC3T3-E1 cells derived from mouse preosteoblasts by enhancing intracellular Ca²⁺ signaling. We examined the effects of treatment with vitamin D receptor (VDR) agonist on the expression and activity of K_Ca3.1 by real-time PCR examination, Western blotting, Ca²⁺ imaging, and patch clamp analyses in mouse MC3T3-E1 cells. Following the downregulation of K_Ca3.1 transcriptional modulators such as Fra-1 and HDAC2, K_Ca3.1 activity was suppressed in MC3T3-E1 cells treated with VDR agonists. Furthermore, application of the K_Ca3.1 activator DCEBIO attenuated the VDR agonist-evoked suppression of cell proliferation rate. These findings suggest that a decrease in K_Ca3.1 activity is involved in the suppression of cell proliferation rate in VDR agonist-treated preosteoblasts. Therefore, K_Ca3.1 plays an important role in bone formation by promoting osteoblastic proliferation under physiological conditions.

VDR in salivary gland homeostasis and cancer

The vitamin D receptor (VDR) and its ligand 1,25(OH)₂D₃ (1,25D) impact differentiation and exert anti-tumor effects in many tissues, but its role in salivary gland has yet to be defined. Using immunohistochemistry (IHC), we have detected strong VDR expression in murine and human salivary gland ducts. Compared to normal gland, VDR protein expression was retained in differentiated human pleomorphic adenoma (PA) but was undetectable in undifferentiated PA and in carcinomas, suggesting deregulation of VDR during salivary cancer progression. To gain insight into the potential role of VDR in salivary cancer, we assessed the effects of vitamin D in vivo and in vitro. Despite the presence of VDR in salivary gland, chronic dietary vitamin D restriction did not alter morphology of the salivary epithelium in C57/Bl6 mice. The localization of VDR in ductal epithelium prompted us to examine the effects of 1,25D in an established cell line (mSGc) derived from normal murine submandibular gland (SMG). This previously characterized cell line consists of multiple stem, progenitor and differentiated cell types as determined by mutually exclusive cellular expression of basal, ductal and myoepithelial markers. We demonstrated VDR expression and regulation of VDR target genes Vdr and Postn by 1,25D in mSGc, indicating functional ligand-mediated transcriptional activity. The effect of VDR signaling on epithelial differentiation markers was assessed by qPCR and IHC in mSGc cells treated with 1,25D. We found that 1,25D reduced mRNA expression of the basal cell progenitor marker keratin 5 (K5) and increased expression of the differentiated ductal cell marker keratin 7 (K7). Further, we found that 1,25D significantly decreased the number of proliferating cells, including proliferating K5⁺ cells. Characterization of cell cycle by Muse cytometry indicated 1,25D treatment decreased cells in S, G₂, and M phase. The inhibition of K5⁺ cell proliferation by 1,25D is of particular interest because K5⁺ basal cells contribute to a wide variety of salivary tumor types. Our studies suggest that 1,25D alters cancer-relevant progenitor and differentiation markers in the salivary gland.

New aspects of vitamin D metabolism and action - addressing the skin as source and target

Vitamin D has a key role in stimulating calcium absorption from the gut and promoting skeletal health, as well as many other important physiological functions. Vitamin D is produced in the skin. It is subsequently metabolized to its hormonally active form, 1,25-dihydroxyvitamin D (1,25(OH)₂D), by the 1-hydroxylase and catabolized by the 24-hydroxylase. In this Review, we pay special attention to the effect of mutations in these enzymes and their clinical manifestations. We then discuss the role of vitamin D binding protein in transporting vitamin D and its metabolites from their source to their targets, the free hormone hypothesis for cell entry and HSP70 for intracellular transport. This is followed by discussion of the vitamin D receptor (VDR) that mediates the cellular actions of 1,25(OH)₂D. Cell-specific recruitment of co-regulatory complexes by liganded VDR leads to changes in gene expression that result in distinct physiological actions by 1,25(OH)₂D, which are disrupted by mutations in the VDR. We then discuss the epidermis and hair follicle, to provide a non-skeletal example of a tissue that expresses VDR that not only makes vitamin D but also can metabolize it to its hormonally active form. This enables vitamin D to regulate epidermal differentiation and hair follicle cycling and, in so doing, to promote barrier function, wound healing and hair growth, while limiting cancer development.

Vitamin D: Newer Concepts of Its Metabolism and Function at the Basic and Clinical Level

The interest in vitamin D continues unabated with thousands of publications contributing to a vast and growing literature each year. It is widely recognized that the vitamin D receptor (VDR) and the enzymes that metabolize vitamin D are found in many cells, not just those involved with calcium and phosphate homeostasis. In this mini review I have focused primarily on recent studies that provide new insights into vitamin D metabolism, mechanisms of action, and clinical applications. In particular, I examine how mutations in vitamin D metabolizing enzymes—and new information on their regulation—links vitamin D metabolism into areas such as metabolism and diseases outside that of the musculoskeletal system. New information regarding the mechanisms governing the function of the VDR elucidates how this molecule can be so multifunctional in a cell-specific fashion. Clinically, the difficulty in determining vitamin D sufficiency for all groups is addressed, including a discussion of whether the standard measure of vitamin D sufficiency, total 25OHD (25 hydroxyvitamin) levels, may not be the best measure—at least by itself. Finally, several recent large clinical trials exploring the role of vitamin D supplementation in nonskeletal diseases are briefly reviewed, with an eye toward what questions they answered and what new questions they raised.

Genetic Predisposition to Colon and Rectal Adenocarcinoma Is Mediated by a Super-enhancer Polymorphism Coactivating CD9 and PLEKHG6

BACKGROUND

Genome-wide association studies (GWAS) have identified dozens of loci associated with colon and rectal adenocarcinoma risk. As tissue-specific super-enhancers (SE) play important roles in tumorigenesis, we systematically investigate SEs and inner variants in established GWAS loci to decipher the underlying biological mechanisms.

METHODS

Through a comprehensive bioinformatics analysis on multi-omics data, we screen potential single-nucleotide polymorphisms (SNP) in cancer-specific SEs, and then subject them to a two-stage case-control study containing 4,929 cases and 7,083 controls from the Chinese population. A series of functional assays, including reporter gene assays, electrophoretic mobility shift assays (EMSA), CRISPR-Cas9 genome editing, chromosome conformation capture (3C) assays, and cell proliferation experiments, are performed to characterize the variant's molecular consequence and target genes.

RESULTS

The SNP rs11064124 in 12p13.31 is found significantly associated with the risk of colon and rectal adenocarcinoma with an odds ratio (OR) of 0.87 [95% confidence interval (CI), 0.82-0.92, P = 8.67E-06]. The protective rs11064124-G weakens the binding affinity with vitamin D receptor (VDR) and increases the enhancer's activity and interactions with two target genes' promoters, thus coactivating the transcription of CD9 and PLEKHG6, which are both putative tumor suppressor genes for colon and rectal adenocarcinoma.

CONCLUSIONS

Our integrative study highlights an SE polymorphism rs11064124 and two susceptibility genes CD9 and PLEKHG6 in 12p13.31 for colon and rectal adenocarcinoma.

IMPACT

These findings suggest a novel insight for genetic pathogenesis of colon and rectal adenocarcinoma, involving transcriptional coactivation of diverse susceptibility genes via the SE element as a gene regulation hub.

DKK1 Induced by 1,25D3 Is Required for the Mineralization of Osteoblasts

1α,25-dihydroxyvitamin D3 (1,25D3), the most popular drug for osteoporosis treatment, drives osteoblast differentiation and bone mineralization. Wnt/β-catenin signaling is involved in commitment and differentiation of osteoblasts, but the role of the Dickkopf-related protein 1 (DKK1), a Wnt antagonist, in osteoblasts remains unknown. Here, we demonstrate the molecular mechanism of DKK1 induction by 1,25D3 and its physiological role during osteoblast differentiation. 1,25D3 markedly promoted the expression of both CCAAT/enhancer binding protein beta (C/EBPβ) and DKK1 at day 7 during osteoblast differentiation. Interestingly, mRNA and protein levels of C/EBPβ and DKK1 in osteoblasts were elevated by 1,25D3. We also found that C/EBPβ, in response to 1,25D3, directly binds to the human DKK1 promoter. Knockdown of C/EBPβ downregulated the expression of DKK1 in osteoblasts, which was partially reversed by 1,25D3. In contrast, overexpression of C/EBPβ upregulated DKK1 expression in osteoblasts, which was enhanced by 1,25D3. Furthermore, 1,25D3 treatment in osteoblasts stimulated secretion of DKK1 protein within the endoplasmic reticulum to extracellular. Intriguingly, blocking DKK1 attenuated calcified nodule formation in mineralized osteoblasts, but not ALP activity or collagen synthesis. Taken together, these observations suggest that 1,25D3 promotes the mineralization of osteoblasts through activation of DKK1 followed by an increase of C/EBPβ.

CCAAT/enhancer-binding protein beta (C/EBPβ) is an important mediator of 1,25 dihydroxyvitamin D3 (1,25D3)-induced receptor activator of nuclear factor kappa-B ligand (RANKL) expression in osteoblasts

Receptor activator of nuclear factor kappa B ligand (RANKL) expression in osteoblasts is regulated by 1,25-dihydroxyvitamin D3 (1,25D3). CCAAT/enhancer-binding protein beta (C/EBPβ) has been proposed to function as a transcription factor and upregulate RANKL expression, but it is still uncertain how C/EBPβ is involved in 1,25D3-induced RANKL expression of osteoblasts. 1,25D3 stimulation increased the expression of RANKL and C/EPBβ genes in osteoblasts and enhanced phosphorylation and stability of these proteins. Moreover, induction of RANKL expression by 1,25D3 in osteoblasts was downregulated upon knockdown of C/EBPβ. In contrast, C/EBPβ overexpression directly upregulated RANKL promoter activity and exhibited a synergistic effect on 1,25D3-induced RANKL expression. In particular, 1,25D3 treatment of osteoblasts increased C/EBPβ protein binding to the RANKL promoter. In conclusion, C/EBPβ is required for induction of RANKL by 1,25D3. [BMB Reports 2019; 52(6): 391-396].

Ezh2-dependent H3K27me3 modification dynamically regulates vitamin D3-dependent epigenetic control of CYP24A1 gene expression in osteoblastic cells

Epigenetic control is critical for the regulation of gene transcription in mammalian cells. Among the most important epigenetic mechanisms are those associated with posttranslational modifications of chromosomal histone proteins, which modulate chromatin structure and increased accessibility of promoter regulatory elements for competency to support transcription. A critical histone mark is trimethylation of histone H3 at lysine residue 27 (H3K27me3), which is mediated by Ezh2, the catalytic subunit of the polycomb group complex PRC2 to repress transcription. Treatment of cells with the active vitamin D metabolite 1,25(OH)₂ D₃ , results in transcriptional activation of the CYP24A1 gene, which encodes a 24-hydroxylase enzyme, that is, essential for physiological control of vitamin D3 levels. We report that the Ezh2-mediated deposition of H3K27me3 at the CYP24A1 gene promoter is a requisite regulatory component during transcriptional silencing of this gene in osteoblastic cells in the absence of 1,25(OH)₂ D₃ . 1,25(OH)₂ D₃ dependent transcriptional activation of the CYP24A1 gene is accompanied by a rapid release of Ezh2 from the promoter, together with the binding of the H3K27me3-specific demethylase Utx/Kdm6a and thereby subsequent erasing of the H3K27me3 mark. Importantly, we find that these changes in H3K27me3 enrichment at the CYP24A1 gene promoter are highly dynamic, as this modification is rapidly reacquired following the withdrawal of 1,25(OH)₂ D₃ .

Identification of a novel melatonin-binding nuclear receptor: Vitamin D receptor

Previous studies confirmed that melatonin regulates Runx2 expression but the mechanism is unclear. There is a direct interaction between Runx2 and the vitamin D receptor (VDR). Herein, we observed a direct interaction between melatonin and the VDR but not Runx2 using isothermal titration calorimetry. Furthermore, this direct binding was detected only in the C-terminal ligand binding domain (LBD) of the VDR but not in the N-terminal DNA-binding domain (DBD) or the hinge region. Spectrophotometry indicated that melatonin and vitamin D3 (VD3) had similar uptake rates, but melatonin's uptake was significantly inhibited by VD3 until the concentration of melatonin was obviously higher than that of VD3 in a preosteoblastic cell line MC3T3-E1. GST pull-down and yeast two-hybrid assay showed that the interactive smallest fragments were on the 319-379 position of Runx2 and the N-terminus 110-amino acid DBD of the VDR. Electrophoretic mobility shift assay (EMSA) demonstrated that Runx2 facilitated the affinity between the VDR and its specific DNA substrate, which was further documented by a fluorescent EMSA assay where Cy3 labeled Runx2 co-localized with the VDR-DNA complex. Another fluorescent EMSA assay confirmed that the binding of the VDR to Runx2 was significantly enhanced with an increasing concentrations of the VDR, especially in the presence of melatonin; it was further documented using a co-immunoprecipitation assay that this direct interaction was markedly enhanced by melatonin treatment in the MC3T3-E1 cells. Thus, the VDR is a novel melatonin-binding nuclear receptor, and melatonin indirectly regulates Runx2 when it directly binds to the LBD and the DBD of the VDR, respectively.

Switches in histone modifications epigenetically control vitamin D3-dependent transcriptional upregulation of the CYP24A1 gene in osteoblastic cells

In bone cells vitamin D dependent regulation of gene expression principally occurs through modulation of gene transcription. Binding of the active vitamin D metabolite, 1,25-dihydroxy vitamin D3 (1,25(OH)₂ D₃ ) to the vitamin D receptor (VDR) induces conformational changes in its C-terminal domain enabling competency for interaction with physiologically relevant coactivators, including SRC-1. Consequently, regulatory complexes can be assembled that support intrinsic enzymatic activities with competency to posttranslationally modify chromatin histones at target genomic sequences to epigenetically alter transcription. Here we examine specific transitions in representation and/or enrichment of epigenetic histone marks during 1,25(OH)₂ D₃ mediated upregulation of CYP24A1 gene expression in osteoblastic cells. This gene encodes the 24-hydroxylase enzyme, essential for biological control of vitamin D levels. We demonstrate that as the CYP24A1 gene promoter remains transcriptionally silent, there is enrichment of H4R3me2s together with its "writing" enzyme PRMT5 and decreased abundance of the istone H3 and H4 acetylation, H3R17me2a, and H4R3me2a marks as well as of their corresponding "writers." Exposure of osteoblastic cells to 1,25(OH)₂ D₃ stimulates the recruitment of a VDR/SRC-1 containing complex to the CYP24A1 promoter to mediate increased H3/H4 acetylation. VDR/SRC-1 binding occurs concomitant with the release of PRMT5 and the recruitment of the arginine methyltransferases CARM1 and PRMT1 to catalyze the deposition of the H3R17me2a and H4R3me2a marks, respectively. Our results indicate that these dynamic transitions of histone marks at the CYP24A1 promoter, provide a "chromatin context" that is transcriptionally competent for activation of the CYP24A1 gene in osteoblastic cells in response to 1,25(OH)₂ D₃ .

A chromatin-based mechanism controls differential regulation of the cytochrome P450 gene Cyp24a1 in renal and non-renal tissues

Cytochrome P450 family 27 subfamily B member 1 (CYP27B1) and CYP24A1 function to maintain physiological levels of 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) in the kidney. Renal Cyp27b1 and Cyp24a1 expression levels are transcriptionally regulated in a highly reciprocal manner by parathyroid hormone (PTH), fibroblast growth factor 23 (FGF23), and 1,25(OH)₂D₃ In contrast, Cyp24a1 regulation in nonrenal target cells (NRTCs) is limited to induction by 1,25(OH)₂D₃ Herein, we used ChIP-Seq analyses of mouse tissues to identify regulatory regions within the Cyp24a1 gene locus. We found an extended region downstream of Cyp24a1 containing a cluster of sites, termed C24-DS1, binding PTH-sensitive cAMP-responsive element-binding protein (CREB) and a cluster termed C24-DS2 binding the vitamin D receptor (VDR). VDR-occupied sites were present in both the kidney and NRTCs, but pCREB sites were occupied only in the kidney. We deleted each segment in the mouse and observed that although the overt phenotypes of both cluster deletions were unremarkable, RNA analysis in the C24-DS1-deleted strain revealed a loss of basal renal Cyp24a1 expression, total resistance to FGF23 and PTH regulation, and secondary suppression of renal Cyp27b1; 1,25(OH)₂D₃ induction remained unaffected in all tissues. In contrast, loss of the VDR cluster in the C24-DS2-deleted strain did not affect 1,25(OH)₂D₃ induction of renal Cyp24a1 expression yet reduced but did not eliminate Cyp24a1 responses in NRTCs. We conclude that a chromatin-based mechanism differentially regulates Cyp24a1 in the kidney and NRTCs and is essential for the specific functions of Cyp24a1 in these two tissue types.

TGF-β1-stimulation of matrix metalloproteinase-13 expression by down-regulation of miR-203a-5p in rat osteoblasts

Transforming growth factor-beta1 (TGF-β1) is a pleiotropic and ubiquitous cytokine involved in bone development and bone remodeling. Matrix metalloproteinase-13 (MMP13) plays a role in the degradation of the extracellular matrix (ECM), and the regulation of this gene is critical in bone remodeling. We previously reported that TGF-β1 stimulates MMP13 expression in rat osteoblasts. Recently, studies have examined the regulation of bone metabolism by microRNAs (miRNAs) to determine their therapeutic potential in osteogenesis. Here, we assessed the effect of TGF-β1 on down-regulation of miRNAs that target MMP13 and stimulation of MMP13 expression in osteoblasts. We used in silico analysis and identified 11 specific miRNAs which directly target rat MMP13. Among these miRNAs, miR-203a-5p expression was significantly decreased by TGF-β1-treatment in rat osteoblasts. Transient transfection of a miR-203a-5p mimic into rat osteoblasts reduced MMP13 expression. A luciferase reporter assay confirmed a direct targeting of miR-miR-203a-5p with the 3' untranslated regions of the MMP13 gene. Hence, we suggest that TGF-β1 stimulated down-regulation of miR-203a-5p, resulting in the stimulation of MMP13 expression in rat osteoblasts. Thus, identification of the role of miR-203a-5p via TGF-β1 and MMP13 in bone remodeling indicated its potential as a biomarker or therapeutic agent for treating bone and bone-related diseases.

The Function of the Vitamin D Receptor and a Possible Role of Enhancer RNA in Epigenomic Regulation of Target Genes: Implications for Bone Metabolism

Vitamin D (VD) is essential for bone health, and VD or its analogues are widely used in clinics to ameliorate bone loss. The targets and mode of VD anti-osteoporotic actions appear to be different from those of other classes of drugs modulating bone remodeling. VD exerts its biological activities through the nuclear VD receptor (VDR)-mediated transcriptional regulation of target mRNA and non-coding RNA genes. VD-induced gene regulation involves epigenetic modifications of chromatin conformation at the target loci as well as reconfiguration of higher-order chromosomal organization through VDR-mediated recruitment of various regulatory factors. Enhancer RNAs (eRNA), a class of non-coding enhancer-derived RNAs, have recently emerged as VDR target gene candidates that act through reorganization of chromatin looping to induce enhancer-promoter interaction in activation of mRNA-encoding genes. This review outlines the molecular mechanisms of VD actions mediated by the VDR and suggests novel function of eRNAs in VDR transactivation.

The impact of the vitamin D-modulated epigenome on VDR target gene regulation

The micronutrient vitamin D significantly modulates the human epigenome via enhancing genome-wide the rate of accessible chromatin and vitamin D receptor (VDR) binding. This study focuses on histone marks of active chromatin at promoter and enhancer regions and investigates, whether these genomic loci are sensitive to vitamin D. The epigenome of THP-1 human monocytes contains nearly 23,000 sites with H3K4me3 histone modifications, 550 of which sites are significantly (p < 0.05) modulated by stimulation with the VDR ligand 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). H3K27ac histone modifications mark active chromatin and 2473 of 45,500 sites are vitamin D sensitive. The two types of ligand-dependent histone marks allow to distinguish promoter and enhancer regulation by vitamin D, respectively. Transcription start site overlap is the prime attribute of ligand-dependent H3K4me3 marks, while VDR co-location is the top ranking parameter describing 1,25(OH)₂D₃-sensitive H3K27ac marks at enhancers. A categorization of 1,25(OH)₂D₃-sensitive histone marks by machine learning algorithms - using the attributes overall peak strength and ligand inducibility - highlights 260 and 287 regions with H3K4me3 and H3K27ac modifications, respectively. These loci are found at the promoter regions of 59 vitamin D target genes and their associated enhancers. In this way, ligand-dependent histone marks provide a link of the effects of 1,25(OH)₂D₃ on the epigenome with previously reported mRNA expression changes of vitamin D target genes. In conclusion, the human epigenome responds also on the level of histone modifications to 1,25(OH)₂D₃ stimulation. This allows a more detailed understanding of vitamin D target gene regulation.

Genomic Effects of the Vitamin D Receptor: Potentially the Link between Vitamin D, Immune Cells, and Multiple Sclerosis

Vitamin D has a plethora of functions that are important for the maintenance of general health and in particular, the functional integrity of the immune system, such as promoting an anti-inflammatory cytokine profile and reducing the Treg/Th17 ratio. Multiple sclerosis (MS) is a chronic, inflammatory, and neurodegenerative central nervous system (CNS) disorder of probable autoimmune origin. MS is characterized by recurring or progressive demyelination and degeneration of the CNS due in part to a misguided immune response to as yet undefined (CNS) antigens, potentially including myelin basic protein and proteolipid protein. MS has also been shown to be associated significantly with environmental factors such as the lack of vitamin D. The role of vitamin D in the pathogenesis and progression of MS is complex. Recent genetic studies have shown that various common MS-associated risk-single-nucleotide polymorphisms (SNPs) are located within or in the vicinity of genes associated with the complex metabolism of vitamin D. The functional aspects of these genetic associations may be explained either by a direct SNP-associated loss- or gain-of-function in a vitamin D-associated gene or due to a change in the regulation of gene expression in certain immune cell types. The development of new genetic tools using next-generation sequencing: e.g., chromatin immunoprecipitation sequencing (ChIP-seq) and the accompanying rapid progress of epigenomics has made it possible to recognize that the association between vitamin D and MS could be based on the extensive and characteristic genomic binding of the vitamin D receptor (VDR). Therefore, it is important to analyze comprehensively the spatiotemporal VDR binding patterns that have been identified using ChIP-seq in multiple immune cell types to reveal an integral profile of genomic VDR interaction. In summary, the aim of this review is to connect genomic effects vitamin D has on immune cells with MS and thus, to contribute to a better understanding of the influence of vitamin D on the etiology and the pathogenesis of this complex autoimmune disease.

The impact of VDR expression and regulation in vivo

The vitamin D receptor (VDR) mediates the pleiotropic biological actions of 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). These actions include orchestration of mineral homeostasis which is coordinated by the kidney, intestine, bone and parathyroid gland wherein the VDR transcriptionally regulates expression of the genes involved in this complex process. Mutations in human VDR (hVDR) cause hereditary vitamin D resistant rickets, a genetic syndrome characterized by hypocalcemia, hyperparathyroidism and rickets resulting from dysregulation of mineral homeostasis. Expression of the VDR is regulated by external stimuli in a tissue-specific manner. However, the mechanisms of this tissue-specificity remain unclear. Studies also suggest that phosphorylation of hVDR at serine 208 impacts the receptor's transcriptional activity. These experiments were conducted in vitro, however, and therefore limited in their conclusions. In this report, we summarize (1) our most recently updated ChIP-seq data from mouse tissues to identify regulatory regions responsible for the tissues-specific regulation of the VDR and (2) our studies to understand the mechanism of hormonal regulation of Vdr expression in bone and kidney in vivo using transgenic mouse strains generated by mouse mini-genes that contain comprehensive genetic information capable of recapitulating endogenous Vdr gene regulation and expression. We also defined the functional human VDR gene locus in vivo by using a human mini-gene comparable to that in the mouse to generate a humanized VDR mouse strain in which the receptor is expressed at normal levels (normal expressor). The present report also shows that a humanized mouse model in which the VDR is expressed at levels about 10-fold lower than the normal expressor mouse rescued the VDR-null phenotype despite its reduced transcriptional activity relative to wildtype expression. We also generated an additional humanized mouse model expressing hVDR bearing a mutation converting serine 208 to alanine (hVDR-S208A). In spite of the mutation, target gene expression induced by the ligand was unchanged relative to a mouse strain expressing comparable levels of wildtype hVDR. Further characterization also showed that serum calcium and parathyroid hormone levels were normal and alopecia was not observed in this hVDR-S208A mouse strain as well. Taken together, our in vivo studies using ChIP-seq analyses and the mini-gene transgenic mice improve our understanding of the tissue-specific regulatory mechanisms of controlling VDR expression and the mechanisms of action of the VDR.

Pathway analysis of transcriptomic data shows immunometabolic effects of vitamin D

Unbiased genomic screening analyses have highlighted novel immunomodulatory properties of the active form of vitamin D, 1,25-dihydroxyvitamin D (1,25(OH)₂D). However, clearer interpretation of the resulting gene expression data is limited by cell model specificity. The aim of the current study was to provide a broader perspective on common gene regulatory pathways associated with innate immune responses to 1,25(OH)₂D, through systematic re-interrogation of existing gene expression databases from multiple related monocyte models (the THP-1 monocytic cell line (THP-1), monocyte-derived dendritic cells (DCs) and monocytes). Vitamin D receptor (VDR) expression is common to multiple immune cell types, and thus, pathway analysis of gene expression using data from multiple related models provides an inclusive perspective on the immunomodulatory impact of vitamin D. A bioinformatic workflow incorporating pathway analysis using PathVisio and WikiPathways was utilized to compare each set of gene expression data based on pathway-level context. Using this strategy, pathways related to the TCA cycle, oxidative phosphorylation and ATP synthesis and metabolism were shown to be significantly regulated by 1,25(OH)₂D in each of the repository models (Z-scores 3.52-8.22). Common regulation by 1,25(OH)₂D was also observed for pathways associated with apoptosis and the regulation of apoptosis (Z-scores 2.49-3.81). In contrast to the primary culture DC and monocyte models, the THP-1 myelomonocytic cell line showed strong regulation of pathways associated with cell proliferation and DNA replication (Z-scores 6.1-12.6). In short, data presented here support a fundamental role for active 1,25(OH)₂D as a pivotal regulator of immunometabolism.

Vitamin D Genomics: From In Vitro to In Vivo

The vitamin D₃ metabolite 1α,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃] is the exclusive high-affinity ligand of the vitamin D receptor (VDR), a transcription factor with direct effects on gene expression. Transcriptome- and epigenome-wide data obtained in THP-1 human monocytes are the basis of the chromatin model of vitamin D signaling. The model describes, how VDR's spatio-temporal binding profile provides key insight into the pleiotropic action of vitamin D. The transcription of some 300 primary target genes is significantly modulated through the action of genomic VDR binding sites in concert with the pioneer transcription factor PU.1 and the chromatin organizer CTCF. In parallel, the short-term vitamin D intervention study VitDbol (NCT02063334) was designed, in order to extrapolate insight into vitamin D signaling from in vitro to in vivo. Before and 24 h after a vitamin D₃ bolus chromatin and RNA were prepared from peripheral blood mononuclear cells for epigenome- and transcriptome-wide analysis. The study subjects showed a personalized response to vitamin D and could be distinguished into high, mid, and low responders. Comparable principles of vitamin D signaling were identified in vivo and in vitro concerning target gene responses as well as changes in chromatin accessibility. In conclusion, short-term vitamin D supplementation studies represent a new type of safe in vivo investigations demonstrating that vitamin D and its metabolites have direct effects on the human epigenome and modulate the response of the transcriptome in a personalized fashion.

Biology and Mechanisms of Action of the Vitamin D Hormone

The central role of hormonal 1,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃] is to regulate calcium and phosphorus homeostasis via actions in intestine, kidney, and bone. These and other actions in many cell types not involved in mineral metabolism are mediated by the vitamin D receptor. Recent studies using genome-wide scale techniques have extended fundamental ideas regarding vitamin D-mediated control of gene expression while simultaneously revealing a series of new concepts. This article summarizes the current view of the biological actions of the vitamin D hormone and focuses on new concepts that drive the understanding of the mechanisms through which vitamin D operates.

Tet-Mediated DNA Demethylation Is Required for SWI/SNF-Dependent Chromatin Remodeling and Histone-Modifying Activities That Trigger Expression of the Sp7 Osteoblast Master Gene during Mesenchymal Lineage Commitment

Here we assess histone modification, chromatin remodeling, and DNA methylation processes that coordinately control the expression of the bone master transcription factor Sp7 (osterix) during mesenchymal lineage commitment in mammalian cells. We find that Sp7 gene silencing is mediated by DNA methyltransferase1/3 (DNMT1/3)-, histone deacetylase 1/2/4 (HDAC1/2/4)-, Setdb1/Suv39h1-, and Ezh1/2-containing complexes. In contrast, Sp7 gene activation involves changes in histone modifications, accompanied by decreased nucleosome enrichment and DNA demethylation mediated by SWI/SNF- and Tet1/Tet2-containing complexes, respectively. Inhibition of DNA methylation triggers changes in the histone modification profile and chromatin-remodeling events leading to Sp7 gene expression. Tet1/Tet2 silencing prevents Sp7 expression during osteoblast differentiation as it impairs DNA demethylation and alters the recruitment of histone methylase (COMPASS)-, histone demethylase (Jmjd2a/Jmjd3)-, and SWI/SNF-containing complexes to the Sp7 promoter. The dissection of these interconnected epigenetic mechanisms that govern Sp7 gene activation reveals a hierarchical process where regulatory components mediating DNA demethylation play a leading role.

Class 3 semaphorins are transcriptionally regulated by 1,25(OH)2D3 in osteoblasts

The vitamin D endocrine system is essential for calcium metabolism and skeletal integrity. 1,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃] regulates bone mineral homeostasis and acts directly on osteoblasts. In the present study we characterized the transcriptional regulation of the class 3 semaphorin (Sema3) gene family by 1,25(OH)₂D₃ in osteoblastic cells. Class 3 semaphorins are secreted proteins that regulate cell growth, morphology and migration, and were recently shown to be involved in bone homeostasis. In ST2, MC3T3-E1 and primary calvarial osteoblast cell cultures we found that all members of the Sema3 gene family were expressed, and that Sema3e and Sema3f were the most strongly induced 1,25(OH)₂D₃ target genes among the studied cell types. In addition, transcription of Sema3b and Sema3c was upregulated, whereas Sema3d and Sema3g was downregulated by 1,25(OH)₂D₃ in different osteoblastic cells. Chromatin immunoprecipitation analysis linked to DNA sequencing (ChIP-seq analysis) revealed the presence of the vitamin D receptor at multiple genomic loci in the proximity of Sema3 genes, demonstrating that the genes are primary 1,25(OH)₂D₃ targets. Furthermore, we showed that recombinant SEMA3E and SEMA3F protein were able to inhibit osteoblast proliferation. However, recombinant SEMA3s did not affect ST2 cell migration. The expression of class 3 semaphorins in osteoblasts together with their regulation by 1,25(OH)₂D₃ suggests that these genes, involved in the regulation of bone homeostasis, are additional mediators for 1,25(OH)₂D₃ signaling in osteoblasts.

Molecular endocrinology of vitamin D on the epigenome level

The molecular endocrinology of vitamin D is based on the facts that i) its metabolite 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) is the high affinity ligand of the nuclear receptor vitamin D receptor (VDR) and ii) the transcription factor VDR is the unique target of 1,25(OH)₂D₃ in the nucleus. Short-term alterations of the epigenome are primarily changes in the post-translational modification status of nucleosome-forming histone proteins, the consequences of which are i) a local increase or decrease in chromatin accessibility and ii) the activation or repression of gene transcription. Vitamin D has via VDR a direct effect on the expression of several hundred primary target genes implying numerous effects on the epigenome. Next-generation sequencing methods, such as ChIP-seq and FAIRE-seq, were applied to cellular model systems of vitamin D signaling, such as THP-1 human monocytes, and provided data for a chromatin model of vitamin D signaling. Key points of this model are that i) in the absence of ligand VDR binds to a limited number of loci within accessible chromatin, ii) a stimulation with ligand increases the number of DNA-bound VDR molecules, iii) VDR's access to genomic DNA is supported by pioneer factors, such as PU.1 in monocytes, iv) VDR binding leads to local opening of chromatin and v) the binding strength of topologically associating domain anchor forming CCCTC-binding factor sites upstream and downstream of prominent VDR binding sites is changing in response to ligand stimulation. This model provides the present basis of the molecular endocrinology of vitamin D and will be in future refined by the integration of vitamin D-sensitive chromatin markers and other genome-wide data, such as the 1,25(OH)₂D₃-sensitive binding of co-factors, chromatin modifying enzymes and chromatin remodeling proteins.

Vitamin D endocrinology of bone mineralization

Bone is a dynamic tissue that is strongly influenced by endocrine factors to restore the balance between bone resorption and bone formation. Bone formation involves the mineralization of the extracellular matrix formed by osteoblasts. In this process the role of vitamin D (1α,25(OH)₂D₃) is both direct and indirect. The direct effects are enabled via the Vitamin D Receptor (VDR); the outcome is dependent on the presence of other factors as well as origin of the osteoblasts, treatment procedures and species differences. Vitamin D stimulates mineralization of human osteoblasts but is often found inhibitory for mineralization of murine osteoblasts. In this review we will overview the current knowledge of the role of the vitamin D endocrine system in controlling the mineralization process in bone.

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.

Comparative analysis of osteoblast gene expression profiles and Runx2 genomic occupancy of mouse and human osteoblasts in vitro

Fast progress of the next generation sequencing (NGS) technology has allowed global transcriptional profiling and genome-wide mapping of transcription factor binding sites in various cellular contexts. However, limited number of replicates and high amount of data processing may weaken the significance of the findings. Comparative analyses of independent data sets acquired in the different laboratories would greatly increase the validity of the data. Runx2 is the key transcription factor regulating osteoblast differentiation and bone formation. We performed a comparative analysis of three published Runx2 data sets of chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) analysis in osteoblasts from mouse and human origin. Moreover, we assessed the similarity of the corresponding transcription data of these studies available online. The ChIP-seq data analysis confirmed general features of Runx2 binding, including location at genic vs intergenic regions and abundant Runx2 binding on promoters of the highly expressed genes. We also found high frequency of Runx2 DNA binding without a consensus Runx2 motif at the binding site. Importantly, mouse and human Runx2 showed moderately similar binding patterns in terms of peak-associated closest genes and their associated genomic ontology (GO) pathways. Accordingly, the gene expression profiles were highly similar and osteoblastic phenotype was prominent in the differentiated stage in both species. In conclusion, ChIP-seq method shows good reproducibility in the context of mature osteoblasts, and mouse and human osteoblast models resemble each other closely in Runx2 binding and in gene expression profiles, supporting the use of these models as adequate tools in studying osteoblast differentiation.

Histone H4 Methyltransferase Suv420h2 Maintains Fidelity of Osteoblast Differentiation

Osteogenic lineage commitment and progression is controlled by multiple signaling pathways (e.g., WNT, BMP, FGF) that converge on bone-related transcription factors. Access of osteogenic transcription factors to chromatin is controlled by epigenetic regulators that generate post-translational modifications of histones ("histone code"), as well as read, edit and/or erase these modifications. Our understanding of the biological role of epigenetic regulators in osteoblast differentiation remains limited. Therefore, we performed next-generation RNA sequencing (RNA-seq) and established which chromatin-related proteins are robustly expressed in mouse bone tissues (e.g., fracture callus, calvarial bone). These studies also revealed that cells with increased osteogenic potential have higher levels of the H4K20 methyl transferase Suv420h2 compared to other methyl transferases (e.g., Suv39h1, Suv39h2, Suv420h1, Ezh1, Ezh2). We find that all six epigenetic regulators are transiently expressed at different stages of osteoblast differentiation in culture, with maximal mRNAs levels of Suv39h1 and Suv39h2 (at day 3) preceding maximal expression of Suv420h1 and Suv420h2 (at day 7) and developmental stages that reflect, respectively, early and later collagen matrix deposition. Loss of function analysis of Suv420h2 by siRNA depletion shows loss of H4K20 methylation and decreased expression of bone biomarkers (e.g., alkaline phosphatase/Alpl) and osteogenic transcription factors (e.g., Sp7/Osterix). Furthermore, Suv420h2 is required for matrix mineralization during osteoblast differentiation. We conclude that Suv420h2 controls the H4K20 methylome of osteoblasts and is critical for normal progression of osteoblastogenesis. J. Cell. Biochem. 118: 1262-1272, 2017. © 2016 Wiley Periodicals, Inc.

Control of the Osteoblast Lineage by Mitogen-Activated Protein Kinase Signaling

PURPOSE OF THE REVIEW

This review will provide a timely assessment of MAP kinase actions in bone development and homeostasis with particular emphasis on transcriptional control of the osteoblast lineage.

RECENT FINDINGS

ERK and p38 MAP kinases function as transducers of signals initiated by the extracellular matrix, mechanical loading, TGF-β, BMPs and FGF2. MAPK signals may also affect and/or interact with other important pathways such as WNT and HIPPO. ERK and p38 MAP kinase pathways phosphorylate specific osteogenic transcription factors including RUNX2, Osterix, ATF4 and DLX5. For RUNX2, phosphorylation at specific serine residues initiates epigenetic changes in chromatin necessary for decondensation and increased transcription. MAPK also suppresses marrow adipogenesis by phosphorylating and inhibiting PPARγ, which may explain the well-known relationship between reduced skeletal loading and marrow fat accumulation.

SUMMARY

MAPKs transduce signals from the extracellular environment to the nucleus allowing bone cells to respond to changes in hormonal/growth factor signaling and mechanical loading thereby optimizing bone structure to meet physiological and mechanical needs of the body.

MAP Kinase-Dependent RUNX2 Phosphorylation Is Necessary for Epigenetic Modification of Chromatin During Osteoblast Differentiation

RUNX2, an essential transcription factor for osteoblast differentiation and bone formation is activated by ERK/MAP kinase-dependent phosphorylation. However, relationship between these early events and specific epigenetic modifications of chromatin during osteoblast differentiation have not been previously examined. Here, we explore these relationships using chromatin immunoprecipitation (ChIP) to detect chromatin modifications in RUNX2-binding regions of Bglap2 and Ibsp. Growth of MC3T3-E1c4 preosteoblast cells in differentiation conditions rapidly induced Bglap2 and lbsp mRNAs. For both genes, osteogenic stimulation increased chromatin-bound P-ERK, P-RUNX2, p300, and RNA polymerase II as well as histone H3K9 and H4K5 acetylation. The level of H3K4 di-methylation, another gene activation-associated histone mark, also increased. In contrast, levels of the gene repressive marks, H3K9 mono-, di-, and tri-methylation in the same regions were reduced. Inhibition of MAP kinase signaling blocked differentiation-dependent chromatin modifications and Bglap2 and Ibsp expression. To evaluate the role of RUNX2 phosphorylation in these responses, RUNX2-deficient C3H10T1/2 cells were transduced with adenovirus encoding wild type or phosphorylation site mutant RUNX2 (RUNX2 S301A/S319A). Wild type RUNX2, but not the non-phosphorylated mutant, increased H3K9 and H4K5 acetylation as well as chromatin-associated P-ERK, p300, and polymerase II. Thus, RUNX2 phosphorylation is necessary for subsequent epigenetic changes required for osteoblast gene expression. Taken together, this study reveals a molecular mechanism through which osteogenic genes are controlled by a MAPK and P-RUNX2-dependent process involving epigenetic modifications of specific promoter regions. J. Cell. Physiol. 232: 2427-2435, 2017. © 2016 Wiley Periodicals, Inc.

Unique Regenerative Mechanism to Replace Bone Lost During Dietary Bone Depletion in Weanling Mice

The present study was undertaken to determine the mechanism whereby calcitropic hormones and mesenchymal stem cell progeny changes are involved in bone repletion, a regenerative bone process that restores the bone lost to calcium deficiency. To initiate depletion, weanling mice with a mixed C57BL/6 (75%) and CD1 (25%) genetic background were fed a calcium-deficient diet (0.01%) for 14 days. For repletion, the mice were fed a control diet containing 1.2% calcium for 14 days. Depletion decreased plasma calcium and increased plasma parathyroid hormone, 1,25(OH)2D (calcitriol), and C-terminal telopeptide of type I collagen. These plasma parameters quickly returned toward normal on repletion. The trabecular bone volume and connectivity decreased drastically during depletion but were completely restored by the end of repletion. This bone repletion process largely resulted from the development of new bone formation. When bromodeoxyuridine (BrdU) was administered in the middle of depletion for 3 days and examined by fluorescence-activated cell sorting at 7 days into repletion, substantial increases in BrdU incorporation were seen in several CD105 subsets of cells of osteoblastic lineage. When BrdU was administered on days 1 to 3 of repletion and examined 11 days later, no increases in BrdU were seen in these subsets. Additionally, osteocytes that stained positively for BrdU were increased during depletion. In conclusion, the results of the present study have established a unique regenerative mechanism to initiate bone repair during the bone insult. Calcium homeostatic mechanisms and the bone repletion mechanism are opposing functions but are simultaneously orchestrated such that both endpoints are optimized. These results have potential clinical relevance for disease entities such as type 2 osteoporosis.

Vitamin D Impacts the Expression of Runx2 Target Genes and Modulates Inflammation, Oxidative Stress and Membrane Vesicle Biogenesis Gene Networks in 143B Osteosarcoma Cells

Osteosarcoma (OS) is an aggressive malignancy of bone affecting children, adolescents and young adults. Understanding vitamin D metabolism and vitamin D regulated genes in OS is an important aspect of vitamin D/cancer paradigm, and in evaluating vitamin D as adjuvant therapy for human OS. Vitamin D treatment of 143B OS cells induced significant and novel changes in the expression of genes that regulate: (a) inflammation and immunity; (b) formation of reactive oxygen species, metabolism of cyclic nucleotides, sterols, vitamins and mineral (calcium), quantity of gap junctions and skeletogenesis; (c) bone mineral density; and (d) cell viability of skeletal cells, aggregation of bone cancer cells and exocytosis of secretory vesicles. Ingenuity pathway analysis revealed significant reduction in Runx2 target genes such as fibroblast growth factor -1, -12 (FGF1 and FGF12), bone morphogenetic factor-1 (BMP1), SWI/SNF related, matrix associated actin dependent regulator of chromatin subfamily a, member 4 (SMARCA4), Matrix extracellular phosphoglycoprotein (MEPE), Integrin, β4 (ITGBP4), Matrix Metalloproteinase -1, -28 (MMP1 and MMP28), and signal transducer and activator of transcription-4 (STAT4) in vitamin D treated 143B OS cells. These genes interact with the inflammation, oxidative stress and membrane vesicle biogenesis gene networks. Vitamin D not only inhibited the expression of Runx2 target genes MMP1, MMP28 and kallikrein related peptidase-7 (KLK7), but also migration and invasion of 143B OS cells. Vitamin D regulated Runx2 target genes or their products represent potential therapeutic targets and laboratory biomarkers for applications in translational oncology.