Noncoding transcription at enhancers: general principles and functional models

Transcription Factor IRF8 Governs Enhancer Landscape Dynamics in Mononuclear Phagocyte Progenitors

Monocytes and dendritic cells (DCs), mononuclear phagocytes essential for immune responses, develop from hematopoietic stem cells via monocyte-DC progenitors (MDPs). The molecular basis of their development remains unclear. Because promoter-distal enhancers are key to cell fate decisions, we analyzed enhancer landscapes during mononuclear phagocyte development in vivo. Monocyte- and DC-specific enhancers were gradually established at progenitor stages before the expression of associated genes. Of the transcription factors predicted to bind to these enhancers, IRF8, essential for monocyte and DC development, was found to be required for the establishment of these enhancers, particularly those common to both monocyte and DC lineages. Although Irf8^-/- mononuclear phagocyte progenitors, including MDPs, displayed grossly normal gene expression patterns, their enhancer landscapes resembled that of an upstream progenitor population. Our results illustrate the dynamic process by which key transcription factors regulate enhancer formation and, therefore, direct future gene expression to achieve mononuclear phagocyte development.

Long Noncoding RNA: Genomics and Relevance to Physiology

The mammalian cell expresses thousands of long noncoding RNAs (lncRNAs) that are longer than 200 nucleotides but do not encode any protein. lncRNAs can change the expression of protein-coding genes through both cis and trans mechanisms, including imprinting and other types of transcriptional regulation, and posttranscriptional regulation including serving as molecular sponges. Deep sequencing, coupled with analysis of sequence characteristics, is the primary method used to identify lncRNAs. Physiological roles of specific lncRNAs can be examined using genetic targeting or knockdown with modified oligonucleotides. Identification of nucleic acids or proteins with which an lncRNA interacts is essential for understanding the molecular mechanism underlying its physiological role. lncRNAs have been reported to contribute to the regulation of physiological functions and disease development in several organ systems, including the cardiovascular, renal, muscular, endocrine, digestive, nervous, respiratory, and reproductive systems. The physiological role of the majority of lncRNAs, many of which are species and tissue specific, remains to be determined. © 2019 American Physiological Society. Compr Physiol 9:933-946, 2019.

LncRNAs and epigenetic regulation of vascular endothelium: genome positioning system and regulators of chromatin modifiers

Epigenetic mechanisms regulate the cell type-specific expression of endothelial-enriched genes. A major question has been how chromatin modifiers without inherent sequence specificity can be targeted to genomic coordinates. Recently, long noncoding RNAs (lncRNAs) have emerged as candidates for specifying genomic positioning for chromatin modifiers. However, lncRNAs function by a number of mechanisms in both the nucleus and the cytoplasm. Recent studies indicate the existence of endothelial-enriched lncRNAs. This review discusses lncRNA regulation in endothelial cells with a focus on four recently described nuclear-enriched lncRNAs: MANTIS, LEENE, STEEL, and GATA6-AS. This emerging work on these lncRNAs contributes to our understanding of epigenetic regulation in the vascular endothelium with links to important themes in endothelial biology, including angiogenesis and shear stress.

Function and Mechanism of Long Noncoding RNAs in Adipocyte Biology

The last two decades have witnessed an explosion of interest in adipocyte biology, coinciding with the upsurge of obesity and metabolic syndrome. Now we have new perspectives on the distinct developmental origins of white, brown, and beige adipocytes and their role in metabolic physiology and disease. Beyond fuel metabolism, adipocytes communicate with the immune system and other tissues by releasing diverse paracrine and endocrine factors to orchestrate adipose tissue remodeling and maintain systemic homeostasis. Significant progress has been made in delineating the regulatory networks that govern different aspects of adipocyte biology. Here we provide an overview on the emerging role of long noncoding RNAs (lncRNAs) in the regulation of adipocyte development and metabolism and discuss the implications of the RNA-protein regulatory interface in metabolic control.

Long Non-Coding RNAs and the Innate Immune Response

Innate immunity provides the initial defence against infection and it is now clear that long non-coding RNAs (lncRNAs) are important regulators of this response. Following activation of the innate response, we commonly see rapid induction of these lncRNAs and this is often mediated via the pro-inflammatory transcription factor, nuclear factor-κB (NF-κB). Knockdown studies have shown that lncRNAs tend to act in trans to regulate the expression of multiple inflammatory mediators and other responses. Mechanistically, many lncRNAs have demonstrated acting through heterogeneous nuclear ribonucleoproteins, complexes that are implicated chromatin re-modelling, transcription process and translation. In addition, these lncRNAs have also been shown to interact with multiple other proteins involved in the regulation of chromatin re-modelling, as well as those proteins involved in intracellular immune signalling, which include NF-κB. In this review, we will describe the evidence that supports this emerging role of lncRNA in the innate immune response.

Sequence Characteristics Distinguish Transcribed Enhancers from Promoters and Predict Their Breadth of Activity

Enhancers and promoters both regulate gene expression by recruiting transcription factors (TFs); however, the degree to which enhancer vs. promoter activity is due to differences in their sequences or to genomic context is the subject of ongoing debate. We examined this question by analyzing the sequences of thousands of transcribed enhancers and promoters from hundreds of cellular contexts previously identified by cap analysis of gene expression. Support vector machine classifiers trained on counts of all possible 6-bp-long sequences (6-mers) were able to accurately distinguish promoters from enhancers and distinguish their breadth of activity across tissues. Classifiers trained to predict enhancer activity also performed well when applied to promoter prediction tasks, but promoter-trained classifiers performed poorly on enhancers. This suggests that the learned sequence patterns predictive of enhancer activity generalize to promoters, but not vice versa. Our classifiers also indicate that there are functionally relevant differences in enhancer and promoter GC content beyond the influence of CpG islands. Furthermore, sequences characteristic of broad promoter or broad enhancer activity matched different TFs, with predicted ETS- and RFX-binding sites indicative of promoters, and AP-1 sites indicative of enhancers. Finally, we evaluated the ability of our models to distinguish enhancers and promoters defined by histone modifications. Separating these classes was substantially more difficult, and this difference may contribute to ongoing debates about the similarity of enhancers and promoters. In summary, our results suggest that high-confidence transcribed enhancers and promoters can largely be distinguished based on biologically relevant sequence properties.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

A critical regulator of Bcl2 revealed by systematic transcript discovery of lncRNAs associated with T-cell differentiation

