Rapid genomic changes by mineralotropic hormones and kinase SIK inhibition drive coordinated renal Cyp27b1 and Cyp24a1 expression via CREB modules.

Vitamin D metabolism centers on kidney regulation of Cyp27b1 by mineralotropic hormones, including induction by parathyroid hormone (PTH), suppression by fibroblast growth factor 23 (FGF23) and 1,25-dihydroxyvitamin D3 (1,25(OH)2D3), and reciprocal regulations for Cyp24a1. This coordinated genomic regulation results in production of endocrine 1,25(OH)2D3, which, together with PTH and FGF23, controls mineral homeostasis. However, how these events are coordinated is unclear. Here, using in vivo chromatin immunoprecipitation sequencing in mouse kidney, we demonstrate that PTH activation rapidly induces increased recruitment of phosphorylated (p-133) CREB (pCREB) and its coactivators, CBP (CREB-binding protein) and CRTC2 (CREB-regulated transcription coactivator 2), to previously defined kidney-specific M1 and M21 enhancers near the Cyp27b1 gene. At distal enhancers of the Cyp24a1 gene, PTH suppression dismisses CBP with only minor changes in pCREB and CRTC2 occupancy, all of which correlate with decreased genomic activity and reduced transcripts. Treatment of mice with salt-inducible kinase inhibitors (YKL-05-099 and SK-124) yields rapid genomic recruitment of CRTC2 to Cyp27b1, limited interaction of CBP, and a transcriptional response for both Cyp27b1 and Cyp24a1 that mirrors the actions of PTH. Surprisingly, we find that 1,25(OH)2D3 suppression increases the occupancy of CRTC2 in the M1 enhancer, a novel observation for CRTC2 and 1,25(OH)2D3 action. Suppressive actions of 1,25(OH)2D3 and FGF23 at the Cyp27b1 gene are associated with reduced CBP recruitment at these CREB-module enhancers that disrupts full PTH induction. Our findings show that CRTC2 contributes to transcription of both Cyp27b1 and Cyp24a1, demonstrate salt-inducible kinase inhibition as a key modulator of vitamin D metabolism, and provide molecular insight into the coordinated mechanistic actions of PTH, FGF23, and 1,25(OH)2D3 in the kidney that regulate mineral homeostasis.

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.

Targeted genomic deletions identify diverse enhancer functions and generate a kidney-specific, endocrine-deficient Cyp27b1 pseudo-null mouse.

Vitamin D3 is terminally bioactivated in the kidney to 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) via cytochrome P450 family 27 subfamily B member 1 (CYP27B1), whose gene is regulated by parathyroid hormone (PTH), fibroblast growth factor 23 (FGF23), and 1,25(OH)2D3 Our recent genomic studies in the mouse have revealed a complex kidney-specific enhancer module within the introns of adjacent methyltransferase-like 1 (Mettl1) and Mettl21b that mediate basal and PTH-induced expression of Cyp27b1 and FGF23- and 1,25(OH)2D3-mediated repression. Gross deletion of these segments in mice has severe effects on Cyp27b1 regulation and skeletal phenotype but does not affect Cyp27b1 expression in nonrenal target cells (NRTCs). Here, we report a bimodal activity in the Mettl1 intronic enhancer with components responsible for PTH-mediated Cyp27b1 induction and 1,25(OH)2D3-mediated repression and additional activities, including FGF23 repression, within the Mettl21b enhancers. Deletion of both submodules eliminated basal Cyp27b1 expression and regulation in the kidney, leading to systemic and skeletal phenotypes similar to those of Cyp27b1-null mice. However, basal expression and lipopolysaccharide-induced regulation of Cyp27b1 in NRTCs was unperturbed. Importantly, dietary normalization of calcium, phosphate, PTH, and FGF23 rescued the skeletal phenotype of this mutant mouse, creating an ideal in vivo model to study nonrenal 1,25(OH)2D3 production in health and disease. Finally, we confirmed a conserved chromatin landscape in human kidney that is similar to that in mouse. These findings define a finely balanced homeostatic mechanism involving PTH and FGF23 together with protection from 1,25(OH)2D3 toxicity that is responsible for both adaptive vitamin D metabolism and mineral regulation.

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 D3 (1,25(OH)2D3) 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)2D3 In contrast, Cyp24a1 regulation in nonrenal target cells (NRTCs) is limited to induction by 1,25(OH)2D3 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)2D3 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)2D3 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.

