Human 5-lipoxygenase regulates transcription by association to euchromatin.

Human 5-lipoxygenase (5-LO) is the key enzyme of leukotriene biosynthesis, mostly expressed in leukocytes and thus a crucial component of the innate immune system. In this study, we show that 5-LO, besides its canonical function as an arachidonic acid metabolizing enzyme, is a regulator of gene expression associated with euchromatin. By Crispr-Cas9-mediated 5-LO knockout (KO) in MonoMac6 (MM6) cells and subsequent RNA-Seq analysis, we identified 5-LO regulated genes which could be clustered to immune/defense response, cell adhesion, transcription and growth/developmental processes. Analysis of differentially expressed genes identified cyclooxygenase-2 (COX2, PTGS2) and kynureninase (KYNU) as strongly regulated 5-LO target genes. 5-LO knockout affected MM6 cell adhesion and tryptophan metabolism via inhibition of the degradation of the immunoregulator kynurenine. By subsequent FAIRE-Seq and 5-LO ChIP-Seq analyses, we found an association of 5-LO with euchromatin, with prominent 5-LO binding to promoter regions in actively transcribed genes. By enrichment analysis of the ChIP-Seq results, we identified potential 5-LO interaction partners. Furthermore, 5-LO ChIP-Seq peaks resemble patterns of H3K27ac histone marks, suggesting that 5-LO recruitment mainly takes place at acetylated histones. In summary, we demonstrate a noncanonical function of 5-LO as transcriptional regulator in monocytic cells.

Nuclear receptor activation shapes spatial genome organization essential for gene expression control: lessons learned from the vitamin D receptor.

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)2D3-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

A hierarchical regulatory network analysis of the vitamin D induced transcriptome reveals novel regulators and complete VDR dependency in monocytes.

The transcription factor vitamin D receptor (VDR) is the high affinity nuclear target of the biologically active form of vitamin D3 (1,25(OH)2D3). In order to identify pure genomic transcriptional effects of 1,25(OH)2D3, we used VDR cistrome, transcriptome and open chromatin data, obtained from the human monocytic cell line THP-1, for a novel hierarchical analysis applying three bioinformatics approaches. We predicted 75.6% of all early 1,25(OH)2D3-responding (2.5 or 4 h) and 57.4% of the late differentially expressed genes (24 h) to be primary VDR target genes. VDR knockout led to a complete loss of 1,25(OH)2D3-induced genome-wide gene regulation. Thus, there was no indication of any VDR-independent non-genomic actions of 1,25(OH)2D3 modulating its transcriptional response. Among the predicted primary VDR target genes, 47 were coding for transcription factors and thus may mediate secondary 1,25(OH)2D3 responses. CEBPA and ETS1 ChIP-seq data and RNA-seq following CEBPA knockdown were used to validate the predicted regulation of secondary vitamin D target genes by both transcription factors. In conclusion, a directional network containing 47 partly novel primary VDR target transcription factors describes secondary responses in a highly complex vitamin D signaling cascade. The central transcription factor VDR is indispensable for all transcriptome-wide effects of the nuclear hormone.

Primary Vitamin D Target Genes of Human Monocytes.

The molecular basis of vitamin D signaling implies that the metabolite 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) of the secosteroid vitamin D3 activates the transcription factor vitamin D receptor (VDR), which in turn modulates the expression of hundreds of primary vitamin D target genes. Since the evolutionary role of nuclear receptors, such as VDR, was the regulation of cellular metabolism, the control of calcium metabolism became the primary function of vitamin D and its receptor. Moreover, the nearly ubiquitous expression of VDR enabled vitamin D to acquire additional physiological functions, such as the support of the innate immune system in its defense against microbes. Monocytes and their differentiated phenotypes, macrophages and dendritic cells, are key cell types of the innate immune system. Vitamin D signaling was most comprehensively investigated in THP-1 cells, which are an established model of human monocytes. This includes the 1,25(OH)2D3-modulated cistromes of VDR, the pioneer transcription factors PU.1 and CEBPA and the chromatin modifier CTCF as well as of the histone markers of promoter and enhancer regions, H3K4me3 and H3K27ac, respectively. These epigenome-wide datasets led to the development of our chromatin model of vitamin D signaling. This review discusses the mechanistic basis of 189 primary vitamin D target genes identified by transcriptome-wide analysis of 1,25(OH)2D3-stimulated THP-1 cells and relates the epigenomic basis of four different regulatory scenarios to the physiological functions of the respective genes.

