Transcription Factor Cooperativity in Early Adipogenic Hotspots and Super-Enhancers

Pan-cancer analyses of the nuclear receptor superfamily

Nuclear receptors (NR) act as an integrated conduit for environmental and hormonal signals to govern genomic responses, which relate to cell fate decisions. We review how their integrated actions with each other, shared co-factors and other transcription factors are disrupted in cancer. Steroid hormone nuclear receptors are oncogenic drivers in breast and prostate cancer and blockade of signaling is a major therapeutic goal. By contrast to blockade of receptors, in other cancers enhanced receptor function is attractive, as illustrated initially with targeting of retinoic acid receptors in leukemia. In the post-genomic era large consortia, such as The Cancer Genome Atlas, have developed a remarkable volume of genomic data with which to examine multiple aspects of nuclear receptor status in a pan-cancer manner. Therefore to extend the review of NR function we have also undertaken bioinformatics analyses of NR expression in over 3000 tumors, spread across six different tumor types (bladder, breast, colon, head and neck, liver and prostate). Specifically, to ask how the NR expression was distorted (altered expression, mutation and CNV) we have applied bootstrapping approaches to simulate data for comparison, and also compared these NR findings to 12 other transcription factor families. Nuclear receptors were uniquely and uniformly downregulated across all six tumor types, more than predicted by chance. These approaches also revealed that each tumor type had a specific NR expression profile but these were most similar between breast and prostate cancer. Some NRs were down-regulated in at least five tumor types (e.g. NR3C2/MR and NR5A2/LRH-1)) whereas others were uniquely down-regulated in one tumor (e.g. NR1B3/RARG). The downregulation was not driven by copy number variation or mutation and epigenetic mechanisms maybe responsible for the altered nuclear receptor expression.

The ubiquitous transcription factor CTCF promotes lineage-specific epigenomic remodeling and establishment of transcriptional networks driving cell differentiation

Cell differentiation relies on tissue-specific transcription factors (TFs) that cooperate to establish unique transcriptomes and phenotypes. However, the role of ubiquitous TFs in these processes remains poorly defined. Recently, we have shown that the CCCTC-binding factor (CTCF) is required for adipocyte differentiation through epigenomic remodelling of adipose tissue-specific enhancers and transcriptional activation of Peroxisome proliferator-activated receptor gamma (PPARG), the main driver of the adipogenic program (PPARG), and its target genes. Here, we discuss how these findings, together with the recent literature, illuminate a functional role for ubiquitous TFs in lineage-determining transcriptional networks.

Super-enhancer lncs to cardiovascular development and disease

Cardiac development, function and pathological remodelling in response to stress depend on the dynamic control of tissue specific gene expression by distant acting transcriptional enhancers. Recently, super-enhancers (SEs), also known as stretch or large enhancer clusters, are emerging as sentinel regulators within the gene regulatory networks that underpin cellular functions. It is becoming increasingly evident that long noncoding RNAs (lncRNAs) associated with these sequences play fundamental roles for enhancer activity and the regulation of the gene programs hardwired by them. Here, we review this emerging landscape, focusing on the roles of SEs and their derived lncRNAs in cardiovascular development and disease. We propose that exploration of this genomic landscape could provide novel therapeutic targets and approaches for the amelioration of cardiovascular disease. Ultimately we envisage a future of ncRNA therapeutics targeting the SE landscape to alleviate cardiovascular disease. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Mechanisms Underlying the Selection and Function of Macrophage-Specific Enhancers

Macrophages populate every tissue of the body and play vital roles in homeostasis, pathogen elimination, and tissue healing. These cells possess the ability to adapt to a multitude of abruptly changing and complex environments. Furthermore, different populations of resident tissue macrophages each show their own defining gene signatures. The enhancer repertoire of these cells underlies both the cellular identity of a given subset of resident macrophage population and their ability to dynamically alter, in an efficient manner, their gene expression programs in response to internal and external signals. Notably, transcription is pervasive at active enhancers and enhancer RNAs, or eRNAs, are tightly correlated to regulated transcription of protein-coding genes. Furthermore, selection and establishment of enhancers is a dynamic and plastic process in which activation of intracellular signaling pathways by factors present in a macrophage's environment play a determining role. Here, we review recent studies providing insights into the distinct mechanisms that contribute to the selection and function of enhancers in macrophages and the relevance of studying these mechanisms to gain a better understanding of complex human diseases.