A kidney-specific genetic control module in mice governs endocrine regulation of the cytochrome P450 gene Cyp27b1 essential for vitamin D3 activation.

The vitamin D endocrine system regulates mineral homeostasis through its activities in the intestine, kidney, and bone. Terminal activation of vitamin D3 to its hormonal form, 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3), occurs in the kidney via the cytochrome P450 enzyme CYP27B1. Despite its importance in vitamin D metabolism, the molecular mechanisms underlying the regulation of the gene for this enzyme, Cyp27b1, are unknown. Here, we identified a kidney-specific control module governed by a renal cell-specific chromatin structure located distal to Cyp27b1 that mediates unique basal and parathyroid hormone (PTH)-, fibroblast growth factor 23 (FGF23)-, and 1,25(OH)2D3-mediated regulation of Cyp27b1 expression. Selective genomic deletion of key components within this module in mice resulted in loss of either PTH induction or FGF23 and 1,25(OH)2D3 suppression of Cyp27b1 gene expression; the former loss caused a debilitating skeletal phenotype, whereas the latter conferred a quasi-normal bone mineral phenotype through compensatory homeostatic mechanisms involving Cyp24a1 We found that Cyp27b1 is also expressed at low levels in non-renal cells, in which transcription was modulated exclusively by inflammatory factors via a process that was unaffected by deletion of the kidney-specific module. These results reveal that differential regulation of Cyp27b1 expression represents a mechanism whereby 1,25(OH)2D3 can fulfill separate functional roles, first in the kidney to control mineral homeostasis and second in extra-renal cells to regulate target genes linked to specific biological responses. Furthermore, we conclude that these mouse models open new avenues for the study of vitamin D metabolism and its involvement in therapeutic strategies for human health and disease.

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.

Selective Distal Enhancer Control of the Mmp13 Gene Identified through Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) Genomic Deletions.

Matrix metalloproteinase 13 (Mmp13, collagenase-3) plays an essential role in bone metabolism and mineral homeostasis. It is regulated by numerous factors, including BMP-2, parathyroid hormone, and 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3), through transcription factors such as Runt-related transcription factor 2 (RUNX2), CCAAT/enhancer-binding protein ß (C/EBPß), OSX, and vitamin D receptor (VDR). During osteoblast maturation, the basal expression of Mmp13 and its sensitivity to 1,25(OH)2D3 are strikingly increased. In this report, ChIP-sequencing analysis in mouse preosteoblasts revealed that the Mmp13 gene was probably regulated by three major enhancers located -10, -20, and -30 kb upstream of the gene promoter, occupied by activated VDR and prebound C/EBPß and RUNX2, respectively. Initially, bacterial artificial chromosome clone recombineering and traditional mutagenesis defined binding sites for VDR and RUNX2. We then employed a CRISPR/Cas9 gene editing approach to delete the -10 and -30 kb Mmp13 enhancers, a region proximal to the promoter, and VDR or RUNX2. VDR-mediated up-regulation of Mmp13 transcription was completely abrogated upon removal of the -10 kb enhancer, resulting in a 1,25(OH)2D3-directed repression of Mmp13. Deletion of either the -30 kb enhancer or RUNX2 resulted in a complete loss of basal transcript activity and a ChIP-identified destabilization of the chromatin enhancer environment and factor binding. Whereas enhancer deletions only affected Mmp13 expression, the RUNX2 deletion led to changes in gene expression, a reduction in cellular proliferation, and an inability to differentiate. We conclude that the Mmp13 gene is regulated via at least three specific distal enhancers that display independent activities yet are able to integrate response from multiple signaling pathways in a model of activation and suppression.