Normal T-cell differentiation requires a complex regulatory network which supports a series of maturation steps, including lineage commitment, T-cell receptor (TCR) gene rearrangement, and thymic positive and negative selection. However, the underlying molecular mechanisms are difficult to assess due to limited T-cell models. Here we explore the use of the pro-T-cell line P5424 to study early T-cell differentiation. Stimulation of P5424 cells by the calcium ionophore ionomycin together with PMA resulted in gene regulation of T-cell differentiation and activation markers, partially mimicking the CD4^-CD8^- double negative (DN) to double positive (DP) transition and some aspects of subsequent T-cell maturation and activation. Global analysis of gene expression, along with kinetic experiments, revealed a significant association between the dynamic expression of coding genes and neighbor lncRNAs including many newly-discovered transcripts, thus suggesting potential co-regulation. CRISPR/Cas9-mediated genetic deletion of Robnr, an inducible lncRNA located downstream of the anti-apoptotic gene Bcl2, demonstrated a critical role of the Robnr locus in the induction of Bcl2. Thus, the pro-T-cell line P5424 is a powerful model system to characterize regulatory networks involved in early T-cell differentiation and maturation.

AP001056.1, A Prognosis-Related Enhancer RNA in Squamous Cell Carcinoma of the Head and Neck

A growing number of long non-coding RNAs (lncRNAs) have been linked to squamous cell carcinoma of the head and neck (SCCHN). A subclass of lncRNAs, termed enhancer RNAs (eRNAs), are derived from enhancer regions and could contribute to enhancer function. In this study, we developed an integrated data analysis approach to identify key eRNAs in SCCHN. Tissue-specific enhancer-derived RNAs and their regulated genes previously predicted using the computational pipeline PreSTIGE, were considered as putative eRNA-target pairs. The interactive web servers, TANRIC (the Atlas of Noncoding RNAs in Cancer) and cBioPortal, were used to explore the RNA levels and clinical data from the Cancer Genome Atlas (TCGA) project. Requiring that key eRNAs should show significant associations with overall survival (Kaplan⁻Meier log-rank test, p < 0.05) and the predicted target (correlation coefficient r > 0.4, p < 0.001), we identified five key eRNA candidates. The most significant survival-associated eRNA was AP001056.1 with ICOSLG encoding an immune checkpoint protein as its regulated target. Another 1640 genes also showed significant correlation with AP001056.1 (r > 0.4, p < 0.001), with the "immune system process" being the most significantly enriched biological process (adjusted p < 0.001). Our results suggest that AP001056.1 is a key immune-related eRNA in SCCHN with a positive impact on clinical outcome.

STAT5-Driven Enhancers Tightly Control Temporal Expression of Mammary-Specific Genes

The de novo formation of milk-secreting mammary epithelium during pregnancy is regulated by prolactin through activation of the transcription factor STAT5, which stimulates the expression of several hundred mammary-specific genes. In addition to its key role in activating gene expression in mammary tissue, STAT5, which is ubiquitously expressed in most cell types, implements T cell-specific programs controlled by interleukins. However, the mechanisms by which STAT5 controls cell-specific genetic programs activated by distinct cytokines remain relatively unknown. Integration of data from genome-wide surveys of chromatin markers and transcription factor binding at regulatory elements may shed light on the mechanisms that drive cell-specific programs. Here, we have illustrated how STAT5 controls cell-specific gene expression through its concentration and an auto-regulatory enhancer supporting its high levels in mammary tissue. The unique genomic features of STAT5-driven enhancers or super-enhancers that regulate mammary-specific genes and their dynamic remodeling in response to pregnancy hormone levels are described. We have further provided biological evidence supporting the in vivo function of a STAT5-driven super-enhancer with the aid of CRISPR/Cas9 genome editing. Finally, we discuss how the functions of mammary-specific super-enhancers are confined by the zinc finger protein, CTCF, to allow exclusive activation of mammary-specific genes without affecting common neighboring genes. This review comprehensively summarizes the molecular pathways underlying differential control of cell-specific gene sets by STAT5 and provides novel insights into STAT5-dependent mammary physiology.

Identification and regulation pattern analysis of long noncoding RNAs in meibomian gland carcinoma

Aim: To identify long noncoding RNAs (lncRNAs) and to elucidate regulation patterns of lncRNAs in meibomian gland carcinoma (MGC). Materials & methods: We used RNA-Seq, gene ontology, ClueGO, Ingenuity Pathway Analysis and co-expression network analyses to profile the expression and regulation patterns of lncRNAs and mRNAs in MGC. Results: We identified 500 lncRNAs and 326 mRNAs as differentially expressed. Co-expression regulatory networks with lncRNAs and mRNAs were constructed. The differentially expressed mRNAs and lncRNAs were enriched by fundamental biological functions that are implicated in the inflammatory signaling pathway and tumor proliferation (IL6 and PTGS2). Conclusion: LncRNAs might play important roles via the competing endogenous RNA regulation pattern in MGC tumorigenesis and contribute to the molecular pathogenesis of MGC.

Endogenous double-stranded Alu RNA elements stimulate IFN-responses in relapsing remitting multiple sclerosis

Various sensors that detect double-stranded RNA, presumably of viral origin, exist in eukaryotic cells and induce IFN-responses. Ongoing IFN-responses have also been documented in a variety of human autoimmune diseases including relapsing-remitting multiple sclerosis (RRMS) but their origins remain obscure. We find increased IFN-responses in leukocytes in relapsing-remitting multiple sclerosis at distinct stages of disease. Moreover, endogenous RNAs isolated from blood cells of these same patients recapitulate this IFN-response if transfected into naïve cells. These endogenous RNAs are double-stranded RNAs, contain Alu and Line elements and are transcribed from leukocyte transcriptional enhancers. Thus, transcribed endogenous retrotransposon elements can co-opt pattern recognition sensors to induce IFN-responses in RRMS.

A Pan-Cancer Analysis of Enhancer Expression in Nearly 9000 Patient Samples

The role of enhancers, a key class of non-coding regulatory DNA elements, in cancer development has increasingly been appreciated. Here, we present the detection and characterization of a large number of expressed enhancers in a genome-wide analysis of 8928 tumor samples across 33 cancer types using TCGA RNA-seq data. Compared with matched normal tissues, global enhancer activation was observed in most cancers. Across cancer types, global enhancer activity was positively associated with aneuploidy, but not mutation load, suggesting a hypothesis centered on "chromatin-state" to explain their interplay. Integrating eQTL, mRNA co-expression, and Hi-C data analysis, we developed a computational method to infer causal enhancer-gene interactions, revealing enhancers of clinically actionable genes. Having identified an enhancer ∼140 kb downstream of PD-L1, a major immunotherapy target, we validated it experimentally. This study provides a systematic view of enhancer activity in diverse tumor contexts and suggests the clinical implications of enhancers.

