Enhancer transcripts mark active estrogen receptor binding sites

Functional annotation of structural ncRNAs within enhancer RNAs in the human genome: implications for human disease

Enhancer RNAs (eRNAs) are a novel class of non-coding RNA (ncRNA) molecules transcribed from the DNA sequences of enhancer regions. Despite extensive efforts devoted to revealing the potential functions and underlying mechanisms of eRNAs, it remains an open question whether eRNAs are mere transcriptional noise or relevant biologically functional species. Here, we identified a catalogue of eRNAs in a broad range of human cell/tissue types and extended our understanding of eRNAs by demonstrating their multi-omic signatures. Gene Ontology (GO) analysis revealed that eRNAs play key roles in human cell identity. Furthermore, we detected numerous known and novel functional RNA structures within eRNA regions. To better characterize the cis-regulatory effects of non-coding variation in these structural ncRNAs, we performed a comprehensive analysis of the genetic variants of structural ncRNAs in eRNA regions that are associated with inflammatory autoimmune diseases. Disease-associated variants of the structural ncRNAs were disproportionately enriched in immune-specific cell types. We also identified riboSNitches in lymphoid eRNAs and investigated the potential pathogenic mechanisms by which eRNAs might function in autoimmune diseases. Collectively, our findings offer valuable insights into the function of eRNAs and suggest that eRNAs might be effective diagnostic and therapeutic targets for human diseases.

Regulation of angiotensin II actions by enhancers and super-enhancers in vascular smooth muscle cells

Angiotensin II (AngII) promotes hypertension and atherosclerosis by activating growth-promoting and pro-inflammatory gene expression in vascular smooth muscle cells (VSMCs). Enhancers and super-enhancers (SEs) play critical roles in driving disease-associated gene expression. However, enhancers/SEs mediating VSMC dysfunction remain uncharacterized. Here, we show that AngII alters vascular enhancer and SE repertoires in cultured VSMCs in vitro, ex vivo, and in AngII-infused mice aortas in vivo. AngII-induced enhancers/SEs are enriched in binding sites for signal-dependent transcription factors and dependent on key signaling kinases. Moreover, CRISPR-Cas9-mediated deletion of candidate enhancers/SEs, targeting SEs with the bromodomain and extra-terminal domain inhibitor JQ1, or knockdown of overlapping long noncoding RNAs (lncRNAs) blocks AngII-induced genes associated with growth-factor signaling and atherosclerosis. Furthermore, JQ1 ameliorates AngII-induced hypertension, medial hypertrophy and inflammation in vivo in mice. These results demonstrate AngII-induced signals integrate enhancers/SEs and lncRNAs to increase expression of genes involved in VSMC dysfunction, and could uncover novel therapies.

DNA methylation at enhancers identifies distinct breast cancer lineages

Breast cancers exhibit genome-wide aberrant DNA methylation patterns. To investigate how these affect the transcriptome and which changes are linked to transformation or progression, we apply genome-wide expression-methylation quantitative trait loci (emQTL) analysis between DNA methylation and gene expression. On a whole genome scale, in cis and in trans, DNA methylation and gene expression have remarkably and reproducibly conserved patterns of association in three breast cancer cohorts (n = 104, n = 253 and n = 277). The expression-methylation quantitative trait loci associations form two main clusters; one relates to tumor infiltrating immune cell signatures and the other to estrogen receptor signaling. In the estrogen related cluster, using ChromHMM segmentation and transcription factor chromatin immunoprecipitation sequencing data, we identify transcriptional networks regulated in a cell lineage-specific manner by DNA methylation at enhancers. These networks are strongly dominated by ERα, FOXA1 or GATA3 and their targets were functionally validated using knockdown by small interfering RNA or GRO-seq analysis after transcriptional stimulation with estrogen.

Super enhancers - new analyses and perspectives on the low hanging fruit

Significant attention has recently been given to a class of enhancers termed "super enhancers", while implying that "typical enhancers" are less important. In this report, we examine criteria for identification of super enhancers and address the need to evaluate the differences between BRD4-occupied "typical" and "super" enhancers.

Nascent RNA sequencing reveals a dynamic global transcriptional response at genes and enhancers to the natural medicinal compound celastrol

Most studies of responses to transcriptional stimuli measure changes in cellular mRNA concentrations. By sequencing nascent RNA instead, it is possible to detect changes in transcription in minutes rather than hours and thereby distinguish primary from secondary responses to regulatory signals. Here, we describe the use of PRO-seq to characterize the immediate transcriptional response in human cells to celastrol, a compound derived from traditional Chinese medicine that has potent anti-inflammatory, tumor-inhibitory, and obesity-controlling effects. Celastrol is known to elicit a cellular stress response resembling the response to heat shock, but the transcriptional basis of this response remains unclear. Our analysis of PRO-seq data for K562 cells reveals dramatic transcriptional effects soon after celastrol treatment at a broad collection of both coding and noncoding transcription units. This transcriptional response occurred in two major waves, one within 10 min, and a second 40-60 min after treatment. Transcriptional activity was generally repressed by celastrol, but one distinct group of genes, enriched for roles in the heat shock response, displayed strong activation. Using a regression approach, we identified key transcription factors that appear to drive these transcriptional responses, including members of the E2F and RFX families. We also found sequence-based evidence that particular transcription factors drive the activation of enhancers. We observed increased polymerase pausing at both genes and enhancers, suggesting that pause release may be widely inhibited during the celastrol response. Our study demonstrates that a careful analysis of PRO-seq time-course data can disentangle key aspects of a complex transcriptional response, and it provides new insights into the activity of a powerful pharmacological agent.

Defining the location of promoter-associated R-loops at near-nucleotide resolution using bisDRIP-seq

R-loops are features of chromatin consisting of a strand of DNA hybridized to RNA, as well as the expelled complementary DNA strand. R-loops are enriched at promoters where they have recently been shown to have important roles in modifying gene expression. However, the location of promoter-associated R-loops and the genomic domains they perturb to modify gene expression remain unclear. To resolve this issue, we developed a bisulfite-based approach, bisDRIP-seq, to map R-loops across the genome at near-nucleotide resolution in MCF-7 cells. We found the location of promoter-associated R-loops is dependent on the presence of introns. In intron-containing genes, R-loops are bounded between the transcription start site and the first exon-intron junction. In intronless genes, the 3' boundary displays gene-specific heterogeneity. Moreover, intronless genes are often associated with promoter-associated R-loop formation. Together, these studies provide a high-resolution map of R-loops and identify gene structure as a critical determinant of R-loop formation.