Mechanisms of Enhancer-mediated Hormonal Control of Vitamin D Receptor Gene Expression in Target Cells.

The biological actions of 1,25-dihydroxyvitamin D3 (1,25(OH)2D3) are mediated by the vitamin D receptor (VDR), whose expression in bone cells is regulated positively by 1,25(OH)2D3, retinoic acid, and parathyroid hormone through both intergenic and intronic enhancers. In this report, we used ChIP-sequencing analysis to confirm the presence of these Vdr gene enhancers in mesenchyme-derived bone cells and to describe the epigenetic histone landscape that spans the Vdr locus. Using bacterial artificial chromosome-minigene stable cell lines, CRISPR/Cas9 enhancer-deleted daughter cell lines, transient transfection/mutagenesis analyses, and transgenic mice, we confirmed the functionality of these bone cell enhancers in vivo as well as in vitro. We also identified VDR-binding sites across the Vdr gene locus in kidney and intestine using ChIP-sequencing analysis, revealing that only one of the bone cell-type enhancers bound VDR in kidney tissue, and none were occupied by the VDR in the intestine, consistent with weak or absent regulation by the 1,25(OH)2D3 hormone in these tissues, respectively. However, a number of additional sites of VDR binding unique to either kidney or intestine were present further upstream of the Vdr gene, suggesting the potential for alternative regulatory loci. Importantly, virtually all of these regions retained histone signatures consistent with those of enhancers and exhibited unique DNase I hypersensitivity profiles that reflected the potential for chromatin access. These studies define mechanisms associated with hormonal regulation of the Vdr and hint at the differential nature of VDR binding activity at the Vdr gene in different primary target tissues in vivo.

1,25-Dihydroxyvitamin D3 Controls a Cohort of Vitamin D Receptor Target Genes in the Proximal Intestine That Is Enriched for Calcium-regulating Components.

1,25-Dihydroxyvitamin D3 (1,25(OH)2D3) plays an integral role in calcium homeostasis in higher organisms through its actions in the intestine, kidney, and skeleton. Interestingly, although several intestinal genes are known to play a contributory role in calcium homeostasis, the entire caste of key components remains to be identified. To examine this issue, Cyp27b1 null mice on either a normal or a high calcium/phosphate-containing rescue diet were treated with vehicle or 1,25(OH)2D3 and evaluated 6 h later. RNA samples from the duodena were then subjected to RNA sequence analysis, and the data were analyzed bioinformatically. 1,25(OH)2D3 altered expression of large collections of genes in animals under either dietary condition. 45 genes were found common to both 1,25(OH)2D3-treated groups and were composed of genes previously linked to intestinal calcium uptake, including S100g, Trpv6, Atp2b1, and Cldn2 as well as others. An additional distinct network of 56 genes was regulated exclusively by diet. We then conducted a ChIP sequence analysis of binding sites for the vitamin D receptor (VDR) across the proximal intestine in vitamin D-sufficient normal mice treated with vehicle or 1,25(OH)2D3. The residual VDR cistrome was composed of 4617 sites, which was increased almost 4-fold following hormone treatment. Interestingly, the majority of the genes regulated by 1,25(OH)2D3 in each diet group as well as those found in common in both groups contained frequent VDR sites that likely regulated their expression. This study revealed a global network of genes in the intestine that both represent direct targets of vitamin D action in mice and are involved in calcium absorption.

1,25-Dihydroxyvitamin D regulates expression of the tryptophan hydroxylase 2 and leptin genes: implication for behavioral influences of vitamin D.

