Enhancer Function: Mechanistic and Genome-Wide Insights Come Together

Enhancer Function: Mechanistic and Genome-Wide Insights Come Together

Molecular Cell Review Enhancer Function: Mechanistic and Genome-Wide Insights Come Together Jennifer L. Plank1 and Ann Dean1,* 1Laboratory of Cellula...

1MB Sizes 0 Downloads 31 Views

Molecular Cell

Review Enhancer Function: Mechanistic and Genome-Wide Insights Come Together Jennifer L. Plank1 and Ann Dean1,* 1Laboratory of Cellular and Developmental Biology, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA *Correspondence: [email protected] http://dx.doi.org/10.1016/j.molcel.2014.06.015

Enhancers establish spatial or temporal patterns of gene expression that are critical for development, yet our understanding of how these DNA cis-regulatory elements function from a distance to increase transcription of their target genes and shape the cellular transcriptome has been gleaned primarily from studies of individual genes or gene families. High-throughput sequencing studies place enhancer-gene interactions within the 3D context of chromosome folding, inviting a new look at enhancer function and stimulating provocative new questions. Here, we integrate these whole-genome studies with recent mechanistic studies to illuminate how enhancers physically interact with target genes, how enhancer activity is regulated during development, and the role of noncoding RNAs transcribed from enhancers in their function. Introduction Enhancers are cis-acting DNA regulatory elements that increase the transcriptional output of target genes to influence the destiny of cells during development and differentiation. Enhancers may reside far from their in vivo targets, raising key mechanistic questions about how they communicate with promoters. Models such as enhancer looping to distant targets, linking by large protein complexes, or tracking along intervening chromatin were animatedly discussed but failed to be definitely resolved (Bulger and Groudine, 2011). The development of chromosome conformation capture (3C) technology (Dekker et al., 2002) allowed physical interaction frequencies between specific enhancers and target genes to be determined (Tolhuis et al., 2002). This significant advance affirmed that enhancers establish proximity with the genes they activate, although the original models are not mutually exclusive and may yet be found to contribute to enhancer function. Enhancer DNA sequences are replete with clusters of binding sites for transcription factors whose occupancy confers upon them tissue specificity. Thus, it had long been proposed that factors binding to enhancers and genes could stabilize chromatin loops between them through homotypic or heterotypic interaction. Indeed, tissue-specific proteins that are critical for looping have been identified in a limited number of mammalian model systems such as the b-globin locus and the human IFNg and TH2 cytokine loci, among others (Krivega and Dean, 2012). In addition, looping interactions of chromosome-organizing proteins, such as the insulator binding protein CTCF and its frequent partner cohesin, have direct and indirect roles in facilitating enhancer-gene chromatin contacts (Ong and Corces, 2014). Furthermore, lineage-specific activators and CTCF/cohesin engage in interactions with components of the RNA polymerase II (Pol II) machinery and with the Mediator transcriptional coactivator complex to tie enhancer loops directly to the transcription apparatus (Maston et al., 2012). Whole-genome studies now place enhancer interactions in a 3D nuclear context. These advances have depended on re-

finements in 3C-related approaches, including 5C, Hi-C, and ChIA-PET, capable of detecting chromatin loops at multiple levels (de Laat and Dekker, 2012). Studies across developmental stages detail the dynamic regulation of enhancers and enhancer looping during development (de Laat and Duboule, 2013). Moreover, recent results suggest a function for enhancer transcription into eRNAs as part of the mechanism of gene activation and possibly looping (Ørom and Shiekhattar, 2013). The new data, coming in a deluge of publications over a very recent period, paint a picture of multiple levels of long-range genome interactions, of which enhancer-gene contacts are a part. Descriptive by nature, the studies nevertheless allow mechanistic insights and raise new questions such as how chromosome folding occurs, what drives the changes in enhancers during development, and how enhancer looping and transcription activation are integrated physically and spatially in the nucleus to achieve a unique gene expression pattern. Excellent recent reviews have covered genome-wide discovery of new enhancers and gene targets, activator and coactivator recruitment to enhancers, and specific features of enhancers such as an open chromatin structure, the H3K4me1 histone mark, and histone acetyltransferase p300 occupancy (Spitz and Furlong, 2012; Maston et al., 2012; Calo and Wysocka, 2013). Here, we specifically focus on mechanistic insights into enhancer function from individual gene and very recent genome-wide studies. We begin by considering the basis for enhancer-gene loops in individual loci and across chromosomes. We then discuss how enhancer activity is regulated during development to modulate gene expression. We end with an account of our current understanding of enhancer RNAs and their roles in enhancer activity and gene regulation. Enhancers Engage in Long-Range Interactions with Target Genes through Lineage-Specific Factors and Ubiquitous Architectural Proteins Enhancers are typically occupied by clusters of transcription factors that exclude nucleosomes and contribute to their DNase I Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc. 5

Molecular Cell

Review A

B

C

D

Figure 1. Enhancers and Promoters Communicate by Chromatin Looping (A) Left: Two lineage-specific genes and an enhancer are depicted along unfolded chromatin with neither gene being transcribed. Right: Lineagespecific transcription factors mediate long-range interaction between the enhancer and one of the genes through homotypic and/or heterotypic protein interaction. The gene in contact with the enhancer is activated; the other gene (inactive) is looped away from the elements that are in proximity. (B) Left: A CTCF binding site and an enhancer are depicted with an inactive gene along unfolded chromatin. Right: The gene is activated by lineage-specific activators that co-opt CTCF into long-range interaction with the gene. (C) Left: A noninteracting enhancer and gene. Right: The enhancer is bridged to the gene promoter by Mediator and cohesin with participation of lineagespecific factors, activating the gene. (D) Left: A locus containing a gene and enhancer reside in an unfolded and inactive state. Center and right: Enhancer-gene looping is depicted as being mediated by lineage-specific activators before accumulation of Pol II and the appearance of a transcription factory and transcription.

hypersensitive character. At select loci, 3C and RNAi studies have shown that specific enhancer binding proteins are required for enhancer-gene looping (Figure 1A). For example, a complex including GATA1 and cofactor FOG1 along with TAL1, LMO2, and LDB1 is required for b-globin locus control region (LCR) looping to globin genes and for transcription activation in mature erythroid cells (Vakoc et al., 2005; Song et al., 2007; Yun et al., 2014). The dimerization domain of LDB1 underlies the enhancer-gene proximity (Krivega et al., 2014). In fact, the dimer6 Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc.