Modulation of vitamin D signaling by the pioneer factor CEBPA.

The myeloid master regulator CCAAT enhancer-binding protein alpha (CEBPA) is known as a pioneer factor. In this study, we report the CEBPA cistrome of THP-1 human monocytes after stimulation with the vitamin D receptor (VDR) ligand 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) for 2, 8 and 24 h. About a third of the genomic VDR binding sites co-located with those of CEBPA. In parallel, the binding strength of 5% of the CEBPA cistrome, i.e. some 1500 sites, is significantly (p < 0.001) affected by 1,25(OH)2D3. Transcriptome-wide analysis after CEBPA silencing indicated that the pioneer factor enhances both the basal expression and ligand inducibility of 70 vitamin D target genes largely involved in lipid signaling and metabolism. In contrast, CEBPA suppresses 82 vitamin D target genes many of which are related to the modulation of T cell activity by monocytes. The inducibility of the promoter-specific histone marker H3K4me3 distinguishes the former class of genes from the latter. Moreover, prominent occupancy of the myeloid pioneer factor PU.1 on 1,25(OH)2D3-sensitive CEBPA enhancers mechanistically explains the dichotomy of vitamin D target genes. In conclusion, CEBPA supports vitamin D signaling concerning actions of the innate immune system, but uses the antagonism with PU.1 for suppressing possible overreactions of adaptive immunity.

In vivo transcriptome changes of human white blood cells in response to vitamin D.

In the vitamin D intervention study VitDbol (NCT02063334) blood samples were drawn directly before an oral bolus (2000 µg vitamin D3) and 24 h later. The focus of phase II of VitDbol was the transcriptome-wide analysis of the effects of vitamin D gene expression in human peripheral blood mononuclear cells (PBMCs). All five participants responded in an individual fashion to the bolus by increases in serum levels of the vitamin D metabolites 25-hydroxyvitamin D3 (25(OH)D3) and 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). RNA sequencing identified 15.040 commonly expressed genes in PBMCs, 702 (4,7%) of which were significantly (p < 0,05) affected by the vitamin D3 bolus. KEGG pathway analysis suggested that these genes are involved in general protein translation, monocyte differentiation and cellular growth control. Previously published transcriptome-wide studies in comparable cell systems confirmed 234 of the 702 vitamin D target genes, leaving many genes, such as HLA-A and HLA-C, as novel discoveries. Interestingly, in vivo stimulated PBMCs of this study showed a larger number of common vitamin D target genes with the monocytic cell line THP-1 than with in vitro stimulated PBMCs. The expression pattern of vitamin D target genes differed significantly between individuals and the average expression change can serve as a marker for vitamin D responsiveness. In conclusion, this study demonstrates that under in vivo conditions changes in 25(OH)D3 and 1,25(OH)2D3 serum concentrations alter the expression of more than 700 vitamin D target genes in human leukocytes.

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 D3 (1,25(OH)2D3). 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)2D3-sensitive H3K27ac marks at enhancers. A categorization of 1,25(OH)2D3-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)2D3 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)2D3 stimulation. This allows a more detailed understanding of vitamin D target gene regulation.

In vivo response of the human epigenome to vitamin D: A Proof-of-principle study.