Chromatin accessibility underlies synthetic lethality of SWI/SNF subunits in ARID1A-mutant cancers

ARID1A, a subunit of the SWI/SNF chromatin remodeling complex, is frequently mutated in cancer. Deficiency in its homolog ARID1B is synthetically lethal with ARID1A mutation. However, the functional relationship between these homologs has not been explored. Here, we use ATAC-seq, genome-wide histone modification mapping, and expression analysis to examine colorectal cancer cells lacking one or both ARID proteins. We find that ARID1A has a dominant role in maintaining chromatin accessibility at enhancers, while the contribution of ARID1B is evident only in the context of ARID1A mutation. Changes in accessibility are predictive of changes in expression and correlate with loss of H3K4me and H3K27ac marks, nucleosome spacing, and transcription factor binding, particularly at growth pathway genes including MET. We find that ARID1B knockdown in ARID1A mutant ovarian cancer cells causes similar loss of enhancer architecture, suggesting that this is a conserved function underlying the synthetic lethality between ARID1A and ARID1B.

DNA Sequence Constraints Define Functionally Active Steroid Nuclear Receptor Binding Sites in Chromatin

Gene regulatory programs are encoded in the sequence of the DNA. Since the completion of the Human Genome Project, millions of gene regulatory elements have been identified in the human genome. Understanding how each of those sites functionally contributes to gene regulation, however, remains a challenge for nearly every field of biology. Transcription factors influence cell function by interpreting information contained within cis-regulatory elements in chromatin. Whereas chromatin immunoprecipitation-sequencing has been used to identify and map transcription factor-DNA interactions, it has been difficult to assign functionality to the binding sites identified. Thus, in this study, we probed the transcriptional activity, DNA-binding competence, and functional activity of select nuclear receptor mutants in cellular and animal model systems and used this information to define the sequence constraints of functional steroid nuclear receptor cis-regulatory elements. Analysis of the architecture within sNR chromatin interacting sites revealed that only a small fraction of all sNR chromatin-interacting events is associated with transcriptional output and that this functionality is restricted to elements that vary from the consensus palindromic elements by one or two nucleotides. These findings define the transcriptional grammar necessary to predict functionality from regulatory sequences, with a multitude of future implications.

Dynamic assembly and activation of estrogen receptor α enhancers through coregulator switching

Although many features of active transcriptional enhancers have been defined by genomic assays, we lack a clear understanding of the order of events leading to enhancer formation and activation as well as the dynamics of coregulator interactions within the enhancer complex. Here, we used selective loss- or gain-of-function mutants of estrogen receptor α (ERα) to define two distinct phases of ligand-dependent enhancer formation. In the first phase (0-20 min), p300 is recruited to ERα by Mediator as well as p300's acetylhistone-binding bromodomain to promote initial enhancer formation, which is not competent for sustained activation. In the second phase (20-45 min), p300 is recruited to ERα by steroid receptor coregulators (SRCs) for enhancer maturation and maintenance. Successful transition between these two phases ("coregulator switching") is required for proper enhancer function. Failure to recruit p300 during either phase leads to abortive enhancer formation and a lack of target gene expression. Our results reveal an ordered and cooperative assembly of ERα enhancers requiring functional interplay among p300, Mediator, and SRCs, which has implications for hormone-dependent gene regulation in breast cancers. More broadly, our results demonstrate the unexpectedly dynamic nature of coregulator interactions within enhancer complexes, which are likely to be a defining feature of all enhancers.

Glucocorticoid Receptor:MegaTrans Switching Mediates the Repression of an ERα-Regulated Transcriptional Program

The molecular mechanisms underlying the opposing functions of glucocorticoid receptors (GRs) and estrogen receptor α (ERα) in breast cancer development remain poorly understood. Here we report that, in breast cancer cells, liganded GR represses a large ERα-activated transcriptional program by binding, in trans, to ERα-occupied enhancers. This abolishes effective activation of these enhancers and their cognate target genes, and it leads to the inhibition of ERα-dependent binding of components of the MegaTrans complex. Consistent with the effects of SUMOylation on other classes of nuclear receptors, dexamethasone (Dex)-induced trans-repression of the estrogen E₂ program appears to depend on GR SUMOylation, which leads to stable trans-recruitment of the GR-N-CoR/SMRT-HDAC3 corepressor complex on these enhancers. Together, these results uncover a mechanism by which competitive recruitment of DNA-binding nuclear receptors/transcription factors in trans to hot spot enhancers serves as an effective biological strategy for trans-repression, with clear implications for breast cancer and other diseases.

An Annotation Agnostic Algorithm for Detecting Nascent RNA Transcripts in GRO-Seq

We present a fast and simple algorithm to detect nascent RNA transcription in global nuclear run-on sequencing (GRO-seq). GRO-seq is a relatively new protocol that captures nascent transcripts from actively engaged polymerase, providing a direct read-out on bona fide transcription. Most traditional assays, such as RNA-seq, measure steady state RNA levels which are affected by transcription, post-transcriptional processing, and RNA stability. GRO-seq data, however, presents unique analysis challenges that are only beginning to be addressed. Here, we describe a new algorithm, Fast Read Stitcher (FStitch), that takes advantage of two popular machine-learning techniques, hidden Markov models and logistic regression, to classify which regions of the genome are transcribed. Given a small user-defined training set, our algorithm is accurate, robust to varying read depth, annotation agnostic, and fast. Analysis of GRO-seq data without a priori need for annotation uncovers surprising new insights into several aspects of the transcription process.

BRD4 promotes p63 and GRHL3 expression downstream of FOXO in mammary epithelial cells

