Noncoding transcription at enhancers: general principles and functional models

Downregulation of enhancer RNA AC003092.1 is associated with poor prognosis in kidney renal clear cell carcinoma

Kidney renal clear cell carcinoma (KIRC) is the most common histological type of renal cancer, enhancer RNA plays a significant role in tumor growth, however, it has been less studied in renal cancer. The aim of this study was to investigate the role of eRNA AC003092.1 in KIRC. Clinical and RNA expression data were downloaded from a TCGA database, and performed bioinformatics analysis, including expression level analysis, survival analysis, clinical correlation analysis, immune correlation analysis. We further confirmed the expression level of AC003092.1 between normal and tumor cell, predicted the biological role of AC003092.1 in KIRC, and performed cell proliferation and wound healing assays, followed by GSEA enrichment analysis and western blot to detect the proteins of the enriched pathway. Bioinformatics results showed that AC003092.1 expression was elevated in tumor tissues, and knockdown of AC003092.1 expression inhibited cell proliferation and migration. GSEA and western blot results showed that knockdown AC003092.1 expression alleviated the extracellular matrix (ECM) process in KIRC cell lines. Our study provides evidence that AC003092.1 play an important role in KIRC, and AC003092.1 may promote tumor cell progression by affecting the ECM process during tumor development.

Porcine cis-acting lnc-CAST positively regulates CXCL8 expression through histone H3K27ac

The chemokine CXCL8, also known as the neutrophil chemotactic factor, plays a crucial role in mediating inflammatory responses and managing cellular immune reactions during viral infections. Porcine reproductive and respiratory syndrome virus (PRRSV) primarily infects pulmonary alveolar macrophages (PAMs), leading to acute pulmonary infections. In this study, we explored a novel long non-coding RNA (lncRNA), termed lnc-CAST, situated within the Cxcl8 gene locus. This lncRNA was found to be highly expressed in porcine macrophages. We observed that both lnc-CAST and CXCL8 were significantly upregulated in PAMs following PRRSV infection, and after treatments with lipopolysaccharide (LPS) or lipoteichoic acid (LTA). Furthermore, we noticed a concurrent upregulation of lnc-CAST and CXCL8 expression in lungs of PRRSV-infected pigs. We then determined that lnc-CAST positively influenced CXCL8 expression in PAMs. Overexpression of lnc-CAST led to an increase in CXCL8 production, which in turn enhanced the migration of epithelial cells and the recruitment of neutrophils. Conversely, inhibiting lnc-CAST expression resulted in reduced CXCL8 production in PAMs, leading to decreased migration levels of epithelial cells and neutrophils. From a mechanistic perspective, we found that lnc-CAST, localized in the nucleus, facilitated the enrichment of histone H3K27ac in CXCL8 promoter region, thereby stimulating CXCL8 transcription in a cis-regulatory manner. In conclusion, our study underscores the pivotal critical role of lnc-CAST in regulating CXCL8 production, offering valuable insights into chemokine regulation and lung damage during PRRSV infection.

A novel piwi-interacting RNA associates with type 2-high asthma phenotypes

BACKGROUND

Piwi-interacting RNAs (piRNAs), comprising the largest noncoding RNA group, regulate transcriptional processes. Whether piRNAs are associated with type 2 (T2)-high asthma is unknown.

OBJECTIVE

We sought to investigate the association between piRNAs and T2-high asthma in childhood asthma.

METHODS

We sequenced plasma samples from 462 subjects in the Childhood Asthma Management Program (CAMP) as the discovery cohort and 1165 subjects in the Genetics of Asthma in Costa Rica Study (GACRS) as a replication cohort. Sequencing reads were filtered first, and piRNA reads were annotated and normalized. Linear regression was used for the association analysis of piRNAs and peripheral blood eosinophil count, total serum IgE level, and long-term asthma exacerbation in children with asthma. Mediation analysis was performed to investigate the effect direction. We then ascertained if the circulating piRNAs were present in asthmatic airway epithelial cells in a Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo) public data set.

RESULTS

Fifteen piRNAs were significantly associated with eosinophil count in CAMP (P ≤ .05), and 3 were successfully replicated in GACRS. Eleven piRNAs were associated with total IgE in CAMP, and one of these was replicated in GACRS. All 22 significant piRNAs were identified in epithelial cells in vitro, and 6 of these were differentially expressed between subjects with asthma and healthy controls. Fourteen piRNAs were associated with long-term asthma exacerbation, and effect of piRNAs on long-term asthma exacerbation are mediated through eosinophil count and serum IgE level.

CONCLUSION

piRNAs are associated with peripheral blood eosinophils and total serum IgE in childhood asthma and may play important roles in T2-high asthma.

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

The Role of Non-coding RNAs in Alzheimer's Disease: Pathogenesis, Novel Biomarkers, and Potential Therapeutic Targets

Long non-coding RNAs (IncRNAs) are regulatory RNA transcripts that have recently been associated with the onset of many neurodegenerative illnesses, including Alzheimer's disease (AD). Several IncRNAs have been found to be associated with AD pathophysiology, each with a distinct mechanism. In this review, we focused on the role of IncRNAs in the pathogenesis of AD and their potential as novel biomarkers and therapeutic targets. Searching for relevant articles was done using the PubMed and Cochrane library databases. Studies had to be published in full text in English in order to be considered. Some IncRNAs were found to be upregulated, while others were downregulated. Dysregulation of IncRNAs expression may contribute to AD pathogenesis. Their effects manifest as the synthesis of beta-amyloid (Aβ) plaques increases, thereby altering neuronal plasticity, inducing inflammation, and promoting apoptosis. Despite the need for more investigations, IncRNAs could potentially increase the sensitivity of early detection of AD. Until now, there has been no effective treatment for AD. Hence, InRNAs are promising molecules and may serve as potential therapeutic targets. Although several dysregulated AD-associated lncRNAs have been discovered, the functional characterization of most lncRNAs is still lacking.

Chromatin structure and 3D architecture define differential functions of PU.1 cis regulatory elements in human blood cell lineages

