A certain degree of chromatin openness is necessary for the activity of transcription-regulating regions within the genome, facilitating accessibility to RNA polymerases and subsequent synthesis of regulatory element RNAs (regRNAs) from these regions. The rapidly increasing number of studies underscores the significance of regRNAs across diverse cellular processes and diseases, challenging the paradigm that these transcripts are non-functional transcriptional noise. This review explores the multifaceted roles of regRNAs in human cells, encompassing rather well-studied entities such as promoter RNAs and enhancer RNAs (eRNAs), while also providing insights into overshadowed silencer RNAs and insulator RNAs. Furthermore, we assess notable examples of shorter regRNAs, like miRNAs, snRNAs, and snoRNAs, playing important roles. Expanding our discourse, we deliberate on the potential usage of regRNAs as biomarkers and novel targets for cancer and other human diseases.

Recent evidence has highlighted the functions of enhancers in modulating transcriptional machinery and affecting the development of human diseases including rheumatoid arthritis (RA). Enhancer RNAs (eRNAs) are RNA molecules transcribed from active enhancer regions. This study investigates the specific function of eRNA in gene transcription and osteoclastogenesis in RA. Regulator of G protein signaling 1 (RGS1)-associated eRNA was highly activated in osteoclasts according to bioinformatics prediction. RGS1 mRNA was increased in mice with collagen-induced arthritis as well as in M-CSF/soluble RANKL-stimulated macrophages (derived from monocytes). This was ascribed to increased RGS1 eRNA activity. Silencing of 5\'-eRNA blocked the binding between forkhead box J3 (FOXJ3) and the RGS1 promoter, thus suppressing RGS1 transcription. RGS1 accelerated osteoclastogenesis through PLC-IP3R-dependent Ca2+ response. Knockdown of either FOXJ3 or RGS1 ameliorated arthritis severity, improved pathological changes, and reduced osteoclastogenesis and bone erosion in vivo and in vitro. However, the effects of FOXJ3 silencing were negated by RGS1 overexpression. In conclusion, this study demonstrates that the RGS1 eRNA-driven transcriptional activation of the FOXJ3/RGS1 axis accelerates osteoclastogenesis through PLC-IP3R dependent Ca2+ response in RA. The finding may offer novel insights into the role of eRNA in gene transcription and osteoclastogenesis in RA.

The Myb proto-oncogene encodes the transcription factor c-MYB, which is critical for hematopoiesis. Distant enhancers of Myb form a hub of interactions with the Myb promoter. We identified a long non-coding RNA (Myrlin) originating from the -81-kb murine Myb enhancer. Myrlin and Myb are coordinately regulated during erythroid differentiation. Myrlin TSS deletion using CRISPR-Cas9 reduced Myrlin and Myb expression and LDB1 complex occupancy at the Myb enhancers, compromising enhancer contacts and reducing RNA Pol II occupancy in the locus. In contrast, CRISPRi silencing of Myrlin left LDB1 and the Myb enhancer hub unperturbed, although Myrlin and Myb expressions were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the KMT2A/MLL1 complex. Myrlin CRISPRi compromised KMT2A occupancy in the Myb locus, decreasing CDK9 and RNA Pol II binding and resulting in Pol II pausing in the Myb first exon/intron. Thus, Myrlin directly participates in activating Myb transcription by recruiting KMT2A.

MYB is a key regulator of hematopoiesis and erythropoiesis, and dysregulation of MYB is closely involved in the development of leukemia, however the mechanism of MYB regulation remains still unclear so far. Our previous study identified a long noncoding RNA (lncRNA) derived from the -34 kb enhancer of the MYB locus, which can promote MYB expression, the proliferation and migration of human leukemia cells, and is therefore termed MY34UE-AS. Then the interacting partner proteins of MY34UE-AS were identified and studied in the present study. hnRNPA0 was identified as a binding partner of MY34UE-AS through RNA pulldown assay, which was further validated through RNA immunoprecipitation (RIP). hnRNPA0 interacted with MY34UE-AS mainly through its RRM2 domain. hnRNPA0 overexpression upregulated MYB and increased the proliferation and migration of K562 cells, whereas hnRNPA0 knockdown showed opposite effects. Rescue experiments showed MY34UE-AS was required for above mentioned functions of hnRNPA0. These results reveal that hnRNPA0 is involved in leukemia through upregulating MYB expression by interacting with MY34UE-AS, suggesting that the hnRNPA0/MY34UE-AS axis could serve as a potential target for leukemia treatment.