Bromodomain-containing protein 4 (BRD4) is a member of the bromo- and extraterminal (BET) domain-containing family of epigenetic readers which is under intensive investigation as a target for anti-tumor therapy. BRD4 plays a central role in promoting the expression of select subsets of genes including many driven by oncogenic transcription factors and signaling pathways. However, the role of BRD4 and the effects of BET inhibitors in non-transformed cells remain mostly unclear. We demonstrate that BRD4 is required for the maintenance of a basal epithelial phenotype by regulating the expression of epithelial-specific genes including TP63 and Grainy Head-like transcription factor-3 (GRHL3) in non-transformed basal-like mammary epithelial cells. Moreover, BRD4 occupancy correlates with enhancer activity and enhancer RNA (eRNA) transcription. Motif analyses of cell context-specific BRD4-enriched regions predicted the involvement of FOXO transcription factors. Consistently, activation of FOXO1 function via inhibition of EGFR-AKT signaling promoted the expression of TP63 and GRHL3. Moreover, activation of Src kinase signaling and FOXO1 inhibition decreased the expression of FOXO/BRD4 target genes. Together, our findings support a function for BRD4 in promoting basal mammary cell epithelial differentiation, at least in part, by regulating FOXO factor function on enhancers to activate TP63 and GRHL3 expression.

Dissecting the genomic activity of a transcriptional regulator by the integrative analysis of omics data

In the study of genomic regulation, strategies to integrate the data produced by Next Generation Sequencing (NGS)-based technologies in a meaningful ensemble are eagerly awaited and must continuously evolve. Here, we describe an integrative strategy for the analysis of data generated by chromatin immunoprecipitation followed by NGS which combines algorithms for data overlap, normalization and epigenetic state analysis. The performance of our strategy is illustrated by presenting the analysis of data relative to the transcriptional regulator Estrogen Receptor alpha (ERα) in MCF-7 breast cancer cells and of Glucocorticoid Receptor (GR) in A549 lung cancer cells. We went through the definition of reference cistromes for different experimental contexts, the integration of data relative to co-regulators and the overlay of chromatin states as defined by epigenetic marks in MCF-7 cells. With our strategy, we identified novel features of estrogen-independent ERα activity, including FoxM1 interaction, eRNAs transcription and a peculiar ontology of connected genes.

Transcriptional response to stress is pre-wired by promoter and enhancer architecture

Programs of gene expression are executed by a battery of transcription factors that coordinate divergent transcription from a pair of tightly linked core initiation regions of promoters and enhancers. Here, to investigate how divergent transcription is reprogrammed upon stress, we measured nascent RNA synthesis at nucleotide-resolution, and profiled histone H4 acetylation in human cells. Our results globally show that the release of promoter-proximal paused RNA polymerase into elongation functions as a critical switch at which a gene's response to stress is determined. Highly transcribed and highly inducible genes display strong transcriptional directionality and selective assembly of general transcription factors on the core sense promoter. Heat-induced transcription at enhancers, instead, correlates with prior binding of cell-type, sequence-specific transcription factors. Activated Heat Shock Factor 1 (HSF1) binds to transcription-primed promoters and enhancers, and CTCF-occupied, non-transcribed chromatin. These results reveal chromatin architectural features that orient transcription at divergent regulatory elements and prime transcriptional responses genome-wide.Heat Shock Factor 1 (HSF1) is a regulator of stress-induced transcription. Here, the authors investigate changes to transcription and chromatin organization upon stress and find that activated HSF1 binds to transcription-primed promoters and enhancers, and to CTCF occupied, untranscribed chromatin.

The interplay of epigenetic marks during stem cell differentiation and development

Chromatin, the template for epigenetic regulation, is a highly dynamic entity that is constantly reshaped during early development and differentiation. Epigenetic modification of chromatin provides the necessary plasticity for cells to respond to environmental and positional cues, and enables the maintenance of acquired information without changing the DNA sequence. The mechanisms involve, among others, chemical modifications of chromatin, changes in chromatin constituents and reconfiguration of chromatin interactions and 3D structure. New advances in genome-wide technologies have paved the way towards an integrative view of epigenome dynamics during cell state transitions, and recent findings in embryonic stem cells highlight how the interplay between different epigenetic layers reshapes the transcriptional landscape.

Single-cell profiling reveals that eRNA accumulation at enhancer-promoter loops is not required to sustain transcription

Enhancers are intergenic DNA elements that regulate the transcription of target genes in response to signaling pathways by interacting with promoters over large genomic distances. Recent studies have revealed that enhancers are bi-directionally transcribed into enhancer RNAs (eRNAs). Using single-molecule fluorescence in situ hybridization (smFISH), we investigated the eRNA-mediated regulation of transcription during estrogen induction in MCF-7 cells. We demonstrate that eRNAs are localized exclusively in the nucleus and are induced with similar kinetics as target mRNAs. However, eRNAs are mostly nascent at enhancers and their steady-state levels remain lower than those of their cognate mRNAs. Surprisingly, at the single-allele level, eRNAs are rarely co-expressed with their target loci, demonstrating that active gene transcription does not require the continuous transcription of eRNAs or their accumulation at enhancers. When co-expressed, sub-diffraction distance measurements between nascent mRNA and eRNA signals reveal that co-transcription of eRNAs and mRNAs rarely occurs within closed enhancer–promoter loops. Lastly, basal eRNA transcription at enhancers, but not E2-induced transcription, is maintained upon depletion of MLL1 and ERα, suggesting some degree of chromatin accessibility prior to signal-dependent activation of transcription. Together, our findings suggest that eRNA accumulation at enhancer–promoter loops is not required to sustain target gene transcription.

Oncogenic long noncoding RNA landscape in breast cancer

Background

Few long noncoding RNAs (lncRNAs) that act as oncogenic genes in breast cancer have been identified.

Methods

Oncogenic lncRNAs associated with tumourigenesis and worse survival outcomes were examined and validated in Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA), respectively. Then, the potential biological functions and expression regulation of these lncRNAs were studied via bioinformatics and genome data analysis. Moreover, progressive breast cancer subtype-specific lncRNAs were investigated via high-throughput sequencing in our cohort and TCGA validation. To elucidate the mechanisms of the regulation of these lncRNAs, genomic alterations from the TCGA, Broad, Sanger and BCCRC data, as well as epigenetic modifications from GEO data, were then applied and examined to meet this objective. Finally, cell proliferation assays, flow cytometry analyses and TUNEL assays were applied to validate the oncogenic roles of these lncRNAs in vitro.

Results