The precise spatio-temporal expression of the hematopoietic ETS transcription factor PU.1 that determines the hematopoietic cell fates is tightly regulated at the chromatin level. However, it remains elusive as to how chromatin signatures are linked to this dynamic expression pattern of PU.1 across blood cell lineages. Here we performed an unbiased and in-depth analysis of the relationship between human PU.1 expression, the presence of trans-acting factors, and 3D architecture at various cis-regulatory elements (CRE) proximal to the PU.1 locus. We identified multiple novel CREs at the upstream region of the gene following an integrative inspection for conserved DNA elements at the chromatin-accessible regions in primary human blood lineages. We showed that a subset of CREs localize within a 10 kb-wide cluster that exhibits that exhibit molecular features of a myeloid-specific super-enhancer involved in mediating PU.1 autoregulation, including open chromatin, unmethylated DNA, histone enhancer marks, transcription of enhancer RNAs, and occupancy of the PU.1 protein itself. Importantly, we revealed the presence of common 35-kb-wide CTCF-bound insulated neighborhood that contains the CRE cluster, forming the chromatin territory for lineage-specific and CRE-mediated chromatin interactions. These include functional CRE-promoter interactions in myeloid and B cells but not in erythroid and T cells. Our findings also provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in defining certain CREs as enhancers or silencers in chromatin regulation of PU.1 expression. The study lays the groundwork for further examination of PU.1 CREs as well as epigenetic regulation in malignant hematopoiesis.

Comprehensive Pan-Cancer Mutation Density Patterns in Enhancer RNA

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.

Noncoding RNA mutations in cancer

Cancer is driven by both germline and somatic genetic changes. Efforts have been devoted to characterizing essential genetic variations in cancer initiation and development. Most attention has been given to mutations in protein-coding genes and associated regulatory elements such as promoters and enhancers. The development of sequencing technologies and in silico and experimental methods has allowed further exploration of cancer predisposition variants and important somatic mutations in noncoding RNAs, mainly for long noncoding RNAs and microRNAs. Association studies including GWAS have revealed hereditary variations including SNPs and indels in lncRNA or miRNA genes and regulatory regions. These mutations altered RNA secondary structures, expression levels, and target recognition and then conferred cancer predisposition to carriers. Whole-exome/genome sequencing comparing cancer and normal tissues has revealed important somatic mutations in noncoding RNA genes. Mutation hotspots and somatic copy number alterations have been identified in various tumor-associated noncoding RNAs. Increasing focus and effort have been devoted to studying the noncoding region of the genome. The complex genetic network of cancer initiation is being unveiled. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Deciphering the Functional Long Non-Coding RNAs Derived from MicroRNA Loci

Albeit the majority of eukaryotic genomes can be pervasively transcribed to a diverse population of lncRNAs and various subtypes of lncRNA are discovered. However, the genome-wide study of miRNA-derived lncRNAs is still lacking. Here, it is reported that over 800 miRNA gene-originated lncRNAs (molncRNAs) are generated from miRNA loci. One of them, molnc-301b from miR-301b and miR-130b, functions as an "RNA decoy" to facilitate dissociation of the chromatin remodeling protein SMARCA5 from chromatin and thereby sequester transcription and mRNA translation. Specifically, molnc-301b attenuates erythropoiesis by mitigating the transcription of erythropoietic and translation-associated genes, such as GATA1 and FOS. In addition, a useful and powerful CRISPR screen platform to characterize the biological functions of molncRNAs at large-scale and single-cell levels is established and 29 functional molncRNAs in hematopoietic cells are identified. Collectively, the focus is on miRNA-derived lncRNAs, deciphering their landscape during normal hematopoiesis, and comprehensively evaluating their potential roles.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

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.

G4access identifies G-quadruplexes and their associations with open chromatin and imprinting control regions

Metazoan promoters are enriched in secondary DNA structure-forming motifs, such as G-quadruplexes (G4s). Here we describe 'G4access', an approach to isolate and sequence G4s associated with open chromatin via nuclease digestion. G4access is antibody- and crosslinking-independent and enriches for computationally predicted G4s (pG4s), most of which are confirmed in vitro. Using G4access in human and mouse cells, we identify cell-type-specific G4 enrichment correlated with nucleosome exclusion and promoter transcription. G4access allows measurement of variations in G4 repertoire usage following G4 ligand treatment, HDAC and G4 helicases inhibitors. Applying G4access to cells from reciprocal hybrid mouse crosses suggests a role for G4s in the control of active imprinting regions. Consistently, we also observed that G4access peaks are unmethylated, while methylation at pG4s correlates with nucleosome repositioning on DNA. Overall, our study provides a new tool for studying G4s in cellular dynamics and highlights their association with open chromatin, transcription and their antagonism to DNA methylation.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

G-quadruplexes from non-coding RNAs

Non-coding RNAs (ncRNAs) are significant regulators of gene expression in a wide range of biological processes, such as transcription, RNA maturation, or translation. ncRNAs interplay with proteins or other RNAs through not only classical sequence-based mechanisms but also unique higher-order structures such as RNA G-quadruplexes (rG4s). rG4s are predictably formed in guanine-rich sequences and are closely related to various human diseases, such as tumors, neurodegenerative diseases, and infections. This review focuses on the vital role of rG4s in ncRNAs, particularly lncRNAs and miRNAs. We outline the dynamic balance between rG4s and RNA stem-loop/hairpin structures and the interplay between ncRNAs and interactors, thereby modulating gene expression and disease progression. A complete understanding of the biological regulatory role and mechanism of rG4s in ncRNAs affirms the critical importance of folding into the appropriate three-dimensional structure in maintaining or modulating the functions of ncRNAs. It makes them novel therapeutic targets for adjusting potential-G4-containing-ncRNAs-associated diseases.

Insights into Enhancer RNAs: Biogenesis and Emerging Role in Brain Diseases

Enhancers are cis-acting elements that control the transcription of target genes and are transcribed into a class of noncoding RNAs (ncRNAs) termed enhancer RNAs (eRNAs). eRNAs have shorter half-lives than mRNAs and long noncoding RNAs; however, the frequency of transcription of eRNAs is close to that of mRNAs. eRNA expression is associated with a high level of histone mark H3K27ac and a low level of H3K27me3. Although eRNAs only account for a small proportion of ncRNAs, their functions are important. eRNAs can not only increase enhancer activity by promoting the formation of enhancer-promoter loops but also regulate transcriptional activation. Increasing numbers of studies have found that eRNAs play an important role in the occurrence and development of brain diseases; however, further research into eRNAs is required. This review discusses the concept, characteristics, classification, function, and potential roles of eRNAs in brain diseases.

An in Silico Approach to Identifying TF Binding Sites: Analysis of the Regulatory Regions of BUSCO Genes from Fungal Species in the Ceratocystidaceae Family