Integrated metabolic modelling reveals cell-type specific epigenetic control points of the macrophage metabolic network

BACKGROUND

The reconstruction of context-specific metabolic models from easily and reliably measurable features such as transcriptomics data will be increasingly important in research and medicine. Current reconstruction methods suffer from high computational effort and arbitrary threshold setting. Moreover, understanding the underlying epigenetic regulation might allow the identification of putative intervention points within metabolic networks. Genes under high regulatory load from multiple enhancers or super-enhancers are known key genes for disease and cell identity. However, their role in regulation of metabolism and their placement within the metabolic networks has not been studied.

METHODS

Here we present FASTCORMICS, a fast and robust workflow for the creation of high-quality metabolic models from transcriptomics data. FASTCORMICS is devoid of arbitrary parameter settings and due to its low computational demand allows cross-validation assays. Applying FASTCORMICS, we have generated models for 63 primary human cell types from microarray data, revealing significant differences in their metabolic networks.

RESULTS

To understand the cell type-specific regulation of the alternative metabolic pathways we built multiple models during differentiation of primary human monocytes to macrophages and performed ChIP-Seq experiments for histone H3 K27 acetylation (H3K27ac) to map the active enhancers in macrophages. Focusing on the metabolic genes under high regulatory load from multiple enhancers or super-enhancers, we found these genes to show the most cell type-restricted and abundant expression profiles within their respective pathways. Importantly, the high regulatory load genes are associated to reactions enriched for transport reactions and other pathway entry points, suggesting that they are critical regulatory control points for cell type-specific metabolism.

CONCLUSIONS

By integrating metabolic modelling and epigenomic analysis we have identified high regulatory load as a common feature of metabolic genes at pathway entry points such as transporters within the macrophage metabolic network. Analysis of these control points through further integration of metabolic and gene regulatory networks in various contexts could be beneficial in multiple fields from identification of disease intervention strategies to cellular reprogramming.

Cell type-selective disease-association of genes under high regulatory load

We previously showed that disease-linked metabolic genes are often under combinatorial regulation. Using the genome-wide ChIP-Seq binding profiles for 93 transcription factors in nine different cell lines, we show that genes under high regulatory load are significantly enriched for disease-association across cell types. We find that transcription factor load correlates with the enhancer load of the genes and thereby allows the identification of genes under high regulatory load by epigenomic mapping of active enhancers. Identification of the high enhancer load genes across 139 samples from 96 different cell and tissue types reveals a consistent enrichment for disease-associated genes in a cell type-selective manner. The underlying genes are not limited to super-enhancer genes and show several types of disease-association evidence beyond genetic variation (such as biomarkers). Interestingly, the high regulatory load genes are involved in more KEGG pathways than expected by chance, exhibit increased betweenness centrality in the interaction network of liver disease genes, and carry longer 3' UTRs with more microRNA (miRNA) binding sites than genes on average, suggesting a role as hubs integrating signals within regulatory networks. In summary, epigenetic mapping of active enhancers presents a promising and unbiased approach for identification of novel disease genes in a cell type-selective manner.

Modulating the Genomic Programming of Adipocytes

The ability to modify the transcriptional program in response to external signals provides a way for mammalian cells to alter their biological fate and properties, thereby adapting to changes in the environment. Adipocytes are excellent examples of differentiated cells that possess a striking transcriptional plasticity when exposed to physiological and metabolic stimuli. In our work, we have focused on understanding the processes responsible for modulating the genomic programming in response to different external signals. Thus, we have shown that browning of human adipocytes with rosiglitazone, an antidiabetic agonist of the key adipocyte transcription factor peroxisome proliferator-activated receptor γ (PPARγ), involves redistribution of PPARγ binding to form browning-selective PPARγ super-enhancers that drive expression of key browning genes. These include genes encoding transcriptional regulators, such as Krüppel-like factor 11 (KLF11) that are essential for modulating the genomic program in white adipocytes to induce browning. Furthermore, we have shown that acute suppression of adipocyte genes by the proinflammatory cytokine, tumor necrosis factor (TNF), involves redistribution of cofactors to enhancers activated by the master inflammatory regulator, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB). Interestingly, this redistribution occurs selectively from enhancers with high-cofactor occupancies, thereby predominantly affecting super-enhancers and their associated genes. We propose that this is a general mechanism contributing to transcriptional repression associated with activation of signal-dependent transcription factors.