Spatiotemporal chromatin dynamics - A telltale of circadian epigenetic gene regulation

Over the course of evolution, nature has forced organisms under selection pressure to hardwire an internal time keeping device that defines 24 h of a daily cycle of physiological and behavioral rhythms, known as circadian rhythms. At the cellular level, the cycle is governed by significant fractions of transcriptomes, which are under the control of transcriptional and translational feedback loop of clock genes. Intriguingly, this feedback loop is regulated at multiple stratums such as at the transcriptional and translational levels, which direct a cell towards producing a robust rhythm by sustaining the repeated stoichiometry of protein products. Moreover, with the advent of state of the art paradigms, epigenetic regulation of circadian rhythms has been becoming more evident at present time. Light-induced recurring fluctuations in chromatin acetylation concurrent with the binding of RNA Pol II and integration of miRNAs monitor the chromatin modifiers or clock genes expression to drive temporal rhythmicity. Furthermore, CLOCK protein intrinsic histone acetyl transferase activity, the interaction of CLOCK-BMAL-1 with HAT enzymes, and the involvement of many histone deacetylases also maintain the rhythmic protein profile. Additionally, the critical role of the rhythmic methylation pattern of clock genes in battery of cancer and metabolic disorders also defines its importance. Therefore, in this review, we focused on accumulating all the present data available on epigenetics and circadian rhythms. Interestingly, we also gathered evidence from the available literature pinpointing towards the dynamic nature of chromatin architecture governed by long and short-range regulatory elements DNA contacts arising daily, that was thought to be steady otherwise.

Transcriptionally induced enhancers in the macrophage immune response to Mycobacterium tuberculosis infection

Background

Tuberculosis is a life-threatening infectious disease caused by Mycobacterium tuberculosis (M.tb). M.tb subverts host immune responses to build a favourable niche and survive inside of host macrophages. Macrophages can control or eliminate the infection, if acquire appropriate functional phenotypes. Transcriptional regulation is a key process that governs the activation and maintenance of these phenotypes. Among the factors orchestrating transcriptional regulation during M.tb infection, transcriptional enhancers still remain unexplored.

Results

We analysed transcribed enhancers in M.tb-infected mouse bone marrow-derived macrophages. We established a link between known M.tb-responsive transcription factors and transcriptional activation of enhancers and their target genes. Our data suggest that enhancers might drive macrophage response via transcriptional activation of key immune genes, such as Tnf, Tnfrsf1b, Irg1, Hilpda, Ccl3, and Ccl4. We report enhancers acquiring transcription de novo upon infection. Finally, we link highly transcriptionally induced enhancers to activation of genes with previously unappreciated roles in M.tb infection, such as Fbxl3, Tapt1, Edn1, and Hivep1.

Conclusions

Our findings suggest the importance of macrophage host transcriptional enhancers during M.tb infection. Our study extends current knowledge of the regulation of macrophage responses to M.tb infection and provides a basis for future functional studies on enhancer-gene interactions in this process.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5450-6) contains supplementary material, which is available to authorized users.

Aberrant enhancer hypomethylation contributes to hepatic carcinogenesis through global transcriptional reprogramming

Hepatocellular carcinomas (HCC) exhibit distinct promoter hypermethylation patterns, but the epigenetic regulation and function of transcriptional enhancers remain unclear. Here, our affinity- and bisulfite-based whole-genome sequencing analyses reveal global enhancer hypomethylation in human HCCs. Integrative epigenomic characterization further pinpoints a recurrent hypomethylated enhancer of CCAAT/enhancer-binding protein-beta (C/EBPβ) which correlates with C/EBPβ over-expression and poorer prognosis of patients. Demethylation of C/EBPβ enhancer reactivates a self-reinforcing enhancer-target loop via direct transcriptional up-regulation of enhancer RNA. Conversely, deletion of this enhancer via CRISPR/Cas9 reduces C/EBPβ expression and its genome-wide co-occupancy with BRD4 at H3K27ac-marked enhancers and super-enhancers, leading to drastic suppression of driver oncogenes and HCC tumorigenicity. Hepatitis B X protein transgenic mouse model of HCC recapitulates this paradigm, as C/ebpβ enhancer hypomethylation associates with oncogenic activation in early tumorigenesis. These results support a causal link between aberrant enhancer hypomethylation and C/EBPβ over-expression, thereby contributing to hepatocarcinogenesis through global transcriptional reprogramming.

There are distinct hypermethylation patterns in gene promoters in hepatocellular carcinomas (HCCs). Here, the authors show that the enhancer of C/EBPβ is recurrently hypomethylated in human HCCs, recapitulating this in a transgenic murine model and linking aberrant enhancer hypomethylation to hepatocarcinogenesis.

Enhancer RNAs: a missing regulatory layer in gene transcription

Enhancers and super-enhancers exert indispensable roles in maintaining cell identity through spatiotemporally regulating gene transcription. Meanwhile, active enhancers and super-enhancers also produce transcripts termed enhancer RNAs (eRNAs) from their DNA elements. Although enhancers have been identified for more than 30 years, widespread transcription from enhancers are just discovered by genome-wide sequencing and considered as the key to understand longstanding questions in gene transcription. RNA-transcribed enhancers are marked by histone modifications such as H3K4m1/2 and H3K27Ac, and enriched with transcription regulatory factors such as LDTFs, P300, CBP, BRD4 and MED1. Those regulatory factors might constitute a Mega-Trans-like complex to potently activate enhancers. Compared to mRNAs, eRNAs are quite unstable and play roles at local. Functionally, it has been shown that eRNAs promote formation of enhancer-promoter loops. Several studies also demonstrated that eRNAs help the binding of RNA polymerase II (RNAPII) or transition of paused RNAPII by de-association of the negative elongation factor (NELF) complex. Nevertheless, these proposed mechanisms are not universally accepted and still under controversy. Here, we comprehensively summarize the reported findings and make perspectives for future exploration. We also believe that super-enhancer derived RNAs (seRNAs) might be informative to understand the nature of super-enhancers.

lncRNA profile study reveals the mRNAs and lncRNAs associated with docetaxel resistance in breast cancer cells

