Vitamin D receptor ChIP-seq in primary CD4+ cells: relationship to serum 25-hydroxyvitamin D levels and autoimmune disease
© Handel et al.; licensee BioMed Central Ltd. 2013
Received: 8 January 2013
Accepted: 20 June 2013
Published: 12 July 2013
Vitamin D insufficiency has been implicated in autoimmunity. ChIP-seq experiments using immune cell lines have shown that vitamin D receptor (VDR) binding sites are enriched near regions of the genome associated with autoimmune diseases. We aimed to investigate VDR binding in primary CD4+ cells from healthy volunteers.
We extracted CD4+ cells from nine healthy volunteers. Each sample underwent VDR ChIP-seq. Our results were analyzed in relation to published ChIP-seq and RNA-seq data in the Genomic HyperBrowser. We used MEMEChIP for de novo motif discovery. 25-Hydroxyvitamin D levels were measured using liquid chromatography–tandem mass spectrometry and samples were divided into vitamin D sufficient (25(OH)D ≥75 nmol/L) and insufficient/deficient (25(OH)D <75 nmol/L) groups.
We found that the amount of VDR binding is correlated with the serum level of 25-hydroxyvitamin D (r = 0.92, P= 0.0005). In vivo VDR binding sites are enriched for autoimmune disease associated loci, especially when 25-hydroxyvitamin D levels (25(OH)D) were sufficient (25(OH)D ≥75: 3.13-fold, P<0.0001; 25(OH)D <75: 2.76-fold, P<0.0001; 25(OH)D ≥75 enrichment versus 25(OH)D <75 enrichment: P= 0.0002). VDR binding was also enriched near genes associated specifically with T-regulatory and T-helper cells in the 25(OH)D ≥75 group. MEME ChIP did not identify any VDR-like motifs underlying our VDR ChIP-seq peaks.
Our results show a direct correlation between in vivo 25-hydroxyvitamin D levels and the number of VDR binding sites, although our sample size is relatively small. Our study further implicates VDR binding as important in gene-environment interactions underlying the development of autoimmunity and provides a biological rationale for 25-hydroxyvitamin D sufficiency being based at 75 nmol/L. Our results also suggest that VDR binding in response to physiological levels of vitamin D occurs predominantly in a VDR motif-independent manner.
KeywordsVitamin D Autoimmune disease ChIP-seq Functional genomics
Vitamin D is a secosteroid produced from 7-dehydrocholesterol by the action of ultraviolet (UV) radiation within the skin and is hydroxylated to its active molecule 1,25-dihydroxyvitamin D (1,25D3) by the liver and kidneys . A role for vitamin D and UV radiation in autoimmune disease was originally suggested by the latitudinal gradient in the prevalence and incidence for many autoimmune disorders . Epidemiological studies have since confirmed the association of low levels of vitamin D with increased susceptibility to autoimmune disease, in some cases when vitamin D levels are measured prior to the clinical onset of disease [3–6]. The ideal dose of vitamin D supplementation to achieve a sufficient level of 25-hydroxyvitamin D is not clear, although it appears to be in excess of 800 international units .
1,25D3 acts intracellularly via the vitamin D receptor (VDR), a nuclear receptor that forms dimers with retinoid X receptors (RXR) to bind DNA and alter gene transcription . Two studies have analyzed genome-wide binding of VDR using chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq); one using a B-lymphoblastic cell line (LCL) and another using a monocytic cell line (MCL) [9, 10]. The methods of vitamin D stimulation used in each study differed markedly and this may contribute to the differences in VDR binding observed in addition to cell-specific differences . Each study determined that the VDR-RXR dimer recognizes a classical motif (DR3) but that this is present only at some of the VDR binding sites detected by ChIP-seq. The LCL ChIP-seq used genetic susceptibility loci drawn from genome-wide association studies to demonstrate significant overlap between autoimmune susceptibility regions and VDR binding sites .
However, in vivo, the situation is likely to be very different, both because DNA accessibility is likely to be altered in cell lines compared with primary immune cells and also because long-term exposure to physiological levels of 1,25D3 is not replicated well by short-term stimulation with high levels of 1,25D3[12–14]. In the present study we, therefore, aimed to use ChIP-seq to study VDR binding in primary CD4+ cells drawn from healthy individuals with measured serum levels of 25-hydroxyvitamin D.
Healthy volunteers were recruited from the general public and nine samples of whole blood obtained (1_VDR, 2_VDR, 3_VDR, 4_VDR, 5_VDR, HB, PD, SP and SR). CD4+ lymphocytes were separated from whole blood using magnetic-activated cell sorting (MACS) as described in . This project was approved by the Mid and South Buckinghamshire Research Ethics Committee (REC Reference # 09/H0607/7).
25-hydroxyvitamin D measurements
25-Hydroxyvitamin D was measured using liquid chromatography–tandem mass spectrometry.
This was performed as in . Briefly, CD4+ cells were fixed with 1% formaldehyde for 15 minutes then quenched with 0.125 M glycerine. Lysis buffer was added to isolate chromatin and the samples were disrupted with a Douce homogenizer. Sonication was used to sheer the resultant protein-DNA complex into 300 to 500 base pair fragments (Misonix, Farmindale, NY 11735, USA). DNA was quantified using a Nanodrop (Wilmington, DE 19810, USA) spectrophotometer.
Aliquots containing 50 μg of chromatin were precleared with protein A agarose beads (Invitrogen, Paisley PA4 9RF, UK). Genomic regions bound by VDR were precipitated out using anti-VDR rabbit antibody (Santa Cruz Biotechnology, sc-1008, Dallas, Texas 75220, USA) and isolated with protein A agarose beads. This was incubated at 4°C overnight, then washed and antibody-bound fragments eluted from the beads with SDS buffer. Samples were treated with proteinase K and RNase. Crosslinks were reversed by incubation overnight at 65°C. ChIP-DNA was purified by subsequent phenol-chloroform extraction and ethanol precipitation.
The purified product was then prepared for sequencing as per the Illumina ChIP-seq library generation protocol. The resultant DNA libraries were sent to Vanderbilt Microarray Shared Resource where they were sequenced on a Genome Analyzer II. Sequence reads (35 bases; 20 to 30 million quality filtered reads/sample) were aligned to the human genome (National Center for Biotechnology Information Build 37) using bowtie (0.10.1, , options ‘-n 2 -a —best —strata -m 1 -p 4’).
ChIP-seq peak calling and artefact filtering
VDR ChIP-seq peaks were called using Zinba (zero-inflated negative binomial algorithm,refine peaks, extension = 200) with the false discovery rate set as <0.1% . We removed peaks that showed overlap with regions known to give false positive ChIP-seq peaks by merging Terry’s blacklist and the list of ultra-high signal artefact regions . ChIP-seq peaks are detailed in the (Additional file 1) dataset. We also called peaks separately using model-based analysis of ChIP-Seq (MACS) for further motif analysis .
MEME-ChIP , Weeder  and ChIPmunk  were used to identify de novo motifs from VDR ChIP-seq peaks from groups of samples with 25-hydroxyvitamin D <75 nM and ≥75 nM, intervals overlapping with LCL/MCL VDR ChIP-seq peaks and intervals overlapping with RXR ChIP-seq peaks from NB4 cells [20, 23]. ChIP-seq peaks were also scanned for known VDR recognition motifs using RSAT  and Fimo .
GREAT gene ontology analysis
25(OH)D ≥75 and 25(OH)D <75 VDR binding sites were input into the Genomic Regions Enrichment of Annotations Tool (GREAT) using the GRCh37 (UCSC hg19, February 2009) assembly and 5 kb proximal and 1 kb distal gene windows .
Overlap and hierarchal clustering analysis
The Genomic HyperBrowser was used to determine overlap and hierarchal clustering between different datasets [27, 28]. Autoimmune disease-associated regions were determined as those 100 kb either side of a SNP associated with an autoimmune disease in the Genome Wide Association Study Catalogue with a P-value ≤1×10-7 (downloaded 13 June 2012). Samples were combined into 25(OH)D ≥75 and 25(OH)D <75 by merging all binding sites from samples with 25-hydroxyvitamin D ≥75 nM (n = 5) and <75 nM (n = 4). Overlap was determined using segment-segment analysis with either 1,000 or 10,000 Monte-Carlo randomizations maintaining the empiric distribution of segment and inter-segment lengths, but randomizing positions. Controlling for gene or immune gene position (obtained from the Gene Ontology project ) used an intensity track created based on the proximity of (pooled) VDR regions to their nearest genes or immune genes, respectively. VDR regions were represented as points (midpoints of VDR binding peaks) and a point-segment analysis using 1,000 Monte-Carlo randomizations with points sampled according to the intensity track, were used to compute P-values (auto-immune regions represented as segments as before). Immune gene-controlled overlap omitted chromosome Y as no immune genes were located there. Comparisons between 25(OH)D <75 and 25(OH)D ≥75 for overlap were performed using case-control tracks generated by the Genomic HyperBrowser and analyzed using valued segment–segment preferential overlap analysis with 10,000 Monte-Carlo randomizations, keeping the location of segments of both tracks constant, while randomly permuting case-control values of the first track in the null model. Heirarchical clustering analysis was performed in the Genomic HyperBroswer by obtaining pairwise overlap-enrichment values for each of the samples and computing distance between samples as the inverse of these values. Th1 DNase I hypersensitivity peaks were obtained from the University of California at Santa Cruz (UCSC) Table Browser and were generated by the Duke group . ChIP-seq peaks for VDR in LCLs and MCLs were obtained from previously published studies using VDR binding intervals after stimulation with calcitriol [9, 10], and co-factor ChIP-seq peaks were obtained from the Encyclopedia of DNA Elements (ENCODE) and Cistrome, using ChIP-seq data from hematopoietic cell lines (GM121878, K562 and NB4) [23, 31–33]. ChIP-seq data on chromatin states (H3K27Ac, H2A.Z, H3K4me1, H3K4me2, H3K4me3, H3K9Ac and H3K9me3) in GM12878 cells and chromatin looping 5C data were obtained from ENCODE [34, 35]. Gene expression data from CD4+ cells was obtained from data published by Birzele and colleagues . Gene expression data from LCLs in response to 1,25D3 treatment was obtained from Ramagopalan and colleagues .
VDR binding sites in CD4+ cells
Number of VDR binding sites
Up and downstream
High 25-hydroxyvitamin D
Low 25-hydroxyvitamin D
We performed hierarchical clustering analysis using pairwise overlap-enrichment of VDR binding sites and this revealed far closer similarity between samples within each group (25(OH)D ≥75 and 25(OH)D <75) than when comparing samples between groups [see Additional file 2: Figure S1]. Binding sites were also frequently shared between samples, but 66.0% of binding sites were unique to a single sample.
VDR binding and gene ontology
VDR binding sites were assessed for overlap with known gene ontology biological pathways in GREAT [See Additional file 3: Table S1] . In 25(OH)D ≥75 samples, binding sites were maximally enriched for pathways involved in RNA processing, gene expression, protein folding and T cell activation or differentiation. In contrast, the top pathways enriched for 25(OH)D <75 VDR binding were involved in RNA splicing, translation and histone modification.
VDR binding motifs
We found that there was no significant enrichment of binding sites containing DR3-like motifs either when searching de novo using MEME-ChIP , CentriMo , Weeder  or ChIPmunk  and analyzing all binding sites, binding sites grouped by high or low vitamin D, binding sites overlapping with the previous LCL or MCL VDR ChIP-seq studies, binding sites common between multiple samples or binding sites overlapping with previous ChIP-seq studies of RXR in NB4 cells . Neither were DR3-like motifs found when each sample was analyzed independently. The top consensus binding sites are shown in Additional file 4: Figure S2 for each analysis approach. Our methods were, however, able to detect the reported DR3 sites in previous VDR ChIP-seq studies [9, 10]. We were also unable to detect VDR-like motifs when restricting our search to only those parts of ChIP-seq intervals common to all samples in the 25(OH)D ≥75 or 25(OH)D <75 groups.
As this was an unexpected finding, we performed an in silico search within the pooled peaks but did not identify an over-representation of known VDR binding motifs using RSAT  and Fimo . The existing RXRA::VDR motif in the Jaspar  and TRANSFAC  databases has been generated from SELEX data, which mainly will represent strong binding without additional co-factors or other context-dependent features. It is, therefore, relevant to search for alternative variants of VDR-like motifs that may be more representative of in vivo binding. Since the CD4+ data set, in particular, shows a lack of centrally enriched binding site motifs, MEME-ChIP and CentriMo are less suitable for this. Therefore, an iterative approach was used, in which the full set of ChIP-Seq regions for LCL, MCL and the merged set of CD4+ regions was searched with MAST and the RXRA::VDR matrix (P-value 0.0001, E-value 100.0) . The significant regions were submitted to MEME for de novo motif discovery. In each data set a VDR-like motif was found. This motif was used as input to MAST again, and the resulting positive set was submitted to MEME, in order to reduce bias from the original RXRA::VDR motif. This process can, in principle, be repeated several times, but in most cases the motifs will start to degenerate after a while into very general motifs with low information content. However, the motifs generated in this case are clearly similar to the classical RXRA::VDR motif, although with distinct differences [See Additional file 5: Figure S3]. They are also similar to the previously published motifs for LCL and MCL. These improved matrices were then used with MAST to make positive and negative subsets for further analysis. Here a slightly higher P-value was used (0.0005) in order to include more borderline motifs, leading to 811 positive sequences (29%) for LCL, 648 (28%) for MCL, and 90 (0.4%) for CD4+. This seems to confirm the lack of VDR-like motifs in the CD4+ set. This was further confirmed using FIMO to search each data set with both the RXRA::VDR matrix and the individually optimized matrices generated above [See Additional file 6: Figure S4]. This showed a clear lack of significant motifs in the CD4+ data, independent of which matrix was used for searching. Analyzing CD4+ binding intervals for other JASPAR motifs showed only a significant overrepresentation of CTCF binding motifs in the 25(OH)D ≥75 but not 25(OH)D <75 group.
We found significant overlap between CD4+ VDR and RXR ChIP-seq peaks drawn from a promyelocytic cell line (NB4; Additional file 7: Table S2) (25(OH)D ≥75 19.77-fold, P= 0.0004; 25(OH)D <75 65.14-fold, P<0.0001 ) and significant overlap between VDR binding sites in CD4+ cells and those observed previously in LCLs (25(OH)D ≥75 70-fold, P<0.0001; 25(OH)D <75 151.7-fold, P<0.0001; 813/2,776 (29.3%) LCL VDR binding sites overlap with VDR binding sites in CD4+ cells) and MCLs (25(OH)D ≥75 28.75-fold, P<0.0001; 25(OH)D <75 37.17-fold, P<0.0001; 353/1,818 (19.4%) MCL VDR binding sites overlap with VDR binding sites in CD4+ cells) making it likely that our data reflect real VDR binding sites.
Motifless binding has been described by the ENCODE project with characteristically greater enrichment of DNase I hypersensitivity than binding sites with classical motifs . We confirmed this in the previous LCL and MCL VDR ChIP-seq datasets by dividing binding sites into those with or without a VDR-like motif as detailed above. Intervals containing the VDR-like motif had less enrichment of DNase I peaks in GM12878 LCLs than those intervals lacking that motif (LCL peaks with a VDR-like motif (LCLmotif), 24.6-fold, P<0.0001; LCL peaks without a VDR-like motif (LCLno motif), 27.8-fold, P<0.0001; LCLmotif versus LCLno motif P= 0.0002; MCLmotif, 13.5-fold, P<0.0001; MCLno motif, 18.0-fold, P<0.0001; MCLmotif versus MCLno motif P= 0.0002). VDR ChIP-seq peaks in the CD4+ cells in this study overlapped more with binding sites in LCLs and MCLs lacking binding motifs than those with motifs (LCLmotif 37.4-fold, P<0.0001; LCLno motif 79.4-fold, P<0.0001; LCLmotif versus LCLno motif P= 0.0002; MCLmotif, 17.7-fold, P<0.0001; MCLno motif, 32.3-fold, P<0.0001; MCLmotif versus MCLno motif P= 0.0002).
VDR co-factors, chromatin state and calcitriol-responsive gene expression
There was significant enrichment of VDR binding within 5 kb of genes responsive to 1,25D3 treatment detected from microarray expression data in LCLs (25(OH)D ≥75 3.86-fold, P<0.0001; 25(OH)D <75 2.98-fold, P= 0.0002; 25(OH)D ≥75 versus 25(OH)D <75 P= 0.004) .
Given the relatively high proportion of intergenic VDR binding sites, we tested for overlap with sites of known chromatin looping in GM12878 cells in pilot ENCODE regions . There was significant but low magnitude overlap of VDR binding and chromatin looping in 25(OH)D ≥75 samples but not 25(OH)D <75 samples (25(OH)D ≥75 1.07-fold, P= 0.002; 25(OH)D <75 0.73-fold, P= 0.83; 25(OH)D ≥75 versus 25(OH)D <75 P= 0.01).
VDR binding sites and autoimmune disease
There was no significant enrichment for genomic regions associated with control conditions (ones in which CD4+ cells would not be expected to play a dominant role), such as coronary heart disease, atopic dermatitis and type 2 diabetes mellitus (P>0.05 for all). Also, in support of separate biochemical pathways for autoimmunity and metabolic effects of vitamin D, VDR binding was not enriched for genomic regions associated with bone mineral density.
The previous study of LCLs had shown VDR enrichment near regions associated with chronic lymphocytic leukemia. However, no significant enrichment was seen for these regions in primary CD4+ cells (25(OH)D ≥75 1.62-fold, P= 0.37; 25(OH)D <75 2.44-fold, P= 0.27; LCLs 20.7-fold, P<0.0001), suggesting that VDR binding in cell lines differs considerably from that seen in primary immune cells.
Although 100 kb was chosen to encompass the likely extent of linkage disequilibrium,both groups showed increased enrichment when the size of the region assessed for overlap decreased. 25(OH)D ≥75 showed consistently greater enrichment for autoimmune regions than 25(OH)D <75 [See Additional file 10: Figure S5].
Several disease-associated SNPs were located within VDR ChIP-seq binding intervals [See Additional file 11: Table S4]. We analyzed these SNPs in Regulome DB and found that several were likely to affect gene expression and/or transcription factor binding .
VDR binding and gene expression in CD4+ cells
We assessed enrichment in VDR binding near genes expressed in different types of CD4+ cells measured by RNA-seq . VDR binding was significantly enriched within 5 kb of genes expressed either specifically in T-regulatory cells or T-helper cells and genes expressed which were common to all CD4+ cells. Enrichment was particularly high for genes associated specifically with T-regulatory and T-helper cells in the 25(OH)D ≥75 group (RNA-seq Treg: 25(OH)D ≥75 4.07-fold, P<0.0001; 25(OH)D <75 2.96-fold, P<0.0001; 25(OH)D ≥75 versus 25(OH)D <75 P= 0.0002; RNA-seq Thelper: 25(OH)D ≥75 3.87-fold, P<0.0001; 25(OH)D <75 2.76-fold, P<0.0001; 25(OH)D ≥75 versus 25(OH)D <75 P= 0.0002; RNA-seq CD4 + common: 25(OH)D ≥75 5.27-fold, P<0.0001; 25(OH)D <75 5.13-fold, P<0.0001; 25(OH)D ≥75 versus 25(OH)D <75 P= 0.0002).
The most arresting finding in this study is that the number of VDR binding sites in primary CD4+ cells is strongly correlated with 25-hydroxyvitamin D levels. The previous VDR ChIP-seq experiments using MCLs and LCLs found an increase in VDR binding site occupancy following treatment with supraphysiological levels of calcitriol [9, 10]. Our finding of a far greater number of VDR binding sites in sufficient vitamin D samples than insufficient samples suggests that this effect also occurs with different in vivo levels of vitamin D. In vivo levels of 25-hydroxyvitamin D are directly associated with the number of VDR binding sites.
VDR binding sites are enriched for markers of active transcription and open chromatin; 25(OH)D ≥75 samples seemed to be less enriched for these markers than 25(OH)D <75, perhaps reflecting binding to open chromatin state in 25(OH)D <75 samples.
We have confirmed that the observation of significant overlap between VDR binding and genomic regions implicated in autoimmune diseases in LCLs is also seen in primary CD4+ cells [9, 10]. Gene ontology analysis suggests that VDR binding in conditions of 25-hydroxyvitamin D sufficiency may be more directly related to immune cell function. This is supported by the observed higher levels of VDR binding near genes expressed specifically in T-regulatory and T-helper cells in 25(OH)D ≥75 but not 25(OH)D <75 samples.
We found a lack of classical VDR binding motifs within the VDR ChIP-seq peaks. In the ChIP-seq studies in MCLs and LCLs the authors identified classical DR3 motifs at differing proportions of sites (32% in MCLs, 67% in LCLs) with SP1-like and ETS-like non-classical peaks identified in the MCL ChIP-seq study (23% and 12% respectively) [9, 10]. We found enrichment of CTCF motifs in several of our samples but were unable to identify any previously described VDR motifs. One possibility is that in vivo VDR binding is modulated by protein-protein interactions with co-factors: SP1 and ETS1 are known to modulate VDR binding, and there is some evidence that interactions between SP1 and VDR may enable modulation of genes that lack a classical VDR recognition motif [43, 44]. Several other proteins are known to bind in association with VDR, including NR4A1 and c-MYC [45, 46]. CTCF is known to modulate DNA binding via protein-protein interactions with other nuclear receptors [47–49]. However, it is unlikely that protein-protein interactions with transcription factors with specific recognition sequences can explain most of these motifless binding sites since one would have expected to find that motif through MEME-ChIP analysis. It may be that in response to physiological levels of 25-hydroxyvitamin D most VDR binding occurs at motifless binding sites similar to those identified by ENCODE , supported by the increased overlap with DNase I peaks. Another possibility is that the lack of motifs may reflect the fact that these CD4+ cells were not stimulated with 1,25D3, as the previous LCL ChIP-seq did not find classical motifs prior to stimulation . Alternatively, current motif-finding methods may be insufficient to locate true VDR binding motifs. Further research will be needed in more lymphocyte subsets to delineate further the role of non-classical binding sites in VDR binding. It would also be useful to obtain 1,25D3, parathyroid hormone and calcium measurements for future study.
The overlap between genomic regions associated with many autoimmune diseases and VDR binding in primary CD4+ cells strongly suggests a role for vitamin D in many of these diseases, as already seen for MCLs and LCLs [9, 10]. This is strengthened by the observation that this effect tends to be stronger in individuals sufficient for 25-hydroxyvitamin D. Interestingly, the magnitude of enrichment for autoimmunity increased as the flanks of the region surrounding implicated SNPs was reduced. This further suggests that this is not a chance finding and that VDR binding may have a functional role in modulating adaptive immunity in autoimmune diseases. We also controlled for genomic architectural features that could bias our results and observed that the results were not substantively altered. Future functional work should focus on the effects of VDR binding on nearby gene expression and targeted sequencing in patients with autoimmune conditions to identify possible rare variants affecting VDR binding.
The role of vitamin D in bone health has long been established. The involvement of this vitamin in autoimmune disease is however heavily debated. We provide here an in vivo mechanism as to how vitamin D deficiency may influence autoimmune disease risk, by directly interacting with disease associated genes. Vitamin D sufficiency has been suggested to have a threshold of approximately 75 nmol/L; we provide here biological evidence in support of this, with significant public health implications.
All subjects gave written informed consent for their samples to be used in this study.
25D3: 1,25-dihydroxyvitamin D
- 25-OH D:
Chromatin immunoprecipitation and massively parallel sequencing
Encyclopedia of DNA Elements
Lymphoblastoid cell line
Magnetic activated cell sorting
Monocytic cell line
Retinoid X receptors
Single nucleotide polymorphism
Vitamin D receptor
- 25(OH)D ≥ 75:
Samples with 25-hydroxyvitamin D ≥75 nM
This work was supported by the Medical Research Council and Queen Mary University of London. AEH was supported by an NIHR Academic Clinical Fellowship.
- Lehmann B, Meurer M: Vitamin D metabolism. Dermatol Ther. 2010, 23: 2-12. 10.1111/j.1529-8019.2009.01286.x.View ArticlePubMedGoogle Scholar
- Handel AE, Handunnetthi L, Ebers GC, Ramagopalan SV: Type 1 diabetes mellitus and multiple sclerosis: common etiological features. Nat Rev Endocrinol. 2009, 5: 655-664. 10.1038/nrendo.2009.216.View ArticlePubMedGoogle Scholar
- Kriegel MA, Manson JE, Costenbader KH: Does vitamin D affect risk of developing autoimmune disease?: a systematic review. Semin Arthritis Rheum. 2011, 40: 512-531. 10.1016/j.semarthrit.2010.07.009.View ArticlePubMedGoogle Scholar
- Munger KL, Levin LI, Hollis BW, Howard NS, Ascherio A: Serum 25-hydroxyvitamin D levels and risk of multiple sclerosis. JAMA. 2006, 296: 2832-2838. 10.1001/jama.296.23.2832.View ArticlePubMedGoogle Scholar
- Orgaz-Molina J, Buendía-Eisman A, Arrabal-Polo MA, Ruiz JC, Arias-Santiago S: Deficiency of serum concentration of 25-hydroxyvitamin D in psoriatic patients: a case-control study. J Am Acad Dermatol. 2012, 67: 931-938. 10.1016/j.jaad.2012.01.040.View ArticlePubMedGoogle Scholar
- Li K, Shi Q, Yang L, Li X, Liu L, Wang L, Li Q, Wang G, Li C-Y, Gao T-W: The association of vitamin D receptor gene polymorphisms and serum 25-hydroxyvitamin D levels with generalized vitiligo. Br J Dermatol. 2012, 167: 815-821. 10.1111/j.1365-2133.2012.11132.x.View ArticlePubMedGoogle Scholar
- McCullough ML, Bostick RM, Daniel CR, Flanders WD, Shaukat A, Davison J, Rangaswamy U, Hollis BW: Vitamin D status and impact of vitamin D3 and/or calcium supplementation in a randomized pilot study in the Southeastern United States. J Am Coll Nutr. 2009, 28: 678-686. 10.1080/07315724.2009.10719801.View ArticlePubMedPubMed CentralGoogle Scholar
- Issa LL, Leong GM, Eisman JA: Molecular mechanism of vitamin D receptor action. Inflamm Res. 1998, 47: 451-475. 10.1007/s000110050360.View ArticlePubMedGoogle Scholar
- Ramagopalan SV, Heger A, Berlanga AJ, Maugeri NJ, Lincoln MR, Burrell A, Handunnetthi L, Handel AE, Disanto G, Orton S-M, Watson CT, Morahan JM, Giovannoni G, Ponting CP, Ebers GC, Knight JC: A ChIP-seq defined genome-wide map of vitamin D receptor binding: associations with disease and evolution. Genome Res. 2010, 20: 1352-1360. 10.1101/gr.107920.110.View ArticlePubMedPubMed CentralGoogle Scholar
- Heikkinen S, Väisänen S, Pehkonen P, Seuter S, Benes V, Carlberg C: Nuclear hormone 1α,25-dihydroxyvitamin D3 elicits a genome-wide shift in the locations of VDR chromatin occupancy. Nucleic Acids Res. 2011, 39: 9181-9193. 10.1093/nar/gkr654.View ArticlePubMedPubMed CentralGoogle Scholar
- Carlberg C, Seuter S, Heikkinen S: The first genome-wide view of vitamin D receptor locations and their mechanistic implications. Anticancer Res. 2012, 32: 271-282.PubMedGoogle Scholar
- Shirley MD, Baugher JD, Stevens EL, Tang Z, Gerry N, Beiswanger CM, Berlin DS, Pevsner J: Chromosomal variation in lymphoblastoid cell lines. Hum Mutat. 2012, 33: 1075-1086. 10.1002/humu.22062.View ArticlePubMedPubMed CentralGoogle Scholar
- Caliskan M, Cusanovich DA, Ober C, Gilad Y: The effects of EBV transformation on gene expression levels and methylation profiles. Hum Mol Genet. 2011, 20: 1643-1652. 10.1093/hmg/ddr041.View ArticlePubMedPubMed CentralGoogle Scholar
- Owen TA, Aronow MS, Barone LM, Bettencourt B, Stein GS, Lian JB: Pleiotropic effects of vitamin D on osteoblast gene expression are related to the proliferative and differentiated state of the bone cell phenotype: dependency upon basal levels of gene expression, duration of exposure, and bone matrix competency in normal rat osteoblast cultures. Endocrinology. 1991, 128: 1496-1504. 10.1210/endo-128-3-1496.View ArticlePubMedGoogle Scholar
- Berglund D, Korsgren O, Lorant T, Schneider K, Tufveson G, Carlsson B: Isolation, expansion and functional assessment of CD4 + CD25 + FoxP3+ regulatory T cells and Tr1 cells from uremic patients awaiting kidney transplantation. Transpl Immunol. 2012, 26: 27-33. 10.1016/j.trim.2011.09.003.View ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.View ArticlePubMedPubMed CentralGoogle Scholar
- Rashid NU, Giresi PG, Ibrahim JG, Sun W, Lieb JD: ZINBA integrates local covariates with DNA-seq data to identify broad and narrow regions of enrichment, even within amplified genomic regions. Genome Biol. 2011, 12: R67-10.1186/gb-2011-12-7-r67.View ArticlePubMedPubMed CentralGoogle Scholar
- Kundaje A: A comprehensive collection of signal artifact blacklist regions in the human genome. ENCODE. [hg19-blacklist-README.doc - EBI]. https://sites.google.com/site/anshulkundaje/projects/blacklists (last accessed 30/06/2013)
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.View ArticlePubMedPubMed CentralGoogle Scholar
- Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, 37: W202-W208. 10.1093/nar/gkp335.View ArticlePubMedPubMed CentralGoogle Scholar
- Pavesi G, Mereghetti P, Mauri G, Pesole G: Weeder Web: discovery of transcription factor binding sites in a set of sequences from co-regulated genes. Nucleic Acids Res. 2004, 32: W199-W203. 10.1093/nar/gkh465.View ArticlePubMedPubMed CentralGoogle Scholar
- Kulakovskiy IV, Boeva VA, Favorov AV, Makeev VJ: Deep and wide digging for binding motifs in ChIP-Seq data. Bioinformatics. 2010, 26: 2622-2623. 10.1093/bioinformatics/btq488.View ArticlePubMedGoogle Scholar
- Martens JHA, Brinkman AB, Simmer F, Francoijs K-J, Nebbioso A, Ferrara F, Altucci L, Stunnenberg HG: PML-RARalpha/RXR alters the epigenetic landscape in acute promyelocytic leukemia. Cancer Cell. 2010, 17: 173-185. 10.1016/j.ccr.2009.12.042.View ArticlePubMedGoogle Scholar
- Thomas-Chollier M, Herrmann C, Defrance M, Sand O, Thieffry D, van Helden J: RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets. Nucleic Acids Res. 2012, 40: e31-10.1093/nar/gkr1104.View ArticlePubMedGoogle Scholar
- Grant CE, Bailey TL, Noble WS: FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011, 27: 1017-1018. 10.1093/bioinformatics/btr064.View ArticlePubMedPubMed CentralGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28: 495-501. 10.1038/nbt.1630.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandve GK, Gundersen S, Rydbeck H, Glad IK, Holden L, Holden M, Liestøl K, Clancy T, Ferkingstad E, Johansen M, Nygaard V, Tøstesen E, Frigessi A, Hovig E: The Genomic HyperBrowser: inferential genomics at the sequence level. Genome Biol. 2010, 11: R121-10.1186/gb-2010-11-12-r121.View ArticlePubMedPubMed CentralGoogle Scholar
- Gundersen S, Kalaš M, Abul O, Frigessi A, Hovig E, Sandve GK: Identifying elemental genomic track types and representing them uniformly. BMC Bioinformatics. 2011, 12: 494-10.1186/1471-2105-12-494.View ArticlePubMedPubMed CentralGoogle Scholar
- Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA. 2009, 106: 9362-9367. 10.1073/pnas.0903103106.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology, The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujita PA, Rhead B, Zweig AS, Hinrichs AS, Karolchik D, Cline MS, Goldman M, Barber GP, Clawson H, Coelho A, Diekhans M, Dreszer TR, Giardine BM, Harte RA, Hillman-Jackson J, Hsu F, Kirkup V, Kuhn RM, Learned K, Li CH, Meyer LR, Pohl A, Raney BJ, Rosenbloom KR, Smith KE, Haussler D, Kent WJ: The UCSC Genome Browser database: update 2011. Nucleic Acids Res. 2011, 39: D876-D882. 10.1093/nar/gkq963.View ArticlePubMedGoogle Scholar
- Liu T, Ortiz JA, Taing L, Meyer CA, Lee B, Zhang Y, Shin H, Wong SS, Ma J, Lei Y, Pape UJ, Poidinger M, Chen Y, Yeung K, Brown M, Turpaz Y, Liu XS: Cistrome: an integrative platform for transcriptional regulation studies. Genome Biol. 2011, 12: R83-10.1186/gb-2011-12-8-r83.View ArticlePubMedPubMed CentralGoogle Scholar
- The ENCODE (ENCyclopedia Of DNA Elements) Project. Science. 2004, 306: 636-640.
- Sanyal A, Lajoie BR, Jain G, Dekker J: The long-range interaction landscape of gene promoters. Nature. 2012, 489: 109-113. 10.1038/nature11279.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, Snyder M: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489: 57-74. 10.1038/nature11247.View ArticleGoogle Scholar
- Birzele F, Fauti T, Stahl H, Lenter MC, Simon E, Knebel D, Weith A, Hildebrandt T, Mennerich D: Next-generation insights into regulatory T cells: expression profiling and FoxP3 occupancy in Human. Nucleic Acids Res. 2011, 39: 7946-7960. 10.1093/nar/gkr444.View ArticlePubMedPubMed CentralGoogle Scholar
- Souberbielle J-C, Body J-J, Lappe JM, Plebani M, Shoenfeld Y, Wang TJ, Bischoff-Ferrari HA, Cavalier E, Ebeling PR, Fardellone P, Gandini S, Gruson D, Guérin AP, Heickendorff L, Hollis BW, Ish-Shalom S, Jean G, von Landenberg P, Largura A, Olsson T, Pierrot-Deseilligny C, Pilz S, Tincani A, Valcour A, Zittermann A: Vitamin D and musculoskeletal health, cardiovascular disease, autoimmunity and cancer: recommendations for clinical practice. Autoimmun Rev. 2010, 9: 709-715. 10.1016/j.autrev.2010.06.009.View ArticlePubMedGoogle Scholar
- Bailey TL, Machanick P: Inferring direct DNA binding from ChIP-seq. Nucleic Acids Res. 2012, 40: e128-10.1093/nar/gks433.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandelin A, Wasserman WW: Constrained binding site diversity within families of transcription factors enhances pattern discovery bioinformatics. J Mol Biol. 2004, 338: 207-215. 10.1016/j.jmb.2004.02.048.View ArticlePubMedGoogle Scholar
- Wingender E, Dietze P, Karas H, Knüppel R: TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res. 1996, 24: 238-241. 10.1093/nar/24.1.238.View ArticlePubMedPubMed CentralGoogle Scholar
- Bailey TL, Gribskov M: Methods and statistics for combining motif match scores. J Comput Biol. 1998, 5: 211-221. 10.1089/cmb.1998.5.211.View ArticlePubMedGoogle Scholar
- Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, Karczewski KJ, Park J, Hitz BC, Weng S, Cherry JM, Snyder M: Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 2012, 22: 1790-1797. 10.1101/gr.137323.112.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang Y-C, Chen J-Y, Hung W-C: Vitamin D3 receptor/Sp1 complex is required for the induction of p27Kip1 expression by vitamin D3. Oncogene. 2004, 23: 4856-4861. 10.1038/sj.onc.1207621.View ArticlePubMedGoogle Scholar
- Tolón RM, Castillo AI, Jiménez-Lara AM, Aranda A: Association with Ets-1 causes ligand- and AF2-independent activation of nuclear receptors. Mol Cell Biol. 2000, 20: 8793-8802. 10.1128/MCB.20.23.8793-8802.2000.View ArticlePubMedPubMed CentralGoogle Scholar
- Marshall PA, Hernandez Z, Kaneko I, Widener T, Tabacaru C, Aguayo I, Jurutka PW: Discovery of novel vitamin D receptor interacting proteins that modulate 1,25-dihydroxyvitamin D3 signaling. J Steroid Biochem Mol Biol. 2012, 132: 147-159. 10.1016/j.jsbmb.2012.05.001.View ArticlePubMedPubMed CentralGoogle Scholar
- Salehi-Tabar R, Nguyen-Yamamoto L, Tavera-Mendoza LE, Quail T, Dimitrov V, An B-S, Glass L, Goltzman D, White JH: Vitamin D receptor as a master regulator of the c-MYC/MXD1 network. Proc Natl Acad Sci U S A. 2012, 109: 18827-18832. 10.1073/pnas.1210037109.View ArticlePubMedPubMed CentralGoogle Scholar
- Weth O, Weth C, Bartkuhn M, Leers J, Uhle F, Renkawitz R: Modular insulators: genome wide search for composite CTCF/thyroid hormone receptor binding sites. PLoS One. 2010, 5: e10119-10.1371/journal.pone.0010119.View ArticlePubMedPubMed CentralGoogle Scholar
- Ross-Innes CS, Brown GD, Carroll JS: A co-ordinated interaction between CTCF and ER in breast cancer cells. BMC Genomics. 2011, 12: 593-10.1186/1471-2164-12-593.View ArticlePubMedPubMed CentralGoogle Scholar
- Lutz M, Burke LJ, LeFevre P, Myers FA, Thorne AW, Crane-Robinson C, Bonifer C, Filippova GN, Lobanenkov V, Renkawitz R: Thyroid hormone-regulated enhancer blocking: cooperation of CTCF and thyroid hormone receptor. EMBO J. 2003, 22: 1579-g1587. 10.1093/emboj/cdg147.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/11/163/prepub