14-3-3ζ coordinates adipogenesis of visceral fat

The proteins that coordinate complex adipogenic transcriptional networks are poorly understood. 14-3-3ζ is a molecular adaptor protein that regulates insulin signalling and transcription factor networks. Here we report that 14-3-3ζ-knockout mice are strikingly lean from birth with specific reductions in visceral fat depots. Conversely, transgenic 14-3-3ζ overexpression potentiates obesity, without exacerbating metabolic complications. Only the 14-3-3ζ isoform is essential for adipogenesis based on isoform-specific RNAi. Mechanistic studies show that 14-3-3ζ depletion promotes autophagy-dependent degradation of C/EBP-δ, preventing induction of the master adipogenic factors, Pparγ and C/EBP-α. Transcriptomic data indicate that 14-3-3ζ acts upstream of hedgehog signalling-dependent upregulation of Cdkn1b/p27(Kip1). Indeed, concomitant knockdown of p27(Kip1) or Gli3 rescues the early block in adipogenesis induced by 14-3-3ζ knockdown in vitro. Adipocyte precursors in 14-3-3ζKO embryos also appear to have greater Gli3 and p27(Kip1) abundance. Together, our in vivo and in vitro findings demonstrate that 14-3-3ζ is a critical upstream driver of adipogenesis.

Global SUMOylation on active chromatin is an acute heat stress response restricting transcription

BACKGROUND

Cells have developed many ways to cope with external stress. One distinctive feature in acute proteotoxic stresses, such as heat shock (HS), is rapid post-translational modification of proteins by SUMOs (small ubiquitin-like modifier proteins; SUMOylation). While many of the SUMO targets are chromatin proteins, there is scarce information on chromatin binding of SUMOylated proteins in HS and the role of chromatin SUMOylation in the regulation of transcription.

RESULTS

We mapped HS-induced genome-wide changes in chromatin occupancy of SUMO-2/3-modified proteins in K562 and VCaP cells using ChIP-seq. Chromatin SUMOylation was further correlated with HS-induced global changes in transcription using GRO-seq and RNA polymerase II (Pol2) ChIP-seq along with ENCODE data for K562 cells. HS induced a rapid and massive rearrangement of chromatin SUMOylation pattern: SUMOylation was gained at active promoters and enhancers associated with multiple transcription factors, including heat shock factor 1. Concomitant loss of SUMOylation occurred at inactive intergenic chromatin regions that were associated with CTCF-cohesin complex and SETDB1 methyltransferase complex. In addition, HS triggered a dynamic chromatin binding of SUMO ligase PIAS1, especially onto promoters. The HS-induced SUMOylation on chromatin was most notable at promoters of transcribed genes where it positively correlated with active transcription and Pol2 promoter-proximal pausing. Furthermore, silencing of SUMOylation machinery either by depletion of UBC9 or PIAS1 enhanced expression of HS-induced genes.

CONCLUSIONS

HS-triggered SUMOylation targets promoters and enhancers of actively transcribed genes where it restricts the transcriptional activity of the HS-induced genes. PIAS1-mediated promoter SUMOylation is likely to regulate Pol2-associated factors in HS.

Increased H3K9me3 drives dedifferentiated phenotype via KLF6 repression in liposarcoma

Liposarcoma (LPS) can be divided into 4 different subtypes, of which well-differentiated LPS (WDLPS) and dedifferentiated LPS (DDLPS) are the most common. WDLPS is typically low grade, whereas DDLPS is high grade, aggressive, and carries a worse prognosis. WDLPS and DDLPS frequently co-occur in patients. However, it is not clear whether DDLPS arises independently from WDLPS, or whether epigenomic alterations underly the histopathological differences of these subtypes. Here, we profiled 9 epigenetic marks in tumor samples from 151 patients with LPS and showed elevated trimethylation of histone H3 at Lys9 (H3K9me3) levels in DDLPS tumors. Integrated ChIP-seq and gene expression analyses of patient-derived cell lines revealed that H3K9me3 mediates differential regulation of genes involved in cellular differentiation and migration. Among these, Kruppel-like factor 6 (KLF6) was reduced in DDLPS, with increased H3K9me3 at associated regulatory regions. Pharmacologic inhibition of H3K9me3 with chaetocin decreased DDLPS proliferation and increased expression of the adipogenesis-associated factors PPARγ, CEBPα, and CEBPβ, suggesting that increased H3K9me3 may mediate DDLPS-associated aggressiveness and dedifferentiation properties. KLF6 overexpression partially phenocopied chaetocin treatment in DDLPS cells and induced phenotypic changes that were consistent with adipocytic differentiation, suggesting that the effects of increased H3K9me3 may be mediated through KLF6. In conclusion, we provide evidence of an epigenetic basis for the transition between WDLPS and DDLPS.