ization domain alone, when linked to the b-globin promoter via a specially designed DNA binding zinc finger protein, is capable of driving loop formation and partially activating transcription in immature erythroid cells where this does not normally occur, arguing for the causality of the loop in the activation (Deng et al., 2012). A large cohort of erythroid genes is regulated by the LDB1 complex, suggesting it can function broadly, likely through chromatin looping, to affect lineage commitment (Li et al., 2013a). Other focused investigations in the IFNg and MYB loci indicate that architectural proteins CTCF/cohesin can participate directly in enhancer-gene looping (Figure 1B) (Sekimata et al., 2009; Stadhouders et al., 2012; Hadjur et al., 2009). For example, enhancers that bind the lineage-specific factor T-BET in TH1 cells are interspersed with CTCF/cohesin binding elements in the IFNg locus (Sekimata et al., 2009). As naive T cells differentiate, CTCF promotes a TH1 cell-specific IFNg locus looped conformation joining both kinds of sites and activating IFNg transcription. CTCF occupancy and looping are dependent on enhancer binding by T-BET, but whether these proteins interact is unknown. These examples illustrate not only how lineage-specific transcription activators mediate enhancer-gene looping but also how they can cooperate with a ubiquitous looping factor (CTCF) to drive a cell-type-specific regional architecture conducive to transcription. Enhancer-gene interactions are dismantled during cell division. Interestingly, recent studies indicate that key lineage factors FOXA1 in hepatoma cells and GATA1 in erythroid cells remain associated with select sites on mitotic chromosomes, some of which have enhancer markings (Caravaca et al., 2013; Kadauke et al., 2012). Moreover, cohesin remains on chromatin during mitosis in colon cancer cells after eviction of clustered transcription factors with which it was associated at enhancerlike sites in interphase (Yan et al., 2013). We speculate that enhancers may be appropriate repositories for mitotic bookmarks to re-establish these long-range interactions pattern after cell division. Enhancer Loops Are Tied to the RNA Polymerase II Transcription Complex and Transcription Activation The enhancer’s main job is to increase transcriptional output. This activity could be manifest at different stages, including transcript initiation, elongation, or termination. Enhancer-promoter interactions can involve components of the basal transcription machinery (Ren et al., 2011; Liu et al., 2011; Koch et al., 2011). Moreover, Mediator occupies enhancers of many ES cell genes together with pluripotency factors OCT4, SOX2, and NANOG (Kagey et al., 2010) and bridges the enhancers to Pol II at target promoters by direct interaction with cohesin (Figure 1C). These studies intimately link enhancer loops to transcription initiation. Other studies suggest a role for enhancers in elongation (Sawado et al., 2003; Deng et al., 2012). New work now shows that certain enhancers specifically function to release Pol II pausing and allow elongation (Liu et al., 2013). These anti-pause enhancers loop to target promoters and permit activation of the P-TEFb complex, which is required for release of Pol II into elongation. Enhancers, thus, have diverse mechanistic functions in transcriptional regulation.

Molecular Cell

Review Further insight emerges from a ChIA-PET study that determined interactions centered on Pol II across the genome in human cell lines (Li et al., 2012). Networks of interacting promoters and promoters/enhancers were revealed. Enhancer-promoter interactions were significantly enriched for cell-type-specific genes. Most genes in multigene complexes were highly transcribed, and gene families regulated by a common transcription factor were significantly overrepresented. A link between these multigene complexes and Pol II foci was elegantly drawn by incorporating 3D DNA-fluorescence in situ hybridization (FISH) and immunostaining at select loci. New work in murine cells of progressive lineage commitment reveals that the gene clustering patterns are different, linking the looping interactions, coordinated transcription, and overall nuclear organization (Zhang et al., 2013). Spatial chromatin connectivity based on Pol II linking of select coregulated genes showed that several pluripotency genes in ES cells are connected within one major hub suggestive of a transcription factory. These studies strongly advocate for enhancer looping as key, and likely required, for transcription activation. This view is supported by data showing formation of new enhancer loops precedes transcription activation in the b-globin locus (Krivega et al., 2014) (Figure 1D). The Pol II chromatin connectivity data (Zhang et al., 2013) are concordant with earlier data showing enhancer looping and colocalization of erythroid genes regulated by a common transcription factor in a shared transcription factory (Schoenfelder et al., 2010). Whether transcription drives this colocalization is unresolved. It seems reasonable to propose that long-range enhancer-gene interactions collectively drive nuclear relocalization and clustering, underlying dynamic focal accumulation of Pol II and the appearance of transcription factories (Buckley and Lis, 2014). Interestingly, not all genes in such foci are active as examples exist where genes with opposite transcriptional responses to a stimulus or differentiation signal were found to colocalize (Hakim et al., 2011; Lin et al., 2012). Enhancer Loops within the Global 3D Organization of the Nucleus High-throughput studies now place enhancer-gene interactions within a genomic context. Thousands of enhancer-gene interactions associated with the unique transcriptional profile of cells have been observed in 5C studies (Sanyal et al., 2012; Thurman et al., 2012). Hi-C revealed that long-range interactions between genomic loci in ES cells and certain differentiated cell types occur predominantly within circumscribed topological domains or topologically associating domains (TADs) (Dixon et al., 2012; Nora et al., 2012). Megabase-sized TADs (0.5–3 Mb) are largely conserved across a range of cell types and during development. TAD borders are enriched for CTCF sites, although CTCF sites also occur within TADs. Long-range interactions within the same TAD are much more frequent than between TADs, and intradomain loops are more variable among cell types, suggesting that the intra-TAD contacts are involved in cell-specific transcription activation. TADs typically have a transcriptionally active or inactive character that can change during development or differentiation when gene activity changes (Denholtz et al., 2013; Lin et al., 2012).

TAD-type organizational structure has been widely observed across multiple mouse tissue and cell types (Shen et al., 2012). Clusters of coregulated enhancers and promoters were revealed that were defined as enhancer-promoter units (EPUs). Hi-C confirmed looping between the identified promoters and enhancers preferentially within and not between EPUs that overlap with defined TADs (Dixon et al., 2012). TADs also encompass ‘‘super enhancers’’ in ES cells (Whyte et al., 2013). These large domains of up to 50 kb contain clusters of individual enhancer elements that are highly occupied by Mediator and pluripotency factors. Enhancers within these regions are consistently associated with genes that encode important regulators of cell identity and interact with these targets with high frequency. Adding to the concept of TADs is the genome-wide description of sub-TADs that vary in a tissue-specific manner (PhillipsCremins et al., 2013). CTCF/cohesin and Mediator play roles at different levels of long-range interactions that contribute to TAD or sub-TAD organization, although a complete description of how these borders differ from one another is still missing. Sub-TAD organization is exemplified by certain HoxA genes and their enhancers that occupy the same TAD but are grouped into distinct topological subdomains in limb buds where the genes are active (Berlivet et al., 2013). 5C experiments suggest that the enhancers and genes interact with each other through contacts between sub-TADs possibly mediated by CTCF and cohesin, although this was not functionally tested. Supporting this idea, enhancer-gene long-range interactions are strengthened by, but do not depend on, enhancer activity. Are TADs functionally relevant for proper long-range enhancer activation of genes? In support, deletion of a TAD border in the Xist locus resulted in new ectopic contacts and long-range transcription misregulation (Nora et al., 2012). Moreover, CTCF (but not cohesin) reduction in HEK293T cells resulted in a gain of inter-TAD interactions, suggesting CTCF helps to maintain TAD borders, although changes in gene expression were modest (Zuin et al., 2014). Cohesin reduction in postmitotic thymocytes had no effect on TAD organization (Seitan et al., 2013), similar to the results of Zuin et al. 2014, but in astrocytes, relaxation of TAD structure and gene misregulation was observed upon cohesin reduction (Sofueva et al., 2013). Differing analysis methodologies may underlie this apparent discrepancy, but overall, the results highlight that much work remains to be done to define the contributions of CTCF, cohesin, Mediator, and other factors to maintenance of TAD boundaries and to understand their function. TAD borders are enriched for housekeeping genes, transfer RNA genes, and SINES (Dixon et al., 2012), so the field is wide open. Is an enhancer able to function outside its normal TAD context? An ectopic b-globin LCR reestablished contact with the b-globin gene and activated transcription, albeit in only a subset of cells (Noordermeer et al., 2011). Moreover, new work suggests that inter-TAD interactions, although less frequent than those within TADs, may be significant for enhancer function. ES cell loci enriched for super enhancers, typically encompassed by TADs (Whyte et al., 2013), colocalize with regions of similar character for transcription (Denholtz et al., 2013). Furthermore, in the HoxD cluster, situated near a TAD border, certain genes are observed to switch enhancer contacts from one TAD to another as gene expression changes across the locus Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc. 7