In vitro cell culture studies showed that the hormonal form of vitamin D3, 1α,25-dihydroxyvitamin D3, significantly (p < 0.05) affects the human epigenome at thousands of genomic loci. Phase II of the VitDbol vitamin D intervention trial (NCT02063334) involved a proof-of-principle study of one individual, who was exposed three times every 28 days to an oral bolus (2000 µg) of vitamin D3. Blood samples were taken directly before each supplementation as well as one and two days after, chromatin was isolated from peripheral blood mononuclear cells without any further in vitro culture and at all nine time points epigenome-wide chromatin accessibility was assessed by applying FAIRE-seq (formaldehyde-assisted isolation of regulatory elements sequencing). The vitamin D3 bolus resulted in an average raise in 25-hydroxyvitamin D3 (25(OH)D3) serum concentration of 11.9 and 19.4 nM within one and two days, respectively. Consistently accessible chromatin was detected at 5205 genomic loci, the 853 most prominent of which a self-organizing map algorithm classified into early, delayed and non-responding genomic regions: 70 loci showed already after one day and 361 sites after two days significant (p < 0.0001) chromatin opening or closing. Interestingly, more than half of these genomic regions overlap with transcription start sites, but the change of chromatin accessibility at these sites has no direct effect on the transcriptome. Some of the vitamin D responsive chromatin sites cluster at specific loci within the human genome, the most prominent of which is the human leukocyte antigen region in chromosome 6. In conclusion, this study demonstrates that under in vivo conditions a rather minor rise in 25(OH)D3 serum levels is sufficient to result in significant changes at hundreds of sites within the epigenome of human leukocytes.

ETS transcription factor family member GABPA contributes to vitamin D receptor target gene regulation.

Binding motifs of the ETS-domain transcription factor GABPA are found with high significance below the summits of the vitamin D receptor (VDR) cistrome. VDR is the nuclear receptor for the biologically most active vitamin D metabolite 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). In this study, we determined the GABPA cistrome in THP-1 human monocytes and found that it is comprised of 3822 genomic loci, some 20% of which were modulated by 1,25(OH)2D3. The GABPA cistrome showed a high overlap rate with accessible chromatin and the pioneer transcription factor PU.1. Interestingly, 23 and 12% of persistent and transient VDR binding sites, respectively, co-localized with GABPA, which is clearly higher than the rate of secondary VDR loci (4%). Some 40% of GABPA binding sites were found at transcription start sites, nearly 100 of which are of 1,25(OH)2D3 target genes. On 593 genomic loci VDR and GABPA co-localized with PU.1, while only 175 VDR sites bound GABPA in the absence of PU.1. In total, VDR sites with GABPA co-localization may control some 450 vitamin D target genes. Those genes that are co-controlled by PU.1 preferentially participate in cellular and immune signaling processes, while the remaining genes are involved in cellular metabolism pathways. In conclusion, GABPA may contribute to differential VDR target gene regulation.

Selective regulation of biological processes by vitamin D based on the spatio-temporal cistrome of its receptor.

The transcription factor vitamin D receptor (VDR) is the exclusive nuclear target of the biologically active form of vitamin D (1,25(OH)2D3). In THP-1 human monocytes we obtained a highly accurate VDR cistrome after 2 and 24h ligand stimulation comprising >11,600 genomic loci, 78% of which were detected exclusively after 24h. In contrast, a group of 510 persistent VDR sites occurred at all conditions and some 2100 VDR loci were only transiently occupied. Machine learning and statistical analysis as well as a comparison with the re-analyzed B cell VDR cistrome indicated a subgroup of 339 highly conserved persistent VDR sites that were suited best for describing vitamin D-triggered gene regulatory scenarios. The 1,25(OH)2D3-dependent transcriptome of THP-1 cells comprised 587 genes, 311 of which were primary targets with main functions in the immune system. More than 97% of the latter genes were located within 1,25(OH)2D3-modulated topologically associated domains (TADs). The number of persistent and transient VDR sites was found to be the main discriminator for sorting these TADs into five classes carrying vitamin D target genes involved in distinct biological processes. In conclusion, specific regulation of biological processes by vitamin D depends on differences in time-dependent VDR binding.

