Lars Grøntved,1,
Hendrik G. Stunnenberg,2 and
Susanne Mandrup1*
Department of Biochemistry and Molecular Biology, University of Southern Denmark, 5230 Odense M, Denmark,1 Department of Molecular Biology, Nijmegen Center for Molecular Life Sciences, Radboud University, 6500 HB Nijmegen, The Netherlands2
Received 25 November 2005/ Returned for modification 2 January 2006/ Accepted 16 May 2006
| ABSTRACT |
|---|
|
|
|---|
2 possesses considerable ligand-dependent as well as independent transactivation potential and that agonists increase the occupancy of PPAR
2/retinoid X receptor at PPAR response elements. Intriguingly, by direct comparison of the PPARs (
,
, and ß/
), we show that the subtypes have very different abilities to gain access to target sites and that in general the genomic occupancy correlates with the ability to activate the corresponding target gene. In addition, the specificity and potency of activation by PPAR subtypes are highly dependent on the cell type. Thus, PPAR subtype-specific activation of genomic target genes involves an intricate interplay between the properties of the subtype- and cell-type-specific settings at the individual target loci. | INTRODUCTION |
|---|
|
|
|---|
, -
, and -ß/
, bind as heterodimers with the retinoid X receptor (RXR) to conserved PPAR response elements (PPREs), i.e., a direct repeat of the hexameric sequence (A/G)GGTCA spaced by one nucleotide (18). The different PPAR subtypes display highly diverse biological functions in vivo. Thus, PPAR
(30) and PPARß/
(11, 67) activate primarily genes encoding enzymes involved in lipid oxidation (catabolic pathways), whereas PPAR
activation leads to induction of genes involved in lipogenesis (anabolic pathways) (44). PPARs are activated by a variety of fatty acids and fatty acid derivatives such as prostaglandins and leukotrienes. In addition, members of the PPAR family are important drug targets for the treatment of insulin resistance, type II diabetes, dyslipidemia, hypertension, atherosclerosis, and some cancers (8, 20).
In spite of the pharmacological and biological importance, the mechanisms by which PPARs induce endogenous target genes in their natural chromatin context has not been thoroughly examined. Approaches such as stable overexpression of the PPAR subtypes or administration of subtype-selective agonists to activate endogenous PPARs in tissues and cell lines have been used to identify novel PPAR target genes (1, 19, 67, 70). However, while such experiments have contributed significantly to the identification of target genes, they are not useful for dissecting the molecular mechanisms by which PPAR subtypes activate endogenous target genes. This is due to the fact that PPARs, as a result of endogenous ligands and/or ligand-independent transactivation, activate target genes in the absence of exogenous ligands. This has been demonstrated in transient transfections (25) as well as in retroviral transductions (7). The relatively high basal expression of endogenous target genes precludes experiments analogous to the ones performed with the steroid receptors, which are kept inactive either by nuclear exclusion and/or by interaction with heat shock protein complexes (28, 38, 56, 57). Thus, in contrast to steroid receptors, the addition of exogenous agonists cannot be used to define clean on/off states of the PPARs with respect to the establishment of a transcriptionally active complex at the promoters.
The different PPAR subtypes show limited specificity in their binding to different PPREs in in vitro mobility shift assays. A preference for PPAR
/RXR heterodimers has been shown for weak PPREs, but for most PPREs no clear subtype specificity has been demonstrated by such assays (27). Similarly, the different PPAR subtypes differ only little in their ability to activate artificial promoter reporter constructs in transient transfections (17, 24, 29, 55). Intriguingly, the limited difference between PPAR subtypes in electrophoretic mobility shift assays and in transient transfections is in sharp contrast to results obtained by in vivo activation of endogenous PPARs using subtype-selective agonists or by experiments with knockout models. The differential biological actions of PPAR subtypes are undoubtedly due in part to the differential expression patterns of the different PPAR subtypes. Thus, PPAR
is mainly expressed in tissues with high fatty acid ß-oxidation rates, such as liver, heart, skeletal muscle, and brown fat, while PPAR
is abundantly expressed in brown and white adipose tissues and macrophages (6). However, adenoviral expression of PPAR
1 in mouse liver leads to induction of several genes which are not readily activated by PPAR
, i.e., genes involved in lipid accumulation and adipogenesis (73). Importantly, results from the few comparative analyses have shown that the PPARs are also biochemically distinct and behave differently when ectopically expressed in the same cell. Thus, PPAR
is by far the most adipogenic of the PPAR subtypes in fibroblasts (7). In addition, we have recently shown that PPAR
and PPAR
induce distinct subsets of genes in pancreatic ß-cells. In these cells ectopic expression of PPAR
results in increased fatty acid uptake and fatty acid oxidation and in promotion of glucose-stimulated insulin secretion, whereas ectopic expression of PPAR
results in increased fatty acid uptake and increased lipid accumulation and inhibition of glucose-stimulated insulin secretion (49). The molecular mechanisms underlying this subtype-specific activation of gene expression in a particular cell are not understood.
Here we show that acute ectopic expression of the PPARs by adenoviral vectors is an ideal approach to investigate PPAR subtype-specific activation of target genes in individual cell types. This system allows us to tightly control the timing and expression level of the PPARs and thereby go from a clean off to a clean on state of PPAR target loci. By combining adenoviral delivery and chromatin immunoprecipitation (ChIP), we investigate the rapid establishment of an active transcriptional complex on PPAR target genes immediately following the ectopic expression of PPARs. We show that the PPARs retain significant subtype specificity in this system and induce target genes in a highly subtype- and cell-type-specific manner. This PPAR subtype-specific activation was shown to be correlated with PPAR occupancy at the given promoters, indicating that PPAR binding to chromatin-embedded target sites is a major determinant of PPAR subtype-specific activation of endogenous target genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Mouse embryonic fibroblasts (MEFs) were isolated from embryonic day 13.5 C57Bl/6J-Bom embryos. The embryos were transferred to cold sterile phosphate-buffered saline (PBS); placenta, membranes, heads, and inner organs were removed, and the rest of the bodies were dissected into small pieces. After centrifugation (1,200 rpm for 5 min at room temperature), the supernatant was removed, and the pellet was resuspended in 25% collagenase (Sigma), 1x PBS, and 20% FBS. Following a 1-h digestion at 37°C, the cells were harvested by centrifugation and resuspended in MEF culture medium (AmnioMax C-100 basal medium [Gibco] supplemented with 7.5% AmnioMax C-100 supplement [Gibco], 7.5% FBS, 1% L-glutamine, 62.5 µg/ml penicillin, and 100 µg/ml streptomycin). Cells were counted and plated on Falcon plates at 37°C with 95% humidity and 7.5% CO2. Medium was changed every day.
Retroviral transduction of NIH 3T3 fibroblasts.
Phoenix cells were transfected using the calcium-phosphate technique with a retroviral LXSN-hCAR
1 vector expressing the truncated coxsackie-adenovirus receptor (CAR
1) (47) at 50% confluence. Two days after transfection, virus supernatant was harvested and centrifuged to remove living phoenix cells. NIH-3T3 cells at 50% confluence were transduced with a 1:1 dilution of virus supernatant and fresh growth medium in the presence of 6 µg/ml Polybrene (Sigma-Aldrich). The following day, the cells were split and subjected to neomycin (G418; 0.7 mg/ml; Bie and Berntsen) selection. Approximately 5 days later, the selected NIH 3T3 (CAR
1) cells were pooled and replated for single colony isolation. Single clones were isolated by replating 80% confluent cells at a 103-fold dilution and culturing them for 3 to 5 days before trypsinization in sterile O-rings. The level of CAR
1 expression in the individual clones was evaluated by fluorescence-activated cell sorter (FACS) analysis. A high CAR
1-expressing clone, termed NIH-CAR, was chosen for further studies.
FACS analysis.
For CAR FACS analysis, 105 cells [NIH 3T3 and NIH 3T3(CAR
1)] were trypsinized, washed in PBS, and incubated for 30 min with 100 µl of PBS containing 0.5% bovine serum albumin (BSA) and CAR antibody (05-644; Upstate) at 1:200 at 4°C. The cells were washed twice in PBS and incubated for 30 min with 150 µl of PBS including 0.5% BSA and fluorescein isothiocyanate-conjugated anti-mouse antibody (12-487; Upstate) at 1:150 at 4°C. The cells were washed with PBS, fixed in 100 µl of PBS containing 1% paraformaldehyde, and analyzed by FACS.
For green fluorescence protein (GFP) FACS analysis, 106 cells [NIH 3T3, NIH 3T3(CAR
1), and NIH-CAR] were transduced with adenoviral GFP (47) (AdGFP; 55 PFU/cell). Cells were washed twice with PBS, trypsinized, fixed in 100 µl of PBS containing 1% paraformaldehyde, and analyzed by FACS.
Adenovirus generation and purification.
Recombinant adenoviruses containing full-length hemagglutinin (HA)-tagged mouse PPAR
, PPAR
2, and PPARß/
were generated using the AdEasy cloning system (Stratagene). The three PPAR subtypes were initially N-terminally HA-tagged and cloned into pShuttle-CMV (Stratagene) (where CMV is cytomegalovirus) using the NotI/XhoI (PPAR
) and KpnI/EcoRV (PPAR
2 and PPARß/
) restriction sites. The CMV-HA-PPAR cassettes were transferred to the AdEasy-1 vector by homologous recombination in Escherichia coli to generate AdHA-PPAR constructs. Plasmids were linearized and transfected into 293-HEK cells using the calcium-phosphate procedure. The adenoviruses were amplified and purified using CsCl gradients, and viral titers were estimated by a plaque assay-based approach as recommended by the AdEasy protocol (Stratagene).
Adenoviral transduction of cell lines.
For adenoviral transduction fresh medium containing AdHA-PPARs (55 PFU/cell) was given to NIH-CAR, MIN6, or AML-12 cells at 80% confluence. After 2 h of viral transduction the virus-containing medium was removed, and new medium containing the following PPAR-specific agonists was added for an additional 2 to 10 h: 1 µM rosiglitazone/BRL49653 (Novo Nordisk), 30 µM Wy14.643 (Calbiochem), 1 µM L-165041 (Merck), and/or PPAR
antagonist 1 µM GW9662 (Sigma-Aldrich).
RNA purification and cDNA synthesis. Adenovirally transduced NIH-CAR, MIN6, and AML-12 cells were harvested in guanidium thiocyanate at the indicated time points (see figures), and RNA was isolated according to a modified Chomczynski-Sacchi protocol (12). RNA preparations were subjected to DNase I (Invitrogen) treatment, and cDNA was synthesized using random deoxynucleic acid hexamers and reverse transcriptase (First-Strand Kit; Invitrogen) as previously described (23). All experiments were performed in duplicates.
Protein analysis.
Adenovirally transduced NIH-CAR, MIN6, and AML-12 cells were harvested in a hypotonic sodium dodecyl sulfate (SDS) sample buffer as described previously (23) and separated by SDS-polyacrylamide gel electrophoresis (PAGE). Proteins were blotted onto polyvinylidene difluoride membranes (Millipore Corp.) and probed with specific antibodies. The following primary antibodies were used: anti-PPAR
(sc-7273; Santa Cruz), anti-TFIIB (sc-225; Santa Cruz) and anti-HA (12CA5 [see reference 40]). The secondary antibodies used were horseradish peroxidase-conjugated goat anti-mouse immunoglobulin G (P0447; DAKO) and swine anti-rabbit immunoglobulin G (P0339; DAKO). Enhanced chemiluminescence (ECL; Amersham Phamacia Biotech) was used for detection.
Transient transfection of AdPPAR-transduced NIH-CAR cells. NIH-CAR cells were transfected at 50 to 60% confluence with a PPRE3-TK-Luc (where TK is thymidine kinase and Luc is luciferase) reporter construct (34) and a pCMV-ß-galactosidase control (Promega) in six-well plates using the calcium-phosphate procedure (10). Following 6 h of incubation with DNA, the medium was changed to DMEM supplemented with 10% resin-charcoal-stripped calf serum containing AdHA-PPARs (55 PFU/cell) for 2 h together with the PPAR-specific ligands 1 µM rosiglitazone/BRL49653 (Novo Nordisk), 30 µM Wy14.643 (Calbiochem), and 1 µM L-165041 (Merck). The transduced cells were subsequently incubated additionally for 10 h in the presence of PPAR-specific ligands before the cells were harvested in a lysis solution (Tropix). Luciferase and ß-galactosidase assays were performed as described previously (24). All experiments were performed in triplicate, and luciferase and ß-galactosidase measurements were done in duplicate.
ChIP.
NIH-CAR (106 to 107) cells were cross-linked in 1% formaldehyde for 20 min. Cross-linking was stopped by the addition of glycine to a final concentration of 0.125 M. Cross-linked cells were removed from the plate and washed once in cold PBS and once in cross-link buffer 1 (0.25% Triton, 10 mM EDTA, 0.5 mM EGTA, 10 mM Tris-HCl, pH 8.0) at 4°C. Finally, the cross-linked cells were washed in cold cross-link buffer 2 (0.2 M NaCl, 1 mM EDTA, 0.5 mM EGTA, 10 mM Tris-HCl, pH 8.0). Cell extract was diluted in 2 ml of chromatin dilution buffer (0.05% SDS, 1 mM EDTA, 0.5 mM EGTA, 50 mM HEPES, pH 8.0, and 1xcomplete protease inhibitor cocktail [Boehringer Mannheim]) and sonicated according to the manufacturers protocol using the Diagenode Bioruptor. The amount of DNA in each sample was determined by the absorbance at 260 nm and diluted to equal amounts. Protein A (Amersham Pharmacia Biotech), protein G (Amersham Pharmacia Biotech), or protein A/G (Santa Cruz) beads were prepared by three washes in immunoprecipitation (IP) dilution buffer (1% Triton, 0.1% SDS, 150 mM NaCl, 2 mM EDTA, 1 mM EGTA, 20 mM Tris, pH 8.0) and subsequently incubated in IP dilution buffer containing 1 µg/µl BSA and 1x complete protease inhibitor mixture. Chromatin was diluted in IP dilution buffer (2% triton, 0.15% SDS, 300 mM NaCl, 4 mM EDTA, 2 mM EGTA, 40 mM Tris, pH 8.0, 2 µg/µl BSA, and 1x complete protease inhibitors) and incubated with prepared protein A, G, or A/G beads at 4°C for 1 h. Beads were removed and supernatant was incubated with antibody over night at 4°C on a rotating wheel, followed by 3 h incubation with protein A, G or A/G beads at 4°C, rotating. Beads were washed at 4°C once with IP wash buffer 1 (1% triton, 0.1% SDS, 150 mM NaCl, 2 mM EDTA, 1 mM EGTA, 20 mM Tris, pH 8.0), twice with IP wash buffer 2 (1% triton, 0.1% SDS, 500 mM NaCl, 2 mM EDTA, 1 mM EGTA, 20 mM Tris, pH 8.0), once with IP wash buffer 3 (0.25 M LiCl, 1% NP-40, 1% deoxycholic acid, 1 mM EDTA, 0.5 mM EGTA, 10 mM Tris, pH 8.0), and finally once with Tris-EDTA buffer, all at 4°C. DNA-protein complexes were eluted with 200 µl of elution buffer (1% SDS and 0.1 M NaHCO3) and decross-linked by adding 0.2 M NaCl and shaking overnight at 65°C. DNA was purified by a PCR purification kit (Amersham Pharmacia Biotech) and dissolved in 50 µl of water. Immunoprecipitated DNA and 1% input DNA were analyzed by real-time PCR. Antibodies used for IP were anti-PPAR
(sc-7273; Santa Cruz), anti-RXR
(sc-774; Santa Cruz), anti-CREB binding protein (anti-CBP; sc-369; Santa Cruz), anti-thyroid hormone receptor-associated protein 220; (anti-TRAP220) (sc-8998; Santa Cruz), anti-TAT a box binding protein (anti-TBP) (SL30 [previously described in reference 52]), anti-RNA polymerase II (anti-Pol II) (8WG16; Covance Research Products), anti-H3-Ac (06-599; Upstate Technologies), anti-H4-Ac (06-866; Upstate Technologies), and anti-HA (12CA5 [previously described in reference 40]).
Real-time PCR. Quantitative three-step real-time PCR was performed on the ABI-7700 PRISM real-time PCR instrument (Applied Biosystems) using 2x SYBR Green Master Mix and Sigma passive reference (Sigma-Aldrich) according to the instructions from the manufacturer. All measurements were performed in duplicate. Primers for real-time PCR (URL: http://www.sdu.dk/Nat/bmb/faculty/pubs/MCB06_supp_data.html) were designed using Primer Express 2.0 (Applied Biosystems), and specificity and efficacy were validated before use.
| RESULTS |
|---|
|
|
|---|
We therefore turned to the adenoviral vector system, which mediates fast and efficient ectopic expression of proteins. However, adenoviral infection of this fibroblastic cell line is known to be very inefficient due to the low expression of the CAR (47). To overcome this problem, we obtained a retroviral vector expressing a truncated version of the CAR protein, CAR
1, which lacks the cytoplasmic tail of CAR. Orlicky and coworkers have previously described and used this elegant system for rendering fibroblasts susceptible to adenoviral infection (47). Using this retroviral vector for expression of CAR
1, the amount of CAR
1 receptors at the cell surface, as determined by FACS analysis, was increased approximately 100-fold compared to parental NIH 3T3 cells (Fig. 1A). The expression of the CAR
1 protein greatly facilitated efficient adenoviral-mediated entry of transgenes into NIH 3T3 (CAR
1) as exemplified by infection with the AdGFP vector (Fig. 1B). To ensure as fast and synchronous uptake of the transgene as possible, we isolated single clones from the pool of CAR
1-transduced NIH 3T3 cells. These clones were screened for their expression of CAR
1 by FACS analysis, and one clone with a high level of CAR
1 expression, termed NIH-CAR, was selected for further use. As illustrated in Fig. 1C, this high-expressing clone (NIH-CAR) exhibited almost 100% efficiency in AdGFP uptake and expression after only 2 h of infection, whereas the mixed NIH 3T3 (CAR
1) population had only 50% efficiency at this early time point.
|
, PPAR
2, and PPARß/
were cloned into adenoviral vectors as previously described (49), except that an HA tag was added to the N terminus of all PPAR subtypes. This tag allowed us to directly compare the expression level as well as recruitment to target sites of all PPAR subtypes by Western blotting and ChIP, respectively.
Initially, the infection ratio of the AdHA-PPAR
2 was adjusted so that the expression of PPAR
after 8 h of infection was equivalent to the expression of PPAR
in differentiated 3T3-L1 adipocytes as determined by Western blotting and probing with a PPAR
antibody (Fig. 2A). Subsequently, using an antibody directed against the HA epitope, the infection ratios of AdHA-PPAR
and AdHA-PPARß/
were adjusted to give protein expression equivalent to that of HA-PPAR
2. Using these infection ratios, elevated PPAR mRNA level could be detected 2 h after adenoviral transduction (Fig. 2B), and prominent PPAR protein expression was observed 4 h after adenoviral transduction (Fig. 2C). Protein and mRNA expression increased between 4 to 8 h and decreased slightly between 8 and 10 h (Fig. 2B and C).
|
was the most potent inducer of transcription, all HA-PPARs were transcriptionally active. Thus, the combination of the NIH-CAR cells and the AdHA-PPAR vectors provides a very useful tool for acute and synchronous ectopic expression of equal amounts of the HA-PPAR subtypes in a background with very low endogenous levels of PPAR. This makes it possible to directly investigate and compare PPAR subtype-specific activation of endogenous PPAR target genes.
Kinetics of AdHA-PPAR activation of target gene expression in NIH-CAR cells.
To investigate the kinetics of PPAR-mediated activation of endogenous target genes in our model system, we chose to focus on PPAR
2-mediated induction of the adipocyte fatty acid binding protein (A-FABP, also called aP2) gene, which is a well-characterized PPAR
target gene (63). Adenoviral expression of PPAR
2 in the presence of the PPAR
-specific agonist rosiglitazone/BRL49653 led to a significant induction of the A-FABP transcript from the endogenous gene after 4 h (Fig. 3A). The fast induction of the A-FABP gene by PPAR
2 is reflected in a significant recruitment of PPAR
2 to the A-FABP enhancer PPREs (Fig. 3B) within 4 h after transduction (Fig. 3C). The coactivators CREB-binding protein (CBP) and TRAP220 are recruited to the enhancer PPREs, and the RNA polymerase II is recruited to the A-FABP promoter within 4 to 5 h. Maximal occupancy is obtained at 5 h following transduction. The activation of the A-FABP gene by PPAR
2 and the subsequent recruitment of coactivators are associated with increased histone H3K9/14di- and H4 tetra-acetylation (Fig. 3C). These data show that the establishment of a transcriptionally active complex and induction of the classical PPAR
target gene, A-FABP, are fast-occurring processes immediately following the expression of PPAR
2 in fibroblasts. This indicates that PPAR
2 is sufficient to trigger activation of the A-FABP gene in NIH-CAR cells.
|
2 displays ligand-dependent as well as ligand-independent transactivation potential.
Up to this point, we have expressed the PPARs in the presence of their selective agonists. A fundamental question in PPAR transactivation that has remained unanswered is the extent to which PPARs display ligand-independent transactivation of endogenous target genes. Thus, it is unknown whether, e.g., the high basal expression of PPAR target genes in adipocytes is due to endogenous PPAR
agonists, ligand-independent activity of PPAR
, or the activity of other transcription factors that activate the expression from these promoters once PPAR
has triggered the induction. In transient transfections all PPARs, and in particular PPAR
, show significant activation of a PPRE reporter construct in the absence of exogenous ligand (Fig. 2D) (25). However, in most experiments it is unclear to what extent this is due to the presence of endogenous ligand or to ligand-independent transactivation. We therefore addressed the question of ligand-independent transactivation by the PPARs using adenoviral delivery of PPAR
2 in NIH-CAR cells in the presence or absence of the specific PPAR
agonist rosiglitazone or the PPAR
antagonist GW9662. The latter binds covalently to PPAR
, thereby permanently abolishing binding of agonists to PPAR
(37).
As seen in Fig. 4A, adenoviral expression of PPAR
2 for 8 h in the presence of the PPAR
-specific agonist rosiglitazone leads to a 1,400-fold induction of A-FABP gene expression relative to nontransduced cells. However, in the absence of exogenous agonist, PPAR
2 in itself leads to a significant (200-fold) induction of the A-FABP gene. The majority of this transactivation appears to be due to ligand-independent transactivation, as the PPAR
antagonist GW9662 decreased PPAR
2-induced A-FABP expression only slightly. In contrast, GW9662 totally abolishes rosiglitazone-dependent PPAR
2 transactivation. Induction of other PPAR
target genes like CD36, perilipin, and lipoprotein lipase (LPL) followed a pattern similar to A-FABP (see Fig. 7A; data not shown). These results clearly demonstrate that PPAR
2 possesses significant ligand-independent as well as ligand-dependent transactivation.
|
|
2 under these conditions, cofactors and basal transcription factors were recruited to the promoter. When the agonist is added, recruitment of PPAR
2 and RXR to the PPRE and recruitment of TBP and Pol II to the promoter are further increased, while the effect on CBP recruitment is more modest. Thus, there is a strong correlation between activation of the PPAR target gene and recruitment of PPAR
2/RXR and auxiliary transcriptional cofactors.
Subtype-specific activation of endogenous PPAR target genes.
To compare the transactivation potential of PPAR subtypes in fibroblasts, NIH-CAR cells were transduced with AdHA-PPAR
, -
, or -ß/
virus in the presence of selective PPAR agonists. RNA was isolated for cDNA synthesis 4, 8, and 12 h posttransduction, and the mRNA expression levels of a large set of PPAR target genes chosen from the literature were evaluated by real-time PCR (Table 1 and Fig. 5A). The ability of the different subtypes to activate endogenous target genes was highly subtype dependent. Based on the previously reported induction by the PPAR subtypes, the selected target genes were divided into different subgroups (Table 1). The first group encompasses characteristic PPAR
target genes investigated: A-FABP, perilipin, and aquaporin-7 (AQP7). These genes were highly activated by ectopic PPAR
2, only modestly activated by PPAR
, and not significantly activated by PPAR
in NIH-CAR cells. The genes reported to be activated by all PPAR subtypes followed (except for adipose differentiation-related protein [ADRP]) a very similar pattern as the PPAR
-type genes; i.e., they were most efficiently activated by PPAR
, less by PPAR
, and not activated by PPARß/
. The reported PPAR
/(
/ß) target genes, by contrast, could be divided into three subgroups. Six target genes, i.e., those encoding hydroxymethyl glutaryl coenzyme A synthetase 2 (HMG-CoA S2), the bifunctional enzyme enoyl-CoA hydrates/L-3hydroxyacyl-CoA dehydrogenase (HD), pyruvate dehydrogenase kinase 4, the long-chain acyl-CoA dehydrogenenase, and carnitine palmitoyltransferase II were most effectively induced by PPAR
. Three of these (HMG-CoA S2, pyruvate dehydrogenase kinase 4, and HD) were induced more than fivefold. Surprisingly, the genes encoding medium chain acyl coenzyme A (acyl-CoA) dehydrogenase, peroxisomal biogenesis factor 11a (PEX11a), and acyl-CoA oxidase 1 (ACOx1) were most effectively induced by PPAR
2 rather than by PPAR
or -
. A large group of known PPAR
/(ß/
) target genes, including, malic enzyme (ME), 3-ketoacyl-CoA thiolase B (thiolase B), cytochrome P450 A10, long chain acyl-CoA dehydrogenenase, and liver-FABP (L-FABP), were not induced by any of the PPARs in NIH-CAR cells. Notably, expression of endogenous PPAR genes was not affected by the acute ectopic expression of PPARs (data not shown). This clearly demonstrates that changes in gene expression cannot be ascribed to activation of endogenous PPARs. Furthermore, most target genes were significantly induced by PPAR
2 (and in some cases also PPAR
) already at the 4-h time point, indicating that PPAR
2 is sufficient for activation of a large number of target genes in NIH-CAR cells.
|
|
and PPAR
2 is slightly less pronounced in MEFs compared to NIH-CAR cells, the PPAR subtypes display very similar induction profiles as observed in NIH-CAR cells. The only major difference is that the L-FABP gene can be activated by PPAR
in MEFs, whereas it is insensitive to PPAR transactivation in NIH-CAR cells. Thus, it appears that during the establishment of the NIH 3T3 cell line from MEFs, the L-FABP gene has been more permanently inactivated.
In summary, despite the equal expression levels of the PPAR subtypes and the lower transactivation potential of PPAR
2 compared to PPAR
as determined by transient transfection (Fig. 2D), PPAR
2 is the most efficient subtype in activating the majority of endogenous chromatin-embedded PPAR target genes in NIH-CAR fibroblasts and in primary MEFs. A small number of genes is more susceptible to PPAR
-mediated transactivation, whereas PPARß/
is a poor activator of the genes investigated (Fig. 5A and B and Table 1), even though PPARß/
is transcriptionally competent in NIH-CAR cells, as shown in transient transfections (Fig. 2D).
PPAR/RXR binding to PPREs is a major determinant of subtype specificity.
The observed inability of the PPAR subtypes to activate certain target genes in the NIH-CAR cells could be due to a reduced potential of the PPARs to gain access to the target sites in the promoters/enhancers or due to the reduced ability of the bound PPAR/RXR heterodimer to induce formation of a transcriptionally active complex. Thus, to investigate how acute activation of PPAR target genes correlates with PPAR occupancy at the corresponding PPREs, we compared subtype-specific recruitment to target sites using ChIP. NIH-CAR cells were transduced with one of the AdHA-PPARs (
,
or ß/
) for 8 h in the presence of specific agonists. Equal expression of the PPAR subtypes at the 8-h time point was verified by Western blotting (Fig. 6A). ChIP was performed with antibodies directed against RXR and the HA epitope of the HA-PPARs, respectively. Input as well as immunoprecipitated DNA was quantified by real-time PCR using primers (available at http://www.sdu.dk/Nat/bmb/faculty/pubs/MCB06_supp_data.html) positioned at the PPREs of a number of different PPAR target genes (Table 2). To be able to compare the different target sites, recoveries were normalized to recoveries for nontransduced NIH-CAR cells (background). Although there is a slightly higher background when the HA antibody is used, there is a good correlation between the ChIP obtained with the HA and the RXR antibodies (Fig. 6B).
|
|
/RXR heterodimer to the PPREs of the A-FABP, perilipin, LPL, AQP7, ACOx1, HMG-CoA S2, HD, and ADRP genes, whereas there was only minor or no recruitment to the L-FABP gene (Fig. 6B). Interestingly, for most genes there was a good correlation between the ability of the subtypes to bind to target sites and their ability to activate transcription of the corresponding genes. The PPAR
/RXR and especially the PPARß/
/RXR heterodimer are much less efficiently bound to PPREs of the PPAR
-specific genes investigated, like A-FABP, perilipin, and AQP7, than is the PPAR
2/RXR heterodimer. This indicates that the failure of PPAR
and PPARß/
to potently activate transcription of these genes in NIH-CAR cells is due in part to poor binding to these target sites in chromatin. By contrast, PPAR
binds most efficiently to the HMG-CoA S2 and HD PPREs and is also the best activator of these genes. In the case of the silent gene L-FABP, which cannot be activated by PPARs in NIH-CAR cells, we did not detect significant PPAR/RXR binding to the corresponding PPRE. Thus, the ability to gain access to the target sites appears to be a major limiting factor in determining whether a gene can be activated by a particular PPAR subtype.
Although the general observation is that binding and activation are correlated, there are clear exceptions. For example, we detect significant binding of PPARß/
to target genes where PPAR
is unable to activate transcription. The most notable example is ADRP, which shows equal binding of all PPAR subtypes, yet the gene can be activated by only PPAR
and PPAR
2. In addition, the already actively transcribed ACOx1 and ME genes respond very modestly (twofold increase) (ACOx1) or not at all (ME) to ectopic PPAR expression (Fig. 6C and Table 1). Yet these genes display significant binding of the PPAR
2/RXR heterodimer to their proximal PPREs (Fig. 6B and C). To investigate whether the failure to further increase expression of the ME and ACO genes was due to the inability of PPAR
2 to increase cofactor occupancy, we determined CBP and TRAP220 recruitment to the PPREs. As seen in Fig. 6C, ectopic expression of PPAR
2 led to a significant increase in cofactor occupancy at both genes. However, in keeping with the high basal transcriptional activity and little PPAR responsiveness of the genes, Pol II occupancy at the proximal promoter was high in the nontransduced cells and only marginally increased by ectopic PPAR expression. This contrasts to the significant PPAR
2-dependent recruitment of Pol II at the A-FABP proximal promoter (Fig. 3C).
Taken together, PPAR
is the PPAR subtype that most efficiently gains access to a large number of PPREs in NIH-CAR cells. Furthermore, although exceptions exist, there is a good correlation between PPAR/RXR binding and target gene activation. Thus, the ability to gain access to the target PPREs appears to be a major, but clearly not the only, determinant in PPAR subtype specific activation of target gene expression.
It is likely that the differential ability of the PPAR subtypes to recruit cofactors to target sites plays an important role in determining PPAR subtype specificity at these sites. To make a first step toward elucidating the role of cofactors, we therefore compared the ability of the subtypes to recruit CBP and TRAP220 to different types of PPAR target genes in the presence of specific agonists. We tested PPAR
-selective target genes (A-FABP and perilipin), a PPAR
-selective target gene (HD), and a gene equally activated by PPAR
2 and PPAR
(ADRP) (Fig. 6D). Notably, recruitment of CBP and TRAP220 to the PPREs correlates to a large extent with binding of the different PPAR subtypes; i.e., subtypes that bind well also efficiently recruit CBP and TRAP220. Interestingly, however, PPARß/
is unable to efficiently recruit TRAP220 to the ADRP PPRE, notwithstanding the efficient binding of the PPARß/
/RXR heterodimer, and this may in part explain why PPARß/
is incapable of activating ADRP (Fig. 5A).
Target genes show different dependencies on agonists for transactivation.
The results presented in Fig. 4 for the A-FABP gene (as well as parallel analyses of the LPL, perilipin, and CD36 genes [data not shown]) clearly demonstrated that PPAR
2 displays significant agonist-independent as well as agonist-dependent binding to target sites and activation of target genes. To investigate to what degree exogenous agonists contribute to PPAR subtype specificity, we compared the activation of target genes by the different PPAR subtypes in the presence and absence of their respective agonists. Selected target genes are shown in Fig. 7A. Interestingly, the relative contribution of the agonists to the transactivation is highly dependent on the target gene and the PPAR subtype. In general, PPAR
transactivation is less dependent on exogenous agonist than PPAR
. Thus, in the absence of exogenous agonists, PPAR
and PPAR
activate perilipin equally well. In addition, the slight preference for PPAR
in the activation of the ADRP gene is reversed in the absence of agonists. However, for genes where the subtype specificity is very pronounced, like the A-FABP gene, subtype specificity is maintained also in the absence of exogenous agonists.
In keeping with the mRNA analysis, agonists had limited effect on binding of PPAR
/RXR and PPARß/
/RXR to the corresponding PPREs under the conditions investigated, while the effect on binding of PPAR
2/RXR is highly dependent on the target gene and to a large extent reflects the effects on gene expression (Fig. 7B). Interestingly, however, we did not observe any effect of agonists on PPAR/RXR binding to the ADRP PPRE. Thus, overall these results are in keeping with our general finding that the ability to gain access to target PPREs is a major, but clearly not the only, determinant in PPAR subtype-specific activation of target gene expression.
PPAR subtype specificity is affected by cell type. To investigate to what extent cell type affects PPAR subtype-specific transactivation of target genes, we compared the ability of HA-PPAR subtypes to induce target genes in NIH-CAR cells with the ability to induce these genes in two other murine cell culture models, the AML-12 hepatocyte cell line (69) and the MIN6 pancreatic ß-cell line (41). Experiments were performed in the presence of specific agonists. In contrast to NIH 3T3 fibroblasts, the AML-12 and MIN6 cells naturally express high levels of the CAR receptor and are therefore efficiently infected with adenoviruses (data not shown). For each cell line the infection ratios of the AdHA-PPARs were adjusted so that the level of HA-PPAR protein after 8 h of infection was similar to the one obtained in NIH-CAR cells (data not shown).
The ability of the PPAR subtypes to induce target genes following 12 h of transduction are summarized in Fig. 8A (for details, see the data at http://www.sdu.dk/Nat/bmb/faculty/pubs/MCB06_supp_data.html). It is evident that PPAR subtype specificity is highly cell-type dependent. In MIN6 and NIH-CAR, a few distinct target genes cannot be activated by any of the PPARs. These genes might be subject to permanent inactivation in the particular cell line (e.g., L-FABP in NIH-CAR and perilipin in MIN6), or they might be actively transcribed but insensitive to further activation by PPARs (e.g., ME in NIH-CAR and MIN6 cells and medium chain acyl-CoA dehydrogenase in MIN6 cells). However, among the target genes that are susceptible to PPAR transactivation, the cell type has great impact on subtype specificity. Thus, whereas PPARß/
is largely unable to directly induce endogenous target gene expression in NIH-CAR cells, it potently induces several target genes, such as ADRP and AQP7, in both AML-12 and MIN-6 cells. PPAR
is also significantly more active in inducing target genes, such as L-FABP and AQP7, in AML-12 and MIN-6 cells compared to NIH-CAR cells. Notably, in all three cell lines PPAR
2 is the only subtype able to potently activate the A-FABP and CD36 genes, suggesting that, irrespective of the subtype, these genes are highly PPAR
specific. By contrast, AQP7 which is a highly PPAR
2-specific target in NIH-CAR cells, is activated by all PPAR subtypes in AML-12 and MIN6 cells. Taken together, our data show that the cell type plays a determining role in rendering a particular target gene susceptible to different PPAR subtypes.
|
-specific genes (exemplified by A-FABP), PPAR
/
target genes (exemplified by HD), and genes activated by all PPAR subtypes (exemplified by ADRP). | DISCUSSION |
|---|
|
|
|---|
Here, we demonstrate that the combination of adenoviral delivery and ChIP is ideal for studying the acute molecular effects of PPAR activation on PPAR target sites and promoters. We show that adenoviral transduction of the HA-PPARs in NIH-CAR fibroblasts results in increased PPAR mRNA levels after 2 h and potent PPAR protein expression within 4 h posttransduction. The HA-PPAR proteins reach maximal levels, similar to that of PPAR
protein in 3T3-L1 adipocytes, 4 to 8 h after transduction (Fig. 2). Binding of PPAR
, RXR, cofactors, and basal transcription factors to the enhancer of the well-established PPAR
target gene, A-FABP, is maximal at 5 h. At this point, the expression from the A-FABP gene is significantly activated. Recruitment of PPAR
to the locus also leads to a significant increase in histone H3 and H4 acetylation at the promoter as well as at the enhancer (Fig. 3). Thus, the adenoviral transgene delivery provides a unique system to tightly control the expression of PPARs and to study the acute and direct effects of PPAR binding to target sites of endogenous genes.
As PPARs are activated by fatty acids and derivatives thereof, the "basal" transcriptional activity (i.e., activity in the absence of exogenous agonists) of PPARs is often ascribed to endogenous ligands or ligands present in the serum. In addition, the ligand-independent activation function of the receptors may contribute to this activity; however, the relative contribution of endogenous ligands and ligand-independent function to "basal" transcriptional activity is unknown. Using the adenoviral delivery system, we show here that PPAR
2 displays significant ligand-dependent as well as -independent transcriptional activity. PPAR
2 leads to profound activation of the A-FABP gene and other target genes like LPL and perilipin in the absence of exogenous ligands. This activation is only slightly decreased by the addition of a specific PPAR
antagonist, GW9662 (Fig. 4A and results not shown). Thus, the considerable transcriptional potential of PPAR
2 in NIH-CAR fibroblasts in the absence of exogenous agonists is predominantly a consequence of the high ligand-independent activity of PPAR
2.
The ligand-independent transcriptional activity of PPAR
is likely to be mediated at least in part by ligand-independent activation function 1 (AF1) in the N-terminal part of the receptor. However, the AF2 of the ligand binding domain (LBD) may also contribute to the ligand-independent transactivation. Helix 12 has been reported to be in an "agonist position" in crystal structures of apo-PPAR
LBD (46), and recent investigations indicated that this is due to the stabilization of helix 12 in the agonist position by specific charged amino acids in the LBD (42). In keeping with this, PPAR
LBD has recently been shown to interact with CBP/p300 (33) as well as p160 CoA (42) in the absence of ligands. Presence of a synthetic ligand further induces stabilization of helix 12, thereby leading to an increased interaction with coactivators (42). Future investigations in our adenoviral system with truncated and mutated PPARs will help clarify whether the ligand-independent transactivation of endogenous genes is due primarily to the AF1 or the AF2 function and whether these effects are promoter specific. Recent data from the Lazar laboratory indicate that there is indeed promoter specificity in the dependency of ligand for cofactor recruitment (21).
Intriguingly, our data show that PPAR
occupancy of the A-FABP PPRE is considerably increased in the presence of the thiazolidinedione-type agonist, rosiglitazone. The increased occupancy is paralleled by increased recruitment of basal transcription factors and of coactivators such as CBP. This ligand-dependent increase in occupancy may in part be explained by a stabilization of the PPAR/RXR/DNA interaction, as suggested by electrophoretic mobility shift assays of PPAR
/RXR (15, 17). In addition, ligand-induced stabilization of receptor-cofactor interactions may lead to a longer residence time of the receptor at the locus. Interestingly, however, the relative importance of agonist for PPAR
2/RXR binding to target sites is highly dependent on the target gene, and for some genes like ADRP, agonists did not appear to affect PPAR
2/RXR binding although it potentiated expression of the gene (Fig. 7A and B).
In general, the effect of the agonist Wy14.643 on PPAR
binding to target sites and activation of target genes is minor (Fig. 7A and B). However, the lack of efficient synthetic PPAR
antagonists for cell culture experiments precludes investigations similar to those shown in Fig. 4 for PPAR
to determine whether the activity in the absence of exogenous agonist is due to agonist-independent transactivation or to the presence of endogenous agonists.
It is clear from numerous in vivo and ex vivo experiments that the different PPAR subtypes serve different physiological functions. However, it is unclear to what extent these different actions are due to intrinsic biochemical properties or to differences in the cellular context in which the PPARs are expressed. Only a few studies have sought to mutually compare PPAR subtype-specific activation of endogenous genes in a particular cell type (7, 60). Although these comparisons have been very useful, they were performed with stable ectopic expression, which makes it difficult to distinguish between direct and indirect target genes. Furthermore, constitutive expression of the receptors hampers investigations of the molecular mechanism involved in target gene activation. The combination of adenoviral delivery and ChIP provides a unique possibility to investigate the molecular mechanisms underlying PPAR subtype specificity by comparing the acute actions of the different PPAR subtypes on endogenous chromatin-embedded target genes. We show that, when adenovirally expressed at equal levels, HA-tagged PPARs have highly different potentials to activate endogenous target genes in the NIH-CAR cells. While PPAR
had more transcriptional activity than PPAR
in transient transfections with a reporter construct consisting of a mulitmerized ACOx1 PPRE, PPAR
was the most potent activator of the majority of the endogenous genes investigated. However, a small number of genes were most potently activated by PPAR
. Surprisingly, however, in NIH-CAR cells as well as in MEFs, PPARß/
was either unable to activate or a poor activator of all genes investigated (Fig. 5A and B and Table 1). The reason for this is unclear, as adenovirally expressed PPARß/
was active in transient transfections of NIH-CAR cells and able to activate endogenous genes in AML-12 and MIN6 cells. Previous experiments with stable retroviral expression of PPARß/
in NIH 3T3 cells have shown that, although PPARß/
is by far the least potent PPAR subtype in these cells, constitutive PPARß/
expression induces low levels of several PPAR target genes (7). However, from these long-term experiments it is impossible to conclude whether the induction is due to direct but slow induction of the target genes or due to indirect activation of the target genes via induction of another transcription factor, e.g., PPAR
. The lack of target gene activation by adenovirally expressed PPARß/
suggests the existence of a global mechanism (e.g., posttranslational modification or lack of a particular cofactor) that renders PPARß/
less transcriptionally active on endogenous genes in the NIH 3T3 fibroblasts.
To address the question how the PPARs retain gene specificity, we used ChIP to investigate if PPAR subtype-specific transcriptional activation is a consequence of differential PPAR binding to promoters or to PPAR subtype-specific recruitment of coactivators. Interestingly, ChIP analysis demonstrated that, in general, the transcriptional potential of the PPAR subtypes is correlated with the level of PPAR/RXR heterodimer binding to the promoters of activated genes (Fig. 6B). There are clear exceptions where a subtype binds without further activating; however, for all target genes investigated, the subtype that binds best is also the subtype that activates most efficiently. This indicates that a major, but clearly not the only, rate-limiting step in establishing a transcriptionally active complex at the target genes investigated is the ability of the receptor subtypes to gain access to the chromatin-embedded binding sites.
Some target genes remained transcriptionally silent in NIH-CAR cells irrespective of the PPAR subtype. As exemplified by the L-FABP gene, we did not observe PPAR/RXR binding to these target sites, and hence these are likely to be embedded in inaccessible chromatin or to require synergy with other cell-type-specific transcription (co)factors to gain access. Other genes (e.g., ME and thiolase B) were transcriptionally active but insensitive to further activation by PPARs. Interestingly, however, PPAR
2/RXR displays significant binding and PPAR
2-dependent cofactor recruitment to the PPREs of the ME gene (Fig. 6C). This suggests that other (co)factors than PPAR/RXR, CBP, and TRAP220 are limiting the expression of this gene. Further detailed ChIP analyses of factor recruitment to this and other loci will be required to understand why receptor binding is unable to further enhance the transcription of such genes.
The differential binding of PPAR subtypes to different PPREs may in part be explained by the sequence composition of the PPREs. In vitro experiments performed in the absence of ligands have shown that PPREs with a relatively low match to the consensus have a preference for PPAR
heterodimers (27). Such "weak" PPREs include the two PPREs of the A-FABP enhancer. However, the LPL PPRE has been shown to bind both PPAR
and -
in vitro and is a recognized PPAR
as well as PPAR
target gene (39, 55). In addition, the ACOx1 PPRE(s) match the consensus very well and is generally regarded as a PPAR
-type PPRE (64, 65). For both the LPL and ACO PPREs, we find that they preferably bind PPAR
heterodimers in the context of chromatin in NIH-CAR cells, indicating that factors other than the affinity for the naked DNA sequence play a role. Furthermore, as discussed below, PPAR subtype specificity differs between the cell lines. Thus, the abilities of the PPAR subtypes to gain access to PPREs in vivo cannot be explained only by their differential affinity for a given PPRE in vitro. It is likely that an important determining factor is the synergy with other transcription factors and the joint ability to effectively recruit the right combination of chromatin remodeling and modifying complexes present in the cell.
In addition to the ability of the PPAR/RXR heterodimer to gain access to target sites, the ability of the PPAR/RXR heterodimer to bind cofactors prior to DNA-binding or to recruit cofactors to target genes once bound to DNA is likely to greatly influence the level of expression of the target genes. An exhaustive analysis of this aspect is beyond the scope of this article; however, we have investigated the differential ability of the PPAR subtypes to recruit TRAP220 and CBP to a small number of selected target genes in NIH-CAR cells (Fig. 6D). Importantly, for most genes the recruitment of these more general cofactors correlates with binding of the PPAR subtypes and therefore does not add to the subtype specificity. An exception may be the ADRP gene to which all PPAR subtypes bind, whereas only PPAR
and PPAR
2 are able to significantly recruit TRAP220. Interestingly, this may in part explain why PPARß/
is unable to activate the ADRP gene. Further insight into the role of cofactors in PPAR subtype specificity for different genes awaits large-scale ChIP analysis of more specific cofactors.
The poor PPAR
and PPARß/
transcriptional activity of endogenous genes in NIH 3T3 fibroblasts prompted us to investigate the ability of PPAR
and PPARß/
to activate endogenous genes in other cell lines such as AML-12 hepatocytes and MIN6 pancreatic ß cells. From these experiments it is evident that the cell type is indeed a very important parameter of PPAR subtype specificity toward several genes (Fig. 8A). PPAR
remains the subtype able to activate the largest number of genes in the cell types investigated. However, for the target genes investigated PPAR
and PPARß/
are on average more transcriptionally active in AML-12 and MIN6 cells than in NIH-CAR cells. Irrespective of the cell line, the A-FABP and CD36 genes remain highly PPAR
specific, but other genes that are PPAR
specific in NIH-CAR cells can be induced by PPAR
and/or PPARß/
in AML-12 and MIN6 cells. Similarly, HMG-CoA S2 and HD, which are PPAR
selective in NIH-CAR cells, are induced by other PPAR subtypes in AML-12 and MIN6 cells. These results show that PPAR subtype specificity is a highly complex phenomenon likely to be dependent on the cell-type-specific setting of the individual target sites, i.e., a combination of other transcription factors, cofactors, and chromatin modifications in a cell. Importantly, however, as observed for the NIH-CAR cells, there was a good correlation between the ability of the PPAR subtypes to bind to target sites and their ability to activate the corresponding gene (Fig. 8B and C).
In conclusion, we demonstrate that the combination of adenoviral delivery and ChIP is an ideal method for investigating the basic molecular events involved in the activation of genomic target genes by PPARs in a broad range of cell lines. Using this approach, we show that PPAR
2 possesses considerable agonist-dependent as well as -independent transactivation potential and that agonist increases the occupancy of PPAR
2/RXR at most but not all PPREs. The ability of a particular PPAR subtype to gain access to a PPRE is highly dependent on the subtype as well as the cell type. However, in general (although there are clear exceptions) there is a good correlation between the extent of PPAR binding and the magnitude of the target gene induction. Future challenges will be to investigate why the PPAR subtypes differ in their ability to gain access to endogenous chromatin-embedded PPREs and how the cell type regulates this ability. In addition, large-scale ChIP analyses will be required to define the role of different cofactors in PPAR subtype-specific transactivation.
| ACKNOWLEDGMENTS |
|---|
We are grateful to D. J. Orlicky for the LXSN-hCAR
1 and AdGFP constructs, R. M. Evans for the PPRE3-TK-LUC reporter, Jun-ichi Mi