BACKGROUND: Dysregulation of enhancer transcription occurs in multiple cancers. Enhancer RNAs (eRNAs) are transcribed products from enhancers that play critical roles in transcriptional control. Characterizing the genetic basis of eRNA expression may elucidate the molecular mechanisms underlying cancers.
METHODS: Initially, a comprehensive analysis of eRNA quantitative trait loci (eRNAQTLs) was performed in The Cancer Genome Atlas (TCGA), and functional features were characterized using multi-omics data. To establish the first eRNAQTL profiles for colorectal cancer (CRC) in China, epigenomic data were used to define active enhancers, which were subsequently integrated with transcription and genotyping data from 154 paired CRC samples. Finally, large-scale case-control studies (34,585 cases and 69,544 controls) were conducted along with multipronged experiments to investigate the potential mechanisms by which candidate eRNAQTLs affect CRC risk.
RESULTS: A total of 300,112 eRNAQTLs were identified across 30 different cancer types, which exert their influence on eRNA transcription by modulating chromatin status, binding affinity to transcription factors and RNA-binding proteins. These eRNAQTLs were found to be significantly enriched in cancer risk loci, explaining a substantial proportion of cancer heritability. Additionally, tumor-specific eRNAQTLs exhibited high responsiveness to the development of cancer. Moreover, the target genes of these eRNAs were associated with dysregulated signaling pathways and immune cell infiltration in cancer, highlighting their potential as therapeutic targets. Furthermore, multiple ethnic population studies have confirmed that an eRNAQTL rs3094296-T variant decreases the risk of CRC in populations from China (OR = 0.91, 95%CI 0.88-0.95, P = 2.92 × 10-7) and Europe (OR = 0.92, 95%CI 0.88-0.95, P = 4.61 × 10-6). Mechanistically, rs3094296 had an allele-specific effect on the transcription of the eRNA ENSR00000155786, which functioned as a transcriptional activator promoting the expression of its target gene SENP7. These two genes synergistically suppressed tumor cell proliferation. Our curated list of variants, genes, and drugs has been made available in CancereRNAQTL ( http://canernaqtl.whu.edu.cn/#/ ) to serve as an informative resource for advancing this field.
CONCLUSIONS: Our findings underscore the significance of eRNAQTLs in transcriptional regulation and disease heritability, pinpointing the potential of eRNA-based therapeutic strategies in cancers.

Medullary thyroid carcinoma (MTC) is a rare but aggressive endocrine malignancy that originates from the parafollicular C cells of the thyroid gland. Enhancer RNAs (eRNAs) are non-coding RNAs transcribed from enhancer regions, which are critical regulators of tumorigenesis. However, the roles and regulatory mechanisms of eRNAs in MTC remain poorly understood. This study aims to identify key eRNAs regulating the malignant phenotype of MTC and to uncover transcription factors involved in the regulation of key eRNAs.
GSE32662 and GSE114068 were used for the identification of differentially expressed genes, eRNAs, enhancers and enhancer-regulated genes in MTC. Metascape and the transcription factor affinity prediction method were used for gene function enrichment and transcription factor prediction, respectively. qRT-PCR was used to detect gene transcription levels. ChIP-qPCR was used to assess the binding of histone H3 lysine 27 acetylation (H3K27ac)-enriched regions to anti- H3K27ac. RIP-qPCR was used to detect the binding between FOXQ1 and LINC00887. CCK8 and Transwell were performed to measure the proliferation and invasion of MTC cells, respectively. Intracellular reactive oxygen species (ROS) levels were quantified using a ROS assay kit.
Four eRNAs (H1FX-AS1, LINC00887, MCM3AP-AS1 and A1BG-AS1) were screened, among which LINC00887 was the key eRNA promoting the proliferation and invasion of MTC cells. A total of 135 genes controlled by LINC00887-regulated enhancers were identified; among them, BCL2, PRDX1, SFTPD, TPO, GSS, RAD52, ZNF580, and ZFP36L1 were significantly enriched in the \"ROS metabolic process\" term. As a transcription factor regulating genes enriched in the \"ROS metabolic process\" term, FOXQ1 could recruit LINC00887. Overexpression of FOXQ1 restored LINC00887 knockdown-induced downregulation of GSS and ZFP36L1 transcription in MTC cells. Additionally, FOXQ1 overexpression counteracted the inhibitory effects of LINC00887 knockdown on the proliferation and invasion of MTC cells and the promotion of intracellular ROS accumulation induced by LINC00887 knockdown.
LINC00887 was identified as a key eRNA promoting the malignant phenotype of MTC cells. The involvement of FOXQ1 was essential for LINC00887 to play a pro-tumorigenic role in MTC. Our findings suggest that the FOXQ1/LINC00887 axis is a potential therapeutic target for MTC.

The nuclear receptor, farnesoid X receptor (FXR/NR1H4), is increasingly recognized as a promising drug target for metabolic diseases, including nonalcoholic steatohepatitis (NASH). Protein-coding genes regulated by FXR are well known, but whether FXR also acts through regulation of long non-coding RNAs (lncRNAs), which vastly outnumber protein-coding genes, remains unknown. Utilizing RNA-seq and global run-on sequencing (GRO-seq) analyses in mouse liver, we found that FXR activation affects the expression of many RNA transcripts from chromatin regions bearing enhancer features. Among these we discovered a previously unannotated liver-enriched enhancer-derived lncRNA (eRNA), termed FXR-induced non-coding RNA (Fincor). We show that Fincor is specifically induced by the hammerhead-type FXR agonists, including GW4064 and tropifexor. CRISPR/Cas9-mediated liver-specific knockdown of Fincor in dietary NASH mice reduced the beneficial effects of tropifexor, an FXR agonist currently in clinical trials for NASH and primary biliary cholangitis (PBC), indicating that amelioration of liver fibrosis and inflammation in NASH treatment by tropifexor is mediated in part by Fincor. Overall, our findings highlight that pharmacological activation of FXR by hammerhead-type agonists induces a novel eRNA, Fincor, contributing to the amelioration of NASH in mice. Fincor may represent a new drug target for addressing metabolic disorders, including NASH.
Non-alcoholic steatohepatitis, also known as NASH, is a severe condition whereby fat deposits around the liver lead to inflammation, swelling, scarring and lasting damage to the organ. Despite being one of the leading causes of liver-related deaths worldwide, the disease has no approved treatment. A protein known as Farnesoid X receptor (or FXR) is increasingly being recognized as a promising drug target for non-alcoholic steatohepatitis. Once activated, FXR helps to regulate the activity of DNA regions which are coding for proteins important for liver health. However, less is known about how FXR may act on non-coding regions, the DNA sequences that do not generate proteins but can be transcribed into RNA molecules with important biological roles. In response, Chen et al. investigated whether FXR activation of non-coding RNAs could be linked to the clinical benefits of hammerhead FXR agonists, a type of synthetic compounds that activates this receptor. To do so, genetic analyses of mouse livers were performed to identify non-coding RNAs generated when FXR was activated by the agonist. These experiments revealed that agonist-activated FXR induced a range of non-coding RNAs transcribed from DNA sequences known as enhancers, which help to regulate gene expression. In particular, hammerhead FXR agonists led to the production of a liver-specific enhancer RNA called Fincor. Additional experiments using tropifexor, a hammerhead FXR agonist currently into clinical trials, showed that this investigational new drug had reduced benefits in a mouse model of non-alcoholic steatohepatitis with low Fincor levels. This suggested that this enhancer RNA may play a key role in mediating the clinical benefits of hammerhead FXR agonists, encouraging further research into its role and therapeutic value.

BACKGROUND: We aimed to comprehensively analyze the clinical value of immune-related eRNAs-driven genes in lung adenocarcinoma (LUAD) and find the potential biomarkers for prognosis and therapeutic response to improve the survival of this malignant disease.
METHODS: Pearson\'s correlation analysis was performed to identify the immune-related eRNAs-driven genes. Cox regression and least absolute shrinkage and selection operator (LASSO) analyses were used to construct this prognostic risk signature. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were used to investigate the underlying molecular mechanism. The single sample gene set enrichment analysis (ssGSEA) algorithm was conducted to evaluate the immune status based on the signature. The quantitative real-time PCR (qRT-PCR) analysis was performed to evaluate the expression value of the signature genes between LUAD tissues and adjacent lung tissues.
RESULTS: Five immune-related eRNAs-driven genes (SHC1, GDF10, CCL14, FYN, and NOD1) were identified to construct a prognostic risk signature with favorable predictive capacity. The patients with high-risk scores based on the signature were significantly associated with the malignant clinical features compared with those with low-risk scores. Kaplan-Meier analysis demonstrated that the sample in the low-risk group had a prolonged survival compared with those in the high-risk group. This risk signature was validated to have a promising predictive capacity and reliability in diverse clinical situations and independent cohorts. The functional enrichment analysis demonstrated that humoral immune response and intestinal immune network for IgA production pathway might be the underlying molecular mechanism related to the signature. The proportion of the vast majority of immune infiltrating cells in the high-risk group was significantly lower than that in the low-risk group, and the immunotherapy response rate in the low-risk group was significantly higher than that in the high-risk group. Moreover, BI-2536, sepantronium bromide, and ULK1 were the potential drugs for the treatment of patients with higher risk scores. Finally, the experiment in vivo and database analysis indicated that CCL14, FYN, NOD1, and GDF10 are the potential LUAD suppressor and SHC1 is a potential treatment target for LUAD.
CONCLUSIONS: Above all, we constructed a prognostic risk signature with favorable predictive capacity in LUAD, which was significantly associated with malignant features, immunosuppressive tumor microenvironment, and immunotherapy response and may provide clinical benefit in clinical decisions.

Single-cell RNA-sequencing (scRNA-seq) provides unprecedented insights into cellular heterogeneity. Although scRNA-seq reads from most prevalent and popular tagged-end protocols are expected to arise from the 3\' end of polyadenylated RNAs, recent studies have shown that \"off-target\" reads can constitute a substantial portion of the read population. In this work, we introduced scCensus, a comprehensive analysis workflow for systematically evaluating and categorizing off-target reads in scRNA-seq. We applied scCensus to seven scRNA-seq datasets. Our analysis of intergenic reads shows that these off-target reads contain information about chromatin structure and can be used to identify similar cells across modalities. Our analysis of antisense reads suggests that these reads can be used to improve gene detection and capture interesting transcriptional activities like antisense transcription. Furthermore, using splice-aware quantification, we find that spliced and unspliced reads provide distinct information about cell clusters and biomarkers, suggesting the utility of integrating signals from reads with different splicing statuses. Overall, our results suggest that off-target scRNA-seq reads contain underappreciated information about various transcriptional activities. These observations about yet-unexploited information in existing scRNA-seq data will help guide and motivate the community to improve current algorithms and analysis methods, and to develop novel approaches that utilize off-target reads to extend the reach and accuracy of single-cell data analysis pipelines.

Immune-related enhancer RNAs (eRNAs) have garnered significant attention in cancer metabolism research, yet their specific roles in ccRCC have remained elusive.
We retrieved eRNA expression profiles from TCGA database and identified immune-related eRNAs (IREs) by assessing their co-expression with immune genes. Utilizing consensus clustering, we organized these IREs into two distinct clusters. The construction of an IREs signature was accomplished through the LASSO and multivariate Cox analysis. Furthermore, we performed Cell Counting Kit-8 and clonogenic assays to assess changes in the proliferative capacity of Caki-1 and 769-P cells.
The existence of two clusters of immune-related eRNAs in ccRCC, each with distinctive prognostic and immunological attributes. Cluster B exhibited immunosuppressive properties and displayed a positive correlation with immunosuppressive cells. Functional enrichment analysis unveiled their involvement in several tumor-promoting pathways, metabolic pathways and immune pathways. The IREs signature demonstrated its potential to accurately predict patient immune and prognostic characteristics. AC003092.1, an eRNA strongly associated with patient survival, emerged as a potential oncogene significantly linked to adverse prognosis and the presence of immunosuppressive cells and checkpoints in ccRCC patients. Notably, AC003092.1 displayed marked upregulation in ccRCC tissues and cell lines, and its knockdown substantially inhibited the proliferation of Caki-1 and 769-P cells.
We established a robust predictive model that played a vital role in determining the prognosis, clinicopathological characteristics and immune cell infiltration patterns of ccRCC patients. IRE, particularly AC003092.1, which was strongly associated with survival, hold promise as novel immunotherapeutic targets for ccRCC.