Epigenomic PU.1-VDR crosstalk modulates vitamin D signaling.

The ETS-domain transcription factor PU.1 acts as a pioneer factor for other transcription factors including nuclear receptors. In this study, we report that in THP-1 human monocytes the PU.1 cistrome comprises 122,319 genomic sites. Interestingly, at 6498 (5.3%) of these loci PU.1 binding was significantly modulated by the vitamin D receptor (VDR) ligand 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). In most cases 1,25(OH)2D3 increased PU.1 association, which correlated strongly with VDR co-location and overlap ratios for canonical DR3-type VDR binding sites. Genome-wide 6488 sites associating both with PU.1 and VDR as well as 5649 non-VDR overlapping, 1,25(OH)2D3-sensitive PU.1 loci represent the PU.1-VDR crosstalk and can be described by four gene regulatory scenarios, each. Chromatin accessibility was the major discriminator between these models. The location of the PU.1 binding loci in open chromatin coincided with a significantly smaller mean distance to the closest 1,25(OH)2D3 target gene. PU.1 knockdown indicated that the pioneer factor is relevant for the transcriptional activation of 1,25(OH)2D3 target genes but its impact differed in magnitude and orientation. In conclusion, PU.1 is an important modulator of VDR signaling in monocytes, including but also exceeding its role as a pioneer factor, but we found no evidence for a direct interaction of both proteins.

Molecular evaluation of vitamin D responsiveness of healthy young adults.

Vitamin D3 has via its metabolites 25-hydroxyvitamin D3 (25(OH)D3) and 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) direct effects on the transcriptome and the epigenome of most human cells. In the VitDbol study we exposed 35 healthy young adults to an oral vitamin D3 dose (2000µg) or placebo and took blood samples directly before the supplementation as well as at days 1, 2 and 30. Within 24h the vitamin D3 intake raised the average serum levels of both 25(OH)D3 and 1,25(OH)2D3 by approximately 20%. However, we observed large inter-individual differences in these serum levels, reflected by the average ratios between 25(OH)D3 and 1,25(OH)2D3 concentrations ranging from 277 to 1365. Interestingly, average serum parathyroid hormone (PTH) levels increased at day 1 by some 10% but then decreased within the following four weeks to levels 5% below baseline. In peripheral blood mononuclear cells (PBMCs) that were isolated at the same time points we determined vitamin D-modulated chromatin accessibility by FAIRE-qPCR at selected genomic loci. This method is well suited to evaluate both short-term and long-term in vivo effects of vitamin D on the epigenome of human subjects. The differential vitamin D responsiveness of the VitDbol study participants was determined via individual changes in their PTH levels or chromatin accessibility in relation to alterations in 25(OH)D3 concentrations. This led to the segregation of the subjects into 14 high, 11 mid and 10 low responders. In summary, the vitamin D responsiveness classification provides additional information compared to a vitamin D status assessment based on single 25(OH)D3 serum measurements. The study was registered at Clinicaltrials.gov (NCT02063334).

Vitamin D-dependent chromatin association of CTCF in human monocytes.

CCCTC-binding factor (CTCF) is a transcription factor being involved in 3D chromatin organization and displays a highly conserved genome-wide binding pattern. In this study, we report the cistrome of CTCF in THP-1 human monocytes and confirm that from the 40,078 CTCF binding sites nearly 85% are identical with those found in K562 monocytes. Quadruplicate chromatin immunoprecipitation sequencing (ChIP-seq) demonstrated that at 2130 loci the association strenght of CTCF with genomic DNA was significantly (p<0.05) modulated by stimulation with the natural vitamin D receptor (VDR) ligand 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). Some 55% of these CTCF sites contribute to DNA looping and mark the anchors of 587 putative topologically associating domains (TADs) containing at least one VDR binding site and one 1,25(OH)2D3 target gene. These TADs can explain the regulatory scenarios of up to 70% of all 1,25(OH)2D3 target genes. A self-organizing map approach subdivided the vitamin D-sensitive CTCF sites into seven classes that can be distinguished by participation in DNA loop formation, binding to open chromatin, carrying binding motifs for CTCF or its relative BORIS, overlap with transcription start site (TSS) regions and binding of VDR. These variant molecular profiles suggest different mechanisms of the 1,25(OH)2D3-dependent action of CTCF. The co-location of VDR and 1,25(OH)2D3-dependent CTCF sites increases in the context of accessible chromatin and TSS regions but does not show any significant correlation with classical DNA binding mechanisms of CTCF. In conclusion, vitamin D-sensitive CTCF sites provide further mechanistic details to the epigenome-wide understanding of 1,25(OH)2D3-mediated gene regulation.