Chromatin states modify network motifs contributing to cell-specific functions

Epigenetic modification can affect many important biological processes, such as cell proliferation and apoptosis. It can alter chromatin conformation and contribute to gene regulation. To investigate how chromatin states associated with network motifs, we assembled chromatin state-modified regulatory networks by combining 269 ChIP-seq data and chromatin states in four cell types. We found that many chromatin states were significantly associated with network motifs, especially for feedforward loops (FFLs). These distinct chromatin state compositions contribute to different expression levels and translational control of targets in FFLs. Strikingly, the chromatin state-modified FFLs were highly cell-specific and, to a large extent, determined cell-selective functions, such as the embryonic stem cell-specific bivalent modification-related FFL with an important role in poising developmentally important genes for expression. Besides, comparisons of chromatin state-modified FFLs between cancerous/stem and primary cell lines revealed specific type of chromatin state alterations that may act together with motif structural changes cooperatively contribute to cell-to-cell functional differences. Combination of these alterations could be helpful in prioritizing candidate genes. Together, this work highlights that a dynamic epigenetic dimension can help network motifs to control cell-specific functions.

Transcriptional and post-transcriptional control of adipocyte differentiation by Jumonji domain-containing protein 6

Jumonji domain-containing protein 6 (JMJD6) is a nuclear protein involved in histone modification, transcription and RNA processing. Although JMJD6 is crucial for tissue development, the link between its molecular functions and its roles in any given differentiation process is unknown. We report that JMJD6 is required for adipogenic gene expression and differentiation in a manner independent of Jumonji C domain catalytic activity. JMJD6 knockdown led to a reduction of C/EBPβ and C/EBPδ protein expression without affecting mRNA levels in the early phase of differentiation. However, ectopic expression of C/EBPβ and C/EBPδ did not rescue differentiation. Further analysis demonstrated that JMJD6 was associated with the Pparγ2 and Cebpα loci and putative enhancers. JMJD6 was previously found associated with bromodomain and extra-terminal domain (BET) proteins, which can be targeted by the bromodomain inhibitor JQ1. JQ1 treatment prevented chromatin binding of JMJD6, Pparγ2 and Cebpα expression, and adipogenic differentiation, yet had no effect on C/EBPβ and C/EBPδ expression or chromatin binding. These results indicate dual roles for JMJD6 in promoting adipogenic gene expression program by post-transcriptional regulation of C/EBPβ and C/EBPδ and direct transcriptional activation of Pparγ2 and Cebpα during adipocyte differentiation.

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.

Acute TNF-induced repression of cell identity genes is mediated by NFκB-directed redistribution of cofactors from super-enhancers

The proinflammatory cytokine tumor necrosis factor (TNF) plays a central role in low-grade adipose tissue inflammation and development of insulin resistance during obesity. In this context, nuclear factor κ-light-chain-enhancer of activated B cells (NFκB) is directly involved and required for the acute activation of the inflammatory gene program. Here, we show that the major transactivating subunit of NFκB, v-rel avian reticuloendotheliosis viral oncogene homolog A (RELA), is also required for acute TNF-induced suppression of adipocyte genes. Notably, this repression does not involve RELA binding to the associated enhancers but rather loss of cofactors and enhancer RNA (eRNA) selectively from high-occupancy sites within super-enhancers. Based on these data, we have developed models that, with high accuracy, predict which enhancers and genes are repressed by TNF in adipocytes. We show that these models are applicable to other cell types where TNF represses genes associated with super-enhancers in a highly cell-type–specific manner. Our results propose a novel paradigm for NFκB-mediated repression, whereby NFκB selectively redistributes cofactors from high-occupancy enhancers, thereby specifically repressing super-enhancer-associated cell identity genes.

"Hit-and-Run" leaves its mark: catalyst transcription factors and chromatin modification