A cluster of oncogenic lncRNAs that was upregulated in breast cancer tissue and was associated with worse survival outcomes was identified. These oncogenic lncRNAs are involved in regulating immune system activation and the TGF-beta and Jak-STAT signalling pathways. Moreover, TINCR, LINC00511, and PPP1R26-AS1 were identified as subtype-specific lncRNAs associated with HER-2, triple-negative and luminal B subtypes of breast cancer, respectively. The up-regulation of these oncogenic lncRNAs is mainly caused by gene amplification in the genome in breast cancer and other solid tumours. Finally, the knockdown of TINCR, DSCAM-AS1 or HOTAIR inhibited breast cancer cell proliferation, increased apoptosis and inhibited cell cycle progression in vitro.

Conclusions

These findings enhance the landscape of known oncogenic lncRNAs in breast cancer and provide insights into their roles. This understanding may potentially aid in the comprehensive management of breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-017-0696-6) contains supplementary material, which is available to authorized users.

Membrane and Nuclear Estrogen Receptor Alpha Actions: From Tissue Specificity to Medical Implications

Estrogen receptor alpha (ERα) has been recognized now for several decades as playing a key role in reproduction and exerting functions in numerous nonreproductive tissues. In this review, we attempt to summarize the in vitro studies that are the basis of our current understanding of the mechanisms of action of ERα as a nuclear receptor and the key roles played by its two activation functions (AFs) in its transcriptional activities. We then depict the consequences of the selective inactivation of these AFs in mouse models, focusing on the prominent roles played by ERα in the reproductive tract and in the vascular system. Evidence has accumulated over the two last decades that ERα is also associated with the plasma membrane and activates non-nuclear signaling from this site. These rapid/nongenomic/membrane-initiated steroid signals (MISS) have been characterized in a variety of cell lines, and in particular in endothelial cells. The development of selective pharmacological tools that specifically activate MISS and the generation of mice expressing an ERα protein impeded for membrane localization have begun to unravel the physiological role of MISS in vivo. Finally, we discuss novel perspectives for the design of tissue-selective ER modulators based on the integration of the physiological and pathophysiological roles of MISS actions of estrogens.

Non-coding Transcripts from Enhancers: New Insights into Enhancer Activity and Gene Expression Regulation

Long non-coding RNAs (lncRNAs) have gained widespread interest in the past decade owing to their enormous amount and surprising functions implicated in a variety of biological processes. Some lncRNAs exert function as enhancers, i.e., activating gene transcription by serving as the cis-regulatory molecules. Furthermore, recent studies have demonstrated that many enhancer elements can be transcribed and produce RNA molecules, which are termed as enhancer RNAs (eRNAs). The eRNAs are not merely the by-product of the enhancer transcription. In fact, many of them directly exert or regulate enhancer activity in gene activation through diverse mechanisms. Here, we provide an overview of enhancer activity, transcription of enhancer itself, characteristics of eRNAs, as well as their roles in regulating enhancer activity and gene expression.

The Diversity of Long Noncoding RNAs and Their Generation

Long noncoding RNAs (lncRNAs) are emerging as potential key regulators in gene expression networks and exhibit a surprising range of shapes and sizes. Several distinct classes of lncRNAs are transcribed from different DNA elements, including promoters, enhancers, and intergenic regions in eukaryotic genomes. Additionally, others are derived from long primary transcripts with noncanonical RNA processing pathways, generating new RNA species with unexpected formats. These lncRNAs can be processed by several mechanisms, including ribonuclease P (RNase P) cleavage to generate mature 3' ends, capping by small nucleolar RNA (snoRNA)-protein (snoRNP) complexes at their ends, or the formation of circular structures. Here we review current knowledge on lncRNAs and highlight the most recent discoveries of the underlying mechanisms related to their formation.

Chromatin loops and causality loops: the influence of RNA upon spatial nuclear architecture

An intrinsic and essential trait exhibited by cells is the properly coordinated and integrated regulation of an astoundingly large number of simultaneous molecular decisions and reactions to maintain biochemical homeostasis. This is especially true inside the cell nucleus, where the recognition of DNA and RNA by a vast range of nucleic acid-interacting proteins organizes gene expression patterns. However, this dynamic system is not regulated by simple "on" or "off" signals. Instead, transcription factor and RNA polymerase recruitment to DNA are influenced by the local chromatin and epigenetic environment, a gene's relative position within the nucleus and the action of noncoding RNAs. In addition, major phase-separated structural features of the nucleus, such as nucleoli and paraspeckles, assemble in direct response to specific transcriptional activities and, in turn, influence global genomic function. Currently, the interpretation of these data is trapped in a causality dilemma reminiscent of the "chicken and the egg" paradox as it is unclear whether changes in nuclear architecture promote RNA function or vice versa. Here, we review recent advances that suggest a complex and interdependent interaction network between gene expression, chromatin topology, and noncoding RNA function. We also discuss the functional links between these essential nuclear processes from the nanoscale (gene looping) to the macroscale (sub-nuclear gene positioning and nuclear body function) and briefly highlight some of the challenges that researchers may encounter when studying these phenomena.

SERPINB2 is regulated by dynamic interactions with pause-release proteins and enhancer RNAs

The SERPINB2 gene is strongly upregulated in inflammatory states. In monocytes, it can constitute up to 1% of total cellular protein. It functions in protection from proteotoxic stress and plays a role in angioedema. The purpose of this study was to define the roles of enhancer RNAs embedded in the SERPIN gene complex. We found that the upstream enhancer RNAs upregulated SERPINB2 and the enhancer RNAs were expressed prior to those of SERPINB2 mRNA. Studies of the SERPINB2 promoter demonstrated the presence of an RNA polymerase II pause-inducing protein, NELF. Stimulation with LPS led to recruitment of the pause-releasing kinase P-TEFb and departure of the pause-inducing protein NELF. RNA immunoprecipitation revealed that NELF and the CDK9 component of P-TEFb bound to the enhancer RNAs after stimulation with distinct kinetics. Knock-down of the enhancer RNAs compromised stimulus induction of promoter and enhancer chromatin changes. Conversely, over-expression was associated with enhanced recruitment of c-JUN and increased expression of SERPINB2 mRNA expression. This study is the first to associate enhancer RNAs with SERPINB2 and is the first demonstration of acquisition of NELF binding by enhancer RNAs on chromatin.

Recognition of long-range enhancer-promoter interactions by adding genomic signatures of segmented regulatory regions

