Enhancer RNAs: mechanisms in transcriptional regulation and functions in diseases

In recent years, increasingly more non-coding RNAs have been detected with the development of high-throughput sequencing technology, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), small nucleolar RNAs (snoRNAs), and piwi-interacting RNA (piRNAs). The discovery of enhancer RNAs (eRNAs) in 2010 has further broadened the range of non-coding RNAs revealed. eRNAs are non-coding RNA molecules produced by the transcription of DNA cis-acting elements, enhancer fragments. Recent studies revealed that the transcription of eRNAs may be a biological marker responding to enhancer activity that can participate in the regulation of coding gene transcription. In this review, we discussed the biological characteristics of eRNAs, their functions in transcriptional regulation, the regulation factors of eRNAs production, and the research progress of eRNAs in different diseases. Video Abstract.

The landscape of enhancer RNA identify prognosis-related molecular subtypes in gastric cancer

BACKGROUND

Enhancer RNAs (eRNAs), the transcriptional products of active enhancers, are of great significance in the initial progression of cancers. However, the biological function and bioinformatics profiles of eRNA in gastric cancer remains largely enigmatic.

METHODS

Firstly, STAD were clustered into three subtypes with the data of eRNA expression from TCeA. Then we explored the difference of the tumor immune microenvironment, transcription levels, and transcription regulation among the three clusters. Finally, samples collected from 12 patients diagnosed with STAD were used to conduct qRT-PCR, verifying the conclusion based on network database.

RESULTS

The three clusters were detected to have different tumor microenvironments: Cluster A has an immune "cold" microenvironment. While cluster B features as more infiltration of immune cells, accompanied with higher expression of immune checkpoints such as PDCD1, LAG3, and TIGIT. Besides, Cluster C shows a higher stromal feature with B lineage, neutrophils, and fibroblasts. Further analyses indicated that CpG island methylation level of Cluster B is different from the other two clusters. Meanwhile, Cluster A and B showed significant enrichment of TP53 and KRAS mutation respectively while Cluster C has higher tumor mutation burden (TMB) and microsatellite instability (MSI). With the elaboration of transcriptional regulation of epigenetic clustering, we detected that Cluster A enriched in epithelial phenotype pathways. Cluster B enriched in cell-cell adhesion. Cluster C enriched in fibroblast proliferation. The clinical cohort show that Cluster B patients have lower interstitial cell characteristics and CAF infiltration.

CONCLUSION

We identified three unique epigenetic clusters of STAD through the differential activation of super-enhancers, and identified Cluster B with a higher immune infiltrating and a better prognosis, which provides a novel understanding of eRNAs and potential clinical applicability of eRNA-based molecular subtypes in gastric cancer.

Evaluate the immune-related eRNA models and signature score to predict the response to immunotherapy in thyroid carcinoma

BACKGROUND

The functional alterations of eRNAs have been reported to be correlated with tumorigenesis. However, the roles of eRNAs in thyroid cancer (THCA) remain still unclear. This study aimed to construct an immune-related eRNA prognostic signature that could effectively predict the survival and prognosis for THCA.

METHODS

The Weighted Gene Co-Expression Network Analysis (WGCNA) was performed to identify THCA-specific immune-related hub genes and immune-related eRNAs were obtained using Pearson correlation analysis. Univariate and least absolute shrinkage and selection operator (LASSO) Cox regression were conducted to construct an immune-related eRNA prognostic signature in training cohort, and the predictive capability was verified in test cohort and entire cohort. Kaplan-Meier analysis, principal component analysis (PCA), receiver operating characteristic (ROC) curves, and nomogram were used to validate the risk signature. Furthermore, CIBERSORT, ESTIMATE and ssGSEA were analyzed to explore the tumor immune microenvironment (TIME) of the risk signature, and the response of potential immunotherapeutic were also discussed.

RESULTS

A total of 125 immune-related eRNAs were obtained and 16 immune-related eRNAs were significantly correlated with overall survival (OS). A 9-immune-related eRNA prognostic signature was constructed, and the risk score was identified as an independent predictor. High-risk groups were associated with a poorer OS. Immune microenvironment analysis indicated that low risk score was correlated with higher immuneScore, high immune cell infiltration, and the better response of immunotherapy. Additionally, we also detected 9 immune-related eRNA expression levels in sixty-two matched tumorous and non-tumorous tissues using qRT-PCR analysis.

CONCLUSION

Our immune-related eRNA risk signature that was an independent prognostic factor was strongly correlated with the immune microenvironment and may be promising for the clinical prediction of prognosis and immunotherapeutic responses in THCA patients.