To investigate vitamin D-related control of brain-expressed genes, candidate vitamin D responsive elements (VDREs) at -7/-10 kb in human tryptophan hydroxylase (TPH)2 were probed. Both VDREs bound the vitamin D receptor (VDR)-retinoid X receptor (RXR) complex and drove reporter gene transcription in response to 1,25-dihydroxyvitamin D3 (1,25D). Brain TPH2 mRNA, encoding the rate-limiting enzyme in serotonin synthesis, was induced 2.2-fold by 10 nM 1,25D in human U87 glioblastoma cells and 47.8-fold in rat serotonergic RN46A-B14 cells. 1,25D regulation of leptin (Lep), encoding a serotoninlike satiety factor, was also examined. In mouse adipocytes, 1,25D repressed leptin mRNA levels by at least 84%, whereas 1,25D induced leptin mRNA 15.1-fold in human glioblastoma cells. Chromatin immunoprecipitation sequencing analysis of the mouse Lep gene in response to 1,25D revealed a cluster of regulatory sites (cis-regulatory module; CRM) at -28 kb that 1,25D-dependently docked VDR, RXR, C/EBPß, and RUNX2. This CRM harbored 3 VDREs and single C/EBPß and RUNX2 sites. Therefore, the expression of human TPH2 and mouse Lep are governed by 1,25D, potentially via respective VDREs located at -7/-10 kb and -28 kb. These results imply that vitamin D affects brain serotonin concentrations, which may be relevant to psychiatric disorders, such as autism, and may control leptin levels and affect eating behavior.

The RUNX2 cistrome in osteoblasts: characterization, down-regulation following differentiation, and relationship to gene expression.

RUNX2 is a transcription factor that is first expressed in early osteoblast-lineage cells and represents a primary determinant of osteoblastogenesis. While numerous target genes are regulated by RUNX2, little is known of sites on the genome occupied by RUNX2 or of the gene networks that are controlled by these sites. To explore this, we conducted a genome-wide analysis of the RUNX2 cistrome in both pre-osteoblastic MC3T3-E1 cells (POB) and their mature osteoblast progeny (OB), characterized the two cistromes and assessed their relationship to changes in gene expression. We found that although RUNX2 was widely bound to the genome in POB cells, this binding profile was reduced upon differentiation to OBs. Numerous sites were lost upon differentiation, new sites were also gained; many sites remained common to both cell states. Additional features were identified as well including location relative to potential target genes, abundance with respect to single genes, the frequent presence of a consensus TGTGGT RUNX2 binding motif, co-occupancy by C/EBPß and the presence of a typical epigenetic histone enhancer signature. This signature was changed quantitatively following differentiation. While RUNX2 binding sites were associated extensively with adjacent genes, the distal nature of the majority of these sites prevented assessment of whether they represented direct targets of RUNX2 action. Changes in gene expression, however, revealed an abundance of genes that contained RUNX2 binding sites and were regulated in concert. These studies establish a basis for further analysis of the role of RUNX2 activity and its function during osteoblast lineage maturation.

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

The biological effects of 1α,25-dihydroxyvitamin D3 (1,25 (OH)2D3) on osteoblast differentiation and function differ significantly depending upon the cellular state of maturation. To explore this phenomenon mechanistically, we examined the impact of 1,25(OH)2D3 on the transcriptomes of both pre-osteoblastic (POBs) and differentiated osteoblastic (OBs) MC3T3-E1 cells, and assessed localization of the vitamin D receptor (VDR) at sites of action on a genome-scale using ChIP sequence analysis. We observed that the 1,25(OH)2D3-induced transcriptomes of POBs and OBs were quantitatively and qualitatively different, supporting not only the altered biology observed but the potential for a change in VDR interaction at the genome as well. This idea was confirmed through discovery that VDR cistromes in POBs and OBs were also strikingly different. Depletion of VDR-binding sites in OBs, due in part to reduced VDR expression, was the likely cause of the loss of VDR-target gene interaction. Continued novel regulation by 1,25(OH)2D3, however, suggested that factors in addition to the VDR might also be involved. Accordingly, we show that transcriptomic modifications are also accompanied by changes in genome binding of the master osteoblast regulator RUNX2 and the chromatin remodeler CCAAT/enhancer-binding protein ß. Importantly, genome occupancy was also highlighted by the presence of epigenetic enhancer signatures that were selectively changed in response to both differentiation and 1,25(OH)2D3. The impact of VDR, RUNX2, and C/EBPß on osteoblast differentiation is exemplified by their actions at the Runx2 and Sp7 gene loci. We conclude that each of these mechanisms may contribute to the diverse actions of 1,25(OH)2D3 on differentiating osteoblasts.

