Sabine Seuter

Nuclear hormone 1α,25-dihydroxyvitamin D3 elicits a genome-wide shift in the locations of VDR chromatin occupancy

A global understanding of the actions of the nuclear hormone 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3) and its vitamin D receptor (VDR) requires a genome-wide analysis of VDR binding sites. In THP-1 human monocytic leukemia cells we identified by ChIP-seq 2340 VDR binding locations, of which 1171 and 520 occurred uniquely with and without 1α,25(OH)2D3 treatment, respectively, while 649 were common. De novo identified direct repeat spaced by 3 nucleotides (DR3)-type response elements (REs) were strongly associated with the ligand-responsiveness of VDR occupation. Only 20% of the VDR peaks diminishing most after ligand treatment have a DR3-type RE, in contrast to 90% for the most growing peaks. Ligand treatment revealed 638 1α,25(OH)2D3 target genes enriched in gene ontology categories associated with immunity and signaling. From the 408 upregulated genes, 72% showed VDR binding within 400u2009kb of their transcription start sites (TSSs), while this applied only for 43% of the 230 downregulated genes. The VDR loci showed considerable variation in gene regulatory scenarios ranging from a single VDR location near the target gene TSS to very complex clusters of multiple VDR locations and target genes. In conclusion, ligand binding shifts the locations of VDR occupation to DR3-type REs that surround its target genes and occur in a large variety of regulatory constellations.

Primary effect of 1α,25(OH)2D3 on IL-10 expression in monocytes is short-term down-regulation

The biologically most active vitamin D compound, 1α,25-dihydroxyvitamin D₃ (1α,25(OH)₂D₃), influences the status of inflammation by modulating the expression of several cytokine genes. In this study, we have examined the mechanism of transcriptional regulation of interleukin 10 (IL-10) by 1α,25(OH)₂D₃ in lipopolysaccharide (LPS)-treated human monocytes (THP-1). Quantitative PCR showed that IL-10 mRNA expression was significantly down-regulated (2.8-fold) during the first 8h of 1α,25(OH)₂D₃ treatment, while after 48 h it was up-regulated (3-fold). Gel shift and quantitative chromatin immunoprecipitation (ChIP) assays showed that the vitamin D receptor (VDR) binds in a cyclical fashion to a promoter region 1500-1700 bp upstream of the IL-10 transcription start site (TSS) containing two conserved VDR binding sites. Targeting of VDR binding sites by enhancer specific duplex RNAs revealed that only the more distal element is functional and chromosome conformation capture analysis suggested that this region loops 1α,25(OH)₂D₃-dependently to the TSS. Quantitative ChIP and micrococcal nuclease assays also revealed 1α,25(OH)₂D₃-dependent cyclical epigenetic changes and nucleosome remodeling at this promoter region. In conclusion, in LPS-treated THP-1 cells the primary effect of 1α,25(OH)₂D₃ on IL-10 expression is down-regulation, which is achieved via a cyclical recruitment of VDR to the promoter.

Chromatin acetylation at transcription start sites and vitamin D receptor binding regions relates to effects of 1α,25-dihydroxyvitamin D3 and histone deacetylase inhibitors on gene expression

The nuclear hormone 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3 or 1,25D) regulates its target genes via activation of the transcription factor vitamin D receptor (VDR) far more specifically than the chromatin modifier trichostatin A (TsA) via its inhibitory action on histone deacetylases. We selected the thrombomodulin gene locus with its complex pattern of five VDR binding sites and multiple histone acetylation and open chromatin regions as an example to investigate together with a number of reference genes, the primary transcriptional responses to 1α,25(OH)2D3 and TsA. Transcriptome-wide, 18.4% of all expressed genes are either up-or down-regulated already after a 90 min TsA treatment; their response pattern to 1α,25(OH)2D3 and TsA sorts them into at least six classes. TsA stimulates a far higher number of genes than 1α,25(OH)2D3 and dominates the outcome of combined treatments. However, 200 TsA target genes can be modulated by 1α,25(OH)2D3 and more than 1000 genes respond only when treated with both compounds. The genomic view on the genes suggests that the degree of acetylation at transcription start sites and VDR binding regions may determine the effect of TsA on mRNA expression and its interference with 1α,25(OH)2D3. Our findings hold true also for other HDAC inhibitors and may have implications on dual therapies using chromatin modifiers and nuclear receptor ligands.

Dynamics of nuclear receptor target gene regulation

Ligand-regulated nuclear receptors, such as estrogen receptors, glucocorticoid receptor, vitamin D receptor, and peroxisome proliferator-activated receptors, belong to the most widely studied and best understood transcription factors. Therefore, the dynamic nature of transcriptional regulation was observed first with different members of the nuclear receptor superfamily, but is now also extended to other transcription factors, such as nuclear factor κB. Dynamic and in part cyclical processes were observed on the level of translocation into the nucleus, association with genomic binding sites, exchange of co-regulators and chromatin modifiers, occurrence of chromatin marks, and activities of RNA polymerase II resulting in mRNA synthesis. In this review, we summarize recent findings on the dynamic regulation of nuclear receptor target genes in the chromatin context.

Primary Vitamin D Target Genes Allow a Categorization of Possible Benefits of Vitamin D3 Supplementation