Molecular Cell

Review A

B

C

during limb development, indicating a certain fluidity of the border (Andrey et al., 2013). The picture coming into focus is of genome folding into conserved topological domains. Within TADs, enhancers cluster together with relevant gene targets. EPUs, super enhancers, and sub-TADs may all be manifestations of the same enhancerdependent clustering phenomenon. Although descriptive, this insight is mechanistically important, because it means enhancers and promoters need only scan a limited nuclear neighborhood to successfully and efficiently encounter each other (Gibcus and Dekker, 2013). Interruption of a contact site within a cluster of coregulated genes affected transcription of other interacting genes, supporting the causality of the contacts for the transcription (Fanucchi et al., 2013). This result is consistent with a model in which enhancer-gene clustering within TADs and association of TADs of similar transcriptional character nucleate transcription factories and argues against the possibility that transcription precipitates clustering. If true, this would emphasize the critical role of enhancers and looping interactions in the spatial organization of transcription within the nucleus (Denholtz et al., 2013). Temporal and Spatial Regulation of Enhancer Activity High-throughput studies across development show that enhancers are progressively modified to activate specific transcrip8 Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc.

Figure 2. Hallmarks of Enhancers at Different Developmental Stages (A) Left: Inactive ES cell enhancers have compact nucleosomes but are occupied by ELL3. Center: Poised enhancers are marked by H3K4me1, H3K27me3, PRC2, ELL3, BRG1, and p300. Right: At active enhancers, H3K27ac replaces H3K27me3 upon loss of PRC2. (B) Left: Latent macrophage enhancers lack transcription factor occupancy and enhancer-like histone modifications and contain compact nucleosomes. Center: Inactive macrophage enhancers differ from latent enhancers as they have acquired H3K4me1 and macrophage transcription factors such as PU.1. Right: Active enhancers are further marked by H3K27ac and p300 occupancy. (C) Center: In ES cells, neuroectoderm enhancers have unmethylated CpGs (white circles) and are occupied by NANOG, SOX2, and OCT4. Right: In the ectoderm, NANOG, SOX2, and OCT4 occupancy is diminished, the region remains unmethylated, and transcription of neural genes ensues. Left: Upon differentiation of ES cells to endoderm or mesoderm, NANOG, SOX2, and OCT4 occupancy is lost, and the region becomes more compact and highly methylated (black circles), inactivating the enhancers.

tional programs, particularly by acquisition of the H3K27ac mark (Rada-Iglesias et al., 2011; Creyghton et al., 2010). For example, ‘‘poised’’ developmental enhancers in ES cells are characterized by p300 and BRG1 occupancy, H3K4me1and low nucleosome density, together with Polycomb Repressive Complex 2 (PRC2) and the associated H3K27me3 mark (Figure 2A). At later developmental stages, these enhancers lose PRC2 and H3K27me3 and acquire H3K27ac and the ability to activate gene expression (Rada-Iglesias et al., 2011; Ferrari et al., 2014). New work now reveals an additional ‘‘inactive’’ state of ES cell enhancers with high nucleosome density but occupancy by the Pol II elongation factor ELL3 (Figure 2A) (Lin et al., 2013). These inactive enhancers must harbor ELL3 or they will be unable to activate target genes upon differentiation. Thus, ELL3 may mark enhancers for later activation of genes by Pol II promoter recruitment, which, interestingly, requires cohesin, consistent with looping. As ES cells differentiate, enhancers of pluripotency genes are inactivated to silence their target genes. Several mechanisms have been proposed. One study invoked the activity of the LSD1 histone H3K4/9 demethylase to remove the H3K4me1 enhancer mark (Whyte et al., 2012). In ES cells lacking LSD1, enhancers of pluripotency genes fail to undergo histone demethylation associated with differentiation, and their target genes are not repressed. Another study demonstrated that PRC2 deposits the H3K27me2 mark at enhancers in ES cells, thereby precluding the active H3K27ac mark and inappropriate enhancer activation (Ferrari et al., 2014). Differentiating cells also require the function of previously inactive enhancers. During the transition from naive to primed stem cells, the transcription factor OTX2 pioneers new enhancer sites allowing OCT4 to relocalize and activate a new

Molecular Cell

Review set of target genes that permit cells to exit the naive state (Buecker et al., 2014). Together, these data provide examples of both negative and positive enhancer regulation to allow orderly developmental progression. Recent studies have identified changes that occur genome wide within the regulatory landscape as cells differentiate and gene expression patterns change. Epigenetic profiling of the H3K27ac active enhancer mark in different developing mouse tissues identified approximately 90,000 enhancers that exhibited precise tissue- and stage-specific windows of predicted activity (Nord et al., 2013). A similar conclusion was reached when DNase I hypersensitive sites were used as the enhancer proxy (Stergachis et al., 2013). Enhancer usage varied widely among cell lineages in an additional study and was linked to recruitment of lineage-determining transcription factors (Kieffer-Kwon et al., 2013). Intriguingly, enhancer usage by commonly expressed genes such as Myc and Pim1 also varies among cell types, suggesting that such genes make use of different tissue specific enhancers and transcription factors as cells progress through development. In differentiated cells, enhancers can respond to external stimuli to activate appropriate genes. For example, TGFb signaling in pro-B cells or myotubes results in SMAD2/3 recruitment to silent enhancers occupied by PU.1 or MYOD, respectively, and target genes are activated (Mullen et al., 2011). A similar mechanism is employed in macrophages. Upon treatment with lipopolysaccharide endotoxin (LPS), ‘‘inactive’’ enhancers marked by H3K4me1 and tissue-specific transcription factors, including PU.1, are rapidly occupied by p300 and activate inflammatory response target genes (Ghisletti et al., 2010) (Figure 2B). Interestingly, a set of ‘‘latent’’ enhancers in macrophages are not sensitive to nucleases and lack both H3K4me1 and PU.1 (Figure 2B) (Ostuni et al., 2013). Upon exposure to LPS, latent enhancers acquire attributes of active enhancers and initiate gene transcription. Although transcription factors can mark inactive enhancers in unstimulated differentiated cells that will be required for later activation of target genes, it is not clear how latent enhancers can be recognized in order to be ‘‘unveiled.’’ Together, these studies reveal a dynamic enhancer landscape at distinct developmental stages linked to cell-specific gene expression, but how modifications ‘‘poise’’ or ‘‘activate’’ enhancers mechanistically and functionally is much less clear. It is tempting to propose that activation of genes by enhancers during development depends on establishment of loops to target promoters. However, in the HoxD locus, select enhancer-gene contacts are detectable before gene activation (Andrey et al., 2013). Moreover, both active and poised enhancers in IMR90 cells are equally likely to be engaged in looping before gene activation by signaling to TNF-a-responsive enhancers (Jin et al., 2013). Therefore, it will be important that future work focus on the precise temporal relationship of enhancer recognition, poising, and activation with enhancer looping and transcriptional activation and to determine the activities required to sequentially specify these outcomes. DNA Methylation as a Modulator of Enhancer Activity Consistent with the idea that DNA methylation interferes with transcription factor binding (Wiench et al., 2011), genome-wide