Transcriptional regulation of the human TNFSF11 gene in T cells via a cell type-selective set of distal enhancers.

In addition to osteoblast lineage cells, the TNF-like factor receptor activator of NF-κB ligand (RANKL) is expressed in both B and T cells and may play a role in bone resorption. Rankl gene (Tnfsf11) expression in mouse T cells is mediated through multiple distal elements marked by increased transcription factor occupancy, histone tail acetylation, and RNA polymerase II recruitment. Little is known, however, of the regulation of human TNFSF11 in T cells. Accordingly, we examined the consequence of T cell activation on the expression of this factor both in Jurkat cells and in primary human T cells. We then explored the mechanism of this regulation by scanning over 400 kb of DNA surrounding the TNFSF11 locus for regulatory enhancers using ChIP-chip analysis. Histone H3/H4 acetylation enrichment identified putative regulatory regions located between -170 and -220 kb upstream of the human TNFSF11 TSS that we designated the human T cell control region (hTCCR). This region showed high sequence conservation with the mouse TCCR. Inhibition of MEK1/2 by U0126 resulted in decreased RANKL expression suggesting that stimulation through MEK1/2 was a prerequisite. ChIP-chip analysis also revealed that c-FOS was recruited to the hTCCR as well. Importantly, both the human TNFSF11 D5a/b (RLD5a/b) enhancer and segments of the hTCCR mediated robust inducible reporter activity following TCR activation. Finally, SNPs implicated in diseases characterized by dysregulated BMD co-localized to the hTCCR region. We conclude that the hTCCR region contains a cell-selective set of enhancers that plays an integral role in the transcriptional regulation of the TNFSF11 gene in human T cells.

Spatial detection and consequences of nonrenal calcitriol production as assessed by targeted mass spectrometry imaging.

The immune benefits of vitamin D3 supplementation beyond calcium and phosphate maintenance are highly clinically debated. Kidney expression of CYP27B1 is the source of endocrine, circulating 1,25(OH)2D3 (active form of vitamin D) that maintains serum calcium and phosphate. 1,25(OH)2D3 may also be made by the CYP27B1 enzyme in non-renal cells, like immune cells, in a process driven by cellular availability of 25(OH)D3 and inflammation. Due to the endocrine nature of 1,25(OH)2D3 in circulation, it is difficult to discern between these two sources. We recently created a regulatory deletion model of Cyp27b1 (M1/M21-DIKO) where mice have normal inflammatory-regulated Cyp27b1 expression in non-renal tissues (unlike global Cyp27b1-KO), but no expression within kidney. Here, utilizing on-tissue chemical derivatization and Matrix Assisted Laser Desorption Ionization-Mass Spectrometry Imaging (MALDI-MSI), we investigated the distribution of 1,25(OH)2D3 and 25(OH)D3 in the kidney, liver, spleen, and thymus. MALDI-MSI demonstrated increased 1,25(OH)2D3 in non-renal tissues such as the spleen after vitamin D3 supplementation in M1/M21-DIKO mice. Additionally, from this we found increased Il4 and decreased Tnfa in the spleen after vitamin D3 supplementation. Taken together, these data demonstrate non-renal production of 1,25(OH)2D3 in vivo and provide a consequence of vitamin D3 supplementation and non-renal 1,25(OH)2D3 production in cytokine changes.

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.

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.