Enhancer-promoter interaction (EPI) is an important cis-regulatory mechanism in the regulation of tissue-specific gene expression. However, it still has limitation to precisely identity these interactions so far. In this paper, using diverse genomic features for various regulatory regions, we presented a computational approach to predict EPIs with improved accuracies. Meanwhile, we comprehensively studied more potential regulatory factors that are important to EPIs prediction, such as nucleosome occupancy, enhancer RNA; and found the cell line-specificity and region-specificity of the contributions of diverse regulatory signatures. By adding genomic signatures of segmented regulatory regions, our best accuracies of cross-validation test were about 11%-16% higher than the previous results, indicating the location-specificity of genomic signatures in a regulatory region for predicting EPIs. Additionally, more training samples and related features can provide reliable performances in new cell lines. Consequently, our study provided additional insights into the roles of diverse signature features for predicting long-range EPIs.

A polymorphism in miR-1262 regulatory region confers the risk of lung cancer in Chinese population

It has been proposed that the majority of disease-associated loci identified by genome-wide association studies (GWAS) are enriched in non-coding regions, such as the promoter, enhancer or non-coding RNA genes. Thus, we performed a two-stage case-control study to systematically evaluate the association of genetic variants in miRNA regulatory regions (promoter and enhancer) with lung cancer risk in 7,763 subjects (discovery stage: 2,331 cases and 3,077 controls; validation stage: 1,065 cases and 1,290 controls). As a result, we identified that rs12740674 (C > T) in miR-1262 enhancer was significantly associated with the increased risk of lung cancer (additive model in discovery stage: adjusted OR = 1.31, 95%CI = 1.13-1.53, p = 3.846 × 10^-4 in Nanjing GWAS; adjusted OR = 1.20, 95%CI = 1.00-1.44, p = 0.041 in Beijing GWAS; validation stage: adjusted OR = 1.20, 95%CI = 1.03-1.41, p = 0.024). In meta-analysis, the p value for the association between rs12740674 and lung cancer risk reached 6.204 × 10^-6 (adjusted OR = 1.24, 95%CI = 1.13-1.36). Using 3DSNP database, The Cancer Genome Atlas (TCGA) data and functional assays, we observed that the risk T allele of rs12740674 reduced the expression level of miR-1262 in lung tissue through chromosomal looping, and overexpression of miR-1262 inhibited lung cancer cell proliferation probably through targeting the expression levels of ULK1 and RAB3D. Our findings confirmed the important role that genetic variants of noncoding sequence play in lung cancer susceptibility and indicated that rs12740674 in miR-1262 may be biologically relevant to lung carcinogenesis.

Enhancer-derived RNA: A Primer

Enhancer-derived RNAs (eRNAs) are a group of RNAs transcribed by RNA polymerase II from the domain of transcription enhancers, a major type of cis-regulatory elements in the genome. The correlation between eRNA production and enhancer activity has stimulated studies on the potential role of eRNAs in transcriptional regulation. Additionally, eRNA has also served as a marker for global identification of enhancers. Here I review the brief history and fascinating properties of eRNAs.

Genomic effects of glucocorticoids

Glucocorticoids and their receptor (GR) have been an important area of research because of their pleiotropic physiological functions and extensive use in the clinic. In addition, the association between GR and glucocorticoids, which is highly specific, leads to rapid nuclear translocation where GR associates with chromatin to regulate gene transcription. This simplified model system has been instrumental for studying the complexity of transcription regulation processes occurring at chromatin. In this review we discuss our current understanding of GR action that has been enhanced by recent developments in genome wide measurements of chromatin accessibility, histone marks, chromatin remodeling and 3D chromatin structure in various cell types responding to glucocorticoids.

Enhancers and chromatin structures: regulatory hubs in gene expression and diseases

Gene expression requires successful communication between enhancer and promoter regions, whose activities are regulated by a variety of factors and associated with distinct chromatin structures; in addition, functionally related genes and their regulatory repertoire tend to be arranged in the same subchromosomal regulatory domains. In this review, we discuss the importance of enhancers, especially clusters of enhancers (such as super-enhancers), as key regulatory hubs to integrate environmental cues and encode spatiotemporal instructions for genome expression, which are critical for a variety of biological processes governing mammalian development. Furthermore, we emphasize that the enhancer–promoter interaction landscape provides a critical context to understand the aetiologies and mechanisms behind numerous complex human diseases and provides new avenues for effective transcription-based interventions.

Multiplexed Engineering and Analysis of Combinatorial Enhancer Activity in Single Cells

The study of enhancers has been hampered by the scarcity of methods to systematically quantify their endogenous activity. We develop Mosaic-seq to systematically perturb enhancers and measure their endogenous activities at single-cell resolution. Mosaic-seq uses a CRISPR barcoding system to jointly measure a cell's transcriptome and its sgRNA modulators, thus quantifying the effects of dCas9-KRAB-mediated enhancer repression in single cells. Applying Mosaic-seq to 71 constituent enhancers from 15 super-enhancers, our analysis of 51,448 sgRNA-induced transcriptomes finds that only a small number of constituents are major effectors of target gene expression. Binding of p300 and RNAPII are key features of these constituents. We determine two key parameters of enhancer activity in single cells: their penetrance in a population and their contribution to expression in these cells. Through combinatorial interrogation, we find that simultaneous repression of multiple weak constituents can alter super-enhancer activity in a manner greatly exceeding repression of individual constituents.

Re-engineering the Pancreas Tumor Microenvironment: A "Regenerative Program" Hacked

The "hallmarks" of pancreatic ductal adenocarcinoma (PDAC) include proliferative, invasive, and metastatic tumor cells and an associated dense desmoplasia comprised of fibroblasts, pancreatic stellate cells, extracellular matrix, and immune cells. The oncogenically activated pancreatic epithelium and its associated stroma are obligatorily interdependent, with the resulting inflammatory and immunosuppressive microenvironment contributing greatly to the evolution and maintenance of PDAC. The peculiar pancreas-specific tumor phenotype is a consequence of oncogenes hacking the resident pancreas regenerative program, a tissue-specific repair mechanism regulated by discrete super enhancer networks. Defined as genomic regions containing clusters of multiple enhancers, super enhancers play pivotal roles in cell/tissue specification, identity, and maintenance. Hence, interfering with such super enhancer-driven repair networks should exert a disproportionately disruptive effect on tumor versus normal pancreatic tissue. Novel drugs that directly or indirectly inhibit processes regulating epigenetic status and integrity, including those driven by histone deacetylases, histone methyltransferase and hydroxylases, DNA methyltransferases, various metabolic enzymes, and bromodomain and extraterminal motif proteins, have shown the feasibility of disrupting super enhancer-dependent transcription in treating multiple tumor types, including PDAC. The idea that pancreatic adenocarcinomas rely on embedded super enhancer transcriptional mechanisms suggests a vulnerability that can be potentially targeted as novel therapies for this intractable disease. Clin Cancer Res; 23(7); 1647-55. ©2017 AACRSee all articles in this CCR Focus section, "Pancreatic Cancer: Challenge and Inspiration."