studies revealed that enhancers are negatively impacted by DNA methylation. A study of the murine methylome in ES and neural progenitor cells established a significant association of DNA hypomethylation at regions distal to promoter CpG islands with transcription factor binding, enhancer histone marks, and the activity of enhancers (Stadler et al., 2011). Subsequent studies in human cells affirmed that DNA hypomethylation at enhancers is associated with an open chromatin state, p300 occupancy, and target gene expression (Thurman et al., 2012; Varley et al., 2013). New studies looked broadly at the dynamics of methylation at regions outside of genes that are likely to be enhancers. While only a small fraction of autosomal CpGs are differentially methylated across numerous human cell and tissue types during development, 26% colocalize with enhancer-like regions (Ziller et al., 2013). Many insights can be gleaned from a large study showing lineage-specific gene silencing through asymmetric loss of enhancers by methylation (Gifford et al., 2013). For example, in human ES cells, enhancers destined to regulate genes required for neural specification are poised, unmethylated, marked by H3K4me1, and occupied by pluripotency factors (Figure 2C). Upon differentiation to neuroectoderm, these enhancers become activated, lose pluripotency factors, and (presumably) acquire neural lineage transcription factors. In contrast, when ES cells differentiate into mesoderm or endoderm, pluripotency factors are lost, and the enhancers become highly methylated. A consistent theme appears to be loss of DNA methylation and transition to H3K4me1 and either H3K27me3 or H3K27ac to poise or activate enhancers, respectively (Gifford et al., 2013; Sheaffer et al., 2014). In a different epigenetic analysis of human ES cells and their differentiated progeny, it was observed that DNA methylation is broadly present at enhancers at all stages and negatively associated with their activity (Xie et al., 2013). In some progeny, but not in others, early enhancers acquired methylation when they became inactive. A further study from this group introduced the term vestigial enhancers to describe early enhancers that retain a hypomethylated landscape in adult tissues, possibly as a memory mark. These enhancers exist in closed chromatin and are heavily marked by repressive H3K27me3, which, it is proposed, may be inhibitory to methylation of local DNA residues (Hon et al., 2013). These studies show that DNA hypomethylation is a common hallmark of active enhancers but highlight that hypomethylation is insufficient for activity of enhancers. Contrarily, recent work suggests the opposite, at least at some loci. Loss of LSH, a modulator of DNA methylation, results in decreased enhancer methylation, increased H3K4me1 occupancy, and increased target gene expression (Yu et al., 2014). More experimentation is needed to clarify this issue. How is DNA methylation at enhancers regulated? At promoters, TET1 has been implicated in regulating DNA demethylation through conversion of 5-methylcytosine to 5-hydroxymethylcytosine (5hmC) (Williams et al., 2012). Enrichment of 5hmC and TET1 occupancy is also associated with enhancers, suggesting that this mechanism might be conserved (Stroud et al., 2011; Pulakanti et al., 2013). Is enhancer DNA methylation the cause of inactivation or a step following other silencing mechanisms? Xie et al. interpreted their data to support the latter Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc. 9

Molecular Cell

Review A

eRNA chroma

looping

B

eRNA chroma

remodeling

Figure 3. Mechanism of eRNA Function (A) To activate target genes, eRNAs can interact with Mediator and cohesin to form long-range interactions between the enhancer and target genes. Loss of eRNAs or Mediator results in decreased target gene expression. (B) eRNAs can mediate target gene expression through influencing chromatin remodeling at target promoters. It is unclear if this mechanism occurs independently of chromatin looping.

?

RNA Pol II

gene

RNA

enhancer

Mediator

cohesin

possibility, consistent with earlier reports (Thurman et al., 2012; Xie et al., 2013; Stadler et al., 2011). Further studies of the dynamics of enhancer DNA methylation-demethylation and function during development may resolve this fundamental question. Regulated Production of Transcripts from Enhancers It had long been known that transcription occurs at enhancer elements, but it was widely assumed to be ancillary to transcription activation of target genes. Enhancer transcripts of various types are now known to occur on a genome scale and can be related to the transcription of select neighboring regulated genes (De Santa et al., 2010). Transcription of an eRNA upstream of the Arc gene requires the intact promoter and gene hinting that eRNA transcription and enhancer-promoter communication are linked (Kim et al., 2010). Several different kinds of stimuli can affect eRNA transcription. Androgen receptor, p53, MYOD, and MYOG positively regulate eRNAs by directly binding to enhancers (Wang et al., 2011; Melo et al., 2013; Mousavi et al., 2013). In contrast, REV-ERB signaling can negatively regulate target eRNAs, as evidenced by reduced eRNA transcription following Rev-Erba/Rev-Erbb overexpression (Lam et al., 2013). Enhancer DNA methylation may also play a role in eRNA transcription (Pulakanti et al., 2013). Demethylation of mature B cell enhancers accompanies their transcription into eRNAs and their activation (Schlesinger et al., 2013). In each of these examples, regulation of eRNA transcription was associated with the appropriate positive or 10 Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc.

negative effect on mRNA transcription of target genes supporting a mechanistic link. Mechanism by which eRNA May Function to Regulate Target Gene Expression High-throughput studies established a significant association between eRNAs, enhancer-promoter interactions, and target gene expression (Sanyal et al., 2012). Does the eRNA transcript per se have a role? Although the mechanism of nuclear targeting remains unclear, siRNAs and locked nucleic acids have been used to successfully deplete nuclear transcripts (Berteaux et al., 2008; Sarma et al., 2010). Using these approaches, when eRNAs transcribed from estrogen receptor-alpha (ERa)-sensitive enhancers were reduced in breast cancer cells, target gene expression was reduced (Li et al., 2013b). The eRNAs were associated with cohesin, and both looping to the genes and cohesin occupancy were reduced after eRNA inhibition, suggesting that eRNAs and cohesin stabilize the contacts (Figure 3A). eRNAs can also interact with the Mediator complex (Lai et al., 2013). Knockdown of MED1 and MED12 compromised eRNA activation of target genes and decreased loop formation between the enhancers and genes. Consistent with this, an eRNA transcribed from an androgen response element is required for MED1 occupancy, chromatin looping, and target gene expression (Hsieh et al., 2014). These studies suggest a direct role of eRNAs in formation or stabilization of loops between enhancers and target genes. In contrast, other studies questioned the role of eRNAs in looping. One study focusing on eRNAs transcribed at p53-dependent enhancers found that the transcribed enhancers physically interacted with multiple target genes whose expression was reduced after eRNA depletion (Melo et al., 2013). However, the loops were not diminished after p53 or eRNA reduction. Similarly, chemical inhibition of transcription elongation of eRNAs activated by ERa binding did not reduce looping (Hah et al., 2013). In these studies, enhancer-gene loops are maintained after reduction of eRNAs required for transcription of target genes, but they do not rule out a requirement for eRNAs to initially form loops or, alternatively, for another function related to target gene activation. In fact, work in which chromatin looping was not directly examined implicates eRNAs in chromatin remodeling at

Molecular Cell