The vitamin D-dependent transcriptome of human monocytes.

Monocytes are important cells of the innate immune system that can differentiate into macrophages and dendritic cells. The biologically active form of vitamin D, 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3), serves as a ligand of the nuclear receptor vitamin D receptor (VDR). A key physiological function of 1,25(OH)2D3 is the defense against pathogens, such as those causing tuberculosis, that involves the modulation of the monocyte transcriptome. THP-1 cells are an established model of human monocytes, for which the at present largest set of 1,25(OH)2D3-affected genome-wide data are available. Here we summarize the insight obtained from the recent transcriptome of 1,25(OH)2D3-stimulated THP-1 cells, that was determined by triplicate RNA sequencing (RNA-seq). Primary and secondary vitamin D target genes being up- and down-regulated were related to changes in the epigenome of THP-1 cells, such as 1,25(OH)2D3-dependent chromatin opening and modulation of the genome-wide association of the transcription factors VDR and CCCTC-binding factor (CTCF) with their respective genomic binding sites. The anti-microbial response is the top-ranking early physiological function represented by 1,25(OH)2D3-stimulated genomic regions and genes, but also other immunity-related pathways, such as IL10 signaling, are activated. Taken together, the epigenomic and transcriptomic responses of THP-1 cells to 1,25(OH)2D3 represent a master example of the impact of vitamin D on human physiology.

Epigenome-wide effects of vitamin D and their impact on the transcriptome of human monocytes involve CTCF.

The physiological functions of vitamin D are mediated by its metabolite 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) activating the transcription factor vitamin D receptor (VDR). In THP-1 human monocytes we demonstrated epigenome-wide effects of 1,25(OH)2D3 at 8979 loci with significantly modulated chromatin accessibility. Maximal chromatin opening was observed after 24 h, while after 48 h most sites closed again. The chromatin-organizing protein CTCF bound to 14% of the 1,25(OH)2D3-sensitive chromatin regions. Interestingly, 1,25(OH)2D3 affected the chromatin association of CTCF providing an additional mechanism for the epigenome-wide effects of the VDR ligand. The 1,25(OH)2D3-modulated transcriptome of THP-1 cells comprised 1284 genes, 77.5% of which responded only 24 h after stimulation. During the 1,25(OH)2D3 stimulation time course the proportion of down-regulated genes increased from 0% to 44.9% and the top-ranking physiological function of the respective genes shifted from anti-microbial response to connective tissue disorders. The integration of epigenomic and transcriptomic data identified 165 physiologically important 1,25(OH)2D3 target genes, including HTT and NOD2, whose expression can be predicted primarily from epigenomic data of their genomic loci. Taken together, a large number of 1,25(OH)2D3-triggered epigenome-wide events precede and accompany the transcriptional activation of target genes of the nuclear hormone.

Relevance of vitamin D receptor target genes for monitoring the vitamin D responsiveness of primary human cells.