Nucleo-cytoplasmic transport of estrogen receptor alpha in breast cancer cells

Approximately 70% cases of breast cancers exhibit high expression and activity levels of estrogen receptor alpha (ERα), a transcription regulator that induces the expression of genes associated with cellular proliferation and survival. These nuclear functions of the receptor are associated with the development of breast cancer. However, ERα localization is not static, but rather, dynamic with continuous shuttling between the nucleus and the cytoplasm. Interestingly, both the nuclear import and export of ERα are modulated by several stimuli that include estradiol, antiestrogens, and growth factors. As ERα nuclear accumulation is critical to the regulation of gene expression, nuclear export of this receptor modulates the intensity and duration of its transcriptional activity. Thus, the subcellular spatial distribution of ERα ensures tight modulation of its concentration in cellular compartments, as well as of its nuclear and extranuclear functions. In this review, we will discuss current findings regarding the biological importance of molecular mechanisms of, and proteins responsible for, the nuclear import and export of ERα in breast cancer cells.

Transcriptional landscape of the human cell cycle

Steady-state gene expression across the cell cycle has been studied extensively. However, transcriptional gene regulation and the dynamics of histone modification at different cell-cycle stages are largely unknown. By applying a combination of global nuclear run-on sequencing (GRO-seq), RNA sequencing (RNA-seq), and histone-modification Chip sequencing (ChIP-seq), we depicted a comprehensive transcriptional landscape at the G0/G1, G1/S, and M phases of breast cancer MCF-7 cells. Importantly, GRO-seq and RNA-seq analysis identified different cell-cycle-regulated genes, suggesting a lag between transcription and steady-state expression during the cell cycle. Interestingly, we identified genes actively transcribed at early M phase that are longer in length and have low expression and are accompanied by a global increase in active histone 3 lysine 4 methylation (H3K4me2) and histone 3 lysine 27 acetylation (H3K27ac) modifications. In addition, we identified 2,440 cell-cycle-regulated enhancer RNAs (eRNAs) that are strongly associated with differential active transcription but not with stable expression levels across the cell cycle. Motif analysis of dynamic eRNAs predicted Kruppel-like factor 4 (KLF4) as a key regulator of G1/S transition, and this identification was validated experimentally. Taken together, our combined analysis characterized the transcriptional and histone-modification profile of the human cell cycle and identified dynamic transcriptional signatures across the cell cycle.

TT-seq captures enhancer landscapes immediately after T-cell stimulation

To monitor transcriptional regulation in human cells, rapid changes in enhancer and promoter activity must be captured with high sensitivity and temporal resolution. Here, we show that the recently established protocol TT-seq ("transient transcriptome sequencing") can monitor rapid changes in transcription from enhancers and promoters during the immediate response of T cells to ionomycin and phorbol 12-myristate 13-acetate (PMA). TT-seq maps eRNAs and mRNAs every 5 min after T-cell stimulation with high sensitivity and identifies many new primary response genes. TT-seq reveals that the synthesis of 1,601 eRNAs and 650 mRNAs changes significantly within only 15 min after stimulation, when standard RNA-seq does not detect differentially expressed genes. Transcription of enhancers that are primed for activation by nucleosome depletion can occur immediately and simultaneously with transcription of target gene promoters. Our results indicate that enhancer transcription is a good proxy for enhancer regulatory activity in target gene activation, and establish TT-seq as a tool for monitoring the dynamics of enhancer landscapes and transcription programs during cellular responses and differentiation.

Functional Enhancers As Master Regulators of Tissue-Specific Gene Regulation and Cancer Development

Tissue-specific transcription is critical for normal development, and abnormalities causing undesirable gene expression may lead to diseases such as cancer. Such highly organized transcription is controlled by enhancers with specific DNA sequences recognized by transcription factors. Enhancers are associated with chromatin modifications that are distinct epigenetic features in a tissue-specific manner. Recently, super-enhancers comprising enhancer clusters co-occupied by lineage-specific factors have been identified in diverse cell types such as adipocytes, hair follicle stem cells, and mammary epithelial cells. In addition, noncoding RNAs, named eRNAs, are synthesized at super-enhancer regions before their target genes are transcribed. Many functional studies revealed that super-enhancers and eRNAs are essential for the regulation of tissue-specific gene expression. In this review, we summarize recent findings concerning enhancer function in tissue-specific gene regulation and cancer development.

A Transcription Factor Pulse Can Prime Chromatin for Heritable Transcriptional Memory

Short-term and long-term transcriptional memory is the phenomenon whereby the kinetics or magnitude of gene induction is enhanced following a prior induction period. Short-term memory persists within one cell generation or in postmitotic cells, while long-term memory can survive multiple rounds of cell division. We have developed a tissue culture model to study the epigenetic basis for long-term transcriptional memory (LTTM) and subsequently used this model to better understand the epigenetic mechanisms that enable heritable memory of temporary stimuli. We find that a pulse of transcription factor CCAAT/enhancer-binding protein alpha (C/EBPα) induces LTTM on a subset of target genes that survives nine cell divisions. The chromatin landscape at genes that acquire LTTM is more repressed than at those genes that do not exhibit memory, akin to a latent state. We show through chromatin immunoprecipitation (ChIP) and chemical inhibitor studies that RNA polymerase II (Pol II) elongation is important for establishing memory in this model but that Pol II itself is not retained as part of the memory mechanism. More generally, our work reveals that a transcription factor involved in lineage specification can induce LTTM and that failure to rerepress chromatin is one epigenetic mechanism underlying transcriptional memory.

High fat diet-induced changes of mouse hepatic transcription and enhancer activity can be reversed by subsequent weight loss