Transcriptional regulation controls gene expression through regulatory promoter regions that contain conserved sequence motifs. These motifs, also known as regulatory elements, are critically important to expression, which is driving research efforts to identify and characterize them. Yeasts have been the focus of such studies in fungi, including in several in silico approaches. This study aimed to determine whether in silico approaches could be used to identify motifs in the Ceratocystidaceae family, and if present, to evaluate whether these correspond to known transcription factors. This study targeted the 1000 base-pair region upstream of the start codon of 20 single-copy genes from the BUSCO dataset for motif discovery. Using the MEME and Tomtom analysis tools, conserved motifs at the family level were identified. The results show that such in silico approaches could identify known regulatory motifs in the Ceratocystidaceae and other unrelated species. This study provides support to ongoing efforts to use in silico analyses for motif discovery.

Behind the scenes: How RNA orchestrates the epigenetic regulation of gene expression

Non-coding DNA accounts for approximately 98.5% of the human genome. Once labeled as "junk DNA", this portion of the genome has undergone a progressive re-evaluation and it is now clear that some of its transcriptional products, belonging to the non-coding RNAs (ncRNAs), are key players in cell regulatory networks. A growing body of evidence demonstrates the crucial impact of regulatory ncRNAs on mammalian gene expression. Here, we focus on the defined relationship between chromatin-interacting RNAs, particularly long non-coding RNA (lncRNA), enhancer RNA (eRNA), non-coding natural antisense transcript (ncNAT), and circular RNA (circRNA) and epigenome, a common ground where both protein and RNA species converge to regulate cellular functions. Through several examples, this review provides an overview of the variety of targets, interactors, and mechanisms involved in the RNA-mediated modulation of loci-specific epigenetic states, a fundamental evolutive strategy to orchestrate mammalian gene expression in a timely and reversible manner. We will discuss how RNA-mediated epigenetic regulation impacts development and tissue homeostasis and how its alteration contributes to the onset and progression of many different human diseases, particularly cancer.

Low input capture Hi-C (liCHi-C) identifies promoter-enhancer interactions at high-resolution

Long-range interactions between regulatory elements and promoters are key in gene transcriptional control; however, their study requires large amounts of starting material, which is not compatible with clinical scenarios nor the study of rare cell populations. Here we introduce low input capture Hi-C (liCHi-C) as a cost-effective, flexible method to map and robustly compare promoter interactomes at high resolution. As proof of its broad applicability, we implement liCHi-C to study normal and malignant human hematopoietic hierarchy in clinical samples. We demonstrate that the dynamic promoter architecture identifies developmental trajectories and orchestrates transcriptional transitions during cell-state commitment. Moreover, liCHi-C enables the identification of disease-relevant cell types, genes and pathways potentially deregulated by non-coding alterations at distal regulatory elements. Finally, we show that liCHi-C can be harnessed to uncover genome-wide structural variants, resolve their breakpoints and infer their pathogenic effects. Collectively, our optimized liCHi-C method expands the study of 3D chromatin organization to unique, low-abundance cell populations, and offers an opportunity to uncover factors and regulatory networks involved in disease pathogenesis.

The epigenetic state of IL-4-polarized macrophages enables inflammatory cistromic expansion and extended synergistic response to TLR ligands

Prior exposure to microenvironmental signals could fundamentally change the response of macrophages to subsequent stimuli. It is believed that T helper-2 (Th2)-cell-type cytokine interleukin-4 (IL-4) and Toll-like receptor (TLR) ligand-activated transcriptional programs mutually antagonize each other, and no remarkable convergence has been identified between them. In contrast, here, we show that IL-4-polarized macrophages established a hyperinflammatory gene expression program upon lipopolysaccharide (LPS) exposure. This phenomenon, which we termed extended synergy, was supported by IL-4-directed epigenomic remodeling, LPS-activated NF-κB-p65 cistrome expansion, and increased enhancer activity. The EGR2 transcription factor contributed to the extended synergy in a macrophage-subtype-specific manner. Consequently, the previously alternatively polarized macrophages produced increased amounts of immune-modulatory factors both in vitro and in vivo in a murine Th2 cell-type airway inflammation model upon LPS exposure. Our findings establish that IL-4-induced epigenetic reprogramming is responsible for the development of inflammatory hyperresponsiveness to TLR activation and contributes to lung pathologies.

A predictive signature based on enhancer RNA associates with immune infiltration and aids treatment decision in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a prevalent urinary malignancy. Despite the recent development of better diagnostic tools and therapy, the five-year survival rate for individuals with advanced and metastatic ccRCC remains dismal. Unfortunately, ccRCC is less susceptible to radiation and chemotherapy. Consequently, targeted therapy and immunotherapy play a crucial role in the treatment of ccRCC. Enhancer RNAs (eRNAs) are noncoding RNAs transcribed by enhancers. Extensive research has shown that eRNAs are implicated in a variety of cancer signaling pathways. However, the biological functions of eRNAs have not been systematically investigated in ccRCC. In this study, we conducted a comprehensive investigation of the role of eRNAs in the onset and management of ccRCC. Patient prognosis-influencing eRNAs and target genes were chosen to construct a predictive signature. On the basis of the median riskscore, ccRCC patients were split into high- and low-risk subgroups. The prediction efficiency was assessed in several cohorts, and multi-omics analysis was carried out to investigate the differences and underlying mechanisms between the high- and low-risk groups. In addition, we investigated its potential to facilitate clinical treatment choices. The riskscore might be used to forecast a patient’s response to immunotherapy and targeted therapy, giving a revolutionary method for selecting treatment regimens with pinpoint accuracy.

Targeting phase separation on enhancers induced by transcription factor complex formations as a new strategy for treating drug-resistant cancers

The limited options for treating patients with drug-resistant cancers have emphasized the need to identify alternative treatment targets. Tumor cells have large super-enhancers (SEs) in the vicinity of important oncogenes for activation. The physical process of liquid-liquid phase separation (LLPS) contributes to the assembly of several membrane-less organelles in mammalian cells. Intrinsically disordered regions (IDRs) of proteins induce LLPS formation by developing condensates. It was discovered that key transcription factors (TFs) undergo LLPS in SEs. In addition, TFs play critical roles in the epigenetic and genetic regulation of cancer progression. Recently, we revealed the essential role of disease-specific TF collaboration changes in advanced prostate cancer (PC). OCT4 confers epigenetic changes by promoting complex formation with TFs, such as Forkhead box protein A1 (FOXA1), androgen receptor (AR) and Nuclear respiratory factor 1 (NRF1), inducing PC progression. It was demonstrated that TF collaboration through LLPS underlying transcriptional activation contributes to cancer aggressiveness and drug resistance. Moreover, the disruption of TF-mediated LLPS inhibited treatment-resistant PC tumor growth. Therefore, we propose that repression of TF collaborations involved in the LLPS of SEs could be a promising strategy for advanced cancer therapy. In this article, we summarize recent evidence highlighting the formation of LLPS on enhancers as a potent therapeutic target in advanced cancers.