Review target genes (Mousavi et al., 2013). Reduction of eRNAs regulating Myod1 or Myog expression results in decreased promoter Pol II occupancy and DNase I hypersensitivity, while eRNA overexpression suggests the transcript could function from an ectopic location (Figure 3B). Together, these data suggest that eRNA transcripts themselves may be functional and directly involved in target gene activation, at least in some cases. However, caution is required in drawing conclusions, and new approaches are needed to strengthen or refute current findings regarding eRNA functions. Moreover, underlying mechanisms of eRNA function on target gene transcription are still unclear. It is of considerable interest that eRNAs may function in chromatin looping via interaction with cohesin and/or the Mediator complex (Li et al., 2013b; Lai et al., 2013; Hsieh et al., 2014), and this issue warrants further study. Conclusions Genome-wide studies have changed the perspective in which we view enhancers and opened new possibilities. One remaining goal is to identify active enhancers in a high-throughput fashion. STARR-seq in Drosophila, provides a promising approach (Arnold et al., 2013). However, of primary importance is the ability to identify natural targets of enhancers, as nearby genes are frequently skipped in favor of more distant ones (Li et al., 2012; Sanyal et al., 2012; Zhang et al., 2013). Computational improvements in Hi-C and greater sequencing depth in ChIA-PET bring us closer to this goal by identifying enhancer-gene pairs with greater resolution (Jin et al., 2013; Zhang et al., 2013; KiefferKwon et al., 2013). Ultimately, however, candidate enhancers must be functionally tested. Genome editing with TALE nucleases (TALENs) has been deployed to delete or inactivate enhancers (Kieffer-Kwon et al., 2013; Mendenhall et al., 2013). TALEN or CRISPR/Cas9 (Jinek et al., 2013; Miller et al., 2011) methodologies appear destined to become the sine qua non to define enhancer activity in an endogenous environment. However, these approaches are still low-throughput, and considering that one gene may have more than one enhancer, and vice versa, they may be painstaking. Do enhancers control gene expression by looping? Single-cell studies capable of visualizing proximity of genes and enhancers at high resolution will be critical to reveal the repertoire and dynamics of enhancer-gene looping and its relationship to transcription. Hi-C can now be performed on single cells revealing considerable variation in the interactome among cells of the same type (Nagano et al., 2013). Approaches like this may allow us to understand the relationship between biochemically determined enhancer-gene loops and proximity determined by DNAFISH, which are not always concordant (Belmont, 2014). How do the patterns of enhancer-gene proximity relate to the firing of nascent transcripts? What is the relationship of physical clustering of genes in transcription factories observed by microscopy and topological features of chromatin, including enhancer looping and clustering in TADs, gleaned from populations of cells by 3Crelated methodologies? These fundamental questions should become approachable—and in live cells (Ghamari et al., 2013). Finally, it is becoming widely appreciated that sequence variation or other disruption of distant enhancers is a critical aspect

of misregulation of gene expression in disease (Maurano et al., 2012). Both sequence and methylation polymorphisms in enhancers are likely to be important, but recent work indicates a better correlation with altered target gene expression of the latter compared to the former (Aran and Hellman, 2013). A more indepth understanding of enhancer function is needed to point the way to disease mechanisms and potential interventions. Improvements to existing technologies are emerging at a fast and furious pace. We can expect our understanding of enhancers to broaden and deepen rapidly. ACKNOWLEDGMENTS We apologize to colleagues whose work we were unable to mention or cite dues to space limitations. We thank Vittorio Sartorelli, Elissa Lei, Judith Kassis, and Ivan Krivega for critical comments on the manuscript. This work was supported by the Intramural Program, NIDDK, NIH. REFERENCES Andrey, G., Montavon, T., Mascrez, B., Gonzalez, F., Noordermeer, D., Leleu, M., Trono, D., Spitz, F., and Duboule, D. (2013). A switch between topological domains underlies HoxD genes collinearity in mouse limbs. Science 340, 1234167. Aran, D., and Hellman, A. (2013). DNA methylation of transcriptional enhancers and cancer predisposition. Cell 154, 11–13. , L.M., Rath, M., and Stark, A. Arnold, C.D., Gerlach, D., Stelzer, C., Boryn (2013). Genome-wide quantitative enhancer activity maps identified by STARR-seq. Science 339, 1074–1077. Belmont, A.S. (2014). Large-scale chromatin organization: the good, the surprising, and the still perplexing. Curr. Opin. Cell Biol. 26, 69–78. Berlivet, S., Paquette, D., Dumouchel, A., Langlais, D., Dostie, J., and Kmita, M. (2013). Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs. PLoS Genet. 9, e1004018. Berteaux, N., Aptel, N., Cathala, G., Genton, C., Coll, J., Daccache, A., Spruyt, N., Hondermarck, H., Dugimont, T., Curgy, J.J., et al. (2008). A novel H19 antisense RNA overexpressed in breast cancer contributes to paternal IGF2 expression. Mol. Cell. Biol. 28, 6731–6745. Buckley, M.S., and Lis, J.T. (2014). Imaging RNA Polymerase II transcription sites in living cells. Curr. Opin. Genet. Dev. 25C, 126–130. Buecker, C., Srinivasan, R., Wu, Z., Calo, E., Acampora, D., Faial, T., Simeone, A., Tan, M., Swigut, T., and Wysocka, J. (2014). Reorganization of enhancer patterns in transition from naive to primed pluripotency. Cell Stem Cell 14, 838–853. Bulger, M., and Groudine, M. (2011). Functional and mechanistic diversity of distal transcription enhancers. Cell 144, 327–339. Calo, E., and Wysocka, J. (2013). Modification of enhancer chromatin: what, how, and why? Mol. Cell 49, 825–837. Caravaca, J.M., Donahue, G., Becker, J.S., He, X., Vinson, C., and Zaret, K.S. (2013). Bookmarking by specific and nonspecific binding of FoxA1 pioneer factor to mitotic chromosomes. Genes Dev. 27, 251–260. Creyghton, M.P., Cheng, A.W., Welstead, G.G., Kooistra, T., Carey, B.W., Steine, E.J., Hanna, J., Lodato, M.A., Frampton, G.M., Sharp, P.A., et al. (2010). Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc. Natl. Acad. Sci. USA 107, 21931–21936. de Laat, W., and Dekker, J. (2012). 3C-based technologies to study the shape of the genome. Methods 58, 189–191. de Laat, W., and Duboule, D. (2013). Topology of mammalian developmental enhancers and their regulatory landscapes. Nature 502, 499–506. De Santa, F., Barozzi, I., Mietton, F., Ghisletti, S., Polletti, S., Tusi, B.K., Muller, H., Ragoussis, J., Wei, C.L., and Natoli, G. (2010). A large fraction of extragenic RNA pol II transcription sites overlap enhancers. PLoS Biol. 8, e1000384.

Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc. 11

Molecular Cell

Review Dekker, J., Rippe, K., Dekker, M., and Kleckner, N. (2002). Capturing chromosome conformation. Science 295, 1306–1311. Deng, W., Lee, J., Wang, H., Miller, J., Reik, A., Gregory, P.D., Dean, A., and Blobel, G.A. (2012). Controlling long-range genomic interactions at a native locus by targeted tethering of a looping factor. Cell 149, 1233–1244. Denholtz, M., Bonora, G., Chronis, C., Splinter, E., de Laat, W., Ernst, J., Pellegrini, M., and Plath, K. (2013). Long-range chromatin contacts in embryonic stem cells reveal a role for pluripotency factors and polycomb proteins in genome organization. Cell Stem Cell 13, 602–616. Dixon, J.R., Selvaraj, S., Yue, F., Kim, A., Li, Y., Shen, Y., Hu, M., Liu, J.S., and Ren, B. (2012). Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature 485, 376–380. Fanucchi, S., Shibayama, Y., Burd, S., Weinberg, M.S., and Mhlanga, M.M. (2013). Chromosomal contact permits transcription between coregulated genes. Cell 155, 606–620.