Ago1 controls myogenic differentiation by regulating eRNA-mediated CBP-guided epigenome reprogramming

The role of chromatin-associated RNAi components in the nucleus of mammalian cells and in particular in the context of developmental programs remains to be elucidated. Here, we investigate the function of nuclear Argonaute 1 (Ago1) in gene expression regulation during skeletal muscle differentiation. We show that Ago1 is required for activation of the myogenic program by supporting chromatin modification mediated by developmental enhancer activation. Mechanistically, we demonstrate that Ago1 directly controls global H3K27 acetylation (H3K27ac) by regulating enhancer RNA (eRNA)-CREB-binding protein (CBP) acetyltransferase interaction, a key step in enhancer-driven gene activation. In particular, we show that Ago1 is specifically required for myogenic differentiation 1 (MyoD) and downstream myogenic gene activation, whereas its depletion leads to failure of CBP acetyltransferase activation and blocking of the myogenic program. Our work establishes a role of the mammalian enhancer-associated RNAi component Ago1 in epigenome regulation and activation of developmental programs.

Enhancer RNA SLIT2 Inhibits Bone Metastasis of Breast Cancer Through Regulating P38 MAPK/c-Fos Signaling Pathway

BACKGROUND

Breast cancer (BRCA) is the most common cancer in women, while the bones are one of the most common sites of metastasis. Although new diagnostic methods or radiation or chemotherapies and targeted therapies have made huge advances, the occurrence of bone metastasis is also linked with poorer survival. Enhancer RNAs (eRNAs) have been demonstrated to participate in the progression of tumorigenesis and metastasis. However, the role of eRNAs in BRCA bone metastasis remains largely unclear.

METHOD

Gene expression profiling of 1,211 primary BRCA and 17 bone metastases samples were retrieved from The Cancer Genome Atlas (TCGA) database, and the significant prognostic eRNAs were identified by Cox regression and least absolute shrinkage and selection operator (LASSO) regression. The acceptable accuracy and discrimination of the nomogram were indicated by the receiver operating characteristic (ROC) and the calibration curves. Then target genes of eRNA, immune cell percentage by CIBERSORT analysis, immune genes by single-sample gene set enrichment analysis (ssGSEA), hallmark of cancer signaling pathway by gene set variation analysis (GSVA), and reverse phase protein array (RPPA) protein chip were used to build a co-expression regulation network and identified the key eRNAs in bone metastasis of BRCA. Finally, Cell Counting Kit-8 (CCK8) assay, cell cycle assay, and transwell assay were used to study changes in cell proliferation, migration, and invasiveness. Immunoprecipitation assay and Western blotting were used to test the interaction and the regulation signaling pathways.

RESULTS

The 27 hub eRNAs were selected, and a survival-related linear risk assessment model with a relatively high accuracy (area under curve (AUC): 0.726) was constructed. In addition, seven immune-related eRNAs (SLIT2, CLEC3B, LBPL1, FRY, RASGEF1B, DST, and ITIH5) as prognostic signatures for bone metastasis of BRCA were further confirmed by LASSO and multivariate Cox regression and CIBERSORT analysis. Finally, in vitro assay demonstrated that overexpression of SLIT2 reduced proliferation and metastasis in BRCA cells. Using high-throughput co-expression regulation network, we identified that SLIT2 may regulating P38 MAPK/c-Fos signaling pathway to promote the effects of metastasis.

CONCLUSION

Based on the co-expression network for bone metastasis of BRCA, we screened key eRNAs to explore a prognostic model in predicting the bone metastasis by bioinformatics analysis. Besides, we identified the potential regulatory signaling pathway of SLIT2 in BRCA bone metastasis, which provides a promising therapeutic strategy for metastasis of BRCA.

A Comprehensive Toolbox to Analyze Enhancer-Promoter Functions

Knowledge in gene transcription and chromatin regulation has been intensely studied for decades, but thanks to next-generation sequencing (NGS) techniques there has been a major leap forward in the last few years. Historically, identification of specific enhancer elements has led to the identification of master transcription factors (TFs) in the 1990s. Genetic and biochemical experiments have identified the key regulators controlling RNA polymerase II (RNAPII) transcription and structurally analyses have elucidated detailed mechanisms. NGS and the development of chromatin immunoprecipitation (ChIP) have accelerated the gain of knowledge in the recent years. By now, we have a dazzling wealth of techniques that are currently used to put gene expression into a genome-wide context. This book is an attempt to assemble useful protocols for many researchers within and nearby research areas. In general, these innovative techniques focus on enhancer and promoter studies. The techniques should also be of interest for related fields such as DNA repair and replication.