Emerging Role of Enhancer RNAs as Potential Diagnostic and Prognostic Biomarkers in Cancer

Enhancers are distal cis-acting elements that are commonly recognized to regulate gene expression via cooperation with promoters. Along with regulating gene expression, enhancers can be transcribed and generate a class of non-coding RNAs called enhancer RNAs (eRNAs). The current discovery of abundant tissue-specific transcription of enhancers in various diseases such as cancers raises questions about the potential role of eRNAs in disease diagnosis and therapy. This review aimed to demonstrate the current understanding of eRNAs in cancer research with a focus on the potential roles of eRNAs as prognostic and diagnostic biomarkers in cancers.

RNA modifications in aging-associated cardiovascular diseases

Cardiovascular disease (CVD) is a leading cause of morbidity and mortality worldwide that bears an enormous healthcare burden and aging is a major contributing factor to CVDs. Functional gene expression network during aging is regulated by mRNAs transcriptionally and by non-coding RNAs epi-transcriptionally. RNA modifications alter the stability and function of both mRNAs and non-coding RNAs and are involved in differentiation, development, and diseases. Here we review major chemical RNA modifications on mRNAs and non-coding RNAs, including N6-adenosine methylation, N1-adenosine methylation, 5-methylcytidine, pseudouridylation, 2′ -O-ribose-methylation, and N7-methylguanosine, in the aging process with an emphasis on cardiovascular aging. We also summarize the currently available methods to detect RNA modifications and the bioinformatic tools to study RNA modifications. More importantly, we discussed the specific implication of the RNA modifications on mRNAs and non-coding RNAs in the pathogenesis of aging-associated CVDs, including atherosclerosis, hypertension, coronary heart diseases, congestive heart failure, atrial fibrillation, peripheral artery disease, venous insufficiency, and stroke.

Comprehensive transcriptomic analysis of immune-related eRNAs associated with prognosis and immune microenvironment in melanoma

Background

Recent evidence suggests that enhancer RNAs (eRNAs) play key roles in cancers. Identification of immune-related eRNAs (ireRNAs) in melanoma can provide novel insights into the mechanisms underlying its genesis and progression, along with potential therapeutic targets.

Aim

To establish an ireRNA-related prognostic signature for melanoma and identify potential drug candidates.

Methods

The ireRNAs associated with the overall survival (OS-ireRNAs) of melanoma patients were screened using data from The Cancer Genome Atlas (TCGA) via WGCNA and univariate Cox analysis. A prognostic signature based on these OS-ireRNAs was then constructed by performing the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. The immune landscape associated with the prognostic model was evaluated by the ESTIMATE algorithm and CIBERSORT method. Finally, the potential drug candidates for melanoma were screened through the cMap database.

Results

A total of 24 OS-ireRNAs were obtained, of which 7 ireRNAs were used to construct a prognostic signature. The ireRNAs-related signature performed well in predicting the overall survival (OS) of melanoma patients. The risk score of the established signature was further verified as an independent risk factor, and was associated with the unique tumor microenvironment in melanoma. We also identified several potential anti-cancer drugs for melanoma, of which corticosterone ranked first.

Conclusions

The ireRNA-related signature is an effective prognostic predictor and provides reliable information to better understand the mechanism of ireRNAs in the progression of melanoma.

Typical Enhancers, Super-Enhancers, and Cancers

Simple Summary

The cancer genome has been exhaustively studied upon the advent of Next-Generation Sequencing technologies. Coding and non-coding sequences have been defined as hotspots of genomic variations that affect the naïve gene expression programs established in normal cells, thus working as endogenous drivers of carcinogenesis. In this review, we comprehensively summarize fundamental aspects of gene expression regulation, with emphasis on the impact of sequence and structural variations mapped across non-coding cis-acting elements of genes encoding for tumor-related transcription factors. Chromatin architecture, epigenome reprogramming, transcriptional enhancers and Super-enhancers, oncogene regulation, cutting-edge technologies, and pharmacological treatment are substantially highlighted.

Abstract

Non-coding segments of the human genome are enriched in cis-regulatory modules that constitute functional elements, such as transcriptional enhancers and Super-enhancers. A hallmark of cancer pathogenesis is the dramatic dysregulation of the “archetype” gene expression profiles of normal human cells. Genomic variations can promote such deficiencies when occurring across enhancers and Super-enhancers, since they affect their mechanistic principles, their functional capacity and specificity, and the epigenomic features of the chromatin microenvironment across which these regulatory elements reside. Here, we comprehensively describe: fundamental mechanisms of gene expression dysregulation in cancers that involve genomic abnormalities within enhancers’ and Super-enhancers’ (SEs) sequences, which alter the expression of oncogenic transcription factors (TFs); cutting-edge technologies applied for the analysis of variation-enriched hotspots of the cancer genome; and pharmacological approaches for the treatment of Super-enhancers’ aberrant function. Finally, we provide an intratumor meta-analysis, which highlights that genomic variations in transcription-factor-driven tumors are accompanied overexpression of genes, a portion of which encodes for additional cancer-related transcription factors.

Identification of prognostic and immunotherapy-related eRNA ID2-AS1 in bladder cancer

AbstractEnhancer RNAs (eRNAs) can participate in enhancer regulation and target gene transcription, thus affecting the occurrence and development of tumors. In this study, we identified eRNAs closely related to bladder cancer (BLCA). Gene expression profiles and clinical information from The Cancer Genome Atlas (TCGA) database were used in this study. The Atlas of Noncoding RNAs in Cancer (TANRIC) co-expression data was also studied to evaluate correlations between the inferred levels of eRNA and its predicted target genes. Moreover, we evaluated differences in tumor microenvironment between high and low ID2-AS1 expression groups, and predicted the response of high- and low-expression groups to immune checkpoint inhibitor (ICI) treatment. Finally, we analyzed the prognostic value of ID2-AS1 in different tumors. ID2-AS1 and ID2 were identified as eRNAs and target genes related to the prognosis of BLCA. Low ID2-AS1 levels were associated with advanced age, low overall survival, high histological grade, and late BLCA staging. ID2-AS1 appeared to regulate epithelial mesenchymal transition, mitotic spindle assembly, and angiogenesis, thereby affecting BLCA progression. The ID2-AS1 high-expression group had better ICI treatment response. In addition, ID2-AS1 also had prognostic value in other cancers. ID2-AS1 helps predict prognostic and immunotherapeutic effects in BLCA.