Understanding how transcription factor (TF) binding is related to gene regulation is a moving target. We recently uncovered genome‐wide evidence for a “Hit‐and‐Run” model of transcription. In this model, a master TF “hits” a target promoter to initiate a rapid response to a signal. As the “hit” is transient, the model invokes recruitment of partner TFs to sustain transcription over time. Following the “run”, the master TF “hits” other targets to propagate the response genome‐wide. As such, a TF may act as a “catalyst” to mount a broad and acute response in cells that first sense the signal, while the recruited TF partners promote long‐term adaptive behavior in the whole organism. This “Hit‐and‐Run” model likely has broad relevance, as TF perturbation studies across eukaryotes show small overlaps between TF‐regulated and TF‐bound genes, implicating transient TF‐target binding. Here, we explore this “Hit‐and‐Run” model to suggest molecular mechanisms and its biological relevance.

The E3 ubiquitin ligase TRIM23 regulates adipocyte differentiation via stabilization of the adipogenic activator PPARγ

Adipocyte differentiation is a strictly controlled process regulated by a series of transcriptional activators. Adipogenic signals activate early adipogenic activators and facilitate the transient formation of early enhanceosomes at target genes. These enhancer regions are subsequently inherited by late enhanceosomes. PPARγ is one of the late adipogenic activators and is known as a master regulator of adipogenesis. However, the factors that regulate PPARγ expression remain to be elucidated. Here, we show that a novel ubiquitin E3 ligase, tripartite motif protein 23 (TRIM23), stabilizes PPARγ protein and mediates atypical polyubiquitin conjugation. TRIM23 knockdown caused a marked decrease in PPARγ protein abundance during preadipocyte differentiation, resulting in a severe defect in late adipogenic differentiation, whereas it did not affect the formation of early enhanceosomes. Our results suggest that TRIM23 plays a critical role in the switching from early to late adipogenic enhanceosomes by stabilizing PPARγ protein possibly via atypical polyubiquitin conjugation.

DOI:http://dx.doi.org/10.7554/eLife.05615.001

The world is facing a global epidemic of obesity, which also increases the risk for diabetes and heart disease. Obesity is caused when excess fat is stored in fat cells, and overweight individuals have larger fat cells compared to healthy weight people. Therefore understanding how fat cells are created in the body can provide new ways to combat obesity.

Fat cells, also known as adipocytes, arise from precursor cells via a process called adipogenesis. This requires the activity of proteins called transcription factors that bind to DNA and switch on the expression of genes. PPARγ is an important transcription factor that drives the expression of the genes that are needed to convert a precursor cell to a mature adipocyte.

For adipogenesis to proceed, cells have to maintain the appropriate levels of PPARγ. If the amount of PPARγ bound to DNA is too low, then it is unable to activate gene expression. However, the mechanisms by which cells maintain the correct levels of PPARγ activity remain poorly understood. Watanabe et al. analyzed this process in mouse cells and identified a protein called TRIM23 that is produced in precursor cells. Cells in which the levels of TRIM23 were artificially lowered failed to mature into fat cells; this suggests that this protein is necessary for adipogenesis. Furthermore, in the absence of TRIM23, the amount of PPARγ that occupied regions of DNA was also markedly reduced. A direct consequence of this was a decline in the expression of several genes that are required for the later steps in the adipogenesis process.

Watanabe et al. next analyzed the mechanism through which TRIM23 had an effect on the levels of PPARγ. It is known from previous work that TRIM23 belongs to a family of enzymes that attach a small molecular tag called ubiquitin onto other proteins. This ubiquitin tag typically marks these proteins for rapid destruction by a large molecular machine called the proteasome. Watanabe et al. found that TRIM23 also modified PPARγ with ubiquitin, but that it did so in an unusual manner that instead prevented the proteasome from recognizing PPARγ and destroying it. As such, TRIM23 stabilizes the levels of PPARγ in cells.

By providing new insights into how adipogenesis is regulated, these findings suggest that TRIM23 may be a potential therapeutic target in the treatment of diabetes and disorders related to obesity.

DOI:http://dx.doi.org/10.7554/eLife.05615.002

Convergence of developmental and oncogenic signaling pathways at transcriptional super-enhancers