Vitamin D3 has transcriptome- and genome-wide effects and activates, via the binding of its metabolite 1α,25-dihydroxyvitamin D3 to the transcription factor vitamin D receptor (VDR), several hundred target genes. Using samples from a 5-month vitamin D3 intervention study (VitDmet), we recently reported that the expression of 12 VDR target genes in peripheral blood mononuclear cells (PBMCs) as well as 12 biochemical and clinical parameters of the study participants are significantly triggered by vitamin D3. In this study, we performed a more focused selection of further 12 VDR target genes and demonstrated that changes of their mRNA expression in PBMCs of VitDmet subjects significantly correlate with alterations of 25-hydroxyvitamin D3 serum levels. Network and self-organizing map analysis of these datasets together with that of the other 24 parameters was followed by relevance calculations and identified changes in parathyroid hormone serum levels and the expression of the newly selected genes STS, BCL6, ITGAM, LRRC25, LPGAT1 and TREM1 as well as of the previously reported genes DUSP10 and CD14 as the most relevant parameters for describing vitamin D responsiveness in vivo. Moreover, parameter relevance ranking allowed the segregation of study subjects into high and low responders. Due to the long intervention period the vitamin D response was not too prominent on the level of transcriptional activation. Therefore, we performed in the separate VitDbol trial a short-term but high dose stimulation with a vitamin D3 bolus. In PBMCs of VitDbol subjects we observed direct transcriptional effects on the selected VDR target genes, such as an up to 2.1-fold increase already one day after supplementation onset. In conclusion, both long-term and short-term vitamin D3 supplementation studies allow monitoring the vitamin D responsiveness of human individuals and represent new types of human in vivo vitamin D3 investigations.

The transcriptional regulator BCL6 participates in the secondary gene regulatory response to vitamin D.

The vitamin D metabolite 1,25-dihydroxyvitamin D3 (1,25(OH)2D3) is the high affinity ligand of the transcription factor vitamin D receptor (VDR) and therefore a direct regulator of transcription. Transcriptome-wide analysis of THP-1 human monocytes had indicated more than 600 genes to be significantly (p<0.05) stimulated after 4h incubation with 1,25(OH)2D3, but only 67 of them where more than 1.5-fold up-regulated. These include the genes encoding for the transcription factors BCL6, NFE2, POU4F2 and ELF4, which are controlled by one or two VDR binding sites within their chromosomal domains. The latter are defined via DNA loop formation mediated by the transcription factor CTCF that is highly conserved in its genome-wide loci. We found BCL6 being most responsive to 1,25(OH)2D3 and selected it for further analysis. An incubation of THP-1 cells with 1,25(OH)2D3 for 24 h resulted in a significant (p<0.001) change in the mRNA expression of more than 1600 genes, of which 132 were at least 2-fold up-regulated. About half of the latter genes are secondary 1,25(OH)2D3 targets, since they do not carry any VDR binding site within their chromosomal domain. Chromatin immunoprecipitation sequencing datasets indicated that the majority of these domains contain a BCL6 binding site. We followed the secondary transcriptional response to 1,25(OH)2D3 for eight representative gene examples and confirmed the binding of CTCF and BCL6 to their respective chromosomal domains. In conclusion, our study indicated that in monocytes most of the physiological responses to 1,25(OH)2D3 involve the action of the transcription factor BCL6.

Primary vitamin D receptor target genes as biomarkers for the vitamin D3 status in the hematopoietic system.