Resistance to adjuvant systemic treatment, including taxanes (docetaxel and paclitaxel) is a major clinical problem for breast cancer patients. lncRNAs (long non-coding RNAs) are non-coding transcripts, which have recently emerged as important players in a variety of biological processes, including cancer development and chemotherapy resistance. However, the contribution of lncRNAs to docetaxel resistance in breast cancer and the relationship between lncRNAs and taxane-resistance genes are still unclear. Here, we performed comprehensive RNA sequencing and analyses on two docetaxel-resistant breast cancer cell lines (MCF7-RES and MDA-RES) and their docetaxel-sensitive parental cell lines. We identified protein coding genes and pathways that may contribute to docetaxel resistance. More importantly, we identified lncRNAs that were consistently up-regulated or down-regulated in both the MCF7-RES and MDA-RES cells. The co-expression network and location analyses pinpointed four overexpressed lncRNAs located within or near the ABCB1 (ATP-binding cassette subfamily B member 1) locus, which might up-regulate the expression of ABCB1. We also identified the lncRNA EPB41L4A-AS2 (EPB41L4A Antisense RNA 2) as a potential biomarker for docetaxel sensitivity. These findings have improved our understanding of the mechanisms underlying docetaxel resistance in breast cancer and have provided potential biomarkers to predict the response to docetaxel in breast cancer patients.

Tyrosine-1 of RNA Polymerase II CTD Controls Global Termination of Gene Transcription in Mammals

The carboxy-terminal domain (CTD) of RNA polymerase (Pol) II is composed of a repetition of YSPTSPS heptads and functions as a loading platform for protein complexes that regulate transcription, splicing, and maturation of RNAs. Here, we studied mammalian CTD mutants to analyze the function of tyrosine1 residues in the transcription cycle. Mutation of 3/4 of the tyrosine residues (YFFF mutant) resulted in a massive read-through transcription phenotype in the antisense direction of promoters as well as in the 3' direction several hundred kilobases downstream of genes. The YFFF mutant shows reduced Pol II at promoter-proximal pause sites, a loss of interaction with the Mediator and Integrator complexes, and impaired recruitment of these complexes to chromatin. Consistent with these observations, Pol II loading at enhancers and maturation of snRNAs are altered in the YFFF context genome-wide. We conclude that tyrosine1 residues of the CTD control termination of transcription by Pol II.

The enhancer RNA lnc-SLC4A1-1 epigenetically regulates unexplained recurrent pregnancy loss (URPL) by activating CXCL8 and NF-kB pathway

Background

Enhancer RNAs (eRNAs) are a group of lncRNAs transcribed from enhancers, whose regulatory effects on gene expression are an emerging area of interest. However, the role of eRNAs in regulating trophoblast cells and unexplained recurrent pregnancy loss (URPL) remains elusive.

Methods

We profiled eRNAs in villi from URPL patients and matched controls by RNA-seq. Functions of URPL-related eRNAs were further investigated in vitro.

Results

We identified lnc-SLC4A1-1, which was transcribed from an active enhancer marked with H3K27ac and H3K4me1 and so-called eRNA, highly expressed in URPL patients. Gain-of-function experiments indicated that lnc-SLC4A1-1 facilitated trophoblast cell migration and apoptosis. Mechanistically, as an eRNA, lnc-SLC4A1-1 was retained in the nuclei and recruited transcription factor NF-κB to bind to CXCL8, resulting in increased H3K27ac in the CXCL8 promoter and subsequent elevation of CXCL8 expression. Activation of CXCL8 exacerbated inflammatory reactions in trophoblast cells by inducing TNF-α and IL-1β, which could be blocked by an antagonist of lnc-SLC4A1-1.

Interpretation

These findings indicate that an eRNA, lnc-SLC4A1-1, alters trophoblast function via activation of immune responses and by regulating the NF-κB/CXCL8 axis. Our study provides new insights in understanding lncRNA/eRNA function in pathological pregnancy, potentially informing on therapeutic strategies for URPL.

Fund

National Natural Science Foundation of China, Natural Science Foundation of Jiangsu Province, National Key Research and Development Program, the Priority Academic Program for the Development of Jiangsu Higher Education Institutions.

Causal Transcription Regulatory Network Inference Using Enhancer Activity as a Causal Anchor

Transcription control plays a crucial role in establishing a unique gene expression signature for each of the hundreds of mammalian cell types. Though gene expression data have been widely used to infer cellular regulatory networks, existing methods mainly infer correlations rather than causality. We developed statistical models and likelihood-ratio tests to infer causal gene regulatory networks using enhancer RNA (eRNA) expression information as a causal anchor and applied the framework to eRNA and transcript expression data from the FANTOM Consortium. Predicted causal targets of transcription factors (TFs) in mouse embryonic stem cells, macrophages and erythroblastic leukaemia overlapped significantly with experimentally-validated targets from ChIP-seq and perturbation data. We further improved the model by taking into account that some TFs might act in a quantitative, dosage-dependent manner, whereas others might act predominantly in a binary on/off fashion. We predicted TF targets from concerted variation of eRNA and TF and target promoter expression levels within a single cell type, as well as across multiple cell types. Importantly, TFs with high-confidence predictions were largely different between these two analyses, demonstrating that variability within a cell type is highly relevant for target prediction of cell type-specific factors. Finally, we generated a compendium of high-confidence TF targets across diverse human cell and tissue types.

Production of Spliced Long Noncoding RNAs Specifies Regions with Increased Enhancer Activity

Active enhancers in mammals produce enhancer RNAs (eRNAs) that are bidirectionally transcribed, unspliced, and unstable. Enhancer regions are also enriched with long noncoding RNA (lncRNA) transcripts, which are typically spliced and substantially more stable. In order to explore the relationship between these two classes of RNAs, we analyzed DNase hypersensitive sites with evidence of bidirectional transcription, which we termed eRNA-producing centers (EPCs). EPCs found very close to transcription start sites of lncRNAs exhibit attributes of both enhancers and promoters, including distinctive DNA motifs and a characteristic chromatin landscape. These EPCs are associated with higher enhancer activity, driven at least in part by the presence of conserved, directional splicing signals that promote lncRNA production, pointing at a causal role of lncRNA processing in enhancer activity. Together, our results suggest that the conserved ability of some enhancers to produce lncRNAs augments their activity in a manner likely mediated through lncRNA maturation.

The Transcription Factor STAT6 Mediates Direct Repression of Inflammatory Enhancers and Limits Activation of Alternatively Polarized Macrophages