Super-enhancers and stretch enhancers (SEs) drive expression of genes that play prominent roles in normal and disease cells, but the functional importance of these clustered enhancer elements is poorly understood, so it is not clear why genes key to cell identity have evolved regulation by such elements. Here, we show that SEs consist of functional constituent units that concentrate multiple developmental signaling pathways at key pluripotency genes in embryonic stem cells and confer enhanced responsiveness to signaling of their associated genes. Cancer cells frequently acquire SEs at genes that promote tumorigenesis, and we show that these genes are especially sensitive to perturbation of oncogenic signaling pathways. Super-enhancers thus provide a platform for signaling pathways to regulate genes that control cell identity during development and tumorigenesis.

Pioneer factors govern super-enhancer dynamics in stem cell plasticity and lineage choice

Adult stem cells (SCs) reside in niches which balance self-renewal with lineage selection and progression during tissue homeostasis. Following injury, culture or transplantation, SCs outside their niche often display fate flexibility^1-4. Here we show that super-enhancers⁵ underlie the identity, lineage commitment and plasticity of adult SCs in vivo. Using hair follicle (HF) as model, we map the global chromatin domains of HFSCs and their committed progenitors in their native microenvironments. We show that super-enhancers and their dense clusters (‘epicenters’) of transcription factor (TF) binding sites change upon lineage progression. New fate is acquired by decommissioning old and establishing new super-enhancers and/or epicenters, an auto-regulatory process that abates one master regulator subset while enhancing another. We further show that when outside their niche, either in vitro or in wound-repair, HFSCs dynamically remodel super-enhancers in response to changes in their microenvironment. Intriguingly, some key super-enhancers shift epicenters, enabling them to remain active and maintain a transitional state in an ever-changing transcriptional landscape. Finally, we identify SOX9 as a crucial chromatin rheostat of HFSC super-enhancers, and provide functional evidence that super-enhancers are dynamic, dense TF-binding platforms which are acutely sensitive to pioneer master regulators whose levels define not only spatial and temporal features of lineage-status, but also stemness, plasticity in transitional states and differentiation.

What are super-enhancers?

The term 'super-enhancer' has been used to describe groups of putative enhancers in close genomic proximity with unusually high levels of Mediator binding, as measured by chromatin immunoprecipitation and sequencing (ChIP-seq). Here we review the identification and composition of super-enhancers, describe links between super-enhancers, gene regulation and disease, and discuss the functional significance of enhancer clustering. We also provide our perspective regarding the proposition that super-enhancers are a regulatory entity conceptually distinct from what was known before the introduction of the term. Our opinion is that there is not yet strong evidence that super-enhancers are a novel paradigm in gene regulation and that use of the term in this context is not currently justified. However, the term likely identifies strong enhancers that exhibit behaviors consistent with previous models and concepts of transcriptional regulation. In this respect, the super-enhancer definition is useful in identifying regulatory elements likely to control genes important for cell type specification.

Enhancer activation requires trans-recruitment of a mega transcription factor complex

Enhancers provide critical information directing cell-type-specific transcriptional programs, regulated by binding of signal-dependent transcription factors and their associated cofactors. Here, we report that the most strongly activated estrogen (E2)-responsive enhancers are characterized by trans-recruitment and in situ assembly of a large 1-2 MDa complex of diverse DNA-binding transcription factors by ERα at ERE-containing enhancers. We refer to enhancers recruiting these factors as mega transcription factor-bound in trans (MegaTrans) enhancers. The MegaTrans complex is a signature of the most potent functional enhancers and is required for activation of enhancer RNA transcription and recruitment of coactivators, including p300 and Med1. The MegaTrans complex functions, in part, by recruiting specific enzymatic machinery, exemplified by DNA-dependent protein kinase. Thus, MegaTrans-containing enhancers represent a cohort of functional enhancers that mediate a broad and important transcriptional program and provide a molecular explanation for transcription factor clustering and hotspots noted in the genome.

The selection and function of cell type-specific enhancers

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

Inflammation-sensitive super enhancers form domains of coordinately regulated enhancer RNAs