Analysis of estrogen-regulated enhancer RNAs identifies a functional motif required for enhancer assembly and gene expression

To better understand the functions of non-coding enhancer RNAs (eRNAs), we annotated the estrogen-regulated eRNA transcriptome in estrogen receptor α (ERα)-positive breast cancer cells using PRO-cap and RNA sequencing. We then cloned a subset of the eRNAs identified, fused them to single guide RNAs, and targeted them to their ERα enhancers of origin using CRISPR/dCas9. Some of the eRNAs tested modulated the expression of cognate, but not heterologous, target genes after estrogen treatment by increasing ERα recruitment and stimulating p300-catalyzed H3K27 acetylation at the enhancer. We identified a ∼40 nucleotide functional eRNA regulatory motif (FERM) present in many eRNAs that was necessary and sufficient to modulate gene expression, but not the specificity of activation, after estrogen treatment. The FERM interacted with BCAS2, an RNA-binding protein amplified in breast cancers. The ectopic expression of a targeted eRNA controlling the expression of an oncogene resulted in increased cell proliferation, demonstrating the regulatory potential of eRNAs in breast cancer.

Multiple sclerosis genetic and non-genetic factors interact through the transient transcriptome

A clinically actionable understanding of multiple sclerosis (MS) etiology goes through GWAS interpretation, prompting research on new gene regulatory models. Our previous investigations suggested heterogeneity in etiology components and stochasticity in the interaction between genetic and non-genetic factors. To find a unifying model for this evidence, we focused on the recently mapped transient transcriptome (TT), that is mostly coded by intergenic and intronic regions, with half-life of minutes. Through a colocalization analysis, here we demonstrate that genomic regions coding for the TT are significantly enriched for MS-associated GWAS variants and DNA binding sites for molecular transducers mediating putative, non-genetic, determinants of MS (vitamin D deficiency, Epstein Barr virus latent infection, B cell dysfunction), indicating TT-coding regions as MS etiopathogenetic hotspots. Future research comparing cell-specific transient and stable transcriptomes may clarify the interplay between genetic variability and non-genetic factors causing MS. To this purpose, our colocalization analysis provides a freely available data resource at www.mscoloc.com.

Comprehensive Transcriptional Profiling and Mouse Phenotyping Reveals Dispensable Role for Adipose Tissue Selective Long Noncoding RNA Gm15551

Cold and nutrient-activated brown adipose tissue (BAT) is capable of increasing systemic energy expenditure via the uncoupled respiration and secretion of endocrine factors, thereby protecting mice against diet-induced obesity and improving insulin response and glucose tolerance in men. Long non-coding RNAs (lncRNAs) have recently been identified as fine-tuning regulators of cellular function. While certain lncRNAs have been functionally characterised in adipose tissue, their overall contribution in the activation of BAT remains elusive. We identified lncRNAs correlating to interscapular brown adipose tissue (iBAT) function in a high fat diet (HFD) and cold stressed mice. We focused on Gm15551, which has an adipose tissue specific expression profile, is highly upregulated during adipogenesis, and downregulated by β-adrenergic activation in mature adipocytes. Although we performed comprehensive transcriptional and adipocyte physiology profiling in vitro and in vivo, we could not detect an effect of gain or loss of function of Gm15551.

Identification of Enhancer RNA CDK6-AS1 as a Potential Novel Prognostic Biomarker in Gastric Cancer

Background: This study aimed to confirm the role of enhancer RNAs (eRNAs) in gastric cancer and their clinical utility.

Methods: We used Cox survival and relevance analysis to identify the candidate eRNAs in gastric cancer and performed Gene Ontology and Reactome pathway enrichment to determine the potential functions of eRNAs. Correlation between eRNA, tumor-infiltrating immune cells, and drug sensitivity was then analyzed.

Results:CDK6-AS1, a long non-coding RNA cyclin-dependent kinase 6, may serve as a poor potential prognostic biomarker candidate in gastric cancer with a positive correlation with its target gene CDK6. The low CDK6-AS1 expression group showed more frequent mutated driver genes than the high expression group. Moreover, CDK6-AS1 is involved in a key oncogenic pathway of the cell cycle and RNA transcription. CDK6-AS1 also shows dysregulations and associations with prognosis at the pan-cancer level. This eRNA may also be associated with immune cell infiltration and drug sensitivity.

Conclusion:CDK6-AS1 may be a potential prognostic biomarker and chemotherapeutic drug sensitivity predictor in gastric cancer.

Enhancer RNA: What we know and what we can achieve

Enhancers are important cis-acting elements that can regulate gene transcription and cell fate alongside promoters. In fact, many human cancers and diseases are associated with the malfunction of enhancers. Recent studies have shown that enhancers can produce enhancer RNAs (eRNAs) by RNA polymerase II. In this review, we discuss eRNA production, characteristics, functions and mechanics. eRNAs can determine chromatin accessibility, histone modification and gene expression by constructing a 'chromatin loop', thereby bringing enhancers to their target gene. eRNA can also be involved in the phase separation with enhancers and other proteins. eRNAs are abundant, and importantly, tissue-specific in tumours, various diseases and stem cells; thus, eRNAs can be a potential target for disease diagnosis and treatment. As eRNA is produced from the active transcription of enhancers and is involved in the regulation of cell fate, its manipulation will influence cell function, and therefore, it can be a new target for biological therapy.

The activity of human enhancers is modulated by the splicing of their associated lncRNAs

Pervasive enhancer transcription is at the origin of more than half of all long noncoding RNAs in humans. Transcription of enhancer-associated long noncoding RNAs (elncRNA) contribute to their cognate enhancer activity and gene expression regulation in cis. Recently, splicing of elncRNAs was shown to be associated with elevated enhancer activity. However, whether splicing of elncRNA transcripts is a mere consequence of accessibility at highly active enhancers or if elncRNA splicing directly impacts enhancer function, remains unanswered. We analysed genetically driven changes in elncRNA splicing, in humans, to address this outstanding question. We showed that splicing related motifs within multi-exonic elncRNAs evolved under selective constraints during human evolution, suggesting the processing of these transcripts is unlikely to have resulted from transcription across spurious splice sites. Using a genome-wide and unbiased approach, we used nucleotide variants as independent genetic factors to directly assess the causal relationship that underpin elncRNA splicing and their cognate enhancer activity. We found that the splicing of most elncRNAs is associated with changes in chromatin signatures at cognate enhancers and target mRNA expression. We provide evidence that efficient and conserved processing of enhancer-associated elncRNAs contributes to enhancer activity.