Epigenetic factors have been suggested to play an important role in metabolic memory by trapping and maintaining initial metabolic changes within the transcriptional regulatory machinery. In this study we fed mice a high fat diet (HFD) for seven weeks followed by additional five weeks of chow, to identify HFD-mediated changes to the hepatic transcriptional program that may persist after weight loss. Mice fed a HFD displayed increased fasting insulin levels, hepatosteatosis and major changes in hepatic gene transcription associated with modulation of H3K27Ac at enhancers, but no significant changes in chromatin accessibility, indicating that HFD-regulated gene transcription is primarily controlled by modulating the activity of pre-established enhancers. After return to the same body weight as chow fed control mice, the fasting insulin, glucose, and hepatic triglyceride levels were fully restored to normal levels. Moreover, HFD-regulated H3K27Ac and mRNA levels returned to similar levels as control mice. These data demonstrates that the transcription regulatory landscape in the liver induced by HFD is highly dynamic and can be reversed by weight loss. This provides hope for efficient treatment of early obesity-associated changes to hepatic complications by simple weight loss intervention without persistent reprograming of the liver transcriptome.

Computational Approaches for Mining GRO-Seq Data to Identify and Characterize Active Enhancers

Transcriptional enhancers are DNA regulatory elements that are bound by transcription factors and act to positively regulate the expression of nearby or distally located target genes. Enhancers have many features that have been discovered using genomic analyses. Recent studies have shown that active enhancers recruit RNA polymerase II (Pol II) and are transcribed, producing enhancer RNAs (eRNAs). GRO-seq, a method for identifying the location and orientation of all actively transcribing RNA polymerases across the genome, is a powerful approach for monitoring nascent enhancer transcription. Furthermore, the unique pattern of enhancer transcription can be used to identify enhancers in the absence of any information about the underlying transcription factors. Here, we describe the computational approaches required to identify and analyze active enhancers using GRO-seq data, including data pre-processing, alignment, and transcript calling. In addition, we describe protocols and computational pipelines for mining GRO-seq data to identify active enhancers, as well as known transcription factor binding sites that are transcribed. Furthermore, we discuss approaches for integrating GRO-seq-based enhancer data with other genomic data, including target gene expression and function. Finally, we describe molecular biology assays that can be used to confirm and explore further the function of enhancers that have been identified using genomic assays. Together, these approaches should allow the user to identify and explore the features and biological functions of new cell type-specific enhancers.

GRO-seq, A Tool for Identification of Transcripts Regulating Gene Expression

The advent of next-generation sequencing (NGS) technologies has revolutionized the way we do research on gene expression. High-throughput transcriptomics became possible with the development of microarray technology, but its widespread application only occurred after the emergence of massive parallel sequencing. Especially, RNA sequencing (RNA-seq) has greatly increased our knowledge about the genome and led to the identification and annotation of novel classes of RNAs in different species. However, RNA-seq measures the steady-state level of a given RNA, which is the equilibrium between transcription, processing, and degradation. In recent years, a number of dedicated RNA-seq technologies were developed to measure specifically transcription events. Global run-on sequencing (GRO-seq) is the most widely used method to measure nascent RNA, and in recent years, it has been applied successfully to study the function and mechanism of action of noncoding RNAs. Here, we describe a detailed protocol of GRO-seq that can be readily applied to investigate different aspects of RNA biology in human cells.

Genome-wide mapping of autonomous promoter activity in human cells

Previous methods to systematically characterize sequence-intrinsic activity of promoters have been limited by relatively low throughput and the length of the sequences that could be tested. Here we present 'survey of regulatory elements' (SuRE), a method that assays more than 10⁸ DNA fragments, each 0.2-2 kb in size, for their ability to drive transcription autonomously. In SuRE, a plasmid library of random genomic fragments upstream of a 20-bp barcode is constructed, and decoded by paired-end sequencing. This library is used to transfect cells, and barcodes in transcribed RNA are quantified by high-throughput sequencing. When applied to the human genome, we achieve 55-fold genome coverage, allowing us to map autonomous promoter activity genome-wide in K562 cells. By computational modeling we delineate subregions within promoters that are relevant for their activity. We show that antisense promoter transcription is generally dependent on the sense core promoter sequences, and that most enhancers and several families of repetitive elements act as autonomous transcription initiation sites.

Searching for convergent pathways in autism spectrum disorders: insights from human brain transcriptome studies

Autism spectrum disorder (ASD) is one of the most heritable neuropsychiatric conditions. The complex genetic landscape of the disorder includes both common and rare variants at hundreds of genetic loci. This marked heterogeneity has thus far hampered efforts to develop genetic diagnostic panels and targeted pharmacological therapies. Here, we give an overview of the current literature on the genetic basis of ASD, and review recent human brain transcriptome studies and their role in identifying convergent pathways downstream of the heterogeneous genetic variants. We also discuss emerging evidence on the involvement of non-coding genomic regions and non-coding RNAs in ASD.

Regulatory Enhancer-Core-Promoter Communication via Transcription Factors and Cofactors

Gene expression is regulated by genomic enhancers that recruit transcription factors and cofactors to activate transcription from target core promoters. Over the past years, thousands of enhancers and core promoters in animal genomes have been annotated, and we have learned much about the domain structure in which regulatory genomes are organized in animals. Enhancer-core-promoter targeting occurs at several levels, including regulatory domains, DNA accessibility, and sequence-encoded core-promoter specificities that are likely mediated by different regulatory proteins. We review here current knowledge about enhancer-core-promoter targeting, regulatory communication between enhancers and core promoters, and the protein factors involved. We conclude with an outlook on open questions that we find particularly interesting and that will likely lead to additional insights in the upcoming years.

Role of non-coding RNA transcription around gene regulatory elements in transcription factor recruitment

Eukaryotic cells produce a variety of non-coding RNAs (ncRNAs), many of which have been shown to play pivotal roles in biological processes such as differentiation, maintenance of pluripotency of stem cells, and cellular response to various stresses. Genome-wide analyses have revealed that many ncRNAs are transcribed around regulatory DNA elements located proximal or distal to gene promoters, but their biological functions are largely unknown. Recently, it has been demonstrated in yeast and mouse that ncRNA transcription around gene promoters and enhancers facilitates DNA binding of transcription factors to their target sites. These results suggest universal roles of promoter/enhancer-associated ncRNAs in the recruitment of transcription factors to their binding sites.