Enhancers are critical genomic elements that define cellular and functional identity through the spatial and temporal regulation of gene expression. Recent studies suggest that key genes regulating cell type-specific functions reside in enhancer-dense genomic regions (i.e., super enhancers, stretch enhancers). Here we report that enhancer RNAs (eRNAs) identified by global nuclear run-on sequencing are extensively transcribed within super enhancers and are dynamically regulated in response to cellular signaling. Using Toll-like receptor 4 (TLR4) signaling in macrophages as a model system, we find that transcription of super enhancer-associated eRNAs is dynamically induced at most of the key genes driving innate immunity and inflammation. Unexpectedly, genes repressed by TLR4 signaling are also associated with super enhancer domains and accompanied by massive repression of eRNA transcription. Furthermore, we find each super enhancer acts as a single regulatory unit within which eRNA and genic transcripts are coordinately regulated. The key regulatory activity of these domains is further supported by the finding that super enhancer-associated transcription factor binding is twice as likely to be conserved between human and mouse than typical enhancer sites. Our study suggests that transcriptional activities at super enhancers are critical components to understand the dynamic gene regulatory network.

A hyper-dynamic nature of bivalent promoter states underlies coordinated developmental gene expression modules

Background

Chromatin remodeling is crucial for proper programing of developmental gene expression. Recent work provides a dynamic view of post-translational histone modifications during differentiation; however there is little insight on the evolution of combinatorial genome-wide patterns of chromatin marks, excluding an essential aspect of developmental gene regulation.

Results

We report here a 15-chromatin state Hidden Markov Model which describes changes in chromatin signatures in relation to transcription profiles during differentiation of human pre-adipocytes into adipocytes. We identify nineteen modules of gene expression reflecting multiple waves of transcriptional up- and down-regulation which characterize adipogenic differentiation. From our model, we developed chromatin state matrices fitting each of these transcription modules to show how the complexity and dynamic nature of chromatin signatures relate to expression patterns. Spatial relationships between chromatin states underlie a high-order chromatin organization in differentiating adipocytes. We show the importance of gene expression level in generating diversity in chromatin signatures, and show that the hyper-dynamic nature of H3K4me2/H3K27me3-marked ‘bivalent’ promoter states underlies many of the gene expression patterns associated with adipogenic differentiation.

Conclusions

Our results reveal the highly dynamic nature of bivalent promoter states within the adipogenic lineage. The data constitute a valuable resource enabling the assessment of possibilities to alter the adipogenic program.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1186) contains supplementary material, which is available to authorized users.

Oncogene regulation. An oncogenic super-enhancer formed through somatic mutation of a noncoding intergenic element

In certain human cancers, the expression of critical oncogenes is driven from large regulatory elements, called super-enhancers, that recruit much of the cell's transcriptional apparatus and are defined by extensive acetylation of histone H3 lysine 27 (H3K27ac). In a subset of T-cell acute lymphoblastic leukemia (T-ALL) cases, we found that heterozygous somatic mutations are acquired that introduce binding motifs for the MYB transcription factor in a precise noncoding site, which creates a super-enhancer upstream of the TAL1 oncogene. MYB binds to this new site and recruits its H3K27 acetylase-binding partner CBP, as well as core components of a major leukemogenic transcriptional complex that contains RUNX1, GATA-3, and TAL1 itself. Additionally, most endogenous super-enhancers found in T-ALL cells are occupied by MYB and CBP, which suggests a general role for MYB in super-enhancer initiation. Thus, this study identifies a genetic mechanism responsible for the generation of oncogenic super-enhancers in malignant cells.

Browning of human adipocytes requires KLF11 and reprogramming of PPARγ superenhancers

Long-term exposure to peroxisome proliferator-activated receptor γ (PPARγ) agonists such as rosiglitazone induces browning of rodent and human adipocytes; however, the transcriptional mechanisms governing this phenotypic switch in adipocytes are largely unknown. Here we show that rosiglitazone-induced browning of human adipocytes activates a comprehensive gene program that leads to increased mitochondrial oxidative capacity. Once induced, this gene program and oxidative capacity are maintained independently of rosiglitazone, suggesting that additional browning factors are activated. Browning triggers reprogramming of PPARγ binding, leading to the formation of PPARγ "superenhancers" that are selective for brown-in-white (brite) adipocytes. These are highly associated with key brite-selective genes. Based on such an association, we identified an evolutionarily conserved metabolic regulator, Kruppel-like factor 11 (KLF11), as a novel browning transcription factor in human adipocytes that is required for rosiglitazone-induced browning, including the increase in mitochondrial oxidative capacity. KLF11 is directly induced by PPARγ and appears to cooperate with PPARγ in a feed-forward manner to activate and maintain the brite-selective gene program.

Identification of the transcription factor ZEB1 as a central component of the adipogenic gene regulatory network

