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.

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.

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.

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.

BACKGROUND: Recently, enhancer RNAs (eRNAs) have garnered attention as pivotal biomarkers for the onset and progression of cancer. However, the landscape of eRNAs and the implications of eRNA-based molecular subtypes in stage II/III colorectal cancer (CRC) remain largely unexplored.
METHODS: Comprehensive profiling of eRNAs was conducted on a public stage II/III CRC cohort with total RNA-seq data. We used unsupervised clustering of prognostic eRNAs to establish an eRNA-based subtyping system. Further evaluations included molecular characteristics, immune infiltration, clinical outcomes, and drug responses. Finally, we validated the eRNA-based subtyping system in The Cancer Genome Atlas (TCGA) CRC cohort.
RESULTS: We identified a total of 6453 expressed eRNAs, among which 237 were prognostic. A global upregulation of eRNAs was observed in microsatellite-stable (MSS) CRCs when compared to microsatellite instability-high (MSI-H) CRCs. Through consensus clustering, two novel molecular subtypes, termed Cluster 1(C1) and Cluster 2(C2), were further identified. C1, associated with the activation of epithelial-mesenchymal transition (EMT), hypoxia, and KRAS signaling pathways, showed poorer prognosis. C2, correlated with the canonical CRC subtype, exhibited superior survival outcomes. In addition, C1 showed enrichment with immune infiltration and more sensitivity to immune checkpoint inhibitors.
CONCLUSIONS: Our study unravels the molecular heterogeneity of stage II/III CRC at the eRNA level and highlights the potential applications of the novel eRNA-based subtyping system in predicting prognosis and guiding immunotherapy.

Abundant evidence suggests that enhancer RNA (eRNA) is closely related to tumorigenesis, and the role of eRNA transcription in promoting genomic instability in cancers is gradually unveiled. However, research on the evaluation of the prognostic value and molecular mechanisms of genomic instability associated eRNAs in breast cancer is long overdue. Here, we integratively analyzed eRNA expression and somatic mutation profiles in breast cancer genome. We identified genomic instability associated eRNAs and developed a prognostic signature based on these eRNAs with the area under the curve (AUC) around 0.8 at 9-year survival. We further found the prognostic value of this signature is independent of common clinical factors and is better than TP53 status. Higher expression of genomic instability associated genes in the high-risk group was observed, suggesting that this eRNA signature may serve as an indicator of genomic instability in breast cancer. We found prognostic eRNA co-expressed genes are mainly enriched in Gene set \'Breast Cancer 8P12-P11 Amplicon\', Gene set \'Metabolism of lipids\' and GO process \'Ubiquitin protein ligase binding\'. Furthermore, 11 eRNA-signature regulated genes are identified by assessing promoter-enhancer interaction. Among these genes, F11R, BHLHE40, and NECTIN4 are previously reported oncogenes and EGOT is a tumor suppressor gene, indicating the direct roles of eRNAs in tumorigenesis.

KRAS is the most commonly mutated oncogene in human cancer, and mutant KRAS is responsible for over 90% of pancreatic ductal adenocarcinoma (PDAC), the most lethal cancer. Here, we show that RNA polymerase II-associated factor 1 complex (PAF1C) is specifically required for survival of PDAC but not normal adult pancreatic cells. We show that PAF1C maintains cancer cell genomic stability by restraining overaccumulation of enhancer RNAs (eRNAs) and promoter upstream transcripts (PROMPTs) driven by mutant Kras. Loss of PAF1C leads to cancer-specific lengthening and accumulation of pervasive transcripts on chromatin and concomitant aberrant R-loop formation and DNA damage, which, in turn, trigger cell death. We go on to demonstrate that the global transcriptional hyperactivation driven by Kras signaling during tumorigenesis underlies the specific demand for PAF1C by cancer cells. Our work provides insights into how enhancer transcription hyperactivation causes general transcription factor addiction during tumorigenesis.

OBJECTIVE: Colon cancer (CC) is a cancer of the large intestine with high prevalence and poor prognosis. enhancer RNAs. Therefore, valuable tools or biomarkers for predicting patient status, directing clinical practice, and reducing overtreatment are needed. Enhancer RNAs (eRNAs), a class of noncoding RNAs transcribed from enhancers, have been shown to function as regulators of oncogene or tumor suppressor gene expression. The aim of our study was to explore the potential roles of eRNAs and their target enhancer-related genes (ERGs) in the prognosis of CC.
METHODS: Selected CC cases (stage I-III) from The Cancer Genome Atlas database were used as a training set, and cases from the Gene Expression Omnibus were used as the validation set. ERGs associated with prognosis were screened through three steps: potential, candidate, and prognosis ERGs. Multivariate Cox proportional hazards analysis was used to identify independent prognostic factors, and a nomogram was created. Calibration curves were drawn by comparing predicted and observed survival probability. For validation, the calibration curves and ROC analysis were also applied to two external validation sets. The biological significance and clinical application of the genes obtained were investigated.
RESULTS: Based on the multiple tiers of strict screening, 11 prognostic ERGs were obtained, which were combined to obtain a prognosis signature. A compound nomogram integrating age, TNM classification, and the prognostic signature was constructed. The model was reliable in distinguishing the risk of patients with stage I-III CC, with AUCs of 0.78 and 0.70 at 5 and 7 years, respectively. There was good reproducibility in calibration curves. The prognostic model also yielded good prediction capability in the validation sets.
CONCLUSIONS: In this study, the usefulness and specificity of the ERGs in prognosis were described, which should be considered a key feature in the clinical guidance of CC patients with early stage. We concluded that the major implications of the eRNAs and ERGs should be valued, which would be an emerging hallmark in the prognosis of cancer.