Genome-wide analyses of gene expression profile identify key genes and pathways involved in skeletal response to phosphate and 1,25-dihydroxyvitamin D3 in vivo.

Phosphate (P) is an essential element involved in various biological actions, such as bone integrity, energy production, cell signaling and molecular component. P homeostasis is modulated by 4 main tissues; intestine, kidney, bone, and parathyroid gland, where 1,25-dihydroxyvitamin D3 (1,25(OH)2D3), parathyroid hormone and fibroblast growth factor 23 (FGF23) are produced and/or have an influence. In bone, serum P level modulates the production of FGF23 which then controls not only P excretion but also vitamin D metabolism in kidney in an endocrine manner. The hormonally active form of vitamin D, 1,25(OH)2D3, also has a significant effect on skeletal cells via its receptor, the vitamin D receptor, to control gene expression which mediates bone metabolism as well as mineral homeostasis. In this study, we adopted RNA-seq analysis to understand genome-wide skeletal gene expression regulation in response to P and 1,25(OH)2D3. We examined lumbar 5 vertebrae from the mice that were fed P deficient diet for a week followed by an acute high P diet for 3, 6, and 24 h as well as mice treated with 1,25(OH)2D3 intraperitoneally for 6 h. Further identification and exploration of the genes regulated by P and 1,25(OH)2D3 showed that P dynamically modulates the expression of skeletal genes involved in various biological processes while 1,25(OH)2D3 regulates genes highly related to bone metabolism. Our in vivo data were then compared with in vitro data that we previously obtained, which suggests that the gene expression profiles presented in this report mainly represent those of osteocytes. Interestingly, it was found that even though the skeletal response to P is distinguished from that to 1,25(OH)2D3, both factors have an effect on Wnt signaling pathway to modulate bone homeostasis. Taken together, this report presents genome-wide data that provide a foundation to understand molecular mechanisms by which skeletal cells respond to P and 1,25(OH)2D3.

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.

Vitamin D receptor cross-talk with p63 signaling promotes epidermal cell fate.

The vitamin D receptor with its ligand 1,25 dihydroxy vitamin D3 (1,25D3) regulates epidermal stem cell fate, such that VDR removal from Krt14 expressing keratinocytes delays re-epithelialization of epidermis after wound injury in mice. In this study we deleted Vdr from Lrig1 expressing stem cells in the isthmus of the hair follicle then used lineage tracing to evaluate the impact on re-epithelialization following injury. We showed that Vdr deletion from these cells prevents their migration to and regeneration of the interfollicular epidermis without impairing their ability to repopulate the sebaceous gland. To pursue the molecular basis for these effects of VDR, we performed genome wide transcriptional analysis of keratinocytes from Vdr cKO and control littermate mice. Ingenuity Pathway analysis (IPA) pointed us to the TP53 family including p63 as a partner with VDR, a transcriptional factor that is essential for proliferation and differentiation of epidermal keratinocytes. Epigenetic studies on epidermal keratinocytes derived from interfollicular epidermis showed that VDR is colocalized with p63 within the specific regulatory region of MED1 containing super-enhancers of epidermal fate driven transcription factor genes such as Fos and Jun. Gene ontology analysis further implicated that Vdr and p63 associated genomic regions regulate genes involving stem cell fate and epidermal differentiation. To demonstrate the functional interaction between VDR and p63, we evaluated the response to 1,25(OH)2D3 of keratinocytes lacking p63 and noted a reduction in epidermal cell fate determining transcription factors such as Fos, Jun. We conclude that VDR is required for the epidermal stem cell fate orientation towards interfollicular epidermis. We propose that this role of VDR involves cross-talk with the epidermal master regulator p63 through super-enhancer mediated epigenetic dynamics.

Highlights from the 24th workshop on vitamin D in Austin, September 2022.

The 24th Workshop on Vitamin D was held September 7-9, 2022 in Austin, Texas and covered a wide diversity of research in the vitamin D field from across the globe. Here, we summarize the meeting, individual sessions, awards and presentations given.