Carsten Carlberg

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 400 kb 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.

Vitamin D receptor signaling mechanisms: Integrated actions of a well-defined transcription factor

The main physiological actions of the biologically most active metabolite of vitamin D, 1α,25-dihydroxyvitamin D(3) (1α,25(OH)(2)D(3)), are calcium and phosphorus uptake and transport and thereby controlling bone formation. Other emergent areas of 1α,25(OH)(2)D(3) action are in the control of immune functions, cellular growth and differentiation. All genomic actions of 1α,25(OH)(2)D(3) are mediated by the transcription factor vitamin D receptor (VDR) that has been the subject of intense study since the 1980s. Thus, vitamin D signaling primarily implies the molecular actions of the VDR. In this review, we present different perspectives on the VDR that incorporate its role as transcription factor and member of the nuclear receptor superfamily, its dynamic changes in genome-wide locations and DNA binding modes, its interaction with chromatin components and its primary protein-coding and non-protein coding target genes and finally how these aspects are united in regulatory networks. By comparing the actions of the VDR, a relatively well-understood and characterized protein, with those of other transcription factors, we aim to build a realistic positioning of vitamin D signaling in the context of other intracellular signaling systems.

Meta-analysis of primary target genes of peroxisome proliferator-activated receptors

BackgroundPeroxisome proliferator-activated receptors (PPARs) are known for their critical role in the development of diseases, such as obesity, cardiovascular disease, type 2 diabetes and cancer. Here, an in silico screening method is presented, which incorporates experiment- and informatics-derived evidence, such as DNA-binding data of PPAR subtypes to a panel of PPAR response elements (PPREs), PPRE location relative to the transcription start site (TSS) and PPRE conservation across multiple species, for more reliable prediction of PPREs.ResultsIn vitro binding and in vivo functionality evidence agrees with in silico predictions, validating the approach. The experimental analysis of 30 putative PPREs in eight validated PPAR target genes indicates that each gene contains at least one functional, strong PPRE that occurs without positional bias relative to the TSS. An extended analysis of the cross-species conservation of PPREs reveals limited conservation of PPRE patterns, although PPAR target genes typically contain strong or multiple medium strength PPREs. Human chromosome 19 was screened using this method, with validation of six novel PPAR target genes.ConclusionAn in silico screening approach is presented, which allows increased sensitivity of PPAR binding site and target gene detection.

Current Status of Vitamin D Signaling and Its Therapeutic Applications

Vitamin D and in particular its biologically most active metabolite, 1α,25-dihydroxyvitamin D₃ (1α,25(OH)₂D₃), are central endocrine molecules that influence many aspects of human physiology, which are not only the well-known calcium and phosphorus up-take and transport controlling bone formation, but also the control of immune functions and of cellular growth and differentiation. Basically all actions of 1α,25(OH)₂D₃ are mediated by the transcription factor vitamin D receptor (VDR). The crystal structure of the VDR and detailed knowledge on its molecular interactions with the ligand provide significant insight into the mechanisms of vitamin D signaling. This applies also on the action of the huge number of synthetic 1α,25(OH)₂D₃ analogues, which have been developed with the goal of a therapeutic application in hyper-proliferative diseases, such as psoriasis, benign prostate hyperplasia and different types of cancer, in immune functions, such as autoimmune diseases and microbial infections, or in bone disorders, such as osteoporosis. Moreover, detailed investigations on many VDR target genes and in particular the recently available genome-wide view on vitamin D signaling allows a more complete view on the potential of the nuclear hormone. In this review we discuss the latest insight into vitamin D signaling in context with the most prominent 1α,25(OH)₂D₃ analogues.

Genome-wide (over)view on the actions of vitamin D

For a global understanding of the physiological impact of the nuclear hormone 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) the analysis of the genome-wide locations of its high affinity receptor, the transcription factor vitamin D receptor (VDR), is essential. Chromatin immunoprecipitation sequencing (ChIP-seq) in GM10855 and GM10861 lymphoblastoid cells, undifferentiated and lipopolysaccharide-differentiated THP-1 monocytes, LS180 colorectal cancer cells and LX2 hepatic stellate cells revealed between 1000 and 13,000 VDR-specific genomic binding sites. The harmonized analysis of these ChIP-seq datasets indicates that the mechanistic basis for the action of the VDR is independent of the cell type. Formaldehyde-assisted isolation of regulatory elements sequencing (FAIRE-seq) data highlight accessible chromatin regions, which are under control of 1,25(OH)2D3. In addition, public data, such as from the ENCODE project, allow to relate the genome-wide actions of VDR and 1,25(OH)2D3 to those of other proteins within the nucleus. For example, locations of the insulator protein CTCF suggest a segregation of the human genome into chromatin domains, of which more than 1000 contain at least one VDR binding site. The integration of all these genome-wide data facilitates the identification of the most important VDR binding sites and associated primary 1,25(OH)2D3 target genes. Expression changes of these key genes can serve as biomarkers for the actions of vitamin D3 and its metabolites in different tissues and cell types of human individuals. Analysis of primary tissues obtained from vitamin D3 intervention studies using such markers indicated a large inter-individual variation for the efficiency of vitamin D3 supplementation. In conclusion, a genome-wide (over)view on the genomic locations of VDR provides a broader basis for addressing vitamin Ds role in health and disease.

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.