Kieffer-Kwon, K.R., Tang, Z., Mathe, E., Qian, J., Sung, M.H., Li, G., Resch, W., Baek, S., Pruett, N., Grøntved, L., et al. (2013). Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation. Cell 155, 1507–1520. Kim, T.K., Hemberg, M., Gray, J.M., Costa, A.M., Bear, D.M., Wu, J., Harmin, D.A., Laptewicz, M., Barbara-Haley, K., Kuersten, S., et al. (2010). Widespread transcription at neuronal activity-regulated enhancers. Nature 465, 182–187. Koch, F., Fenouil, R., Gut, M., Cauchy, P., Albert, T.K., Zacarias-Cabeza, J., Spicuglia, S., de la Chapelle, A.L., Heidemann, M., Hintermair, C., et al. (2011). Transcription initiation platforms and GTF recruitment at tissue-specific enhancers and promoters. Nat. Struct. Mol. Biol. 18, 956–963. Krivega, I., and Dean, A. (2012). Enhancer and promoter interactions-long distance calls. Curr. Opin. Genet. Dev. 22, 79–85. Krivega, I., Dale, R.K., and Dean, A. (2014). Role of LDB1 in the transition from chromatin looping to transcription activation. Genes Dev. 28, 1278–1290.

Ferrari, K.J., Scelfo, A., Jammula, S., Cuomo, A., Barozzi, I., Stu¨tzer, A., Fischle, W., Bonaldi, T., and Pasini, D. (2014). Polycomb-dependent H3K27me1 and H3K27me2 regulate active transcription and enhancer fidelity. Mol. Cell 53, 49–62.

Lai, F., Orom, U.A., Cesaroni, M., Beringer, M., Taatjes, D.J., Blobel, G.A., and Shiekhattar, R. (2013). Activating RNAs associate with Mediator to enhance chromatin architecture and transcription. Nature 494, 497–501.

Ghamari, A., van de Corput, M.P., Thongjuea, S., van Cappellen, W.A., van Ijcken, W., van Haren, J., Soler, E., Eick, D., Lenhard, B., and Grosveld, F.G. (2013). In vivo live imaging of RNA polymerase II transcription factories in primary cells. Genes Dev. 27, 767–777.

Lam, M.T., Cho, H., Lesch, H.P., Gosselin, D., Heinz, S., Tanaka-Oishi, Y., Benner, C., Kaikkonen, M.U., Kim, A.S., Kosaka, M., et al. (2013). Rev-Erbs repress macrophage gene expression by inhibiting enhancer-directed transcription. Nature 498, 511–515.

Ghisletti, S., Barozzi, I., Mietton, F., Polletti, S., De Santa, F., Venturini, E., Gregory, L., Lonie, L., Chew, A., Wei, C.L., et al. (2010). Identification and characterization of enhancers controlling the inflammatory gene expression program in macrophages. Immunity 32, 317–328.

Li, G., Ruan, X., Auerbach, R.K., Sandhu, K.S., Zheng, M., Wang, P., Poh, H.M., Goh, Y., Lim, J., Zhang, J., et al. (2012). Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation. Cell 148, 84–98.

Gibcus, J.H., and Dekker, J. (2013). The hierarchy of the 3D genome. Mol. Cell 49, 773–782.

Li, L., Freudenberg, J., Cui, K., Dale, R., Song, S.H., Dean, A., Zhao, K., Jothi, R., and Love, P.E. (2013a). Ldb1-nucleated transcription complexes function as primary mediators of global erythroid gene activation. Blood 121, 4575– 4585.

Gifford, C.A., Ziller, M.J., Gu, H., Trapnell, C., Donaghey, J., Tsankov, A., Shalek, A.K., Kelley, D.R., Shishkin, A.A., Issner, R., et al. (2013). Transcriptional and epigenetic dynamics during specification of human embryonic stem cells. Cell 153, 1149–1163. Hadjur, S., Williams, L.M., Ryan, N.K., Cobb, B.S., Sexton, T., Fraser, P., Fisher, A.G., and Merkenschlager, M. (2009). Cohesins form chromosomal cis-interactions at the developmentally regulated IFNG locus. Nature 460, 410–413. Hah, N., Murakami, S., Nagari, A., Danko, C.G., and Kraus, W.L. (2013). Enhancer transcripts mark active estrogen receptor binding sites. Genome Res. 23, 1210–1223. Hakim, O., Sung, M.H., Voss, T.C., Splinter, E., John, S., Sabo, P.J., Thurman, R.E., Stamatoyannopoulos, J.A., de Laat, W., and Hager, G.L. (2011). Diverse gene reprogramming events occur in the same spatial clusters of distal regulatory elements. Genome Res. 21, 697–706. Hon, G.C., Rajagopal, N., Shen, Y., McCleary, D.F., Yue, F., Dang, M.D., and Ren, B. (2013). Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues. Nat. Genet. 45, 1198–1206. Hsieh, C.L., Fei, T., Chen, Y., Li, T., Gao, Y., Wang, X., Sun, T., Sweeney, C.J., Lee, G.S., Chen, S., et al. (2014). Enhancer RNAs participate in androgen receptor-driven looping that selectively enhances gene activation. Proc. Natl. Acad. Sci. USA 111, 7319–7324. Jin, F., Li, Y., Dixon, J.R., Selvaraj, S., Ye, Z., Lee, A.Y., Yen, C.A., Schmitt, A.D., Espinoza, C.A., and Ren, B. (2013). A high-resolution map of the threedimensional chromatin interactome in human cells. Nature 503, 290–294. Jinek, M., East, A., Cheng, A., Lin, S., Ma, E., and Doudna, J. (2013). RNA-programmed genome editing in human cells. Elife 2, e00471. Kadauke, S., Udugama, M.I., Pawlicki, J.M., Achtman, J.C., Jain, D.P., Cheng, Y., Hardison, R.C., and Blobel, G.A. (2012). Tissue-specific mitotic bookmarking by hematopoietic transcription factor GATA1. Cell 150, 725–737. Kagey, M.H., Newman, J.J., Bilodeau, S., Zhan, Y., Orlando, D.A., van Berkum, N.L., Ebmeier, C.C., Goossens, J., Rahl, P.B., Levine, S.S., et al. (2010). Mediator and cohesin connect gene expression and chromatin architecture. Nature 467, 430–435.

12 Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc.