Adipose tissue is a key determinant of whole body metabolism and energy homeostasis. Unraveling the regulatory mechanisms underlying adipogenesis is therefore highly relevant from a biomedical perspective. Our current understanding of fat cell differentiation is centered on the transcriptional cascades driven by the C/EBP protein family and the master regulator PPARγ. To elucidate further components of the adipogenic gene regulatory network, we performed a large-scale transcription factor (TF) screen overexpressing 734 TFs in mouse pre-adipocytes and probed their effect on differentiation. We identified 22 novel pro-adipogenic TFs and characterized the top ranking TF, ZEB1, as being essential for adipogenesis both in vitro and in vivo. Moreover, its expression levels correlate with fat cell differentiation potential in humans. Genomic profiling further revealed that this TF directly targets and controls the expression of most early and late adipogenic regulators, identifying ZEB1 as a central transcriptional component of fat cell differentiation.

DOI:http://dx.doi.org/10.7554/eLife.03346.001

The growing rates of obesity and related metabolic diseases are a major public health concern worldwide. People who are overweight or obese are at increased risk for a range of diseases including diabetes and heart disease, which may reduce their quality of life and shorten their lifespans. Obese people have more, larger fat cells than individuals of healthy weight, and so understanding how the body creates fat cells may provide new insights into ways to prevent or treat obesity.

Fat cells arise from a population of stem cells that can also give rise to bone cells and cartilage. Some of the proteins—called transcription factors—that work together to turn on the expression of genes needed for a stem cell to become a fat cell have been identified. However, the exact regulatory processes that cause an unspecialized cell to develop into a fat cell remain unclear.

Gubelmann et al. set out to identify more of the transcription factors that cause stem cells to become fat cells. A high-throughput, automated process was used to test the effect of over-expressing each of 734 transcription factors in mouse cells that are primed to become fat cells. Twenty-six transcription factors were found to increase the number of these primed cells that became mature fat cells—most of which had not previously been shown to affect how fat cells develop. The most powerful driver of fat cell development was ZEB1: a transcription factor that has previously been implicated in many other biological processes. Most notably, ZEB1 was linked to a transition during development that allows cells to migrate to the correct location in the embryo, but also to a mechanism that allows cancerous cells to spread to new tissues.

Using studies of mouse cells and live mice, computational analyses, and biopsies from obese patients, Gubelmann et al. show that ZEB1 regulates numerous other transcription factors that promote the development of fat cells. These include factors that initially set an unspecialized cell on the path to becoming a fat cell and those that guide the changes that occur as the fat cell matures. Further studies will be required to show whether the ZEB1 protein itself is needed to prime stem cells to start becoming fat cells.

DOI:http://dx.doi.org/10.7554/eLife.03346.002

PPARγ and the global map of adipogenesis and beyond

Peroxisome proliferator-activated receptor γ (PPARγ) is a member of the nuclear receptor (NR) superfamily of ligand-dependent transcription factors (TFs) and function as a master regulator of adipocyte differentiation and metabolism. We review recent breakthroughs in the understanding of PPARγ gene regulation and function in the chromatin context. It is now clear that multiple TFs team up to induce PPARγ during adipogenesis, and that other TFs cooperate with PPARγ to ensure adipocyte-specific genomic binding and function. We discuss how this differs in other PPARγ-expressing cells such as macrophages and how these genome-wide mechanisms are preserved across species despite modest conservation of specific binding sites. These emerging considerations inform our understanding of PPARγ function as well as of adipocyte development and physiology.

Molecular architecture of transcription factor hotspots in early adipogenesis

Transcription factors have recently been shown to colocalize in hotspot regions of the genome, which are further clustered into super-enhancers. However, the detailed molecular organization of transcription factors at hotspot regions is poorly defined. Here, we have used digital genomic footprinting to precisely define factor localization at a genome-wide level during the early phase of 3T3-L1 adipocyte differentiation, which allows us to obtain detailed molecular insight into how transcription factors target hotspots. We demonstrate the formation of ATF-C/EBP heterodimers at a composite motif on chromatin, and we suggest that this may be a general mechanism for integrating external signals on chromatin. Furthermore, we find evidence of extensive recruitment of transcription factors to hotspots through alternative mechanisms not involving their known motifs and demonstrate that these alternative binding events are functionally important for hotspot formation and activity. Taken together, these findings provide a framework for understanding transcription factor cooperativity in hotspots.