Isoform-specific functions of PPARγ in gene regulation and metabolism.

Peroxisome proliferator-activated receptor γ (PPARγ) is a nuclear receptor that is a vital regulator of adipogenesis, insulin sensitivity, and lipid metabolism. Activation of PPARγ by antidiabetic thiazolidinediones (TZD) reverses insulin resistance but also leads to weight gain that limits the use of these drugs. There are two main PPARγ isoforms, but the specific functions of each are not established. Here we generated mouse lines in which endogenous PPARγ1 and PPARγ2 were epitope-tagged to interrogate isoform-specific genomic binding, and mice deficient in either PPARγ1 or PPARγ2 to assess isoform-specific gene regulation. Strikingly, although PPARγ1 and PPARγ2 contain identical DNA binding domains, we uncovered isoform-specific genomic binding sites in addition to shared sites. Moreover, PPARγ1 and PPARγ2 regulated a different set of genes in adipose tissue depots, suggesting distinct roles in adipocyte biology. Indeed, mice with selective deficiency of PPARγ1 maintained body temperature better than wild-type or PPARγ2-deficient mice. Most remarkably, although TZD treatment improved glucose tolerance in mice lacking either PPARγ1 or PPARγ2, the PPARγ1-deficient mice were protected from TZD-induced body weight gain compared with PPARγ2-deficient mice. Thus, PPARγ isoforms have specific and separable metabolic functions that may be targeted to improve therapy for insulin resistance and diabetes.

Cell-Intrinsic Tumorigenic Functions of PPARγ in Bladder Urothelial Carcinoma.

The role of PPAR gamma (PPARγ) has been well characterized in the developmental process of adipogenesis, yet its aberrant expression patterns and functions in cancer subtypes are less understood. Although PPARγ has been recently demonstrated to play non-cell-autonomous roles in promoting bladder urothelial carcinoma (UC) progression, underlying mechanisms of the cell-intrinsic oncogenic activity remain unknown. Here, we report robust expression and nuclear accumulation of PPARγ in 47% of samples of patients with UC, exceeding mRNA expression patterns published by The Cancer Genome Atlas. In vitro assays revealed for the first time that treatment of UC cells with PPARγ inverse agonist or PPARG knockout by CRISPR-Cas9 reduces proliferation, migration, and invasion of multiple established UC cell lines, most strongly in those characterized by PPARG genomic amplification or activating mutations of RXRA, the obligate heterodimer of PPARγ. Through genome-wide approaches including chromatin immunoprecipitation sequencing and RNA sequencing, we define a novel set of PPARγ-regulated genes in UC, including Sonic Hedgehog (SHH). Similar to PPARγ, genetic inhibition of SHH reduces proliferation and motility. Finally, we demonstrate the PPARγ dependency of UC tumors in vivo by genetic and pharmacologic PPARγ inhibition in subcutaneous xenografts. Collectively, our data indicate that PPARγ promotes UC progression in a subset of patients, at least in part, through cell-autonomous mechanisms linked to SHH signaling. IMPLICATIONS: Genome-wide analysis of DNA-binding sites for oncogenic factor PPARγ revealed SHH as a novel downstream target involved in UC progression, providing important insight into the tumorigenic nature and molecular mechanism of PPARγ signaling in UC.

Shared PPARα/γ Target Genes Regulate Brown Adipocyte Thermogenic Function.

Brown adipose tissue (BAT) generates heat to maintain body temperature and suppress obesity. Agonists for nuclear receptors PPARα and PPARγ both affect brown adipocyte function, yet the interplay between these factors in BAT is uncertain. Here, we report that PPARα shares most genomic binding sites with PPARγ, and these common binding sites are more related to BAT function than PPARγ-selective sites without PPARα. Integrating PPARα and PPARγ genomic occupancy with cold-responsive BAT transcriptomes identifies a subset of 16 genes with potential relevance to BAT function. Among these, we focused on the lysosomal protease cathepsin Z (CTSZ) and showed it is necessary for mitochondrial respiration in both mouse and human brown adipocytes. Thus, CTSZ is a shared PPARα/γ target gene in BAT and a regulator of brown adipocyte thermogenic function.