Li, W., Notani, D., Ma, Q., Tanasa, B., Nunez, E., Chen, A.Y., Merkurjev, D., Zhang, J., Ohgi, K., Song, X., et al. (2013b). Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation. Nature 498, 516–520. Lin, Y.C., Benner, C., Mansson, R., Heinz, S., Miyazaki, K., Miyazaki, M., Chandra, V., Bossen, C., Glass, C.K., and Murre, C. (2012). Global changes in the nuclear positioning of genes and intra- and interdomain genomic interactions that orchestrate B cell fate. Nat. Immunol. 13, 1196–1204. Lin, C., Garruss, A.S., Luo, Z., Guo, F., and Shilatifard, A. (2013). The RNA Pol II elongation factor Ell3 marks enhancers in ES cells and primes future gene activation. Cell 152, 144–156. Liu, Z., Scannell, D.R., Eisen, M.B., and Tjian, R. (2011). Control of embryonic stem cell lineage commitment by core promoter factor, TAF3. Cell 146, 720–731. Liu, W., Ma, Q., Wong, K., Li, W., Ohgi, K., Zhang, J., Aggarwal, A.K., and Rosenfeld, M.G. (2013). Brd4 and JMJD6-associated anti-pause enhancers in regulation of transcriptional pause release. Cell 155, 1581–1595. Maston, G.A., Landt, S.G., Snyder, M., and Green, M.R. (2012). Characterization of enhancer function from genome-wide analyses. Annu. Rev. Genomics Hum. Genet. 13, 29–57. Maurano, M.T., Humbert, R., Rynes, E., Thurman, R.E., Haugen, E., Wang, H., Reynolds, A.P., Sandstrom, R., Qu, H., Brody, J., et al. (2012). Systematic localization of common disease-associated variation in regulatory DNA. Science 337, 1190–1195. Melo, C.A., Drost, J., Wijchers, P.J., van de Werken, H., de Wit, E., Oude Vrielink, J.A., Elkon, R., Melo, S.A., Le´veille´, N., Kalluri, R., de Laat, W., and Agami, R. (2013). eRNAs are required for p53-dependent enhancer activity and gene transcription. Mol. Cell 49, 524–535. Mendenhall, E.M., Williamson, K.E., Reyon, D., Zou, J.Y., Ram, O., Joung, J.K., and Bernstein, B.E. (2013). Locus-specific editing of histone modifications at endogenous enhancers. Nat. Biotechnol. 31, 1133–1136.

Molecular Cell

Review Miller, J.C., Tan, S., Qiao, G., Barlow, K.A., Wang, J., Xia, D.F., Meng, X., Paschon, D.E., Leung, E., Hinkley, S.J., et al. (2011). A TALE nuclease architecture for efficient genome editing. Nat. Biotechnol. 29, 143–148. Mousavi, K., Zare, H., Dell’orso, S., Grontved, L., Gutierrez-Cruz, G., Derfoul, A., Hager, G.L., and Sartorelli, V. (2013). eRNAs promote transcription by establishing chromatin accessibility at defined genomic loci. Mol. Cell 51, 606–617. Mullen, A.C., Orlando, D.A., Newman, J.J., Love´n, J., Kumar, R.M., Bilodeau, S., Reddy, J., Guenther, M.G., DeKoter, R.P., and Young, R.A. (2011). Master transcription factors determine cell-type-specific responses to TGFb signaling. Cell 147, 565–576. Nagano, T., Lubling, Y., Stevens, T.J., Schoenfelder, S., Yaffe, E., Dean, W., Laue, E.D., Tanay, A., and Fraser, P. (2013). Single-cell Hi-C reveals cell-tocell variability in chromosome structure. Nature 502, 59–64. Noordermeer, D., de Wit, E., Klous, P., van de Werken, H., Simonis, M., LopezJones, M., Eussen, B., de Klein, A., Singer, R.H., and de Laat, W. (2011). Variegated gene expression caused by cell-specific long-range DNA interactions. Nat. Cell Biol. 13, 944–951. Nora, E.P., Lajoie, B.R., Schulz, E.G., Giorgetti, L., Okamoto, I., Servant, N., Piolot, T., van Berkum, N.L., Meisig, J., Sedat, J., et al. (2012). Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature 485, 381–385. Nord, A.S., Blow, M.J., Attanasio, C., Akiyama, J.A., Holt, A., Hosseini, R., Phouanenavong, S., Plajzer-Frick, I., Shoukry, M., Afzal, V., et al. (2013). Rapid and pervasive changes in genome-wide enhancer usage during mammalian development. Cell 155, 1521–1531. Ong, C.T., and Corces, V.G. (2014). CTCF: an architectural protein bridging genome topology and function. Nat. Rev. Genet. 15, 234–246. Ørom, U.A., and Shiekhattar, R. (2013). Long noncoding RNAs usher in a new era in the biology of enhancers. Cell 154, 1190–1193. Ostuni, R., Piccolo, V., Barozzi, I., Polletti, S., Termanini, A., Bonifacio, S., Curina, A., Prosperini, E., Ghisletti, S., and Natoli, G. (2013). Latent enhancers activated by stimulation in differentiated cells. Cell 152, 157–171. Phillips-Cremins, J.E., Sauria, M.E., Sanyal, A., Gerasimova, T.I., Lajoie, B.R., Bell, J.S., Ong, C.T., Hookway, T.A., Guo, C., Sun, Y., et al. (2013). Architectural protein subclasses shape 3D organization of genomes during lineage commitment. Cell 153, 1281–1295.

Seitan, V.C., Faure, A.J., Zhan, Y., McCord, R.P., Lajoie, B.R., Ing-Simmons, E., Lenhard, B., Giorgetti, L., Heard, E., Fisher, A.G., et al. (2013). Cohesinbased chromatin interactions enable regulated gene expression within preexisting architectural compartments. Genome Res. 23, 2066–2077. Sekimata, M., Pe´rez-Melgosa, M., Miller, S.A., Weinmann, A.S., Sabo, P.J., Sandstrom, R., Dorschner, M.O., Stamatoyannopoulos, J.A., and Wilson, C.B. (2009). CCCTC-binding factor and the transcription factor T-bet orchestrate T helper 1 cell-specific structure and function at the interferon-gamma locus. Immunity 31, 551–564. Sheaffer, K.L., Kim, R., Aoki, R., Elliott, E.N., Schug, J., Burger, L., Schu¨beler, D., and Kaestner, K.H. (2014). DNA methylation is required for the control of stem cell differentiation in the small intestine. Genes Dev. 28, 652–664. Shen, Y., Yue, F., McCleary, D.F., Ye, Z., Edsall, L., Kuan, S., Wagner, U., Dixon, J., Lee, L., Lobanenkov, V.V., and Ren, B. (2012). A map of the cis-regulatory sequences in the mouse genome. Nature 488, 116–120. Sofueva, S., Yaffe, E., Chan, W.C., Georgopoulou, D., Vietri Rudan, M., MiraBontenbal, H., Pollard, S.M., Schroth, G.P., Tanay, A., and Hadjur, S. (2013). Cohesin-mediated interactions organize chromosomal domain architecture. EMBO J. 32, 3119–3129. Song, S.-H., Hou, C., and Dean, A. (2007). A positive role for NLI/Ldb1 in longrange b-globin locus control region function. Mol. Cell 28, 810–822. Spitz, F., and Furlong, E.E. (2012). Transcription factors: from enhancer binding to developmental control. Nat. Rev. Genet. 13, 613–626. Stadhouders, R., Thongjuea, S., Andrieu-Soler, C., Palstra, R.J., Bryne, J.C., van den Heuvel, A., Stevens, M., de Boer, E., Kockx, C., van der Sloot, A., et al. (2012). Dynamic long-range chromatin interactions control Myb protooncogene transcription during erythroid development. EMBO J. 31, 986–999. Stadler, M.B., Murr, R., Burger, L., Ivanek, R., Lienert, F., Scho¨ler, A., van Nimwegen, E., Wirbelauer, C., Oakeley, E.J., Gaidatzis, D., et al. (2011). DNA-binding factors shape the mouse methylome at distal regulatory regions. Nature 480, 490–495. Stergachis, A.B., Neph, S., Reynolds, A., Humbert, R., Miller, B., Paige, S.L., Vernot, B., Cheng, J.B., Thurman, R.E., Sandstrom, R., et al. (2013). Developmental fate and cellular maturity encoded in human regulatory DNA landscapes. Cell 154, 888–903. Stroud, H., Feng, S., Morey Kinney, S., Pradhan, S., and Jacobsen, S.E. (2011). 5-Hydroxymethylcytosine is associated with enhancers and gene bodies in human embryonic stem cells. Genome Biol. 12, R54.

