Previous Article | Next Article ![]()
Molecular and Cellular Biology, September 2006, p. 6469-6486, Vol. 26, No. 17
0270-7306/06/$08.00+0 doi:10.1128/MCB.00353-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
B Ligand Gene Expression by 1,25-Dihydroxyvitamin D3 Is Mediated through Multiple Long-Range Enhancers
Department of Biochemistry, University of WisconsinMadison, Madison, Wisconsin 53706
Received 27 February 2006/ Returned for modification 4 May 2006/ Accepted 12 June 2006
|
|
|---|
|
|
|---|
Although many factors participate in osteoclastogenesis, the molecule that is now considered to be both necessary and sufficient in vivo and in vitro is the receptor activator of NF-
B ligand (RankL). RankL is a tumor necrosis factor (TNF)-like factor that is produced by stromal cells and osteoblasts as well as a variety of other cell types (38). This factor not only actively promotes the process of osteoclast differentiation, but also is required for the cell's bone-resorbing activity and for its survival (22). The interaction of RankL with receptor activator of NF-
B (Rank), an integral receptor protein located on the surface of osteoclast precursors, triggers a number of signaling cascades that include the IKK/IKß/NF-
B transduction pathway and the mitogen-activated protein kinase (MAPK), Src, and phosphatidylinositol 3-kinase (PI3K)/AKT pathways as well (64). A novel calcium oscillation pathway that involves ITAM coreceptors is also involved in the downstream effects of RankL (35). Importantly, stimulation of these pathways culminates in the activation of multiple transcription factors, including c-fos, NF-
B, and NFATc1, all of which play strategic roles in the differentiation process at the genetic level (44, 61). Overall, the timely activation of these transcription factors initiates growth arrest and promotes osteoclast differentiation, fusion, activation, and survival (63). The evidence that supports the essentiality of both RankL and its receptor in osteoclast formation is most strongly supported by the skeletal phenotypes of both RankL- and Rank-null mice, neither of which are capable of producing osteoclasts in vivo and thus are phenotypically osteopetrotic (20, 37).
RankL is synthesized and expressed on the surface of regulatory cells in response to a myriad of both local and systemic factors, many of which are essential to physiologic bone turnover. These include the two hormones integral to calcium homeostasis, 1,25-dihydroxyvitamin D3 [1,25(OH)2D3] (60) and parathyroid hormone (PTH) (36, 41). RankL expression can also be influenced by the glucocorticoid (GC) stress hormones (18, 54), inflammatory cytokines such as TNF-
and interleukin-1 (IL-1) (23), the gp130-activating cytokines IL-6 and IL-11 (51, 62), certain prostaglandins (66), and transforming growth factor ß (TGF-ß) (28). While these factors can function in a physiologic setting, their activities are often manifested during disease. Thus, RankL overproduction can in a number of circumstances lead to the pathological bone resorption associated with age-related and postmenopausal bone disease, rheumatoid arthritis and osteoarthritis, multiple myelomas, diabetic neuropathy, metastatic cancer, general hypercalcemia of malignancy, and a variety of other syndromes that impact the skeleton (19).
The calciotropic hormones 1,25(OH)2D3 and PTH, as well as certain growth factors, cytokines, and prostaglandins, all regulate the expression of RankL from stromal cells and osteoblasts (63). Interestingly, the molecular mechanisms responsible for activation by these regulators remain largely unknown. With respect to vitamin D, 1,25(OH)2D3 is known to induce RankL upregulation primarily through actions initiated by the vitamin D receptor (VDR) (29). Kitazawa and colleagues (34) reported a 1,25(OH)2D3 response in the mouse RankL gene promoter and mapped this activity to a 16-bp vitamin D response element (VDRE) that was located 935 bp upstream of the transcriptional start site (TSS). Inducibility was modest, however, and subsequent studies by a number of laboratories have failed to confirm this finding (11, 48). More recently, Kabe et al. (24) have explored the ability of a potential VDRE located immediately downstream of the RankL TSS to mediate 1,25(OH)2D3 activity. Although this element is similar to that of a consensus VDRE, it is not conserved within RankL genes of other species and is capable of only a modest activity in the context of the RankL proximal promoter. Thus, it remains unclear at this stage whether or how these two elements contribute to the RankL gene's regulation by 1,25(OH)2D3.
PTH also regulates RankL expression (40), although the mechanism through which this peptide acts has remained equally elusive. Recent studies by O'Brien and coworkers (12) have established that PTH is capable of both stabilizing RankL mRNA and inducing its expression, the latter via stimulation of the protein kinase A (PKA) pathway and activation of the CREB transcription factor. It seems likely that activation of this transcription factor may underlie the ability of the prostaglandin PGE2 to stimulate RankL expression as well (25). Despite these insights, however, the mechanism by which CREB induces transcriptional activation at the level of the RankL gene promoter remains undefined. The accompanying article by Fu, Manolagas, and O'Brien delineates the molecular mechanism whereby PTH stimulates RankL gene expression (13).
Finally, the actions of GCs on the skeleton are highly complex. Enhanced exposure to these hormones, however, generally results in osteoporosis (42). At the level of the osteoblast, GCs reduce both the functional capacity of these cells to form bone and the period of time during which this functionality occurs (43, 67). GCs also influence the production of osteoclasts (18, 54). These actions occur through a direct activity on the osteoblast primarily to suppress the expression of osteoprotegerin, a decoy receptor which blocks RankL activity, but also to increase RankL gene expression. Although a role for the glucocorticoid receptor (GR) in RankL induction is likely, the mechanism whereby the stress hormone induces RankL expression remains unexplored.
As outlined above, the mechanism of action of 1,25(OH)2D3 involves a ligand-initiated interaction between the VDR and the regulatory regions of target genes wherein VDR functions together with its retinoid X receptor (RXR) partner to recruit coregulatory complexes that are essential for modulation of transcriptional output (46, 58). Binding sites for the VDR have been found in a number of genes and are generally, although not always, comprised of two hexanucleotide half-sites separated by a short spacer and frequently located within the first kilobase or two upstream of the transcriptional start site (50). The absence of bona fide target sites for VDR action within the first 8 kb of the RankL gene, as reported by us and others (11, 48), prompted a more expansive approach to delineate RankL regulatory regions. We therefore used contemporary chromatin immunoprecipitation (ChIP)/chip analysis (as described below) in this endeavor. We discovered five sites of action located at increasing distances from the mouse RankL gene TSS, the furthest approximately 76 kb upstream. The latter region, which we termed the RankL distal control region (RL-DCR), was transcriptionally active in transfection studies and contained an unusual VDRE sequence. Our studies, and the results described in the accompanying article by O'Brien and colleagues (13), provide essential clues as to how 1,25(OH)2D3 and PTH regulate RankL induction.
|
|
|---|
-MEM) was purchased from Mediatech (Herndon, VA), and minimum essential medium alpha (MEM-
) was obtained from Invitrogen Corporation (Carlsbad, CA). Oligonucleotide primers were obtained from IDT (Coralville, IA). Anti-VDR (Sc-1008), -RXR (Sc-774), -C/EBPß (Sc-150), and -GR (Sc-1004) antibodies were obtained from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA). Anti-acetyl H4 antibody (06-866) was obtained from Upstate (Charlottesville, VA), and anti-RNA polymerase II antibody (8WG16) was obtained from Berkeley Antibody Company (Richmond, CA). Lipofectamine Plus was obtained from Invitrogen Corporation (Carlsbad, CA). [
-32P]dATP was obtained from NEN Life Science Products, Inc. (Boston, MA). Dexamethasone (D4902), anti-rat immunoglobulin G (IgG; R5128), and RU486 (M8046) were purchased from Sigma Chemical Co. (St. Louis, MO).
Cell culture.
Mouse MC3T3-E1 and ST2 osteoblastic cells were cultured in
-MEM and MEM-
, respectively. Primary calvarial osteoblasts (mOBs) were obtained as previously described (53) and cultured in
-MEM. Human osteosarcoma MG63 cells were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with 1% nonessential amino acids. COS-7 fibroblasts were also cultured in DMEM. Each medium was supplemented with 10% fetal bovine serum (FBS) obtained from HyClone (Logan, UT), 100 U/ml penicillin, and 100 µg/ml streptomycin. All ligands were added in ethanol (0.1% maximum final concentration) or dimethyl sulfoxide.
RNA isolation and analysis. Total RNA was isolated from cells using Triazol reagent obtained from MRC (Cincinnati, OH). The isolated RNA was reverse transcribed using the SuperScript III RNase H reverse transcriptase kit from Invitrogen (Carlsbad, CA) and then subjected to PCR analysis using standard PCR methods. Primers used include those for the mouse ß-actin gene (mß-actin) (forward, TGTTTGAGACCTTCAACACCC; reverse, CGTTGCCAATAGTGATGACCT), mCyp24a1 (forward, GTGCGGATTTCCTTTGTGATA; reverse, GGTAGCGTGTATTCACCCAGA), mVDR (forward, TCACTGATGTCTCCAGAGCTGGGC; reverse, TGGATAGGCGGTCCTGAATGGC), and mOpn (forward, CTAACTACGACCATGAGATTGGCAG; reverse, CTTTAGTTGACCTCAGAAGATGAA) and mRankL (forward, GAATCCTGAGACTCCATGAAAACGC; reverse, CCATGAGCCTTCCATCATAGCTGG). Primers for the human genes used included those for the human ß-actin gene (hß-actin) (forward, TTAGTTGCGTTACACCCTTTC; reverse, GTCACCTTCACCGTTCCAGTT), hCYP24A1 (forward, CTTTGCTTCCTTTTCCCAGAAT; reverse, CGCCGTAGATGTCACCAGTC), and hRANKL (forward, AACAGGCCTTTCAAGGAGCTGTGC; reverse, AAGAGGACAGACTCACTTTATGGG).
siRNA studies. All small interfering RNA (siRNA) duplexes were obtained from Dharmacon RNA Technologies (Lafayette, CO). ST2 cells were seeded into six-well plates at a concentration of 1.5 x 105 cells/well and transfected approximately 24 h later using Lipofectamine Plus in serum and antibiotic-free medium. A 20 nM concentration of mVDR siRNA (D-058923-01), nontargeting siRNA pool (D-001206-13), or cyclophilin B siRNA (D-001136-01) was used for transfection. After transfection, the cells were cultured in medium supplemented with 10% FBS for 48 h before they were treated with a routine concentration of 107 M 1,25(OH)2D3 for 6 h. RNA isolation and standard PCR analysis were carried out using the primers listed above. Western blot analysis confirmed depletion of VDR protein expression in ST2 cells (data not shown).
ChIP assay. Chromatin immunoprecipitation assays were performed as previously described (30). Primer sets used for amplifying mouse and human Cyp24a1, osteopontin (Opn), and RankL gene regions of interest are all listed in Table 1. Densitometric analysis was carried out using Kodak ID Image Analysis (software version 3.5).
|
View this table: [in a new window] |
TABLE 1. Primers used for ChIP analysis
|
500-bp amplicons were then labeled with Cy3 or Cy5 dyes using an indirect labeling protocol. In this method, biotinylated dUTP is first incorporated into the individual amplicons by standard procedures and the conjugated DNA is labeled subsequently with either Cy5- or Cy3-conjugated streptavidin. Cy3- and Cy5-labeled DNA samples are then mixed in the presence of CoT-1 DNA, denatured, and cohybridized to custom oligonucleotide microarrays (Nimblegen Systems Inc., Madison, WI) as described. The microarrays are washed extensively and scanned using an Axon 4000B scanner at the appropriate wavelengths. Custom oligonucleotide arrays were synthesized by Nimblegen Systems, Inc (Madison, WI). The microarray probes consisted of maskless-array, in situ-synthesized 50-mer oligonucleotides at 2-bp intervals representing a screen of over 300 kb of the mouse RANKL gene locus from 200 kb upstream of the gene's TSS to 100 kb downstream of the final 3' noncoding exon. The tiled arrays were synthesized in duplicate in both the forward as well as reverse directions, providing four independent measurements at each site within the gene. In addition, each analysis was carried out using two independently derived ChIP DNA samples. Both the Cyp24a1 and Opn genes (as well as other candidate VDR target genes that are not discussed) were tiled in a similar fashion. A series of comparisons were made between (i) IgG in the presence or absence of hormone, (ii) VDR in the presence or absence of hormone, (iii) VDR in the presence of hormone versus input DNA, (iv) RXR in the presence or absence of hormone, and (v) RXR in the presence of hormone versus input DNA. After sample cohybridization, the logarithmic enrichment ratios of Cy5 to Cy3 hybridization intensities (log2) were plotted as a function of chromosome nucleotide position. While all of the peaks representative of enhanced VDR or RXR binding either in the presence of 1,25(OH)2D3 or as compared to input DNA are presented as the raw data, a peak-finding algorithm was utilized to score the relative levels of binding between the three regions identified (14).
Plasmids.
Full-length hVDR and hRXR
were cloned into the pET-29b vector obtained from Novagen (Darmstadt, Germany) and expressed with C-terminal His6 tags. The pCH110-ß-galactosidase (pCH110-ßgal) reporter plasmid and the pcDNA-hVDR vector or a mutant (hVDR L417A/E420A) version [pcDNA-hVDR(m)] and pRSV-GR
were previously described (70). The parent thymidine kinase (TK) and luciferase (luc) vectors pTK-luc and pGL3-luc were utilized in subsequent cloning efforts. mRLD1 (16.4 to 15.2), mRLD2 (23.1 to 21.5), mRLD3 (60.4 to 59.3), mRLD4 (69.0 to 68.1), and mRLD5 (76045 to 74973) were amplified using primers that contained HindIII, HindIII/SalI, or HindIII/BamHI restriction ends and then cloned into the corresponding sites within the pTK-luc vector. pTK-mRLD5 fragments pTK-mRLD5-5/1 (75724 to 74973), pTK-mRLD5-5/2 (75475 to 74973), pTK-mRLD5-5/3 (75228 to 74973), pTK-mRLD5-3/1 (75724 to 75228), and pTK-mRLD5-3/2 (75724 to 75475) were subcloned similarly into the pTK-luc vector using HindIII/SalI restriction sites. The hRLD5 region of the human RANKL gene (96903 to 95805) was amplified from genomic DNA and cloned into the HindIII/BamHI sites of the TK-luc vector to produce pTK-hRLD5. mRL-VDRE (75620 to 75590) and the hRL-VDRE (96467 to 96431) as well as mRL-VDRE1 and mRL-VDRE2, each containing several overhanging 5' and 3' nucleotides, were synthesized, annealed, and similarly cloned into the HindIII/BamHI sites of pTK-luc. Triplet mutations in the mRL-VDRE half-sites and in the hRL-VDRE half-sites in the context of pTK-mRLD5 and pTK-hRLD5 were created using the Quikchange mutagenesis kit from Stratagene (San Diego, CA). mRL-VDRE segments containing mutations within each half-site were also synthesized and cloned into the HindIII/BamHI sites of pTK-luc. The pmRL(100) vector was prepared by introducing an amplified segment of the mouse RankL gene promoter (101 to +54 relative to the RankL TSS) into the pGL3-luc expression vector at the XhoI/HindIII sites. Each of the mRLD regions, mRLD1 (16.4 to 15.2), mRLD2 (23.1 to 21.5), mRLD3 (60.4 to 59.3), mRLD4 (69.0 to 68.1), and mRLD5 (76045 to 74973), as well as mRL-VDRE (75620 to 75590), was then amplified and cloned upstream of pmRL(100) using the MluI/XhoI restriction sites. All plasmid constructs were sequenced to verify successful cloning.
Transfection assays.
MC3T3-E1 and/or ST2 cells were seeded into 24-well plates at appropriate densities and cultured in
-MEM or MEM-
containing 10% FBS. Cells were transfected 24 h later with Lipofectamine Plus in serum and antibiotic-free medium. Individual wells were transfected with 250 ng of a luciferase reporter vector, 50 ng of pCH110-ßgal, and 50 ng of pcDNA-hVDR (which was routinely transfected with all luciferase reporters unless otherwise indicated) or 50 ng phRSV-GR
(added only when specifically indicated). Nontargeting, cyclophilin B or VDR siRNA pools (50 nM) were also transfected where indicated. After transfection, the cells were cultured first for 48 h in a medium supplemented with 20% FBS and subsequently for an additional 24 h with or without 1,25(OH)2D3. Cells were then harvested, and the lysates were assayed for luciferase and ß-galactosidase activities as previously described (70). Luciferase activity was normalized to ß-galactosidase activity in all cases.
Protein purification.
Human VDR and RXR
proteins were produced using the bacterial expression vectors pET-hVDR and pET-hRXR
in BL21(DE3) codon Plus RIL cells obtained from Stratagene (San Diego, CA). Soluble full-length hVDR and hRXR
proteins were purified to homogeneity using sequential Ni-nitrilotriacetic acid (NTA) and SP-Sepharose column chromatography (70). Two forms of RXR
were present due to redundant start sites.
DNA band shift analysis.
The duplex oligonucleotide probes comprised of the mouse osteopontin VDRE, mRL-VDRE, mRL-VDRE1m (M1), mRL-mVDRE2m (M2), and hRL-VDRE, as documented in Table 2, were end labeled using [
-32P]dATP. Probes were incubated at room temperature with the indicated concentrations of hVDR and hRXR
in 10 mM HEPES, pH 7.4, 1 mM EDTA, 5 mM MgCl2, 10% glycerol, 0.5 mM dithiothreitol, 0.7 mM phenylmethylsulfonyl fluoride, and 50 or 150 mM KCl in the absence or presence of 1,25(OH)2D3 for 30 min. Complexes were resolved on nondenaturing 6% polyacrylamide gels, dried, and then visualized using autoradiography. Densitometric analyses of complex 1 and complex 2 were carried out using Kodak ID Image Analysis (software version 3.5).
|
View this table: [in a new window] |
TABLE 2. DNA sequences for EMSA
|
|
|
|---|
![]() View larger version (31K): [in a new window] |
FIG. 1. Induction of mRankL mRNA by 1,25(OH)2D3 in ST2 cells is enhanced by DEX and mediated by the VDR. (A) Induction of RankL expression levels by 1,25(OH)2D3 and DEX in vitro. ST2 cells were treated for periods of up to 24 h with either 1,25(OH)2D3 (107 M), DEX (107 M), or both (at 107 M). Total RNA was isolated and subjected to reverse transcription-PCR (RT-PCR) analysis using primers specific to mouse Cyp24a1 (30 cycles), osteopontin (Opn) (15 cycles), RankL (30 cycles), or ß-actin (20 cycles) as documented in Materials and Methods. The results are typical of multiple similar experiments. (B) Effects of mVDR siRNA on 1,25(OH)2D3-induced RankL expression levels in ST2 cells. ST2 cells were transfected with 20 nM nontargeting siRNA, cyclophilin B (Cyclo B) siRNA, or mVDR siRNA. After 48 h, the cells were treated for an additional 6 h with either vehicle or 1,25(OH)2D3 (107 M). Total RNA was isolated and subjected to standard RT-PCR analysis using the primers identified in panel A above. The numbers of cycles for amplification of each transcript are as follows: VDR, 20 cycles; Cyp24a1, 25 cycles; Opn, 15 cycles, and ß-actin, 15 cycles. These results are typical of several independent experiments.
|
![]() View larger version (35K): [in a new window] |
FIG. 2. ChIP/chip analysis reveals five VDR/RXR-interacting regions at significant distances upstream of the mRankL gene TSS. (A) ST2 cells were treated with either vehicle or 1,25(OH)2D3 (107 M) for 6 h and then subjected to ChIP analysis using antibodies to VDR, RXR(pan), or control IgG ( VDR, RXR, and IgG, respectively). Immunoprecipitated DNA was isolated and then amplified using ChIP primer sets to either Cyp24a1 or Opn as indicated in Materials and Methods and Table 1. Input DNA was obtained prior to precipitation. (B) Schematic diagram of the mouse RankL gene with its five exons and its position relative to adjacent downstream (AK034132) and upstream (AK129178) genes on chromosome 14. The reverse arrow indicates the direction of transcription on the reverse strand. The nucleotide base pairs indicate nucleotide location on chromosome 14 (December 2004 assembly). (C, upper panel) Individual data tracks representing the enrichment ratio of Cy5-to-Cy3 hybridization intensity (log2) for IgG ± hormone or VDR ± hormone. The nucleotide base pairs on the x axis indicate the position on chromosome 14. (C, lower panel) An expanded view of the Cy5/Cy3 signal enrichment ratios for VDR ± hormone, VDR ± versus input, and RXR plus hormone versus input. The highlighted areas indicate the peaks of interest which are designated mRLD1 to mRLD5.
|
![]() View larger version (25K): [in a new window] |
FIG. 3. ChIP analysis of the upstream region of the mRankL gene reveals complex 1,25(OH)2D3- and DEX-stimulated transactivator DNA binding and RNA pol II recruitment. ST2 cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), or 1,25(OH)2D3 and DEX (107 M) for 6 h and then subjected to ChIP analysis using antibodies to VDR, GR, RNA pol II, or control IgG ( VDR, GR, pol II, and IgG, respectively). The immunoprecipitated DNA was isolated and then amplified using the primer sets whose positions are illustrated in the top panel and whose sequences are documented in Table 1. Amplification utilized 31 cycles for mRLD1 to mRLD5 and mRLIS1 to mRLIS7 and 34 cycles for the TSS. PCR analyses were all performed within the linear range of amplification. These results are typical of several similar studies.
|
![]() View larger version (21K): [in a new window] |
FIG. 4. The mRLD5 region of the mRankL gene mediates transcriptional induction by 1,25(OH)2D3 and DEX. (A) Features of the cloned mRLD1-to-mRLD5 regions of the mRankL gene, including the sizes of the fragments cloned for evaluation, the distances from the RankL TSS, and the boundaries of the fragments on chromosome 14 (February 2006 assembly). (B) Basal and hormone-inducible activities of mRLD1 to mRLD5 in ST2 cells. ST2 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng) or pTK-mRLD1, pTK-mRLD2, pTK-mRLD3, pTK-mRLD4, or pTK-mRLD5. Cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), or DEX (107 M) plus increasing concentrations of 1,25(OH)2D3 (1010 to 107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities as described in Materials and Methods. Each point represents the normalized relative light unit average ± standard error of the mean for a triplicate set of transfections. These data are representative of three or more similar experiments.
|
![]() View larger version (26K): [in a new window] |
FIG. 5. Hormone-inducible activity of the mRLD5 region is mediated by the VDR and GR. (A) The VDR antagonist ZK159222 blocks activation of mRLD5 by 1,25(OH)2D3. ST2 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng) or pTK-mRLD5. Cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), DEX (107 M) plus 1,25(OH)2D3 (107 M), ZK159222 (106 M), or 1,25(OH)2D3 (107 M) plus ZK159222 (106 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities as described in Materials and Methods. (B) The GR antagonist blocks activation of mRLD5 by DEX. ST2 cells were treated as in panel A above with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), DEX (107 M) plus 1,25(OH)2D3 (107 M), RU486 (106 M), RU486 (106 M) plus 1,25(OH)2D3 (107 M), or RU486 (106 M) plus 1,25(OH)2D3 (107 M) plus DEX (107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities as described in Materials and Methods. (C) 1,25(OH)2D3 activates mRLD5 via endogenous VDR. ST2 cells were transfected with pCH110-ßgal (50 ng) and either the pTK control vector (250 ng) or pTK-mRLD5 and then treated with either vehicle, DEX (107 M), DEX (107 M) plus 1,25(OH)2D3 (1010 to 107 M), or 1,25(OH)2D3 (107 M) alone. Luciferase and ß-galactosidase activities were assessed 24 h later as described in Materials and Methods. (D) The ability of VDR siRNA to block activation of mRLD5-mediated transcription by 1,25(OH)2D3 is rescued through the addition of a functional VDR expression vector. ST2 cells were cotransfected with 10 ng pCH110-ßgal, pTK-mRLD5, pcDNA-hVDR(wtp) (wild type), or pcDNA-hVDR(m) (mutant) and 50 nM of either nontargeted siRNA or mVDR siRNA as described in Materials and Methods. Transfected cells were cultured for 48 h and then treated for an additional 24 h with either vehicle or 1,25(OH)2D3 (107 M) and processed as in panel B. Each point represents the normalized relative light unit average ± standard error of the mean for a triplicate set of transfections. All studies were repeated with similar results.
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Deletion analysis of the mRLD5 region reveals the approximate location of transcriptional response to both 1,25(OH)2D3 and DEX. (A) Schematic of the mRankL deletion fragments of the mRLD5 region that were used to map 1,25(OH)2D3 and DEX responses. The numbering at the bottom represents the distance from the RankL TSS (February 2006 assembly). (B) Transcriptional activities of the mRLD5 subfragments in response to 1,25(OH)2D3 and DEX. ST2 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng) or the pTK-mRLD5 deletion constructs indicated. Cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), or DEX (107 M) plus 1,25(OH)2D3 (107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities as described in Materials and Methods. (C) Transcriptional activities of the mRLD5 subfragments in response to DEX. ST2 cells were transfected with pCH110-ßgal (50 ng), pRSV-hGR (50 ng), and either the pTK control vector (250 ng) or the pTK-mRLD5 deletion constructs indicated. Cells were treated with either vehicle or DEX (107 M) and evaluated after 24 h as in panel B. These results were confirmed via at least three separate experiments.
|
proteins. As can be seen in Fig. 7C, the VDR/RXR heterodimer bound to both the Opn- and the mRL-VDREs in a salt-sensitive, hormone-dependent fashion. Similar binding of endogenous VDR and RXR was also observed in the presence of ST2 nuclear extracts (data not shown). Competition studies suggest that the relative binding affinities of the VDR/RXR heterodimer for these two VDREs were similar (see Fig. S6A in the supplemental material). Interestingly, the results in Fig. 7C also show that incubation of VDR/RXR with the mRL-VDRE but not the Opn-VDRE results in the appearance of a second, more slowly migrating species, suggestive of a specific complex comprised of two VDR/RXR heterodimers. Support for this contention is provided by the observation that the higher-order complex fails to form when purified VDR/RXR is incubated with mRL-VDREs that contain mutations in either VDRE1 or VDRE2 (see Fig. S6B in the supplemental material). Taken together, these results suggest that the element located at 75620 to 75590 upstream of the mouse RankL TSS is capable of both binding two VDR/RXR heterodimers and mediating the transactivation potential of this complex in transfected cells.
![]() View larger version (28K): [in a new window] |
FIG. 7. Mapping the mRL-VDRE. (A) DNA sequence of a putative mRL-VDRE revealed by in silico analysis (http://mordor.cgb.ki.se/cgi-bin/CONSITE/consite). The nucleotide numbering represents the boundaries of the mRL-VDRE relative to the RankL TSS (February 2006 assembly). Nucleotide bases above the arrows represent the triplet alterations introduced by site-directed mutagenesis into each of the half-sites (HS1 to HS4) of the two VDREs (VDRE1 and VDRE2) that constitute the mRL-VDRE. Mutants are designated Mu1 to Mu4. (B) Transcriptional activity of wild-type pTK-mRLD5 or pTK-mRLD5 containing mutation 1 (Mu1), Mu2, Mu3, or Mu4. ST2 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng), pTK-mRLD5, or pTK-mRLD5 containing Mu1 to Mu4 as illustrated. Cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), or the combination and evaluated after 24 h for both luciferase and ß-galactosidase activities. These results were repeated with similar findings. (C) Purified VDR and RXR bind to the mRL-VDRE as a pair of heterodimers. Labeled duplex DNA probes from the mouse osteopontin VDRE (mOpn DR3) or the sequence corresponding to the mRL-VDRE depicted in panel A above were incubated with the indicated amounts of purified VDR and RXR in 50 mM or 150 mM KCl without () or with (+) 1,25(OH)2D3 (5 x 109 M) as indicated. Complexes were resolved on 6% nondenaturing polyacrylamide gels, dried, and visualized by autoradiography. Complexes 1 and 2 are indicated. These results were confirmed through at least three separate EMSA analyses.
|
![]() View larger version (16K): [in a new window] |
FIG. 8. The mRL-VDRE confers significant VDR-dependent, 1,25(OH)2D3 response to the heterologous TK promoter. (A) Transcriptional activity of the mRL-VDRE or its two-component VDREs in ST2 cells. ST2 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng), pTK-mRL-VDRE, pTK-mRL-VDRE1, or pTK-mRL-VDRE2. Cells were treated with either vehicle or increasing concentrations of 1,25(OH)2D3 (1010 to 107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities. (B) Effect of half-site mutations within the mRL-VDRE on response to 1,25(OH)2D3. Duplex mRL-VDRE oligonucleotides containing either the wild-type sequence or individual triplet mutations as depicted in Fig. 7A were cloned into the TK vector and transfected into ST2 cells, and their activity was evaluated in response to 1,25(OH)2D3 (107 M) after 24 h. Due to baseline differences, induction (fold) is reported. (C) Ability of mVDR siRNA to block activation of mRL-VDRE-mediated transcription by 1,25(OH)2D3 is rescued through the addition of a functional VDR expression vector. ST2 cells were cotransfected with 10 ng pCH110-ßgal, pTK-mRL-VDRE, pcDNA-hVDR(wtp) (wild type), or pcDNA hVDR(m) (mutant) and 50 nM either nontargeted siRNA or mVDR siRNA as described in the legend to Fig. 5. Transfected cells were cultured for 48 h and then treated for an additional 24 h with either vehicle or 1,25(OH)2D3 (107 M) and processed as in panel B. Each point represents the normalized relative light unit average ± standard error of the mean for a triplicate set of transfections.
|
Activation of the five regulatory regions of the mouse RankL gene results in histone acetylation. While mRLD4 to mRLD1 failed to exhibit inducible transcriptional activity in the context of either the TK promoter or a homologous mRL promoter, the results of the ChIP/chip and direct ChIP assays using VDR, RXR, and RNA pol II strongly support a role for mRLD4 to mRLD1 in mediating the actions of 1,25(OH)2D3 on RankL gene expression. We therefore asked whether 1,25(OH)2D3-induced binding of VDR and RXR to the five RankL gene regions was capable of inducing modifications to histones within those regions of the gene. Previous studies have suggested that this modification is essential to the upregulation of Cyp24a1 (30). ST2 cells were therefore treated with 1,25(OH)2D3, and the cells were subjected to ChIP analysis for increasing times using antibodies to tetra-acetylated histone 4 as well as VDR. 1,25(OH)2D3 induced both a time-dependent accumulation of VDR at the mRLD1-D5 regions and an increase in histone 4 acetylation, as can be seen in Fig. 9 (and see Fig. S9 in the supplemental material). Surprisingly, however, increased H4 acetylation did not occur at mRLD1 to mRLD5 in the RL gene locus (with the possible exception of mRLD3), but rather at sites located between the five enhancers and at the RankL TSS. This finding suggests that while 1,25(OH)2D3 can induce changes in acetylation within the mRankL gene locus, these changes do not appear to correlate directly with regions of VDR/RXR heterodimer binding.
![]() View larger version (20K): [in a new window] |
FIG. 9. 1,25(OH)2D3 induces histone acetylation in the upstream regions of the mRankL gene. ST2 cells were treated with 1,25(OH)2D3 (107 M) for periods of up to 6 h and then subjected to ChIP analysis using antibodies to VDR, tetra-acetylated histone 4, or IgG ( VDR, AcH4, and IgG, respectively). Precipitated DNA was isolated and evaluated by PCR using Cyp24a1 primers (see Materials and Methods) or RankL primers whose locations are depicted in Fig. 3 (upper panel) and whose sequences are delineated in Table 1. Amplifications were carried out as in Fig. 3. Similar results were obtained in at least four separate experiments.
|
![]() View larger version (30K): [in a new window] |
FIG. 10. The highly evolutionarily conserved hRLD5 region mediates the transcriptional activity of 1,25(OH)2D3 in the human RANKL gene. (A) RANKL mRNA is induced by 1,25(OH)2D3 and enhanced by DEX in MG63 cells. MG63 cells were treated for periods up to 24 h with either vehicle, DEX (107 M), 1,25(OH)2D3 (107 M), or both (at 107 M). Total RNA was isolated and subjected to reverse transcription-PCR (RT-PCR) analysis using primers specific to human CYP24A1 (27 cycles), RANKL (30 cycles), or ß-actin (20 cycles) as documented in Materials and Methods. The results are typical of multiple similar experiments. (B) 1,25(OH)2D3 induces VDR binding to the hRLD5 region of the human RANKL gene. MG63 cells were treated with either vehicle, DEX (107 M), 1,25(OH)2D3 (107 M), or 1,25(OH)2D3 and DEX (107 M) and then subjected to ChIP analysis using antibodies to VDR or control IgG ( VDR and IgG, respectively). The immunoprecipitated DNA was isolated and then amplified using the primer sets documented in Table 1. PCR was performed for 31 cycles. Two separate primer sets were used for validation. (C) The hRLD5 region of the hRANKL gene mediates transcriptional induction by 1,25(OH)2D3 and DEX. MG63 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng), pTK-mRLD5, or pTK-hRLD5. Cells were treated with either vehicle, 1,25(OH)2D3 (107 M), DEX (107 M), or DEX (107 M) plus 1,25(OH)2D3 (107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities. Each point represents the normalized relative light unit (RLU) average ± standard error of the mean for a triplicate set of transfections. (D) Transcriptional activities of wild-type pTK-hRLD5 or mutant pTK-hRLD5 containing triplet changes in the hRL-VDRE. MG63 cells were transfected with pCH110-ßgal (50 ng), pcDNA-hVDR (50 ng), and either the pTK control vector (250 ng), pTK-mRLD5, pTK-hRLD5, or pTK-hRLD5 containing triplet mutations in half-site 2 (Mu2) or 4 (Mu4). Cells were treated with either vehicle or increasing concentrations of 1,25(OH)2D3 (109 to 107 M) and evaluated after 24 h for both luciferase and ß-galactosidase activities. Each point represents the normalized RLU average ± standard error of the mean for a triplicate set of transfections. (E) VDR and RXR bind directly to the hRL-VDRE in a salt- and hormone-dependent fashion. Labeled duplex DNA probes comprising the mRL-VDRE or the hRL-VDRE were incubated with the indicated amounts of purified VDR and RXR in 50 mM () or 150 mM (+) KCl without or with 1,25(OH)2D3 (109 M) as indicated. Complexes were resolved on 6% nondenaturing polyacrylamide gels, dried, and visualized using autoradiography. The results of these separate studies were reproduced at least three times.
|
|
|
|---|
Our studies using ChIP/chip analysis, direct ChIP, and more traditional molecular approaches identified five distinct regions located at significant distances upstream of the TSS that were responsible for regulating the expression of the RankL gene. VDR and RXR localize to each of these regions in response to 1,25(OH)2D3 and appear to be accompanied by other factors as well, including GR. This supports the idea that each of these regions may represent enhancer modules integral to RankL gene expression. A typical consequence of transactivation at enhancers is the rapid recruitment of coregulatory enzymes and the subsequent modification of local chromatin structure, in part through acetylation (5, 52, 55, 56). It is believed that these epigenetic modifications function to promote the chromatin decondensation that is often necessary for increased transcription. In our studies, while it is clear that 1,25(OH)2D3 induced histone 4 acetylation at the RankL gene locus, these modifications did not occur at sites of VDR/RXR binding, but rather at sites located in between. Thus, although direct recruitment of histone-modifying enzymes by the VDR/RXR heterodimer may be involved, we believe it more likely that the VDR precipitates events at these sites that lead in turn to the increased acetylation seen between the enhancer regions. The nature of these events remains to be determined. Interestingly, 1,25(OH)2D3 also induced the recruitment of RNA pol II to the mRLD1-to-mRLD5 regions of the RankL gene. This finding suggests that the mRLD1-to-mRLD5 enhancers may also act as recruitment centers for components of the transcriptional apparatus as well (59). Although speculative, it is possible that the recruitment of transcriptional components to the five enhancer regions and the consequences of this recruitment may represent the initiating event for the histone modifications seen above. Regardless, our data suggest that the ChIP/chip and ChIP approaches taken here were strategic in identifying key regions within the RankL locus that are responsible for transcriptional regulation of this gene by 1,25(OH)2D3 and DEX.
While all of the sites appear to contribute to the response initiated by 1,25(OH)2D3, only the mRLD5 region was capable of conferring sensitivity to 1,25(OH)2D3 in a transcriptional reporter plasmid assay. The respective responses to 1,25(OH)2D3, DEX, and the combination are consistent with that seen at the level of RankL mRNA, thereby validating further the relevance of mRLD5 to RankL gene expression. The roles of VDR and RXR in this regulation were further established using both selective antagonists as well as siRNA, which reduced VDR mRNA levels and compromised the capacity of 1,25(OH)2D3 to induce transcription via the mRLD5 segment. Because ZK159222 and RU486 do not prevent receptor DNA binding, but rather alter their ability to recruit coregulators (16, 30), these results provide additional support for cofactor involvement in enhanced RankL gene expression. We defined the boundaries of mRLD5 based upon the VDR and RXR binding activities seen in the ChIP/chip analysis. The fragment was comprised of approximately 1,100 bp of highly conserved sequence present across multiple RankL genes. An additional highly conserved region of 800 bp immediately upstream of the mRLD5 region was also noted. This segment corresponds directly to the region that mediates PTH response, as described in the accompanying article by Fu, Manolagas, and O'Brien (13). Our additional studies using ChIP analysis, transient transcription assays, and mutagenesis fully confirm these observations (S. Kim and J. W. Pike, unpublished observations). Thus, we propose that this entire conserved region be designated the RankL distal control region, or RL-DCR.
The transcriptional activity of the mRLD5 region allowed us both to map the mRL-VDRE and to characterize a specific GC response. The mRL-VDRE retains a unique structure comprised of two separate VDRE sequences linked via a single base pair. This element is clearly capable of the simultaneous binding of two VDR/RXR heterodimers and is likewise responsible for the 1,25(OH)2D3 sensitivity observed within the mRLD5 region. Interestingly, the half-elements alone bind single VDR/RXR heterodimers, but do not elicit a significant transcriptional response when cloned and analyzed independently. We did not map the existing GREs within this region of the RankL gene, however. The absence of activity in two contiguous but nonoverlapping fragments and the DEX sensitivity observed in an mRLD5 fragment which spanned the boundary between these two fragments suggest the presence of several GREs that may function synergistically to mediate DEX activation. Further studies will be necessary to define the putative GREs that we observed using in silico analysis (data not shown). Although DEX activity on mRLD5 appears to be synergistic, the addition of exogenous GR enabled us to detect a direct response to the GC. This effect was also noted when the human RLD5 region (hRLD5) was evaluated in the context of MG63 cells. Regardless of the nature of these overall effects, it seems clear that 1,25(OH)2D3 action is mediated via a single VDRE, whereas DEX activity is likely mediated by several independent GREs.
The appearance of VDR and RXR, histone acetylation, and RNA pol II recruitment all support the idea that mRLD1 to mRLD5 are active in the regulation of the RankL gene by 1,25(OH)2D3. It is therefore curious as to why the mRLD1-to-mRLD4 upstream regions of the gene are incapable of mediating a transcriptional response to 1,25(OH)2D3 in transfection studies. We believe that this resistance to hormonal induction highlights the importance of context, wherein the presence of both positive as well as negative cis elements and their respective transregulators exert significant influence on transcriptional readouts obtained during transient transfection assays. Interestingly, even the activity of a natural promoter in a plasmid context may be deceptive. The large distances that exist between mRLD1 to mRLD5 and the location of the most distance element (76 kb) made it difficult for us to assess the transient activity of a "full-length" RankL promoter by these means. Fu, Manolagas, and O'Brien overcame this problem, however, by using recombineering methods to produce large, bacterial artificial chromosome (BAC) clone-derived RankL constructs which could be stably transfected into host cells (13). Interestingly, despite the differences in the two approaches, both of our groups were able to identify the same highly conserved distal region within the RankL gene that displayed sensitivity to the two hormones 1,25(OH)2D3 and PTH. While the above results clearly demonstrate the overall dominance of the RL-DCR, the ability of 1,25(OH)2D3 to induce a large RankL DNA fragment that no longer contained the mRLD5 (13) provides further evidence that the mRLD1-to-mRLD4 regions are also important contributors to 1,25(OH)2D3 response. Therefore, we believe that ChIP analysis may represent a more reliable method of assessing the presence of functional enhancers than the alternative analysis using plasmid transfection. Studies are ongoing in our laboratory, however, to identify the precise locations of the VDREs lodged within the mRLD1-to-mRLD4 regions that regulate RankL gene expression.
Regardless of the properties of mRLD1 to mRLD4, the studies of both Fu, Manolagas, and O'Brien (13) as well as our own highlight the important role of the mRL-DCR in mediating both PTH and 1,25(OH)2D3 responses. It is therefore noteworthy that in the accompanying article, Fu et al. (13) provide in vivo evidence that genetic deletion of the mRL-DCR in mice leads to a loss of PTH-mediated RankL induction. We would predict that these mice will also display a similar resistance to 1,25(OH)2D3, although whether this resistance will be only partial or complete remains to be determined. Collaborative studies are currently under way to more fully investigate the properties of this altered RankL gene locus in vivo and, more importantly, perhaps, to determine the nature of the skeletal phenotype that results from this genetic change.
The regulatory regions for RankL are widely dispersed across a rather large segment of upstream DNA, the furthest some 76 kb from the TSS. It is therefore clear that an understanding of the function of each of these regions and their individual roles in modulating the expression of the RankL gene will require significant additional work. In that vein, there is increasing evidence that many genes contain distant regulatory elements such that this mechanism of regulation may be more frequent than previously believed (2, 21, 39). The extended pattern and distance of these elements within the RankL gene as well as in other genes almost certainly highlight the crucial impact of chromatin structure and organization on the expression of these genes such that the elements can directly impact activity at the TSS. One might speculate that these regions converge directly on the proximal RankL promoter by virtue of extensive chromatin looping (6, 9, 49, 65). New technologies are currently available or in development to test distance relationships between regulatory regions and their functional promoters, thus making it possible to explore this intriguing hypothesis (7).
A final question that emerges from these studies is why RNA pol II might be recruited to each of the five enhancer domains within the RankL gene. One might imagine that the presence of RNA pol II at these sites might also require the simultaneous recruitment of basal transcription factors such as TF-IIA, TF-IIB and TAF-II250 as well (59). Indeed, we have seen in preliminary studies that 1,25(OH)2D3 can induce the recruitment of TAF-II250 to the mRLD5 region (S. Kim and J. Pike, unpublished data). One frequent consequence of the assembly of such factors at upstream regulatory regions, however, is the production of nascent noncoding mRNAs (3, 4, 10, 26, 27). The role of this transcription and these transcripts is currently unknown, although it has been suggested that their production may be essential to the maintenance of an open chromatin state necessary for gene regulation (15). An alternative proposal is that upstream enhancers may function as centers for the recruitment of basal transcriptional machinery, thereby providing a source of such factors for authentic promoter-driven transcription (59). Now that some of the fundamentals of hormonal regulation have been established within the RankL gene locus, future studies will focus on these issues involving RNA pol II recruitment and on understanding the spatial arrangement that likely exists between the five regulatory regions and the TSS.
In summary, we have shown that 1,25(OH)2D3 and its receptor induce the expression of RankL via five regulatory domains located at significant distances upstream of the TSS. This regulation is facilitated by DEX and the GR, the latter appearing to localize to several but not all of these regions. Transcription factor binding within these regions is associated with the recruitment of cofactors such as RNA pol II and adjacent histone acetylation. Mapping studies of mRLD5, perhaps the dominant control region for the RankL gene, led to the delineation of an unusual VDRE. Support for mRLD5 as a primary enhancer of mouse RankL gene expression was increased by the discovery that an analogous component is located within the human RANKL gene.
This work was supported by National Institutes of Health grant DK-74993 (to J.W.P.).
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
|
|
|---|
B ligand gene expression via a distant transcriptional enhancer. Mol. Cell. Biol. 26:6453-6468.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»