Patient Adipose Stem Cell-Derived Adipocytes Reveal Genetic Variation that Predicts Antidiabetic Drug Response.

Thiazolidinedione drugs (TZDs) target the transcriptional activity of peroxisome proliferator activated receptor γ (PPARγ) to reverse insulin resistance in type 2 diabetes, but side effects limit their clinical use. Here, using human adipose stem cell-derived adipocytes, we demonstrate that SNPs were enriched at sites of patient-specific PPARγ binding, which correlated with the individual-specific effects of the TZD rosiglitazone (rosi) on gene expression. Rosi induction of ABCA1, which regulates cholesterol metabolism, was dependent upon SNP rs4743771, which modulated PPARγ binding by influencing the genomic occupancy of its cooperating factor, NFIA. Conversion of rs4743771 from the inactive SNP allele to the active one by clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated editing rescued PPARγ binding and rosi induction of ABCA1 expression. Moreover, rs4743771 is a major determinant of undesired serum cholesterol increases in rosi-treated diabetics. These data highlight human genetic variation that impacts PPARγ genomic occupancy and patient responses to antidiabetic drugs, with implications for developing personalized therapies for metabolic disorders.

Shared nucleotide flanks confer transcriptional competency to bZip core motifs.

Sequence-specific DNA binding recruits transcription factors (TFs) to the genome to regulate gene expression. Here, we perform high resolution mapping of CEBP proteins to determine how sequence dictates genomic occupancy. We demonstrate a fundamental difference between the sequence repertoire utilized by CEBPs in vivo versus the palindromic sequence preference reported by classical in vitro models, by identifying a palindromic motif at <1% of the genomic binding sites. On the native genome, CEBPs bind a diversity of related 10 bp sequences resulting from the fusion of degenerate and canonical half-sites. Altered DNA specificity of CEBPs in cells occurs through heterodimerization with other bZip TFs, and approximately 40% of CEBP-binding sites in primary human cells harbor motifs characteristic of CEBP heterodimers. In addition, we uncover an important role for sequence bias at core-motif-flanking bases for CEBPs and demonstrate that flanking bases regulate motif function across mammalian bZip TFs. Favorable flanking bases confer efficient TF occupancy and transcriptional activity, and DNA shape may explain how the flanks alter TF binding. Importantly, motif optimization within the 10-mer is strongly correlated with cell-type-independent recruitment of CEBPß, providing key insight into how sequence sub-optimization affects genomic occupancy of widely expressed CEBPs across cell types.

PPARγ is dispensable for clear cell renal cell carcinoma progression.

OBJECTIVE: Clear cell renal cell carcinoma (ccRCC) is a subtype of kidney cancer defined by robust lipid accumulation, which prior studies have indicated plays an important role in tumor progression. We hypothesized that the peroxisome proliferator-activated receptor gamma (PPARγ), detected in both ccRCC tumors and cell lines, promotes lipid storage in ccRCC and contributes to tumorigenesis in this setting. PPARγ transcriptionally regulates a number of genes involved in lipid and glucose metabolism in adipocytes, yet its role in ccRCC has not been described. The objective of this study was to elucidate endogenous PPARγ function in ccRCC cells. METHODS AND RESULTS: Using chromatin immunoprecipitation followed by deep sequencing (ChIP-seq), we found that PPARγ and its heterodimer RXR occupy the canonical DR1 PPAR binding motif at approximately 1000 locations throughout the genome that can be subdivided into adipose-shared and ccRCC-specific sites. CRISPR-Cas9 mediated, loss-of-function studies determined that PPARγ is dispensable for viability, proliferation, and migration of ccRCC cells in vitro and in vivo. Also, surprisingly, PPARγ deletion had little effect on the robust lipid accumulation that typifies the "clear cell" phenotype of kidney cancer. CONCLUSION: Our results suggest that PPARγ plays neither a tumor suppressive nor oncogenic role in advanced ccRCC, and thus single-agent therapeutics targeting PPARγ are unlikely to be effective for the treatment of this disease. The unique cistrome of PPARγ in ccRCC cells demonstrates the importance of cell type in determining the functions of PPARγ.