Most, if not all, active enhancers are transcribed, giving rise to a plethora of transcripts, including enhancer-associated long noncoding RNAs (elncRNAs). Changes in elncRNA levels impacts cognate enhancer activity. Recently splicing of elncRNA has also been found to associate with enhancer activity. Whether this associations reflects a contribution of elncRNA splicing to increased enhancer activity or else is simply the consequence of increased chromatin accessibility that promotes transcriptional elongation and allows for spurious splicing events remains unknown. We show that natural selection has acted, at the species and population level, to preserve DNA elements required for frequent and efficient elncRNA splicing Importantly, using a genome-wide and unbiased statistical population genomics approach, we demonstrate that elncRNA splicing is associated with cognate enhancer function, contributing to chromatin status and enhancer activity. Our results provides strong evidence that efficient elncRNA splicing contributes to enhancer activity genome-wide.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

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.

Kdm6a deficiency restricted to mouse hematopoietic cells causes an age- and sex-dependent myelodysplastic syndrome-like phenotype

Kdm6a/Utx, a gene on the X chromosome, encodes a histone H3K27me3 demethylase that has an orthologue on the Y chromosome (Uty) (Zheng et al. 2018). We previously identified inactivating mutations of Kdm6a in approximately 50% of mouse acute promyelocytic leukemia samples; however, somatic mutations of KDM6A are more rare in human AML samples, ranging in frequency from 2–15% in different series of patients, where their role in pathogenesis is not yet clear. In this study, we show that female Kdm6a^flox/flox mice (with allele inactivation initiated by Vav1-Cre in hematopoietic stem and progenitor cells (HSPCs) have a sex-specific phenotype that emerges with aging, with features resembling a myelodysplastic syndrome (MDS). Female Kdm6a-knockout (KO) mice have an age-dependent expansion of their HSPCs with aberrant self-renewal, but they did not differentiate normally into downstream progeny. These mice became mildly anemic and thrombocytopenic, but did not develop overt leukemia, or die from these cytopenias. ChIP-seq and ATAC-seq studies showed only minor changes in H3K27me3, H3K27ac, H3K4me, H3K4me3 and chromatin accessibility between Kdm6a-WT and Kdm6a-KO mice. Utilizing scRNA-seq, Kdm6a loss was linked to the transcriptional repression of genes that mediate hematopoietic cell fate determination. These data demonstrate that Kdm6a plays an important role in normal hematopoiesis, and that its inactivation may contribute to AML pathogenesis.

The expression of Drosophila melanogaster Hox gene Ultrabithorax is not overtly regulated by the intronic long non-coding RNA lncRNA: PS4 in a wild type genetic background

Long noncoding RNAs (lncRNAs) have been implicated in a variety of processes in development, differentiation, and disease. In Drosophila melanogaster, the bithorax Hox cluster contains three Hox genes [Ultrabithorax (Ubx), abdominal-A, and Abdominal-B], along with a number of lncRNAs, most with unknown functions. Here, we investigated the function of a lncRNA, lncRNA:PS4 that originates in the second intron of Ubx and is transcribed in the antisense orientation to Ubx. The expression pattern of lncRNA:PS4 is complementary to Ubx in the thoracic primordia, and the lncRNA:PS4 coding region overlaps the location of the large insertion that causes the dominant homeotic mutation Contrabithorax-1 (Ubx^Cbx-1), which partially transforms Drosophila wings into halteres via ectopic activation of Ubx. This led us to investigate the potential role of this lncRNA in regulation of Ubx expression. The Ubx^Cbx-1 mutation dramatically changes the pattern of lncRNA:PS4, eliminating the expression of most lncRNA:PS4 sequences from parasegment 4 (where Ubx protein is normally absent) and ectopically activating lncRNA:PS4 at high levels in the abdomen (where Ubx is normally expressed). These changes, however, did not lead to changes in the Ubx embryonic transcription pattern. Targeted deletion of the two promoters of lncRNA:PS4 did not result in the change of Ubx expression in the embryos. In the genetic background of a Ubx^Cbx-1 mutation, the lncRNA:PS4 mutation does slightly enhance the ectopic activation of Ubx protein expression in wing discs and also slightly enhances the wing phenotype seen in Ubx^Cbx-1 heterozygotes.

"Enhancing" mechanosensing: Enhancers and enhancer-derived long non-coding RNAs in endothelial response to flow

Endothelial cells (ECs), uniquely localized and strategically forming the inner lining of vascular wall, constitute the largest cell surface by area in the human body. The dynamic sensing and response of ECs to mechanical cues, especially shear stress, is crucial for maintenance of vascular homeostasis. It is well recognized that different flow patterns associated with atheroprotective vs atheroprone regions in the arterial tree, result in distinct EC functional phenotypes with differential transcriptome profiles. Mounting evidence has demonstrated an integrative and essential regulatory role of non-coding genome in EC biology. In particular, recent studies have begun to reveal the importance of enhancers and enhancer-derived transcripts in flow-regulated EC gene expression and function. In this minireview, we summarize studies in this area and discuss examples in support of the emerging importance of enhancers and enhancer(-derived) long non-coding RNAs (elncRNAs) in EC mechanosensing, with a focus on flow-responsive EC transcription. Finally, we will provide perspective and discuss standing questions to elucidate the role of these novel regulators in EC mechanobiology.

Noncoding sequence variants define a novel regulatory element in the first intron of the N-acetylglutamate synthase gene

N-acetylglutamate synthase deficiency is an autosomal recessive urea cycle disorder caused either by decreased expression of the NAGS gene or defective NAGS enzyme resulting in decreased production of N-acetylglutamate (NAG), an allosteric activator of carbamylphosphate synthetase 1 (CPS1). NAGSD is the only urea cycle disorder that can be effectively treated with a single drug, N-carbamylglutamate (NCG), a stable NAG analog, which activates CPS1 to restore ureagenesis. We describe three patients with NAGSD due to four novel noncoding sequence variants in the NAGS regulatory regions. All three patients had hyperammonemia that resolved upon treatment with NCG. Sequence variants NM_153006.2:c.427-222G>A and NM_153006.2:c.427-218A>C reside in the 547 bp-long first intron of NAGS and define a novel NAGS regulatory element that binds retinoic X receptor α. Sequence variants NC_000017.10:g.42078967A>T (NM_153006.2:c.-3065A>T) and NC_000017.10:g.42078934C>T (NM_153006.2:c.-3098C>T) reside in the NAGS enhancer, within known HNF1 and predicted glucocorticoid receptor binding sites, respectively. Reporter gene assays in HepG2 and HuH-7 cells demonstrated that all four substitutions could result in reduced expression of NAGS. These findings show that analyzing noncoding regions of NAGS and other urea cycle genes can reveal molecular causes of disease and identify novel regulators of ureagenesis.

A Novel Regulatory Player in the Innate Immune System: Long Non-Coding RNAs

Long non-coding RNAs (lncRNAs) represent crucial transcriptional and post-transcriptional gene regulators during antimicrobial responses in the host innate immune system. Studies have shown that lncRNAs are expressed in a highly tissue- and cell-specific- manner and are involved in the differentiation and function of innate immune cells, as well as inflammatory and antiviral processes, through versatile molecular mechanisms. These lncRNAs function via the interactions with DNA, RNA, or protein in either cis or trans pattern, relying on their specific sequences or their transcriptions and processing. The dysregulation of lncRNA function is associated with various human non-infectious diseases, such as inflammatory bowel disease, cardiovascular diseases, and diabetes mellitus. Here, we provide an overview of the regulation and mechanisms of lncRNA function in the development and differentiation of innate immune cells, and during the activation or repression of innate immune responses. These elucidations might be beneficial for the development of therapeutic strategies targeting inflammatory and innate immune-mediated diseases.

Analyses of Promoter , Enhancer, and Nucleosome Organization in Mammalian Cells by MNase-Seq

MNase-Seq is a genome-wide procedure that allows mapping of DNA associated to nucleosomes following micrococcal nuclease digestion. It is a rapid and robust technology useful for the analysis of chromatin properties genome-wide at the resolution of mono-nucleosomes. Here, we describe how to produce high-resolution nucleosome maps of cells grown in suspension or adherent mammalian cells. After only three steps: nuclei or cell preparation, native MNase digestion and DNA purification, libraries for high-throughput sequencing can be prepared. Genome-wide nucleosome maps allow analyzing chromatin opening at promoters or enhancers, nucleosome displacement, or labile nucleosome occupancy depending on the digestion condition used. As presented, MNase-Seq is a versatile tool for investigating chromatin dynamics, regulation, and to define open chromatin regions of regulatory elements in mammalian genomes.

Gene Model Related to m6A Predicts the Prognostic Effect of Immune Infiltration on Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a highly aggressive solid tumor. Because most studies have focused on the intrinsic carcinogenic pathways of tumors, we focused on the relationship between N6-methyladenosine (m6A) and the prognosis of HNSCC in the tumor immune microenvironment. We downloaded RNA-seq data from the TCGA dataset and used univariate Cox regression to screen m6A-related lncRNAs. The expression value of LASSO-screened genes was the sum of LASSO regression coefficients. We then evaluated relationships between the risk score and cellular components or cellular immune response. Differences in immune response under various algorithms were visualized with heat maps. The GSVA package in R was used to analyze GO, BP, KEGG, and hallmark gene sets of immune checkpoint clusters and immune checkpoint scores. The GSEA analysis was performed with the cluster profile package, yielding 21 m6A genes. Related lncRNAs were screened with Pearson's correlations, and the resulting 442 lncRNAs were screened using single-factor analysis. Eight lncRNAs closely related to prognosis were identified through survival random forest. Survival analysis showed that patients with a high risk score had a poor prognosis. Low- and high-risk-score groups differed significantly in m6A gene expression. Prognostic scores from different algorithms were significantly correlated with B cells, T cells, and memory cells in the immune microenvironment. Expression of immune checkpoints and signal pathways differed significantly across risk-score groups, suggesting that m6A could mediate lncRNA-induced immune system dysfunction and affect HNSCC development. A comprehensive study of tumor-cell immune characteristics should provide more insight into the complex immune microenvironment, thus contributing to the development of new immunotherapeutic agents.

Enhancer RNA m6A methylation facilitates transcriptional condensate formation and gene activation

The mechanistic understanding of nascent RNAs in transcriptional control remains limited. Here, by a high sensitivity method methylation-inscribed nascent transcripts sequencing (MINT-seq), we characterized the landscapes of N6-methyladenosine (m6A) on nascent RNAs. We uncover heavy but selective m6A deposition on nascent RNAs produced by transcription regulatory elements, including promoter upstream antisense RNAs and enhancer RNAs (eRNAs), which positively correlates with their length, inclusion of m6A motif, and RNA abundances. m6A-eRNAs mark highly active enhancers, where they recruit nuclear m6A reader YTHDC1 to phase separate into liquid-like condensates, in a manner dependent on its C terminus intrinsically disordered region and arginine residues. The m6A-eRNA/YTHDC1 condensate co-mixes with and facilitates the formation of BRD4 coactivator condensate. Consequently, YTHDC1 depletion diminished BRD4 condensate and its recruitment to enhancers, resulting in inhibited enhancer and gene activation. We propose that chemical modifications of eRNAs together with reader proteins play broad roles in enhancer activation and gene transcriptional control.

Enhancing Gonadotrope Gene Expression Through Regulatory lncRNAs

The world of long non-coding RNAs (lncRNAs) has opened up massive new prospects in understanding the regulation of gene expression. Not only are there seemingly almost infinite numbers of lncRNAs in the mammalian cell, but they have highly diverse mechanisms of action. In the nucleus, some are chromatin-associated, transcribed from transcriptional enhancers (eRNAs) and/or direct changes in the epigenetic landscape with profound effects on gene expression. The pituitary gonadotrope is responsible for activation of reproduction through production and secretion of appropriate levels of the gonadotropic hormones. As such, it exemplifies a cell whose function is defined through changes in developmental and temporal patterns of gene expression, including those that are hormonally induced. Roles for diverse distal regulatory elements and eRNAs in gonadotrope biology have only just begun to emerge. Here, we will present an overview of the different kinds of lncRNAs that alter gene expression, and what is known about their roles in regulating some of the key gonadotrope genes. We will also review various screens that have detected differentially expressed pituitary lncRNAs associated with changes in reproductive state and those whose expression is found to play a role in gonadotrope-derived nonfunctioning pituitary adenomas. We hope to shed light on this exciting new field, emphasize the open questions, and encourage research to illuminate the roles of lncRNAs in various endocrine systems.

TBX5-AS1, an enhancer RNA, is a potential novel prognostic biomarker for lung adenocarcinoma

BACKGROUND

Enhancer RNAs (eRNAs) are demonstrated to be closely associated with tumourigenesis and cancer progression. However, the role of eRNAs in lung adenocarcinoma (LUAD) remains largely unclear. Thus, a comprehensive analysis was constructed to identify the key eRNAs, and to explore the clinical utility of the identified eRNAs in LUAD.

METHODS

First, LUAD expression profile data from the Cancer Genome Atlas (TCGA) dataset and eRNA-relevant information were integrated for Kaplan-Meier survival analysis and Spearman's correlation analysis to filtered the key candidate eRNAs that was associated with survival rate and their target genes in LUAD. Then, the key eRNA was selected for subsequent clinical correlation analysis. KEGG pathway enrichment analyses were undertaken to explore the potential signaling pathways of the key eRNA. Data from the human protein atlas (HPA) database were used to validate the outcomes and the quantitative real time-polymerase chain reaction (qRT-PCR) analysis was conducted to measure eRNA expression levels in tumor tissues and paired normal adjacent tissues from LUAD patients. Finally, the eRNAs were validated in pan-cancer.

RESULTS

As a result, TBX5-AS1 was identified as the key eRNA, which has T-box transcription factor 5 (TBX5) as its regulatory target. KEGG analysis indicated that TBX5-AS1 may exert a vital role via the PI3K/AKT pathway, Ras signaling pathway, etc. Additionally, the qRT-PCR results and the HPA database indicated that TBX5-AS1 and TBX5 were significantly downregulated in tumour samples compared to matched-adjacent pairs. The pan-cancer validation results showed that TBX5-AS1 was associated with survival in four tumors, namely, adrenocortical carcinoma (ACC), LUAD, lung squamous cell carcinoma (LUSC), and uterine corpus endometrial carcinoma (UCEC). Correlations were found between TBX5-AS1 and its target gene, TBX5, in 26 tumor types.

CONCLUSION

Collectively, our results indicated that TBX5-AS1 may be a potential prognostic biomarker for LUAD patients and promote the targeted therapy of LUAD.

The Biological Function Delineated Across Pan-Cancer Levels Through lncRNA-Based Prognostic Risk Assessment Factors for Pancreatic Cancer

Long non-coding RNAs (lncRNAs) play key roles in tumors and function not only as important molecular markers for cancer prognosis, but also as molecular characteristics at the pan-cancer level. Because of the poor prognosis of pancreatic cancer, accurate assessment of prognosis is a key issue in the development of treatment plans for pancreatic cancer. Here we analyzed pancreatic cancer data from The Cancer Genome Atlas and The Genotype Tissue Expression database using Cox regression and lasso regression in analyses using a combination of the two databases as well as only The Cancer Genome Atlas database (Cancer Genome Atlas Research Network et al., 2013). A prognostic risk score model with significant correlation with pancreatic cancer survival was constructed, and two lncRNAs were investigated. Additional analysis of 33 cancers using the two lncRNAs showed that lncRNA TsPOAP1-AS1 was a prognostic marker of seven cancers, among which pancreatic cancer was the most significant, and lncRNA mi600hg was a prognostic marker of ovarian cancer and pancreatic cancer. LncRNA TsPOAP1-AS1 is associated with clinical stage and tumor mutation burden of some cancers as well as a strong degree of immune infiltration in many cancers, while a strong correlation between lncRNA mi600hg and microsatellite instability was observed in several cancers. The results of this study help further our understanding of the different functions of lncRNAs in cancer and may aid in the clinical application of lncRNAs as prognostic factors for cancer.

Chromatin Accessibility Identifies Regulatory Elements Predictive of Gene Expression and Disease Outcome in Multiple Myeloma

PURPOSE

Multiple myeloma is a malignancy of plasma cells. Extensive genetic and transcriptional characterization of myeloma has identified subtypes with prognostic and therapeutic implications. In contrast, relatively little is known about the myeloma epigenome.

EXPERIMENTAL DESIGN

CD138+CD38+ myeloma cells were isolated from fresh bone marrow aspirate or the same aspirate after freezing for 1-6 months. Gene expression and chromatin accessibility were compared between fresh and frozen samples by RNA sequencing (RNA-seq) and assay for transpose accessible chromatin sequencing (ATAC-seq). Chromatin accessible regions were used to identify regulatory RNA expression in more than 700 samples from newly diagnosed patients in the Multiple Myeloma Research Foundation CoMMpass trial (NCT01454297).

RESULTS

Gene expression and chromatin accessibility of cryopreserved myeloma recapitulated that of freshly isolated samples. ATAC-seq performed on a series of biobanked specimens identified thousands of chromatin accessible regions with hundreds being highly coordinated with gene expression. More than 4,700 of these chromatin accessible regions were transcribed in newly diagnosed myelomas from the CoMMpass trial. Regulatory element activity alone recapitulated myeloma gene expression subtypes, and in particular myeloma subtypes with immunoglobulin heavy chain translocations were defined by transcription of distal regulatory elements. Moreover, enhancer activity predicted oncogene expression implicating gene regulatory mechanisms in aggressive myeloma.

CONCLUSIONS

These data demonstrate the feasibility of using biobanked specimens for retrospective studies of the myeloma epigenome and illustrate the unique enhancer landscapes of myeloma subtypes that are coupled to gene expression and disease progression.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

Pausing of transcribing RNA polymerase is regulated and creates opportunities to control gene expression. Research in metazoans has so far mainly focused on RNA polymerase II (Pol II) promoter-proximal pausing leaving the pervasive nature of pausing and its regulatory potential in mammalian cells unclear. Here, we developed a pause detecting algorithm (PDA) for nucleotide-resolution occupancy data and a new native elongating transcript sequencing approach, termed nested NET-seq, that strongly reduces artifactual peaks commonly misinterpreted as pausing sites. Leveraging PDA and nested NET-seq reveal widespread genome-wide Pol II pausing at single-nucleotide resolution in human cells. Notably, the majority of Pol II pauses occur outside of promoter-proximal gene regions primarily along the gene-body of transcribed genes. Sequence analysis combined with machine learning modeling reveals DNA sequence properties underlying widespread transcriptional pausing including a new pause motif. Interestingly, key sequence determinants of RNA polymerase pausing are conserved between human cells and bacteria. These studies indicate pervasive sequence-induced transcriptional pausing in human cells and the knowledge of exact pause locations implies potential functional roles in gene expression.