Patterns of Genome-Wide VDR Locations

The genome-wide analysis of the binding sites of the transcription factor vitamin D receptor (VDR) is essential for a global appreciation the physiological impact of the nuclear hormone 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). Genome-wide analysis of lipopolysaccharide (LPS)-polarized THP-1 human monocytic leukemia cells via chromatin immunoprecipitation sequencing (ChIP-seq) resulted in 1,318 high-confidence VDR binding sites, of which 789 and 364 occurred uniquely with and without 1,25(OH)2D3 stimulation, while only 165 were common. We re-analyzed five public VDR ChIP-seq datasets with identical peak calling settings (MACS, version 2) and found, using a novel consensus summit identification strategy, in total 23,409 non-overlapping VDR binding sites, 75% of which are unique within the six analyzed cellular models. LPS-differentiated THP-1 cells have 22% more genomic VDR locations than undifferentiated cells and both cell types display more overlap in their VDR locations than the other investigated cell types. In general, the intersection of VDR binding profiles of ligand-stimulated cells is higher than those of unstimulated cells. De novo binding site searches and HOMER screening for binding motifs formed by direct repeats spaced by three nucleotides (DR3) suggest for all six VDR ChIP-seq datasets that these sequences are found preferentially at highly ligand responsive VDR loci. Importantly, all VDR ChIP-seq datasets display the same relationship between the VDR occupancy and the percentage of DR3-type sequences below the peak summits. The comparative analysis of six VDR ChIP-seq datasets demonstrated that the mechanistic basis for the action of the VDR is independent of the cell type. Only the minority of genome-wide VDR binding sites contains a DR3-type sequence. Moreover, the total number of identified VDR binding sites in each ligand-stimulated cell line inversely correlates with the percentage of peak summits with DR3 sites.

Distinct HDACs regulate the transcriptional response of human cyclin-dependent kinase inhibitor genes to trichostatin A and 1α,25-dihydroxyvitamin D3

The anti-proliferative effects of histone deacetylase (HDAC) inhibitors and 1α,25-dihydroxyvitamin D3 [1α,25(OH)2D3] converge via the interaction of un-liganded vitamin D receptor (VDR) with co-repressors recruiting multiprotein complexes containing HDACs and via the induction of cyclin-dependent kinase inhibitor (CDKI) genes of the INK4 and Cip/Kip family. We investigated the effects of the HDAC inhibitor Trichostatin A (TSA) and 1α,25(OH)2D3 on the proliferation and CDKI gene expression in malignant and non-malignant mammary epithelial cell lines. TSA induced the INK4-family genes p18 and p19, whereas the Cip/Kip family gene p21 was stimulated by 1α,25(OH)2D3. Chromatin immunoprecipitation and RNA inhibition assays showed that the co-repressor NCoR1 and some HDAC family members complexed un-liganded VDR and repressed the basal level of CDKI genes, but their role in regulating CDKI gene expression by TSA and 1α,25(OH)2D3 were contrary. HDAC3 and HDAC7 attenuated 1α,25(OH)2D3-dependent induction of the p21 gene, for which NCoR1 is essential. In contrast, TSA-mediated induction of the p18 gene was dependent on HDAC3 and HDAC4, but was opposed by NCoR1 and un-liganded VDR. This suggests that the attenuation of the response to TSA by NCoR1 or that to 1α,25(OH)2D3 by HDACs can be overcome by their combined application achieving maximal induction of anti-proliferative target genes.

The Down-regulation of the Human MYC Gene by the Nuclear Hormone 1α,25-dihydroxyvitamin D3 is Associated with Cycling of Corepressors and Histone Deacetylases

MYC is a pleiotropic transcription factor that coordinates the expression of diverse programs that are together necessary for the growth and expansion of somatic cells. The nuclear hormone 1alpha,25-dihydroxyvitamin D(3) down-regulates MYC expression, but the exact mechanism is still elusive. We found in RWPE-1 normal human prostate cells that 1alpha,25-dihydroxyvitamin D(3) down-regulates MYC mRNA with a periodicity of 30-90 min. In silico screening of the MYC gene locus identified six putative binding sites [vitamin D response elements (VDREs)] for the vitamin D receptor (VDR). Two of these VDREs efficiently bound VDR-retinoid X receptor heterodimers in vitro, and their genomic regions associated with VDR in RWPE-1 cells. Gene-specific small inhibitory RNA silencing indicated that basal MYC mRNA expression, as well as its down-regulation, depended on the exchange factor TBL1X (transducer beta-like 1, X-linked), the corepressor silencing mediator for retinoid and thyroid hormone receptor, and histone deacetylases (HDACs) 2, 6, and 11. Assaying the association of these five proteins with the VDRE-containing genomic regions of the MYC gene locus showed characteristic ligand-dependent profiles of TBL1X, silencing mediator for retinoid and thyroid hormone receptor, HDAC6, and HDAC11, in particular on an evolutionarily conserved VDRE. In conclusion, our data suggest that dynamically composed protein complexes that dock via VDR to the two VDREs may explain the repression of the MYC gene.

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.