Rev-erbα dynamically modulates chromatin looping to control circadian gene transcription.

Mammalian physiology exhibits 24-hour cyclicity due to circadian rhythms of gene expression controlled by transcription factors that constitute molecular clocks. Core clock transcription factors bind to the genome at enhancer sequences to regulate circadian gene expression, but not all binding sites are equally functional. We found that in mice, circadian gene expression in the liver is controlled by rhythmic chromatin interactions between enhancers and promoters. Rev-erbα, a core repressive transcription factor of the clock, opposes functional loop formation between Rev-erbα-regulated enhancers and circadian target gene promoters by recruitment of the NCoR-HDAC3 co-repressor complex, histone deacetylation, and eviction of the elongation factor BRD4 and the looping factor MED1. Thus, a repressive arm of the molecular clock operates by rhythmically modulating chromatin loops to control circadian gene transcription.

Histone deacetylase 3 prepares brown adipose tissue for acute thermogenic challenge.

Brown adipose tissue is a thermogenic organ that dissipates chemical energy as heat to protect animals against hypothermia and to counteract metabolic disease. However, the transcriptional mechanisms that determine the thermogenic capacity of brown adipose tissue before environmental cold are unknown. Here we show that histone deacetylase 3 (HDAC3) is required to activate brown adipose tissue enhancers to ensure thermogenic aptitude. Mice with brown adipose tissue-specific genetic ablation of HDAC3 become severely hypothermic and succumb to acute cold exposure. Uncoupling protein 1 (UCP1) is nearly absent in brown adipose tissue lacking HDAC3, and there is also marked downregulation of mitochondrial oxidative phosphorylation genes resulting in diminished mitochondrial respiration. Remarkably, although HDAC3 acts canonically as a transcriptional corepressor, it functions as a coactivator of oestrogen-related receptor α (ERRα) in brown adipose tissue. HDAC3 coactivation of ERRα is mediated by deacetylation of PGC-1α and is required for the transcription of Ucp1, Ppargc1a (encoding PGC-1α), and oxidative phosphorylation genes. Importantly, HDAC3 promotes the basal transcription of these genes independently of adrenergic stimulation. Thus, HDAC3 uniquely primes Ucp1 and the thermogenic transcriptional program to maintain a critical capacity for thermogenesis in brown adipose tissue that can be rapidly engaged upon exposure to dangerously cold temperature.

Nuclear Receptor Function through Genomics: Lessons from the Glucocorticoid Receptor.

Unlocking the therapeutic potential of the glucocorticoid receptor (GR) has motivated a search for small molecules that selectively modulate its ability to activate or repress gene transcription. Recently, breakthrough studies in the field of genomics have reinvigorated debate over longstanding transcriptional models explaining how GR controls tissue-specific gene expression. Here, we highlight these genomic studies with the dual goals of advancing understanding of nuclear receptor-mediated transcription and stimulating thought on the development of anti-inflammatory and immunosuppressive ligands for GR that have reduced harmful effects on metabolism.

Targeting PPARγ in the epigenome rescues genetic metabolic defects in mice.

Obesity causes insulin resistance, and PPARγ ligands such as rosiglitazone are insulin sensitizing, yet the mechanisms remain unclear. In C57BL/6 (B6) mice, obesity induced by a high-fat diet (HFD) has major effects on visceral epididymal adipose tissue (eWAT). Here, we report that HFD-induced obesity in B6 mice also altered the activity of gene regulatory elements and genome-wide occupancy of PPARγ. Rosiglitazone treatment restored insulin sensitivity in obese B6 mice, yet, surprisingly, had little effect on gene expression in eWAT. However, in subcutaneous inguinal fat (iWAT), rosiglitazone markedly induced molecular signatures of brown fat, including the key thermogenic gene Ucp1. Obesity-resistant 129S1/SvImJ mice (129 mice) displayed iWAT browning, even in the absence of rosiglitazone. The 129 Ucp1 locus had increased PPARγ binding and gene expression that were preserved in the iWAT of B6x129 F1-intercrossed mice, with an imbalance favoring the 129-derived alleles, demonstrating a cis-acting genetic difference. Thus, B6 mice have genetically defective Ucp1 expression in iWAT. However, when Ucp1 was activated by rosiglitazone, or by iWAT browning in cold-exposed or young mice, expression of the B6 version of Ucp1 was no longer defective relative to the 129 version, indicating epigenomic rescue. These results provide a framework for understanding how environmental influences like drugs can affect the epigenome and potentially rescue genetically determined disease phenotypes.