The molecular basis of signal-dependent transcriptional activation has been extensively studied in macrophage polarization, but our understanding remains limited regarding the molecular determinants of repression. Here we show that IL-4-activated STAT6 transcription factor is required for the direct transcriptional repression of a large number of genes during in vitro and in vivo alternative macrophage polarization. Repression results in decreased lineage-determining transcription factor, p300, and RNA polymerase II binding followed by reduced enhancer RNA expression, H3K27 acetylation, and chromatin accessibility. The repressor function of STAT6 is HDAC3 dependent on a subset of IL-4-repressed genes. In addition, STAT6-repressed enhancers show extensive overlap with the NF-κB p65 cistrome and exhibit decreased responsiveness to lipopolysaccharide after IL-4 stimulus on a subset of genes. As a consequence, macrophages exhibit diminished inflammasome activation, decreased IL-1β production, and pyroptosis. Thus, the IL-4-STAT6 signaling pathway establishes an alternative polarization-specific epigenenomic signature resulting in dampened macrophage responsiveness to inflammatory stimuli.

•

IL-4-activated STAT6 acts as a transcriptional repressor in macrophages

•

IL-4-STAT6-repressed enhancers associate with reduced LDTF and p300 binding

•

Inflammatory responsiveness of the IL-4-repressed enhancers is attenuated

•

IL-4 limits the LPS-induced inflammasome activation, IL-1β production, and pyroptosis

Moonlighting Phosphatase Activity of Klenow DNA Polymerase in the Presence of RNA

RNA is a key player in the cellular central dogma, including RNA transcription and protein synthesis. However, it is unknown whether RNA can directly interfere with DNA synthesis. Recently, we have found in vitro that while binding to DNA polymerase nonspecifically, RNA can transform DNA polymerase to display a moonlighting activity, dNTP phosphatase, in turn interfering with DNA synthesis. This phosphatase activity removes the γ-phosphate from dNTPs (generating dNDPs) and subsequently removes the β-phosphate from the formed dNDPs (generating dNMPs), confirmed by the noncleavable α,β-CH₂-dGTP and β,γ-CH₂-dGTP analogues. We also found that dGTP is the best substrate for the phosphatase, and the dNTP phosphatase activity is sensitive to the reaction medium. In addition, we have revealed that RNA can tune the activity of closely related proteins and give rise to new catalytic functions with subtle differences. Moreover, we have demonstrated in vitro that at the lower dNTP level, this phosphatase can directly inhibit DNA synthesis by dNTP depletion, though the phosphatase activity is 690-fold slower than the polymerase activity. Our observation in vitro suggests a plausible strategy for RNA to directly interfere with DNA polymerase and DNA synthesis in vivo.

Extensive and functional overlap of the STAT6 and RXR cistromes in the active enhancer repertoire of human CD14+ monocyte derived differentiating macrophages

Macrophages are able to differentiate into classically polarized (M1) or alternatively polarized (M2) states upon encountering pro-inflammatory cytokines such as interferon (IFN) γ or anti-inflammatory cytokines such as interleukin (IL) -4/IL-13, respectively. Moreover, macrophages are known to regulate lipid metabolism via multiple members of the nuclear hormone receptor family, including the retinoid X receptors (RXR). It has been also documented that cytokines are able to modulate macrophage responses to lipid signals but the nature of these interactions and the underlying mechanisms of these processes especially at the level of the chromatinized genome are not well understood. Previous work from our laboratory suggested that STAT6 is a facilitator of nuclear receptor mediated transcriptional activity acting at the genome level. This prompted us to investigate genome-wide DNA binding events and the development of cistromes in human CD14+ monocyte-derived macrophages upon exposure to IL-4. We determined the impact of IL-4 on the PU.1, RXR and STAT6 cistromes within the active enhancer regions marked by H3K27-acetylation using chromatin immunoprecipitation followed by deep sequencing and integrated bioinformatics analyses. We found that about 2/3rd of the IL-4 induced STAT6 peaks co-localized with RXR peaks. These STAT6/RXR co-peaks differed at least in part from the non-overlapping RXR peaks regarding the most enriched de novo transcription factor binding motifs. Interestingly, RXR-binding was not regulated at the STAT6/RXR co-bound enhancers following IL-4 stimulation, but differential enhancer interactions were observed between the IL-4/STAT6 and RXR signaling pathways acting in a gene selective manner. Our results suggest that there is a novel, so far uncharacterized cistromic crosstalk between RXR and STAT6 that is likely to contribute to the formation of the active enhancer repertoire, transcriptome and differential signal-specific gene regulation of polarized macrophages.

Role of Pseudogenes in Tumorigenesis

Functional genomics has provided evidence that the human genome transcribes a large number of non-coding genes in addition to protein-coding genes, including microRNAs and long non-coding RNAs (lncRNAs). Among the group of lncRNAs are pseudogenes that have not been paid attention in the past, compared to other members of lncRNAs. However, increasing evidence points the important role of pseudogenes in diverse cellular functions, and dysregulation of pseudogenes are often associated with various human diseases including cancer. Like other types of lncRNAs, pseudogenes can also function as master regulators for gene expression and thus, they can play a critical role in various aspects of tumorigenesis. In this review we discuss the latest developments in pseudogene research, focusing on how pseudogenes impact tumorigenesis through different gene regulation mechanisms. Given the high sequence homology with the corresponding parent genes, we also discuss challenges for pseudogene research.

Human Enhancers Harboring Specific Sequence Composition, Activity, and Genome Organization Are Linked to the Immune Response

The FANTOM5 consortium recently characterized 65,423 human enhancers from 1829 cell and tissue samples using the Cap Analysis of Gene Expression technology. We showed that the guanine and cytosine content at enhancer regions distinguishes two classes of enhancers harboring distinct DNA structural properties at flanking regions. A functional analysis of their predicted gene targets highlighted one class of enhancers as significantly enriched for associations with immune response genes. Moreover, these enhancers were specifically enriched for regulatory motifs recognized by transcription factors involved in immune response. We observed that enhancers enriched for links to immune response genes were more cell-type specific, preferentially activated upon bacterial infection, and with specific response activity. Looking at chromatin capture data, we found that the two classes of enhancers were lying in distinct topologically associating domains and chromatin loops. Our results suggest that specific nucleotide compositions encode for classes of enhancers that are functionally distinct and specifically organized in the human genome.

Genome-wide Identification and Characterization of Enhancers Across 10 Human Tissues