Crosstalk between androgen and pro-inflammatory signaling remodels androgen receptor and NF-κB cistrome to reprogram the prostate cancer cell transcriptome

Inflammatory processes and androgen signaling are critical for the growth of prostate cancer (PC), the most common cancer among males in Western countries. To understand the importance of potential interplay between pro-inflammatory and androgen signaling for gene regulation, we have interrogated the crosstalk between androgen receptor (AR) and NF-κB, a key transcriptional mediator of inflammatory responses, by utilizing genome-wide chromatin immunoprecipitation sequencing and global run-on sequencing in PC cells. Co-stimulation of LNCaP cells with androgen and pro-inflammatory cytokine TNFα invoked a transcriptome which was very distinct from that induced by either stimulation alone. The altered transcriptome that included gene programs linked to cell migration and invasiveness was orchestrated by significant remodeling of NF-κB and AR cistrome and enhancer landscape. Although androgen multiplied the NF-κB cistrome and TNFα restrained the AR cistrome, there was no general reciprocal tethering of the AR to the NF-κB on chromatin. Instead, redistribution of FOXA1, PIAS1 and PIAS2 contributed to the exposure of latent NF-κB chromatin-binding sites and masking of AR chromatin-binding sites. Taken together, concomitant androgen and pro-inflammatory signaling significantly remodels especially the NF-κB cistrome, reprogramming the PC cell transcriptome in fashion that may contribute to the progression of PC.

A generative model for the behavior of RNA polymerase

MOTIVATION

Transcription by RNA polymerase is a highly dynamic process involving multiple distinct points of regulation. Nascent transcription assays are a relatively new set of high throughput techniques that measure the location of actively engaged RNA polymerase genome wide. Hence, nascent transcription is a rich source of information on the regulation of RNA polymerase activity. To fully dissect this data requires the development of stochastic models that can both deconvolve the stages of polymerase activity and identify significant changes in activity between experiments.

RESULTS

We present a generative, probabilistic model of RNA polymerase that fully describes loading, initiation, elongation and termination. We fit this model genome wide and profile the enzymatic activity of RNA polymerase across various loci and following experimental perturbation. We observe striking correlation of predicted loading events and regulatory chromatin marks. We provide principled statistics that compute probabilities reminiscent of traveler's and divergent ratios. We finish with a systematic comparison of RNA Polymerase activity at promoter versus non-promoter associated loci.

AVAILABILITY AND IMPLEMENTATION

Transcription Fit (Tfit) is a freely available, open source software package written in C/C ++ that requires GNU compilers 4.7.3 or greater. Tfit is available from GitHub (https://github.com/azofeifa/Tfit).

CONTACT

robin.dowell@colorado.eduSupplementary information: Supplementary data are available at Bioinformatics online.

BRD4 localization to lineage-specific enhancers is associated with a distinct transcription factor repertoire

Proper temporal epigenetic regulation of gene expression is essential for cell fate determination and tissue development. The Bromodomain-containing Protein-4 (BRD4) was previously shown to control the transcription of defined subsets of genes in various cell systems. In this study we examined the role of BRD4 in promoting lineage-specific gene expression and show that BRD4 is essential for osteoblast differentiation. Genome-wide analyses demonstrate that BRD4 is recruited to the transcriptional start site of differentiation-induced genes. Unexpectedly, while promoter-proximal BRD4 occupancy correlated with gene expression, genes which displayed moderate expression and promoter-proximal BRD4 occupancy were most highly regulated and sensitive to BRD4 inhibition. Therefore, we examined distal BRD4 occupancy and uncovered a specific co-localization of BRD4 with the transcription factors C/EBPb, TEAD1, FOSL2 and JUND at putative osteoblast-specific enhancers. These findings reveal the intricacies of lineage specification and provide new insight into the context-dependent functions of BRD4.

CTCF modulates Estrogen Receptor function through specific chromatin and nuclear matrix interactions

Enhancer regions and transcription start sites of estrogen-target regulated genes are connected by means of Estrogen Receptor long-range chromatin interactions. Yet, the complete molecular mechanisms controlling the transcriptional output of engaged enhancers and subsequent activation of coding genes remain elusive. Here, we report that CTCF binding to enhancer RNAs is enriched when breast cancer cells are stimulated with estrogen. CTCF binding to enhancer regions results in modulation of estrogen-induced gene transcription by preventing Estrogen Receptor chromatin binding and by hindering the formation of additional enhancer-promoter ER looping. Furthermore, the depletion of CTCF facilitates the expression of target genes associated with cell division and increases the rate of breast cancer cell proliferation. We have also uncovered a genomic network connecting loci enriched in cell cycle regulator genes to nuclear lamina that mediates the CTCF function. The nuclear lamina and chromatin interactions are regulated by estrogen-ER. We have observed that the chromatin loops formed when cells are treated with estrogen establish contacts with the nuclear lamina. Once there, the portion of CTCF associated with the nuclear lamina interacts with enhancer regions, limiting the formation of ER loops and the induction of genes present in the loop. Collectively, our results reveal an important, unanticipated interplay between CTCF and nuclear lamina to control the transcription of ER target genes, which has great implications in the rate of growth of breast cancer cells.

SETDB2 Links Glucocorticoid to Lipid Metabolism through Insig2a Regulation

Transcriptional and chromatin regulations mediate the liver response to nutrient availability. The role of chromatin factors involved in hormonal regulation in response to fasting is not fully understood. We have identified SETDB2, a glucocorticoid-induced putative epigenetic modifier, as a positive regulator of GR-mediated gene activation in liver. Insig2a increases during fasting to limit lipid synthesis, but the mechanism of induction is unknown. We show Insig2a induction is GR-SETDB2 dependent. SETDB2 facilitates GR chromatin enrichment and is key to glucocorticoid-dependent enhancer-promoter interactions. INSIG2 is a negative regulator of SREBP, and acute glucocorticoid treatment decreased active SREBP during refeeding or in livers of Ob/Ob mice, both systems of elevated SREBP-1c-driven lipogenesis. Knockdown of SETDB2 or INSIG2 reversed the inhibition of SREBP processing. Overall, these studies identify a GR-SETDB2 regulatory axis of hepatic transcriptional reprogramming and identify SETDB2 as a potential target for metabolic disorders with aberrant glucocorticoid actions.