The Nuclear Receptor Rev-erbα Regulates Adipose Tissue-specific FGF21 Signaling.

FGF21 is an atypical member of the FGF family that functions as a hormone to regulate carbohydrate and lipid metabolism. Here we demonstrate that the actions of FGF21 in mouse adipose tissue, but not in liver, are modulated by the nuclear receptor Rev-erbα, a potent transcriptional repressor. Interrogation of genes induced in the absence of Rev-erbα for Rev-erbα-binding sites identified ßKlotho, an essential coreceptor for FGF21, as a direct target gene of Rev-erbα in white adipose tissue but not liver. Rev-erbα ablation led to the robust elevated expression of ßKlotho. Consequently, the effects of FGF21 were markedly enhanced in the white adipose tissue of mice lacking Rev-erbα. A major Rev-erbα-controlled enhancer at the Klb locus was also bound by the adipocytic transcription factor peroxisome proliferator-activated receptor (PPAR) γ, which regulates its activity in the opposite direction. These findings establish Rev-erbα as a specific modulator of FGF21 signaling in adipose tissue.

Tribbles-1 regulates hepatic lipogenesis through posttranscriptional regulation of C/EBPα.

Variants near the gene TRIB1 are significantly associated with several plasma lipid traits, circulating liver enzymes, and the development of coronary artery disease in humans; however, it is not clear how its protein product tribbles-1 regulates lipid metabolism. Here, we evaluated mice harboring a liver-specific deletion of Trib1 (Trib1_LSKO) to elucidate the role of tribbles-1 in mammalian hepatic lipid metabolism. These mice exhibited increased hepatic triglyceride (TG) content, lipogenic gene transcription, and de novo lipogenesis. Microarray analysis revealed altered transcription of genes that are downstream of the transcription factor C/EBPα, and Trib1_LSKO mice had increased hepatic C/EBPα protein. Hepatic overexpression of C/EBPα in WT mice phenocopied Trib1_LSKO livers, and hepatic knockout of Cebpa in Trib1_LSKO mice revealed that C/EBPα is required for the increased lipogenesis. Using ChIP-Seq, we found that Trib1_LSKO mice had increased DNA-bound C/EBPα near lipogenic genes and the Trib1 gene, which itself was transcriptionally upregulated by C/EBPα overexpression. Together, our results reveal that tribbles-1 regulates hepatic lipogenesis through posttranscriptional regulation of C/EBPα, which in turn transcriptionally upregulates Trib1. These data suggest an important role for C/EBPα in mediating the lipogenic effects of hepatic Trib1 deletion and provide insight into the association between TRIB1 and plasma lipids, and liver traits in humans.

Genetic Variation Determines PPARγ Function and Anti-diabetic Drug Response In Vivo.

SNPs affecting disease risk often reside in non-coding genomic regions. Here, we show that SNPs are highly enriched at mouse strain-selective adipose tissue binding sites for PPARγ, a nuclear receptor for anti-diabetic drugs. Many such SNPs alter binding motifs for PPARγ or cooperating factors and functionally regulate nearby genes whose expression is strain selective and imbalanced in heterozygous F1 mice. Moreover, genetically determined binding of PPARγ accounts for mouse strain-specific transcriptional effects of TZD drugs, providing proof of concept for personalized medicine related to nuclear receptor genomic occupancy. In human fat, motif-altering SNPs cause differential PPARγ binding, provide a molecular mechanism for some expression quantitative trait loci, and are risk factors for dysmetabolic traits in genome-wide association studies. One PPARγ motif-altering SNP is associated with HDL levels and other metabolic syndrome parameters. Thus, natural genetic variation in PPARγ genomic occupancy determines individual disease risk and drug response.