Vitamin D(3) belongs to the few nutritional compounds that has, via the binding of its metabolite 1α,25-dihydroxyvitamin D(3) (1,25(OH)(2)D(3)) to the transcription factor vitamin D receptor (VDR), a direct effect on gene regulation. The relation of thousands of genomic VDR-binding sites to a few hundred primary 1,25(OH)(2)D(3) target genes is still largely unresolved. We studied chromatin domains containing genes for the adhesion molecules CD97 and LRRC8A, the glucose transporter SLC37A2 and the coactivator NRIP1. These domains vary significantly in size (7.3 to 956 kb) but contain each one major VDR-binding site. In monocytic cells these four sites are associated with open chromatin and occupied by VDR, while in macrophage-like cells only the sites of LRRC8A, SLC37A2 and NRIP1 are accessible and receptor bound. The VDR site of CD97 does, in contrast to the three other loci, not carry any DR3-type binding sequence. CD97, LRRC8A, SLC37A2 and NRIP1 are early responding 1,25(OH)(2)D(3) target genes in monocytic cells, while in macrophage-like cells they respond less and, in part, delayed. In primary human peripheral blood mononuclear cells from 71 prediabetic subjects of a vitamin D(3) intervention study (VitDmet) CD97, LRRC8A, SLC37A2 and NRIP1 can be used as transcriptomic biomarkers for classifying human individuals for their possible benefit from vitamin D(3) supplementation. In particular, NRIP1 exceeds the potential of the previously identified marker CD14 by more than 40% and seems to be a well-suited molecular marker for the vitamin D(3) status in the hematopoietic system.

Characterization of genomic vitamin D receptor binding sites through chromatin looping and opening.

The vitamin D receptor (VDR) is a transcription factor that mediates the genomic effects of 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). Genome-wide there are several thousand binding sites and hundreds of primary 1,25(OH)2D3 target genes, but their functional relation is largely elusive. In this study, we used ChIA-PET data of the transcription factor CTCF in combination with VDR ChIP-seq data, in order to map chromatin domains containing VDR binding sites. In total, we found 1,599 such VDR containing chromatin domains and studied in THP-1 human monocytic leukemia cells four representatives of them. Our combined ChIP-seq and FAIRE-seq time course data showed that each of these four domains contained a master VDR binding site, where an increase of VDR binding pairs with 1,25(OH)2D3-promoted chromatin opening and the presence of a highly significant DR3-type sequence below the peak summit. These sites differed in their relative VDR binding but not in their kinetics, while other loci either had a weaker and delayed VDR association or could not be confirmed at all. All studied chromatin domains contained at least one primary 1,25(OH)2D3 target gene demonstrating a characteristic slope of mRNA increase, while neighboring genes responded delayed, if at all. In conclusion, the observation of ligand-inducible VDR binding and chromatin opening combined with a DR3-type sequence highlighted genome-wide 160 VDR loci that have within their chromatin domain a more than 4-fold increased likelihood to identify a primary 1,25(OH)2D3 target gene than in the vicinity of other genomic VDR binding sites.

The ASAP2 gene is a primary target of 1,25-dihydroxyvitamin D3 in human monocytes and macrophages.

A genome-wide data set on vitamin D receptor (VDR) binding sites in human THP-1 cells (monocytes) led us to the genomic region around the ASAP2 (Arf-GAP with SH3 domain, ankyrin repeat and PH domain 2) gene, whose product is involved in the regulation of vesicular transport, cellular migration and autophagy. Using ENCODE data, we demonstrated that the ASAP2 gene is flanked by conserved binding sites of the insulating transcription factor CTCF. These sites define different chromosomal domains containing the ASAP2 gene, up to six additional genes and three VDR binding sites. In human monocytes (THP-1 cells) the ASAP2 gene is more weakly expressed but more and faster inducible by the biologically active form of vitamin D, 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3), than in M2-type macrophages (phorbol ester-differentiated THP-1 cells). Within the investigated genomic region, the basal mRNA expressions of the neighboring genes are comparably high in both monocytes and macrophages, but the ASAP2 gene is the only primary 1,25(OH)2D3 target. The three VDR binding sites located 54, 436 and 574kb downstream of the ASAP2 transcription start site each carry a sequence formed by a direct repeat with three intervening nucleotides (DR3). Ligand-inducible VDR binding was confirmed to all three genomic sites in monocytes and macrophages. Taken together, the region around the ASAP2 gene is genome-wide highlighted as a special attraction point for the VDR, but the presently sole known functional consequence of the binding of VDR to three sites within this chromosomal region is that ASAP2 is a primary 1,25(OH)2D3 target gene in monocytes and macrophages. This article is part of a Special Issue entitled '16th Vitamin D Workshop'.