Transcriptionally active enhancers in human cancer cells

The growth of human cancer cells is driven by aberrant enhancer and gene transcription activity. Here, we use transient transcriptome sequencing (TT-seq) to map thousands of transcriptionally active putative enhancers in fourteen human cancer cell lines covering seven types of cancer. These enhancers were associated with cell type-specific gene expression, enriched for genetic variants that predispose to cancer, and included functionally verified enhancers. Enhancer-promoter (E-P) pairing by correlation of transcription activity revealed ~ 40,000 putative E-P pairs, which were depleted for housekeeping genes and enriched for transcription factors, cancer-associated genes, and 3D conformational proximity. The cell type specificity and transcription activity of target genes increased with the number of paired putative enhancers. Our results represent a rich resource for future studies of gene regulation by enhancers and their role in driving cancerous cell growth.

PAN RNA: transcriptional exhaust from a viral engine

Kaposi’s sarcoma-associated herpesvirus (KSHV), also designated human herpesvirus 8 (HHV-8), has been linked to Kaposi’s sarcoma, as well as to primary effusion lymphoma (PEL), and a subset of multicentric Castleman’s disease. KSHV genomes are maintained as episomes within infected cells and the virus exhibits a biphasic life cycle consisting of a life-long latent phase during which only a few viral genes are expressed and no viral progeny are produced and a transient lytic reactivation phase, in which a full repertoire of ~ 80 lytic genes are activated in a temporally regulated manner culminating in the release of new virions. Lytic replication is initiated by a single viral protein, K-Rta (ORF50), which activates more than 80 viral genes from multiple resident viral episomes (i.e., viral chromosomes). One of the major targets of K-Rta is a long non-coding nuclear RNA, PAN RNA (polyadenylated nuclear RNA), a lncRNA that accumulates to exceedingly high levels in the nucleus during viral reactivation. K-Rta directly binds to the PAN RNA promoter and robustly activates PAN RNA expression. Although PAN RNA has been known for over 20 years, its role in viral replication is still incompletely understood. In this perspective, we will briefly review the current understanding of PAN RNA and then describe our current working model of this RNA. The model is based on our observations concerning events that occur during KSHV lytic reactivation including (i) a marked accumulation of RNA Pol II at the PAN promoter, (ii) genomic looping emanating from the PAN locus, (iii) interaction of a second viral lytic protein (ORF57) with K-Rta, PAN RNA and RNA Pol II, (iv) the essential requirement for PAN RNA expression in cis for optimal transcriptional execution needed for the entire lytic program, and (v) ORF57 recruitment of RNA Pol II to the PAN genomic locus. Together our results generate a model in which the PAN locus serves as a hub for sequestration/trapping of the cellular transcriptional machinery proximal to viral episomes. Sequestration at the PAN locus facilitates high levels of viral transcription throughout the viral genome during lytic replication. ORF57 acts as a transcription-dependent transactivator at the PAN locus by binding to both Rta and PAN to locally trap RNA Pol II. The resulting accumulation of high levels of nuclear PAN RNA created by this process is an inducible enhancer-derived (eRNA) by-product that litters the infected cell nucleus.

Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs

Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.

Enhancer RNAs: Insights Into Their Biological Role

Enhancers play a central role in the transcriptional regulation of metazoans. Almost a decade ago, the discovery of their pervasive transcription into noncoding RNAs, termed enhancer RNAs (eRNAs), opened a whole new field of study. The presence of eRNAs correlates with enhancer activity; however, whether they act as functional molecules remains controversial. Here we review direct experimental evidence supporting a functional role of eRNAs in transcription and provide a general pipeline that could help in the design of experimental approaches to investigate the function of eRNAs. We propose that induction of transcriptional activity at enhancers promotes an increase in its activity by an RNA-mediated titration of regulatory proteins that can impact different processes like chromatin accessibility or chromatin looping. In a few cases, transcripts originating from enhancers have acquired specific molecular functions to regulate gene expression. We speculate that these transcripts are either nonannotated long noncoding RNAs (lncRNAs) or are evolving toward functional lncRNAs. Further work will be needed to comprehend better the biological activity of these transcripts.

Perturbing Enhancer Activity in Cancer Therapy