ATF4 licenses C/EBPß activity in human mesenchymal stem cells primed for adipogenesis.

A well-established cascade of transcription factor (TF) activity orchestrates adipogenesis in response to chemical cues, yet how cell-intrinsic determinants of differentiation such as cell shape and/or seeding density inform this transcriptional program remain enigmatic. Here, we uncover a novel mechanism licensing transcription in human mesenchymal stem cells (hMSCs) adipogenically primed by confluence. Prior to adipogenesis, confluency promotes heterodimer recruitment of the bZip TFs C/EBPß and ATF4 to a non-canonical C/EBP DNA sequence. ATF4 depletion decreases both cell-density-dependent transcription and adipocyte differentiation. Global profiling in hMSCs and a novel cell-free assay reveals that ATF4 requires C/EBPß for genomic binding at a motif distinct from that bound by the C/EBPß homodimer. Our observations demonstrate that C/EBPß bridges the transcriptional programs in naïve, confluent cells and early differentiating pre-adipocytes. Moreover, they suggest that homo- and heterodimer formation poise C/EBPß to execute diverse and stage-specific transcriptional programs by exploiting an expanded motif repertoire.

GENE REGULATION. Discrete functions of nuclear receptor Rev-erbα couple metabolism to the clock.

Circadian and metabolic physiology are intricately intertwined, as illustrated by Rev-erbα, a transcription factor (TF) that functions both as a core repressive component of the cell-autonomous clock and as a regulator of metabolic genes. Here, we show that Rev-erbα modulates the clock and metabolism by different genomic mechanisms. Clock control requires Rev-erbα to bind directly to the genome at its cognate sites, where it competes with activating ROR TFs. By contrast, Rev-erbα regulates metabolic genes primarily by recruiting the HDAC3 co-repressor to sites to which it is tethered by cell type-specific transcription factors. Thus, direct competition between Rev-erbα and ROR TFs provides a universal mechanism for self-sustained control of the molecular clock across all tissues, whereas Rev-erbα uses lineage-determining factors to convey a tissue-specific epigenomic rhythm that regulates metabolism tailored to the specific need of that tissue.

Genomic redistribution of GR monomers and dimers mediates transcriptional response to exogenous glucocorticoid in vivo.

Glucocorticoids (GCs) are commonly prescribed drugs, but their anti-inflammatory benefits are mitigated by metabolic side effects. Their transcriptional effects, including tissue-specific gene activation and repression, are mediated by the glucocorticoid receptor (GR), which is known to bind as a homodimer to a palindromic DNA sequence. Using ChIP-exo in mouse liver under endogenous corticosterone exposure, we report here that monomeric GR interaction with a half-site motif is more prevalent than homodimer binding. Monomers colocalize with lineage-determining transcription factors in both liver and primary macrophages, and the GR half-site motif drives transcription, suggesting that monomeric binding is fundamental to GR's tissue-specific functions. In response to exogenous GC in vivo, GR dimers assemble on chromatin near ligand-activated genes, concomitant with monomer evacuation of sites near repressed genes. Thus, pharmacological GCs mediate gene expression by favoring GR homodimer occupancy at classic palindromic sites at the expense of monomeric binding. The findings have important implications for improving therapies that target GR.

PRDM16 binds MED1 and controls chromatin architecture to determine a brown fat transcriptional program.

PR (PRD1-BF1-RIZ1 homologous) domain-containing 16 (PRDM16) drives a brown fat differentiation program, but the mechanisms by which PRDM16 activates brown fat-selective genes have been unclear. Through chromatin immunoprecipitation (ChIP) followed by deep sequencing (ChIP-seq) analyses in brown adipose tissue (BAT), we reveal that PRDM16 binding is highly enriched at a broad set of brown fat-selective genes. Importantly, we found that PRDM16 physically binds to MED1, a component of the Mediator complex, and recruits it to superenhancers at brown fat-selective genes. PRDM16 deficiency in BAT reduces MED1 binding at PRDM16 target sites and causes a fundamental change in chromatin architecture at key brown fat-selective genes. Together, these data indicate that PRDM16 controls chromatin architecture and superenhancer activity in BAT.