Pulakanti, K., Pinello, L., Stelloh, C., Blinka, S., Allred, J., Milanovich, S., Kiblawi, S., Peterson, J., Wang, A., Yuan, G.C., and Rao, S. (2013). Enhancer transcribed RNAs arise from hypomethylated, Tet-occupied genomic regions. Epigenetics 8, 1303–1320.

Thurman, R.E., Rynes, E., Humbert, R., Vierstra, J., Maurano, M.T., Haugen, E., Sheffield, N.C., Stergachis, A.B., Wang, H., Vernot, B., et al. (2012). The accessible chromatin landscape of the human genome. Nature 489, 75–82.

Rada-Iglesias, A., Bajpai, R., Swigut, T., Brugmann, S.A., Flynn, R.A., and Wysocka, J. (2011). A unique chromatin signature uncovers early developmental enhancers in humans. Nature 470, 279–283.

Tolhuis, B., Palstra, R.J., Splinter, E., Grosveld, F., and de Laat, W. (2002). Looping and interaction between hypersensitive sites in the active b-globin locus. Mol. Cell 10, 1453–1465.

Ren, X., Siegel, R., Kim, U., and Roeder, R.G. (2011). Direct interactions of OCA-B and TFII-I regulate immunoglobulin heavy-chain gene transcription by facilitating enhancer-promoter communication. Mol. Cell 42, 342–355.

Vakoc, C.R., Letting, D.L., Gheldof, N., Sawado, T., Bender, M.A., Groudine, M., Weiss, M.J., Dekker, J., and Blobel, G.A. (2005). Proximity among distant regulatory elements at the b-globin locus requires GATA-1 and FOG-1. Mol. Cell 17, 453–462.

Sanyal, A., Lajoie, B.R., Jain, G., and Dekker, J. (2012). The long-range interaction landscape of gene promoters. Nature 489, 109–113. Sarma, K., Levasseur, P., Aristarkhov, A., and Lee, J.T. (2010). Locked nucleic acids (LNAs) reveal sequence requirements and kinetics of Xist RNA localization to the X chromosome. Proc. Natl. Acad. Sci. USA 107, 22196–22201.

Varley, K.E., Gertz, J., Bowling, K.M., Parker, S.L., Reddy, T.E., Pauli-Behn, F., Cross, M.K., Williams, B.A., Stamatoyannopoulos, J.A., Crawford, G.E., et al. (2013). Dynamic DNA methylation across diverse human cell lines and tissues. Genome Res. 23, 555–567.

Sawado, T., Halow, J., Bender, M.A., and Groudine, M. (2003). The b -globin locus control region (LCR) functions primarily by enhancing the transition from transcription initiation to elongation. Genes Dev. 17, 1009–1018.

Wang, D., Garcia-Bassets, I., Benner, C., Li, W., Su, X., Zhou, Y., Qiu, J., Liu, W., Kaikkonen, M.U., Ohgi, K.A., et al. (2011). Reprogramming transcription by distinct classes of enhancers functionally defined by eRNA. Nature 474, 390–394.

Schlesinger, F., Smith, A.D., Gingeras, T.R., Hannon, G.J., and Hodges, E. (2013). De novo DNA demethylation and noncoding transcription define active intergenic regulatory elements. Genome Res. 23, 1601–1614.

Whyte, W.A., Bilodeau, S., Orlando, D.A., Hoke, H.A., Frampton, G.M., Foster, C.T., Cowley, S.M., and Young, R.A. (2012). Enhancer decommissioning by LSD1 during embryonic stem cell differentiation. Nature 482, 221–225.

Schoenfelder, S., Sexton, T., Chakalova, L., Cope, N.F., Horton, A., Andrews, S., Kurukuti, S., Mitchell, J.A., Umlauf, D., Dimitrova, D.S., et al. (2010). Preferential associations between co-regulated genes reveal a transcriptional interactome in erythroid cells. Nat. Genet. 42, 53–61.

Whyte, W.A., Orlando, D.A., Hnisz, D., Abraham, B.J., Lin, C.Y., Kagey, M.H., Rahl, P.B., Lee, T.I., and Young, R.A. (2013). Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 153, 307–319.

Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc. 13

Molecular Cell

Review Wiench, M., John, S., Baek, S., Johnson, T.A., Sung, M.H., Escobar, T., Simmons, C.A., Pearce, K.H., Biddie, S.C., Sabo, P.J., et al. (2011). DNA methylation status predicts cell type-specific enhancer activity. EMBO J. 30, 3028–3039. Williams, K., Christensen, J., and Helin, K. (2012). DNA methylation: TET proteins-guardians of CpG islands? EMBO Rep. 13, 28–35. Xie, W., Schultz, M.D., Lister, R., Hou, Z., Rajagopal, N., Ray, P., Whitaker, J.W., Tian, S., Hawkins, R.D., Leung, D., et al. (2013). Epigenomic analysis of multilineage differentiation of human embryonic stem cells. Cell 153, 1134–1148. Yan, J., Enge, M., Whitington, T., Dave, K., Liu, J., Sur, I., Schmierer, B., Jolma, A., Kivioja, T., Taipale, M., and Taipale, J. (2013). Transcription factor binding in human cells occurs in dense clusters formed around cohesin anchor sites. Cell 154, 801–813. Yu, W., Briones, V., Lister, R., McIntosh, C., Han, Y., Lee, E.Y., Ren, J., Terashima, M., Leighty, R.M., Ecker, J.R., and Muegge, K. (2014). CG hypomethylation in Lsh-/- mouse embryonic fibroblasts is associated with de novo

14 Molecular Cell 55, July 3, 2014 ª2014 Elsevier Inc.

H3K4me1 formation and altered cellular plasticity. Proc. Natl. Acad. Sci. USA 111, 5890–5895. Yun, W.J., Kim, Y.W., Kang, Y., Lee, J., Dean, A., and Kim, A. (2014). The hematopoietic regulator TAL1 is required for chromatin looping between the b-globin LCR and human g-globin genes to activate transcription. Nucleic Acids Res. 42, 4283–4293. Zhang, Y., Wong, C.H., Birnbaum, R.Y., Li, G., Favaro, R., Ngan, C.Y., Lim, J., Tai, E., Poh, H.M., Wong, E., et al. (2013). Chromatin connectivity maps reveal dynamic promoter-enhancer long-range associations. Nature 504, 306–310. Ziller, M.J., Gu, H., Mu¨ller, F., Donaghey, J., Tsai, L.T., Kohlbacher, O., De Jager, P.L., Rosen, E.D., Bennett, D.A., Bernstein, B.E., et al. (2013). Charting a dynamic DNA methylation landscape of the human genome. Nature 500, 477–481. Zuin, J., Dixon, J.R., van der Reijden, M.I., Ye, Z., Kolovos, P., Brouwer, R.W., van de Corput, M.P., van de Werken, H.J., Knoch, T.A., van IJcken, W.F., et al. (2014). Cohesin and CTCF differentially affect chromatin architecture and gene expression in human cells. Proc. Natl. Acad. Sci. USA 111, 996–1001.