Tight regulation of gene transcription is essential for normal development, tissue homeostasis, and disease-free survival. Enhancers are distal regulatory elements in the genome that provide specificity to gene expression programs and are frequently misregulated in cancer. Recent studies examined various enhancer-driven malignant dependencies and identified different approaches to specifically target these programs. In this review, we describe numerous features that make enhancers good transcriptional targets in cancer therapy and discuss different approaches to overcome enhancer perturbation. Interestingly, a number of approved therapeutic agents, such as cyclosporine, steroid hormones, and thiazolidinediones, actually function by affecting enhancer landscapes by directly targeting very specific transcription factor programs. More recently, a broader approach to targeting deregulated enhancer programs has been achieved via Bromodomain and Extraterminal (BET) inhibition or perturbation of transcription-related cyclin-dependent kinases (CDK). One challenge to enhancer-targeted therapy is proper patient stratification. We suggest that monitoring of enhancer RNA (eRNA) expression may serve as a unique biomarker of enhancer activity that can help to predict and monitor responsiveness to enhancer-targeted therapies. A more thorough investigation of cancer-specific enhancers and the underlying mechanisms of deregulation will pave the road for an effective utilization of enhancer modulators in a precision oncology approach to cancer treatment.

Production of Spliced Long Noncoding RNAs Specifies Regions with Increased Enhancer Activity. Cell Syst. 2018;7:537-547.e3. [PMID: 30447999 DOI: 10.1016/j.cels.2018.10.009]

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.

Nanog Expression in Embryonic Stem Cells - An Ideal Model System to Dissect Enhancer Function

Embryonic stem cells (ESCs) are derived from the preimplantation embryo and can differentiate into virtually any other cell type (termed pluripotency), which is governed by lineage specific transcriptions factors (TFs) binding to cis regulatory elements (CREs) to mediate changes in gene expression. The reliance on transcriptional regulation to maintain pluripotency makes ESCs a valuable model to study the role of distal CREs such as enhancers in modulating gene expression to affect cell fate decisions. This review will highlight recent advance on transcriptional enhancers, focusing on studies performed in ESCs. In addition, we argue that the Nanog locus, which encodes for an ESC-critical TF, is particularly informative because it contains multiple co-regulated genes and enhancers in close proximity to one another. The unique landscape at Nanog permits the study of ongoing questions including whether multiple enhancers function additively versus synergistically, determinants of gene specificity, and cell-to-cell variability in gene expression.

Super-Enhancers at the Nanog Locus Differentially Regulate Neighboring Pluripotency-Associated Genes

Super-enhancers are tissue-specific cis-regulatory elements that drive expression of genes associated with cell identity and malignancy. A cardinal feature of super-enhancers is that they are transcribed to produce enhancer-derived RNAs (eRNAs). It remains unclear whether super-enhancers robustly activate genes in situ and whether their functions are attributable to eRNAs or the DNA element. CRISPR/Cas9 was used to systematically delete three discrete super-enhancers at the Nanog locus in embryonic stem cells, revealing functional differences in Nanog transcriptional regulation. One distal super-enhancer 45 kb upstream of Nanog (-45 enhancer) regulates both nearest neighbor genes, Nanog and Dppa3. Interestingly, eRNAs produced at the -45 enhancer specifically regulate Dppa3 expression by stabilizing looping of the -45 enhancer and Dppa3. Our work illustrates that genomic editing is required to determine enhancer function and points to a method to selectively target a subset of super-enhancer-regulated genes by depleting eRNAs.

Noncoding Regulatory RNAs in Hematopoiesis

Hematopoiesis is a dynamic process in which blood cells are continuously generated from hematopoietic stem cells (HSCs). The regulatory mechanisms controlling HSC fate have been studied extensively over the past several decades. Although many protein-coding genes have been shown to regulate hematopoietic differentiation, additional levels of HSC regulation are not well studied. Advances in deep sequencing have revealed many new classes of regulatory noncoding RNAs (ncRNAs), such as enhancer RNAs and antisense ncRNAs. Functional analysis of some of these ncRNAs has provided insights into the molecular mechanisms that regulate hematopoietic development and disease. In this review, we summarize recent advances in our understanding of functional regulatory ncRNAs associated with hematopoietic self-renewal and differentiation, as well as those dysregulated ncRNAs involved in hematologic malignancies.

The many faces of long noncoding RNAs

Over the past few years, the field of noncoding RNAs has grown from a niche for geneticists into a prominent domain of mainstream biology. Advances in genomic technologies have provided a more comprehensive view of the mammalian genome, improving our knowledge of regions of the genome devoid of protein-coding potential. A large body of evidence supports the proposal that noncoding RNAs account for a large proportion of the transcriptional output of any given cell and tissue type. This review will delve into the biogenesis and function of long noncoding RNAs. We will discuss our current understanding of these molecules as major chromatin players, and explore future directions in the field.