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.

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.

Significant advances have been achieved in understanding the critical role of enhancer RNAs (eRNAs) in the complex field of gene regulation. However, notable uncertainty remains concerning the biology of eRNAs, highlighting the need for continued research to uncover their exact functions in cellular processes and diseases. We present a comprehensive study to scrutinize mutation density patterns, mutation strand bias, and mutation burden in eRNAs across multiple cancer types. Our findings reveal that eRNAs exhibit mutation strand bias akin to that observed in protein-coding RNAs. We also identified a novel pattern, in which mutation density is notably diminished around the central region of the eRNA, but conspicuously elevated towards both the beginning and end. This pattern can be potentially explained by a mechanism involving heightened transcriptional activity and the activation of transcription-coupled repair. The central regions of the eRNAs appear to be more conserved, hinting at a potential mechanism preserving their structural and functional integrity, while the extremities may be more susceptible to mutations due to increased exposure. The evolutionary trajectory of this mutational pattern suggests a nuanced adaptation in eRNAs, where stability at their core coexists with flexibility at their extremities, potentially facilitating their diverse interactions with other genetic entities.

Gene transcription is controlled and modulated by regulatory regions, including enhancers and promoters. These regions are abundant in unstable, non-coding bidirectional transcription. Using nascent RNA transcription data across hundreds of human samples, we identified over 800,000 regions containing bidirectional transcription. We then identify highly correlated transcription between bidirectional and gene regions. The identified correlated pairs, a bidirectional region and a gene, are enriched for disease associated SNPs and often supported by independent 3D data. We present these resources as an SQL database which serves as a resource for future studies into gene regulation, enhancer associated RNAs, and transcription factors.

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 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 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 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.

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 expression were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the MLL1 complex. Myrlin CRISPRi compromised MLL1 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 MLL1.

BACKGROUND: Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states.
METHODS: We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing.
RESULTS: We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease.
CONCLUSIONS: This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.