Vitamin D deficiency has been associated with an increased risk of developing a number of diseases. Here we investigated samples from 71 pre-diabetic individuals of the VitDmet study, a 5-month high dose vitamin D3 intervention trial during Finnish winter, for their changes in serum 25-hydroxyvitamin D3 (25(OH)D3) concentrations and the expression of primary vitamin D target genes in peripheral blood mononuclear cells and adipose tissue. A negative correlation between serum concentrations of parathyroid hormone and 25(OH)D3 suggested an overall normal physiological vitamin D response among the participants. The genes CD14 and thrombomodulin (THBD) are up-regulated primary vitamin D targets and showed to be suitable gene expression markers for vitamin D signaling in both primary tissues. However, in a ranking of the samples concerning their expected response to vitamin D only the top half showed a positive correlation between the changes of CD14 or THBD mRNA and serum 25(OH)D3 concentrations. Interestingly, this categorization allows unmasking a negative correlation between changes in serum concentrations of 25(OH)D3 and the inflammation marker interleukin 6. We propose the genes CD14 and THBD as transcriptomic biomarkers, from which the effects of a vitamin D3 supplementation can be evaluated. These biomarkers allow the classification of subjects into those, who might benefit from a vitamin D3 supplementation, and others who do not.

Mechanism of 1α,25-dihydroxyvitamin D3-dependent repression of interleukin-12B

Interleukin 12 (IL-12) is a heterodimeric, pro-inflammatory cytokine that plays a central role in activation and differentiation of CD4(+) T cells into interferon-γ secreting T-helper type 1 cells. IL-12B, a gene encoding the larger subunit of active IL-12, has been reported to be down-regulated by the nuclear hormone 1α,25-dihydroxyvitamin D(3) (1α,25(OH)(2)D(3)), but the mechanism of the regulation is unknown. In this study, we have examined the molecular mechanism of transcriptional regulation of the IL-12B gene by 1α,25(OH)(2)D(3) in lipopolysaccharide (LPS)-treated human monocytes (THP-1). Quantitative RT-PCR showed that IL-12B mRNA displays a cyclical expression profile and is down-regulated 2.8-fold during the first 8h and even 12.1-fold 24h after exposure to 1α,25(OH)(2)D(3). Gel shift and quantitative chromatin immunoprecipitation (ChIP) assays demonstrated vitamin D receptor (VDR) binding to genomic regions 480 and 6300bp upstream of the IL-12B transcription start site (TSS). Quantitative ChIP assays also revealed that together with VDR and its partner RXR the above regions recruited the co-repressor NCOR2/SMRT and histone deacetylase 3 leading to a decreased histone 4 acetylation and increased histone 3 trimethylation at the IL-12B promoter and its TSS. We suggest that these repressive epigenetic changes eventually cause down-regulation of IL-12 expression. This article is part of a Special Issue entitled: 11th European Symposium on Calcium.

Trichostatin A and structurally related histone deacetylase inhibitors induce 5-lipoxygenase promoter activity

Abstract 5-Lipoxygenase (5-LO) mRNA expression in Mono Mac 6 cells is induced by the histone deacetylase inhibitor trichostatin A (TsA). In order to study the effects of TsA and several structurally related compounds such as MD85, D237 and M232 on 5-LO promoter activity, we have analyzed the response of a 5-lipoxygenase (5-LO) promoter luciferase reporter gene construct to histone deacetylase inhibitors in transiently transfected Mono Mac 6 and HeLa cells. We show that the activity of 5-LO promoter constructs comprising the sequences -778 to +53 and of several successive deletions of the 5-LO promoter is strongly increased upon TsA treatment. The data suggest a significant involvement of histone deacetylases in the regulation of 5-LO gene transcription. The basal activity of the 5-LO promoter strongly depends on the presence of multiple Sp1-binding sites (GC-boxes), five of which are positioned in tandem. Deletion of the five tandemized GCboxes in the 5-LO reporter gene construct revealed that the induction of 5-LO promoter activity by TsA seems to be independent of these GCboxes. Methylation of 5-LO reporter gene constructs by M.HpaII reduced 5-LO promoter activity but did not prevent induction of promoter activity by TsA, although the activated reporter gene activities were lower compared to the unmethylated plasmid, indicating the dominance of methylation over TsA-sensitive histone deacetylation in silencing of the 5-LO gene. The structure-activity data obtained for histone deacetylase inhibitors suggest that this assay system might serve as a cellular screening tool for the development of HDAC inhibitors.

Dynamics of 1α,25-dihydroxyvitamin D3-dependent chromatin accessibility of early vitamin D receptor target genes

The signaling cascade of the transcription factor vitamin D receptor (VDR) is triggered by its specific ligand 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3). In this study we demonstrate that in THP-1 human monocytic leukemia cells 87.4% of the 1034 most prominent genome-wide VDR binding sites co-localize with loci of open chromatin. At 165 of them 1α,25(OH)2D3 strongly increases chromatin accessibility and has at further 217 sites weaker effects. Interestingly, VDR binding sites in 1α,25(OH)2D3-responsive chromatin regions are far more often composed of direct repeats with 3 intervening nucleotides (DR3s) than those in ligand insensitive regions. DR3-containing VDR sites are enriched in the neighborhood of genes that are involved in controling cellular growth, while non-DR3 VDR binding is often found close to genes related to immunity. At the example of six early VDR target genes we show that the slope of their 1α,25(OH)2D3-induced transcription correlates with the basal chromatin accessibility of their major VDR binding regions. However, the chromatin loci controlling these genes are indistinguishable in their VDR association kinetics. Taken together, ligand responsive chromatin loci represent dynamically regulated contact points of VDR with the genome, from where it controls early 1α,25(OH)2D3 target genes.

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.

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.