Background: Enhancers can act as cis-regulatory elements (CREs) to control development and cellular function by regulating gene expression in a tissue-specific and ubiquitous manner. However, the regulatory network and characteristic of different types of enhancers (e.g., transcribed/non-transcribed enhancers, tissue-specific/ubiquitous enhancers) across multiple tissues are still unclear. Results: Here, a total of 53,924 active enhancers and 10,307 enhancer-associated RNAs (eRNAs) in 10 tissues (adrenal, brain, breast, heart, liver, lung, ovary, placenta, skeletal muscle and kidney) were identified through the integration of histone modifications (H3K4me1, H3K27ac and H3K4me3) and DNase I hypersensitive sites (DHSs) data. Moreover, 40,101 tissue-specific enhancers (TS-Enh), 1,241 ubiquitously expressed enhancers (UE-Enh) as well as transcribed enhancers (T-Enh), including 7,727 unidirectionally transcribed enhancers (1D-Enh) and 1,215 bidirectionally transcribed enhancers (2D-Enh) were defined in 10 tissues. The results show that enhancers exhibited high GC content, genomic variants and transcription factor binding sites (TFBS) enrichment in all tissues. These characteristics were significantly different between TS-Enh and UE-Enh, T-Enh and NT-Enh, 2D-Enh and 1D-Enh. Furt hermore, the results showed that enhancers obviously upregulate the expression of adjacent target genes which were remarkably correlated with the functions of corresponding tissues. Finally, a free user-friendly tissue-specific enhancer database, TiED (http://lcbb.swjtu.edu.cn/TiED), has been built to store, visualize, and confer these results. Conclusion: Genome-wide analysis of the regulatory network and characteristic of various types of enhancers showed that enhancers associated with TFs, eRNAs and target genes appeared in tissue specificity and function across different tissues.

A Muscle-Specific Enhancer RNA Mediates Cohesin Recruitment and Regulates Transcription In trans

The enhancer regions of the myogenic master regulator MyoD give rise to at least two enhancer RNAs. Core enhancer eRNA (CEeRNA) regulates transcription of the adjacent MyoD gene, whereas DRReRNA affects expression of Myogenin in trans. We found that DRReRNA is recruited at the Myogenin locus, where it colocalizes with Myogenin nascent transcripts. DRReRNA associates with the cohesin complex, and this association correlates with its transactivating properties. Despite being expressed in undifferentiated cells, cohesin is not loaded on Myogenin until the cells start expressing DRReRNA, which is then required for cohesin chromatin recruitment and maintenance. Functionally, depletion of either cohesin or DRReRNA reduces chromatin accessibility, prevents Myogenin activation, and hinders muscle cell differentiation. Thus, DRReRNA ensures spatially appropriate cohesin loading in trans to regulate gene expression.

An Enhancer-Derived RNA Muscles In to Regulate Myogenin In trans

Noncoding RNAs (ncRNAs) transcribed from active enhancers are known as enhancer RNAs (eRNAs). eRNAs have generally been shown to contribute to transcriptional activation of target genes in cis. In this issue, Tsai et al. (2018) demonstrate that an eRNA expressed from a distal enhancer of Myogenic differentiation1 (MyoD) (DRReRNA) does not regulate its neighboring MyoD; instead, it promotes myogenic differentiation by activating Myogenin, which is located on a different chromosome.

Long non-coding RNAs in ovarian cancer

Long non-coding RNAs (lncRNAs) refer to functional cellular RNAs molecules longer than 200 nucleotides in length. Unlike microRNAs, which have been widely studied, little is known about the enigmatic role of lncRNAs. However, lncRNAs have motivated extensively attention in the past few years and are emerging as potentially important regulators in pathological processes, including in cancer. We now understand that lncRNAs play role in cancer through their interactions with DNA, protein, and RNA in many instances. Moreover, accumulating evidence has recognized that large classes of lncRNAs are functional for ovarian cancer. Nevertheless, the biological phenomena and molecular mechanisms of lncRNAs in ovarian cancer remain to be better identified. In this review, we outline the dysregulated expression of lncRNAs and their potential clinical implications in ovarian cancer, with a particular emphasis on discussing the well characterized mechanisms underlying lncRNAs in ovarian cancer.

Enhancer RNAs (eRNAs): New Insights into Gene Transcription and Disease Treatment

Enhancers are cis-acting elements that have the ability to increase the expression of target genes. Recent studies have shown that enhancers can act as transcriptional units for the production of enhancer RNAs (eRNAs), which are hallmarks of activity enhancers and are involved in the regulation of gene transcription. The in-depth study of eRNAs is of great significance for us to better understand enhancer function and transcriptional regulation in various diseases. Therefore, eRNAs may be a potential therapeutic target for diseases. Here, we review the current knowledge of the characteristics of eRNAs, the molecular mechanisms of eRNAs action, as well as diseases related to dysregulation of eRNAs.

Long non-coding regulatory RNAs in sponges and insights into the origin of animal multicellularity

How animals evolved from a single-celled ancestor over 700 million years ago is poorly understood. Recent transcriptomic and chromatin analyses in the sponge Amphimedon queenslandica, a morphologically-simple representative of one of the oldest animal phyletic lineages, have shed light on what innovations in the genome and its regulation underlie the emergence of animal multicellularity. Comparisons of the regulatory genome of this sponge with those of more complex bilaterian model species and even simpler unicellular relatives have revealed that fundamental changes in genome regulatory complexity accompanied the evolution of animal multicellularity. Here, we review and discuss the results of these recent investigations by specifically focusing on the contribution of long non-coding RNAs to the evolution of the animal regulatory genome.

Current Advances on the Important Roles of Enhancer RNAs in Gene Regulation and Cancer

Revealing the gene regulation networks governing cancer initiation and development is necessary while it remains uncompleted. In recent years, enhancers have been reported to be widely transcribed, resulting in the generation of enhancer RNAs (eRNAs). Previous studies have reported that eRNAs are a subclass of long noncoding RNAs (lncRNAs), which play a critical role in gene regulation and cancer development. These eRNAs can promote enhancer-promoter (E-P) looping formation by binding to other protein factors or propel expression of downstream protein-coding gene. In this review, we have focused on the characteristics of eRNAs and illustrated the biological function and potential mechanism of eRNAs in regulating gene expression and cancer development.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.

JMJD6 Licenses ERα-Dependent Enhancer and Coding Gene Activation by Modulating the Recruitment of the CARM1/MED12 Co-activator Complex

Whereas the actions of enhancers in gene transcriptional regulation are well established, roles of JmjC-domain-containing proteins in mediating enhancer activation remain poorly understood. Here, we report that recruitment of the JmjC-domain-containing protein 6 (JMJD6) to estrogen receptor alpha (ERα)-bound active enhancers is required for RNA polymerase II recruitment and enhancer RNA production on enhancers, resulting in transcriptional pause release of cognate estrogen target genes. JMJD6 is found to interact with MED12 in the mediator complex to regulate its recruitment. Unexpectedly, JMJD6 is necessary for MED12 to interact with CARM1, which methylates MED12 at multiple arginine sites and regulates its chromatin binding. Consistent with its role in transcriptional activation, JMJD6 is required for estrogen/ERα-induced breast cancer cell growth and tumorigenesis. Our data have uncovered a critical regulator of estrogen/ERα-induced enhancer coding gene activation and breast cancer cell potency, providing a potential therapeutic target of ER-positive breast cancers.

The Super-Enhancer-Derived alncRNA-EC7/Bloodlinc Potentiates Red Blood Cell Development in trans

Enhancer-derived RNAs are thought to act locally by contributing to their parent enhancer function. Whether large domains of clustered enhancers (super-enhancers) also produce cis-acting RNAs, however, remains unclear. Unlike typical enhancers, super-enhancers form large spans of robustly transcribed chromatin, amassing capped and polyadenylated RNAs that are sufficiently abundant to sustain trans functions. Here, we show that one such RNA, alncRNA-EC7/Bloodlinc, is transcribed from a super-enhancer of the erythroid membrane transporter SLC4A1/BAND3 but diffuses beyond this site. Bloodlinc localizes to trans-chromosomal loci encoding critical regulators and effectors of terminal erythropoiesis and directly binds chromatin-organizing and transcription factors, including the chromatin attachment factor HNRNPU. Inhibiting Bloodlinc or Hnrnpu compromises the terminal erythropoiesis gene program, blocking red cell production, whereas expressing Bloodlinc ectopically stimulates this program and can promote erythroblast proliferation and enucleation in the absence of differentiation stimuli. Thus, Bloodlinc is a trans-acting super-enhancer RNA that potentiates red blood cell development.

SP8 Transcriptional Regulation of Cyclin D1 During Mouse Early Corticogenesis

Multiple signals control the balance between proliferation and differentiation of neural progenitor cells during corticogenesis. A key point of this regulation is the control of G1 phase length, which is regulated by the Cyclin/Cdks complexes. Using genome-wide chromatin immunoprecipitation assay and mouse genetics, we have explored the transcriptional regulation of Cyclin D1 (Ccnd1) during the early developmental stages of the mouse cerebral cortex. We found evidence that SP8 binds to the Ccnd1 locus on exon regions. In vitro experiments show SP8 binding activity on Ccnd1 gene 3'-end, and point to a putative role for SP8 in modulating PAX6-mediated repression of Ccnd1 along the dorso-ventral axis of the developing pallium, creating a medialLow-lateralHigh gradient of neuronal differentiation. Activation of Ccnd1 through the promoter/5'-end of the gene does not depend on SP8, but on βcatenin (CTNNB1). Importantly, alteration of the Sp8 level of expression in vivo affects Ccnd1 expression during early corticogenesis. Our results indicate that Ccnd1 regulation is the result of multiple signals and that SP8 is a player in this regulation, revealing an unexpected and potentially novel mechanism of transcriptional activation.

Comprehensive analysis of lncRNAs microarray profile and mRNA-lncRNA co-expression in oncogenic HPV-positive cervical cancer cell lines

Long non-coding RNAs are emerging to be novel regulators in gene expression. In current study, lncRNAs microarray and lncRNA-mRNA co-expression analysis were performed to explore the alternation and function of lncRNAs in cervical cancer cells. We identified that 4750 lncRNAs (15.52%) were differentially expressed in SiHa (HPV-16 positive) (2127 up-regulated and 2623 down-regulated) compared with C-33A (HPV negative), while 5026 lncRNAs (16.43%) were differentially expressed in HeLa (HPV-18 positive) (2218 up-regulated and 2808 down-regulated) respectively. There were 5008 mRNAs differentially expressed in SiHa and 4993 in HeLa, which were all cataloged by GO terms and KEGG pathway. With the help of mRNA-lncRNA co-expression network, we found that ENST00000503812 was significantly negative correlated with RAD51B and IL-28A expression in SiHa, while ENST00000420168, ENST00000564977 and TCONS_00010232 had significant correlation with FOXQ1 and CASP3 expression in HeLa. Up-regulation of ENST00000503812 may inhibit RAD51B and IL-28A expression and result in deficiency of DNA repair pathway and immune responses in HPV-16 positive cervical cancer cell. Up-regulation of ENST00000420168, ENST00000564977 and down-regulation of TCONS_00010232 might stimulate FOXQ1 expression and suppress CASP3 expression in HPV-18 positive cervical cancer cell, which lead to HPV-induced proliferation and deficiency in apoptosis. These results indicate that changes of lncRNAs and related mRNAs might impact on several cellular pathways and involve in HPV-induced proliferation, which enriches our understanding of lncRNAs and coding transcripts anticipated in HPV oncogenesis of cervical cancer.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Individual effects of the copia and gypsy enhancer and insulator on chromatin marks, eRNA synthesis, and binding of insulator proteins in transfected genetic constructs

Enhancers and insulators are involved in the regulation of gene expression, but the basic underlying mechanisms of action of these elements are unknown. We analyzed the individual effects of the enhancer and the insulator from Drosophila mobile elements copia [enh(copia)] and gypsy using transfected genetic constructs in S2 cells. This system excludes the influence of genomic cis regulatory elements. The enhancer-induced synthesis of 350-1050-nt-long enhancer RNAs (eRNAs) and H3K4me3 and H3K18ac marks, mainly in the region located about 300bp downstream of the enhancer. Insertion of the insulator between the enhancer and the promoter reduced these effects. We also observed the binding of dCTCF to the enhancer and to gypsy insulator. Our data indicate that a single gypsy insulator interacts with both the enhancer and the promoter, while two copies of the gypsy insulator preferentially interact with each other. Our results suggest the formation of chromatin loops that are shaped by the enhancer and the insulator.

The degree of enhancer or promoter activity is reflected by the levels and directionality of eRNA transcription

Here, Mikhaylichenko et al. investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. The authors demonstrate that while the timing of enhancer transcription is correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers and conclude that this is likely an inherent sequence property of the elements themselves.

Gene expression is regulated by promoters, which initiate transcription, and enhancers, which control their temporal and spatial activity. However, the discovery that mammalian enhancers also initiate transcription questions the inherent differences between enhancers and promoters. Here, we investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. We show that while the timing of enhancer transcription is generally correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers. To assess how this impacts function, we developed a dual transgenic assay to simultaneously measure enhancer and promoter activities from a single element in the same embryo. Extensive transgenic analysis revealed a relationship between the direction of endogenous transcription and the ability to function as an enhancer or promoter in vivo, although enhancer RNA (eRNA) production and activity are not always strictly coupled. Some enhancers (mainly bidirectional) can act as weak promoters, producing overlapping spatio–temporal expression. Conversely, bidirectional promoters often act as strong enhancers, while unidirectional promoters generally cannot. The balance between enhancer and promoter activity is generally reflected in the levels and directionality of eRNA transcription and is likely an inherent sequence property of the elements themselves.

Knockdown of Nuclear-Located Enhancer RNAs and Long ncRNAs Using Locked Nucleic Acid GapmeRs

The human genome is widely transcribed outside of protein-coding genes, producing thousands of noncoding RNAs from different subfamilies including enhancer RNAs. Functional studies to determine the role of individual genes are challenging with noncoding RNAs appearing to be more difficult to knockdown than mRNAs. One factor that may have hindered progress is that the majority of noncoding RNAs are thought to be located within the nucleus, where the efficiency of traditional RNA interference techniques is debatable. Here we present an alternative RNA interference technique utilizing Locked Nucleic Acids, which is able to efficiently knockdown noncoding RNAs irrespective of intracellular location.

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.

CRISPR-Cas9-mediated functional dissection of 3'-UTRs

Many studies using reporter assays have demonstrated that 3′ untranslated regions (3′-UTRs) regulate gene expression by controlling mRNA stability and translation. Due to intrinsic limitations of heterologous reporter assays, we sought to develop a gene editing approach to investigate the regulatory activity of 3′-UTRs in their native context. We initially used dual-CRISPR (clustered, regularly interspaced, short palindromic repeats)-Cas9 targeting to delete DNA regions corresponding to nine chemokine 3′-UTRs that destabilized mRNA in a reporter assay. Targeting six chemokine 3′-UTRs increased chemokine mRNA levels as expected. However, targeting CXCL1, CXCL6 and CXCL8 3′-UTRs unexpectedly led to substantial mRNA decreases. Metabolic labeling assays showed that targeting these three 3′-UTRs increased mRNA stability, as predicted by the reporter assay, while also markedly decreasing transcription, demonstrating an unexpected role for 3′-UTR sequences in transcriptional regulation. We further show that CRISPR–Cas9 targeting of specific 3′-UTR elements can be used for modulating gene expression and for highly parallel localization of active 3′-UTR elements in the native context. Our work demonstrates the duality and complexity of 3′-UTR sequences in regulation of gene expression and provides a useful approach for modulating gene expression and for functional annotation of 3′-UTRs in the native context.

Intragenic Enhancers Attenuate Host Gene Expression

Eukaryotic gene transcription is regulated at many steps, including RNA polymerase II (Pol II) recruitment, transcription initiation, promoter-proximal Pol II pause release, and transcription termination; however, mechanisms regulating transcription during productive elongation remain poorly understood. Enhancers, which activate gene transcription, themselves undergo Pol II-mediated transcription, but our understanding of enhancer transcription and enhancer RNAs (eRNAs) remains incomplete. Here we show that transcription at intragenic enhancers interferes with and attenuates host gene transcription during productive elongation. While the extent of attenuation correlates positively with nascent eRNA expression, the act of intragenic enhancer transcription alone, but not eRNAs, explains the attenuation. Through CRISPR/Cas9-mediated deletions, we demonstrate a physiological role for intragenic enhancer-mediated transcription attenuation in cell fate determination. We propose that intragenic enhancers not only enhance transcription of one or more genes from a distance but also fine-tune transcription of their host gene through transcription interference, facilitating differential utilization of the same regulatory element for disparate functions.

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints

In the era of genome-wide association studies (GWAS) and personalized medicine, predicting the impact of single nucleotide polymorphisms (SNPs) in regulatory elements is an important goal. Current approaches to determine the potential of regulatory SNPs depend on inadequate knowledge of cell-specific DNA binding motifs. Here, we present Sasquatch, a new computational approach that uses DNase footprint data to estimate and visualize the effects of noncoding variants on transcription factor binding. Sasquatch performs a comprehensive k-mer-based analysis of DNase footprints to determine any k-mer's potential for protein binding in a specific cell type and how this may be changed by sequence variants. Therefore, Sasquatch uses an unbiased approach, independent of known transcription factor binding sites and motifs. Sasquatch only requires a single DNase-seq data set per cell type, from any genotype, and produces consistent predictions from data generated by different experimental procedures and at different sequence depths. Here we demonstrate the effectiveness of Sasquatch using previously validated functional SNPs and benchmark its performance against existing approaches. Sasquatch is available as a versatile webtool incorporating publicly available data, including the human ENCODE collection. Thus, Sasquatch provides a powerful tool and repository for prioritizing likely regulatory SNPs in the noncoding genome.

Mutations in ACTRT1 and its enhancer RNA elements lead to aberrant activation of Hedgehog signaling in inherited and sporadic basal cell carcinomas

Basal cell carcinoma (BCC), the most common human cancer, results from aberrant activation of the Hedgehog signaling pathway. Although most cases of BCC are sporadic, some forms are inherited, such as Bazex-Dupré-Christol syndrome (BDCS)-a cancer-prone genodermatosis with an X-linked, dominant inheritance pattern. We have identified mutations in the ACTRT1 gene, which encodes actin-related protein T1 (ARP-T1), in two of the six families with BDCS that were examined in this study. High-throughput sequencing in the four remaining families identified germline mutations in noncoding sequences surrounding ACTRT1. These mutations were located in transcribed sequences encoding enhancer RNAs (eRNAs) and were shown to impair enhancer activity and ACTRT1 expression. ARP-T1 was found to directly bind to the GLI1 promoter, thus inhibiting GLI1 expression, and loss of ARP-T1 led to activation of the Hedgehog pathway in individuals with BDCS. Moreover, exogenous expression of ACTRT1 reduced the in vitro and in vivo proliferation rates of cell lines with aberrant activation of the Hedgehog signaling pathway. In summary, our study identifies a disease mechanism in BCC involving mutations in regulatory noncoding elements and uncovers the tumor-suppressor properties of ACTRT1.

Gene regulation in the immune system by long noncoding RNAs

Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression in the immune system. Studies have shown that lncRNAs are expressed in a highly lineage-specific manner and control the differentiation and function of innate and adaptive cell types. In this Review, we focus on mechanisms used by lncRNAs to regulate genes encoding products involved in the immune response, including direct interactions with chromatin, RNA and proteins. In addition, we address new areas of lncRNA biology, such as the functions of enhancer RNAs, circular RNAs and chemical modifications to RNA in cellular processes. We emphasize critical gaps in knowledge and future prospects for the roles of lncRNAs in the immune system and autoimmune disease.