Anti-diabetic rosiglitazone remodels the adipocyte transcriptome by redistributing transcription to PPARγ-driven enhancers.

Rosiglitazone (rosi) is a powerful insulin sensitizer, but serious toxicities have curtailed its widespread clinical use. Rosi functions as a high-affinity ligand for peroxisome proliferator-activated receptor γ (PPARγ), the adipocyte-predominant nuclear receptor (NR). The classic model, involving binding of ligand to the NR on DNA, explains positive regulation of gene expression, but ligand-dependent repression is not well understood. We addressed this issue by studying the direct effects of rosi on gene transcription using global run-on sequencing (GRO-seq). Rosi-induced changes in gene body transcription were pronounced after 10 min and correlated with steady-state mRNA levels as well as with transcription at nearby enhancers (enhancer RNAs [eRNAs]). Up-regulated eRNAs occurred almost exclusively at PPARγ-binding sites, to which rosi treatment recruited coactivators, including MED1, p300, and CBP. In contrast, transcriptional repression by rosi involved a loss of coactivators from eRNA sites devoid of PPARγ and enriched for other transcription factors, including AP-1 factors and C/EBPs. Thus, rosi activates and represses transcription by fundamentally different mechanisms that could inform the future development of anti-diabetic drugs.

Integrative genomic analysis of CREB defines a critical role for transcription factor networks in mediating the fed/fasted switch in liver.

BACKGROUND: Metabolic homeostasis in mammals critically depends on the regulation of fasting-induced genes by CREB in the liver. Previous genome-wide analysis has shown that only a small percentage of CREB target genes are induced in response to fasting-associated signaling pathways. The precise molecular mechanisms by which CREB specifically targets these genes in response to alternating hormonal cues remain to be elucidated. RESULTS: We performed chromatin immunoprecipitation coupled to high-throughput sequencing of CREB in livers from both fasted and re-fed mice. In order to quantitatively compare the extent of CREB-DNA interactions genome-wide between these two physiological conditions we developed a novel, robust analysis method, termed the 'single sample independence' (SSI) test that greatly reduced the number of false-positive peaks. We found that CREB remains constitutively bound to its target genes in the liver regardless of the metabolic state. Integration of the CREB cistrome with expression microarrays of fasted and re-fed mouse livers and ChIP-seq data for additional transcription factors revealed that the gene expression switches between the two metabolic states are associated with co-localization of additional transcription factors at CREB sites. CONCLUSIONS: Our results support a model in which CREB is constitutively bound to thousands of target genes, and combinatorial interactions between DNA-binding factors are necessary to achieve the specific transcriptional response of the liver to fasting. Furthermore, our genome-wide analysis identifies thousands of novel CREB target genes in liver, and suggests a previously unknown role for CREB in regulating ER stress genes in response to nutrient influx.

An Evi1-C/EBPß complex controls peroxisome proliferator-activated receptor γ2 gene expression to initiate white fat cell differentiation.

Fibroblastic preadipocyte cells are recruited to differentiate into new adipocytes during the formation and hyperplastic growth of white adipose tissue. Peroxisome proliferator-activated receptor γ (PPARγ), the master regulator of adipogenesis, is expressed at low levels in preadipocytes, and its levels increase dramatically and rapidly during the differentiation process. However, the mechanisms controlling the dynamic and selective expression of PPARγ in the adipocyte lineage remain largely unknown. We show here that the zinc finger protein Evi1 increases in preadipocytes at the onset of differentiation prior to increases in PPARγ levels. Evi1 expression converts nonadipogenic cells into adipocytes via an increase in the predifferentiation levels of PPARγ2, the adipose-selective isoform of PPARγ. Conversely, loss of Evi1 in preadipocytes blocks the induction of PPARγ2 and suppresses adipocyte differentiation. Evi1 binds with C/EBPß to regulatory sites in the Pparγ locus at early stages of adipocyte differentiation, coincident with the induction of Pparγ2 expression. These results indicate that Evi1 is a key regulator of adipogenic competency.