Enhancer transcripts mark active estrogen receptor binding sites

TF Profiler: a transcription factor inference method that broadly measures transcription factor activity and identifies mechanistically distinct networks

TF Profiler is a method of inferring transcription factor (TF) regulatory activity, i.e., when a TF is present and actively participating in the regulation of transcription, directly from nascent sequencing assays such as PRO-seq and GRO-seq. While ChIP assays have measured DNA localization, they fall short of identifying when and where the effector domain of a transcription factor is active. Our method uses RNA polymerase activity to infer TF effector domain activity across hundreds of data sets and transcription factors. TF Profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

A transcription coupling model for how enhancers communicate with their target genes

How enhancers communicate with their target genes to influence transcription is an unresolved question of fundamental importance. Current models of the mechanism of enhancer-target gene or enhancer-promoter (E-P) communication are transcription-factor-centric and underappreciate major findings, including that enhancers are themselves transcribed by RNA polymerase II, which correlates with enhancer activity. In this Perspective, we posit that enhancer transcription and its products, enhancer RNAs, are elementary components of enhancer-gene communication. Specifically, we discuss the possibility that transcription at enhancers and at their cognate genes are linked and that this coupling is at the basis of how enhancers communicate with their targets. This model of transcriptional coupling between enhancers and their target genes is supported by growing experimental evidence and represents a synthesis of recent key discoveries.

The epipliancy journey: Tumor initiation at the mercy of identity crisis and epigenetic drift

Cellular pliancy refers to the unique disposition of different stages of cellular differentiation to transform when exposed to specific oncogenic insults. This concept highlights a strong interconnection between cellular identity and tumorigenesis, and implies overcoming of epigenetic barriers defining cellular states. Emerging evidence suggests that the cell-type-specific response to intrinsic and extrinsic stresses is modulated by accessibility to certain areas of the genome. Understanding the interplay between epigenetic mechanisms, cellular differentiation, and oncogenic insults is crucial for deciphering the complex nature of tumorigenesis and developing targeted therapies. Hence, cellular pliancy relies on a dynamic cooperation between the cellular identity and the cellular context through epigenetic control, including the reactivation of cellular mechanisms, such as epithelial-to-mesenchymal transition (EMT). Such mechanisms and pathways confer plasticity to the cell allowing it to adapt to a hostile environment in a context of tumor initiation, thus changing its cellular identity. Indeed, growing evidence suggests that cancer is a disease of cell identity crisis, whereby differentiated cells lose their defined identity and gain progenitor characteristics. The loss of cell fate commitment is a central feature of tumorigenesis and appears to be a prerequisite for neoplastic transformation. In this context, EMT-inducing transcription factors (EMT-TFs) cooperate with mitogenic oncoproteins to foster malignant transformation. The aberrant activation of EMT-TFs plays an active role in tumor initiation by alleviating key oncosuppressive mechanisms and by endowing cancer cells with stem cell-like properties, including the ability to self-renew, thus changing the course of tumorigenesis. This highly dynamic phenotypic change occurs concomitantly to major epigenome reorganization, a key component of cell differentiation and cancer cell plasticity regulation. The concept of pliancy was initially proposed to address a fundamental question in cancer biology: why are some cells more likely to become cancerous in response to specific oncogenic events at particular developmental stages? We propose the concept of epipliancy, whereby a difference in epigenetic configuration leads to malignant transformation following an oncogenic insult. Here, we present recent studies furthering our understanding of how the epigenetic landscape may impact the modulation of cellular pliancy during early stages of cancer initiation.

A Super Enhancer-Derived Enhancer RNA Acts Together with CTCF/Cohesin in Trans to Regulate Erythropoiesis

Background/Objectives: Enhancer RNAs (eRNAs) function in diverse modes and increasing studies have shown that they play important roles in normal development and disease. However, their role in erythropoiesis is not fully understood. Methods: We analyzed published RNA-seq and Promoter Capture Hi-C data from mouse E14.5 fetal liver cells to identify enhancer RNAs in erythroid cells with long-range interactions. Results: We discovered an erythroid-specific enhancer RNA (CpoxeRNA) transcribed from an enhancer region upstream of Cpox, an enzyme important for heme synthesis. CpoxeRNA is important for erythropoiesis, as the knockdown of CpoxeRNA by shRNA results in impaired enucleation and cell proliferation during terminal differentiation. CpoxeRNA interacts with cohesin and acts both in cis and trans to regulate erythroid genes. Conclusions: we have identified a trans-acting eRNA, CpoxeRNA, as a potential regulator of terminal erythropoiesis.

An Integrated Database for Exploring Alternative Promoters in Animals

Alternative promoter (AP) events, as a major pre-transcriptional mechanism, can initiate different transcription start sites to generate distinct mRNA isoforms and regulate their expression. At present, hundreds of thousands of APs have been identified across human tissues, and a considerable number of APs have been demonstrated to be associated with complex traits and diseases. Recent researches have also proven important effects of APs on animals. However, the landscape of APs in animals has not been fully recognized. In this study, 102,349 AP profiles from 23,077 samples across 12 species were systematically characterized. We further identified tissue-specific APs and investigated trait-related promoters among various species. In addition, we analyzed the associations between APs and enhancer RNAs (eRNA)/transcription factors (TF) as a means of identifying potential regulatory factors. Integrating these findings, we finally developed Animal-APdb, a database for the searching, browsing, and downloading of information related to Animal APs. Animal-APdb is expected to serve as a valuable resource for exploring the functions and mechanisms of APs in animals.

A role for NFIB in SOX2 downregulation and epigenome accessibility changes due to long-term estrogen treatment of breast cancer epithelial cells

Estrogen (E2) regulates the differentiation and proliferation of mammary progenitor cells by modulating the transcription of multiple genes. One of the genes that is downregulated by E2 is SOX2, a transcription factor associated with stem and progenitor cells that is overexpressed during breast tumourigenesis. To elucidate the mechanisms underlying E2-mediated SOX2 repression, we investigated epigenome and transcriptome changes following short- and long-term E2 exposure in breast cancer cells. We found that short-term E2 exposure reduces chromatin accessibility at the downstream SOX2 SRR134 enhancer, decreasing SOX2 expression. In contrast, long-term E2 exposure completely represses SOX2 transcription while maintaining accessibility at the SRR124-134 enhancer cluster, keeping it poised for reactivation. This repression was accompanied by widespread epigenome and transcriptome changes associated with commitment towards a more differentiated and less invasive luminal phenotype. Finally, we identified a role for the transcription factor NFIB in this process, suggesting it collaborates with the estrogen receptor to mediate SOX2 repression and genome-wide epigenome accessibility changes.

An RNA-centric view of transcription and genome organization

Foundational models of transcriptional regulation involve the assembly of protein complexes at DNA elements associated with specific genes. These assemblies, which can include transcription factors, cofactors, RNA polymerase, and various chromatin regulators, form dynamic spatial compartments that contribute to both gene regulation and local genome architecture. This DNA-protein-centric view has been modified with recent evidence that RNA molecules have important roles to play in gene regulation and genome structure. Here, we discuss evidence that gene regulation by RNA occurs at multiple levels that include assembly of transcriptional complexes and genome compartments, feedback regulation of active genes, silencing of genes, and control of protein kinases. We thus provide an RNA-centric view of transcriptional regulation that must reside alongside the more traditional DNA-protein-centric perspectives on gene regulation and genome architecture.

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explored the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild-type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln modulates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation of ERα on serine 118 (S118), as well as by reducing the E2-dependent binding of ERα across the genome. These effects were associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of some genes and enhanced expression of others. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

Cut from the same cloth: RNAs transcribed from regulatory elements

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

Relaxin Modulates the Genomic Actions and Biological Effects of Estrogen in the Myometrium

Estradiol (E2) and relaxin (Rln) are steroid and polypeptide hormones, respectively, with important roles in the female reproductive tract, including myometrium. Some actions of Rln, which are mediated by its membrane receptor RXFP1, require or are augmented by E2 signaling through its cognate nuclear steroid receptor, estrogen receptor alpha (ERα). In contrast, other actions of Rln act in opposition to the effects of E2. Here we explored the molecular and genomic mechanisms that underlie the functional interplay between E2 and Rln in the myometrium. We used both ovariectomized female mice and immortalized human myometrial cells expressing wild-type or mutant ERα (hTERT-HM-ERα cells). Our results indicate that Rln modulates the genomic actions and biological effects of estrogen in the myometrium and myometrial cells by reducing phosphorylation of ERα on serine 118 (S118), as well as by reducing the E2-dependent binding of ERα across the genome. These effects were associated with changes in the hormone-regulated transcriptome, including a decrease in the E2-dependent expression of some genes and enhanced expression of others. The inhibitory effects of Rln cotreatment on the E2-dependent phosphorylation of ERα required the nuclear dual-specificity phosphatases DUSP1 and DUSP5. Moreover, the inhibitory effects of Rln were reflected in a concomitant inhibition of the E2-dependent contraction of myometrial cells. Collectively, our results identify a pathway that integrates Rln/RXFP1 and E2/ERα signaling, resulting in a convergence of membrane and nuclear signaling pathways to control genomic and biological outcomes.

Mechanistic analysis of enhancer sequences in the estrogen receptor transcriptional program

Estrogen Receptor α (ERα) is a major lineage determining transcription factor (TF) in mammary gland development. Dysregulation of ERα-mediated transcriptional program results in cancer. Transcriptomic and epigenomic profiling of breast cancer cell lines has revealed large numbers of enhancers involved in this regulatory program, but how these enhancers encode function in their sequence remains poorly understood. A subset of ERα-bound enhancers are transcribed into short bidirectional RNA (enhancer RNA or eRNA), and this property is believed to be a reliable marker of active enhancers. We therefore analyze thousands of ERα-bound enhancers and build quantitative, mechanism-aware models to discriminate eRNAs from non-transcribing enhancers based on their sequence. Our thermodynamics-based models provide insights into the roles of specific TFs in ERα-mediated transcriptional program, many of which are supported by the literature. We use in silico perturbations to predict TF-enhancer regulatory relationships and integrate these findings with experimentally determined enhancer-promoter interactions to construct a gene regulatory network. We also demonstrate that the model can prioritize breast cancer-related sequence variants while providing mechanistic explanations for their function. Finally, we experimentally validate the model-proposed mechanisms underlying three such variants.

Transcriptomic Profiles of Long Noncoding RNAs and Their Target Protein-Coding Genes Reveals Speciation Adaptation on the Qinghai-Xizang (Tibet) Plateau in Orinus

Long noncoding RNAs (lncRNAs) are RNA molecules longer than 200 nt, which lack the ability to encode proteins and are involved in multifarious growth, development, and regulatory processes in plants and mammals. However, the environmental-regulated expression profiles of lncRNAs in Orinus that may associated with their adaptation on the Qinghai-Xizang (Tibet) Plateau (QTP) have never been characterized. Here, we utilized transcriptomic sequencing data of two Orinus species (O. thoroldii and O. kokonoricus) to identify 1624 lncRNAs, including 1119 intergenic lncRNAs, 200 antisense lncRNAs, five intronic lncRNAs, and 300 sense lncRNAs. In addition, the evolutionary relationships of Orinus lncRNAs showed limited sequence conservation among 39 species, which implied that Orinus-specific lncRNAs contribute to speciation adaptation evolution. Furthermore, considering the cis-regulation mechanism, from 286 differentially expressed lncRNAs (DElncRNAs) and their nearby protein coding genes (PCGs) between O. thoroldii and O. kokonoricus, 128 lncRNA-PCG pairs were obtained in O. thoroldii, whereas 92 lncRNA-PCG pairs were obtained in O. kokonoricus. In addition, a total of 19 lncRNA-PCG pairs in O. thoroldii and 14 lncRNA-PCG pairs in O. kokonoricus were found to participate in different biological processes, indicating that the different expression profiles of DElncRNAs between O. thoroldii and O. kokonoricus were associated with their adaptation at different elevations on the QTP. We also found several pairs of DElncRNA nearby transcription factors (TFs), indicating that these DElncRNAs regulate the expression of TFs to aid O. thoroldii in adapting to the environment. Therefore, this work systematically identified a series of lncRNAs in Orinus, laying the groundwork for further exploration into the biological function of Orinus in environmental adaptation.

Lineage-determining transcription factor-driven promoters regulate cell type-specific macrophage gene expression

Mammalian promoters consist of multifarious elements, which make them unique and support the selection of the proper transcript variants required under diverse conditions in distinct cell types. However, their direct DNA-transcription factor (TF) interactions are mostly unidentified. Murine bone marrow-derived macrophages (BMDMs) are a widely used model for studying gene expression regulation. Thus, this model serves as a rich source of various next-generation sequencing data sets, including a large number of TF cistromes. By processing and integrating the available cistromic, epigenomic and transcriptomic data from BMDMs, we characterized the macrophage-specific direct DNA-TF interactions, with a particular emphasis on those specific for promoters. Whilst active promoters are enriched for certain types of typically methylatable elements, more than half of them contain non-methylatable and prototypically promoter-distal elements. In addition, circa 14% of promoters-including that of Csf1r-are composed exclusively of 'distal' elements that provide cell type-specific gene regulation by specialized TFs. Similar to CG-rich promoters, these also contain methylatable CG sites that are demethylated in a significant portion and show high polymerase activity. We conclude that this unusual class of promoters regulates cell type-specific gene expression in macrophages, and such a mechanism might exist in other cell types too.

Identification of key eRNAs for intervertebral disc degeneration by integrated multinomial bioinformatics analysis

BACKGROUND

Intervertebral disc degeneration (IVDD) is a common degenerative condition leading to abnormal stress distribution under load, causing intervertebral stenosis, facet joint degeneration, and foraminal stenosis. Very little is known about the molecular mechanism of eRNAs in IVDD.

METHODS

Gene expression profiles of 38 annulus disc samples composed of 27 less degenerated discs (LDs) and 11 more degenerated discs (MDs) were retrieved from the GEO database. Then, differentially expressed enhancer RNAs (DEeRNAs), differentially expressed target genes (DETGs), and differentially expressed transcription factors (DETFs), hallmark of cancer signalling pathways according to GSVA; the types and quantity of immune cells according to CIBERSORT; and immune gene sets according to ssGSEA were analysed to construct an IVDD-related eRNA network. Then, multidimensional validation was performed to explore the interactions among DEeRNAs, DETFs and DEGs in space.

RESULTS

A total of 53 components, 14 DETGs, 15 DEeRNAs, 3 DETFs, 5 immune cells, 9 hallmarks, and 7 immune gene sets, were selected to construct the regulatory network. After validation by online multidimensional databases, 21 interactive DEeRNA-DEG-DETF axes related to IVDD exacerbation were identified, among which the C1S-CTNNB1-CHD4 axis was the most significant.

CONCLUSION

Based upon the results of our study, we theorize that the C1S-CTNNB1-CHD4 axis plays a vital role in IVDD exacerbation. Specifically, C1S recruits CTNNB1 and upregulates the expression of CHD4 in IVDD, and subsequently, CHD4 suppresses glycolysis and activates oxidative phosphorylation, thus generating insoluble collagen fibre deposits and leading to the progression of IVDD. Overall, these DEeRNAs could comprise promising therapeutic targets for IVDD due to their high tissue specificity.

Transcription decouples estrogen-dependent changes in enhancer-promoter contact frequencies and spatial proximity

How enhancers regulate their target genes in the context of 3D chromatin organization is extensively studied and models which do not require direct enhancer-promoter contact have recently emerged. Here, we use the activation of estrogen receptor-dependent enhancers in a breast cancer cell line to study enhancer-promoter communication at two loci. This allows high temporal resolution tracking of molecular events from hormone stimulation to efficient gene activation. We examine how both enhancer-promoter spatial proximity assayed by DNA fluorescence in situ hybridization, and contact frequencies resulting from chromatin in situ fragmentation and proximity ligation, change dynamically during enhancer-driven gene activation. These orthogonal methods produce seemingly paradoxical results: upon enhancer activation enhancer-promoter contact frequencies increase while spatial proximity decreases. We explore this apparent discrepancy using different estrogen receptor ligands and transcription inhibitors. Our data demonstrate that enhancer-promoter contact frequencies are transcription independent whereas altered enhancer-promoter proximity depends on transcription. Our results emphasize that the relationship between contact frequencies and physical distance in the nucleus, especially over short genomic distances, is not always a simple one.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

A transcription factor (TF) inference method that broadly measures TF activity and identifies mechanistically distinct TF networks

TF profiler is a method of inferring transcription factor regulatory activity, i.e. when a TF is present and actively regulating transcription, directly directly from nascent sequencing assays such as PRO-seq and GRO-seq. Transcription factors orchestrate transcription and play a critical role in cellular maintenance, identity and response to external stimuli. While ChIP assays have measured DNA localization, they fall short of identifying when and where transcription factors are actively regulating transcription. Our method, on the other hand, uses RNA polymerase activity to infer TF activity across hundreds of data sets and transcription factors. Based on these classifications we identify three distinct classes of transcription factors: ubiquitous factors that play roles in cellular homeostasis, driving basal gene programs across tissues and cell types, tissue specific factors that act almost exclusively at enhancers and are themselves regulated at transcription, and stimulus responsive TFs which are regulated post-transcriptionally but act predominantly at enhancers. TF profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Differential Expression of Noncoding RNAs Revealed Enhancer RNA AC016735.2 as a Potential Pathogenic Marker of Congenital Microtia Patients

PURPOSE

Congenital microtia is a complex maxillofacial malformation with various risk factors. This study aimed to find potential pathogenic noncoding RNAs for congenital microtia patients.

METHODS

We collected 3 pairs of residual ear cartilage samples and corresponding normal ear cartilage samples from nonsyndromic congenital microtia patients for microarray experiments. The differentially expressed RNAs were screened, and enrichment analysis and correlation expression analysis were performed to elucidate the function of the differentially expressed genes (DEGs). We further investigated the most significantly differentially expressed long noncoding RNA (lncRNA), AC016735.2, through follow-up analyses including RT-qPCR and Western blotting, to validate its differential expression in residual ear cartilage compared with normal ear cartilage. SiRNA was designed to study the regulatory role of AC016735.2, and cell proliferation experiments were conducted to explore its impact on residual ear chondrocytes.

RESULTS

Analysis of the microarray data revealed a total of 1079 differentially expressed RNAs, including 305 mRNAs and x lncRNAs, using a threshold of FC>1.5 and P<0.05 for mRNA, and FC>1.0 and P<0.05 for lncRNA. Enrichment analysis indicated that these genes are mainly involved in extracellular matrix regulation and embryonic development. AC016735.2 showed the highest differential expression among the eRNAs, being upregulated in residual ear cartilage. It acts in cis to regulate the nearby coding gene ZFP36L2, indirectly affecting downstream genes such as BMP4, TWSG1, COL2A1, and COL9A2.

CONCLUSION

Significant differences were observed in the expression of lncRNAs and mRNAs between residual ear cartilage and normal auricular cartilage tissues in the same genetic background of congenital microtia. These differentially expressed lncRNAs and mRNAs may play crucial roles in the occurrence and development of microtia through pathways associated with extracellular matrix regulation and gastrulation. Particularly, AC016735.2, an eRNA acting in cis, could serve as a potential pathogenic noncoding gene.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Topological regulation of the estrogen transcriptional response by ZATT-mediated inhibition of TOP2B activity

Human type-II topoisomerases, TOP2A and TOP2B, remove transcription associated DNA supercoiling, thereby affecting gene-expression programs, and have recently been associated with 3D genome architecture. Here, we study the regulatory roles of TOP2 paralogs in response to estrogen, which triggers an acute transcriptional induction that involves rewiring of genome organization. We find that, whereas TOP2A facilitates transcription, as expected for a topoisomerase, TOP2B limits the estrogen response. Consistent with this, TOP2B activity is locally downregulated upon estrogen treatment to favor the establishment and stabilization of regulatory chromatin contacts, likely through an accumulation of DNA supercoiling. We show that estrogen-mediated inhibition of TOP2B requires estrogen receptor α (ERα), a non-catalytic function of TOP2A, and the action of the atypical SUMO-ligase ZATT. This mechanism of topological transcriptional-control, which may be shared by additional gene-expression circuits, highlights the relevance of DNA topoisomerases as central actors of genome dynamics.

Genome-wide identification of transcriptional enhancers during human placental development and association with function, differentiation, and disease†

The placenta is a dynamic organ that must perform a remarkable variety of functions during its relatively short existence in order to support a developing fetus. These functions include nutrient delivery, gas exchange, waste removal, hormone production, and immune barrier protection. Proper placenta development and function are critical for healthy pregnancy outcomes, but the underlying genomic regulatory events that control this process remain largely unknown. We hypothesized that mapping sites of transcriptional enhancer activity and associated changes in gene expression across gestation in human placenta tissue would identify genomic loci and predicted transcription factor activity related to critical placenta functions. We used a suite of genomic assays [i.e., RNA-sequencing (RNA-seq), Precision run-on-sequencing (PRO-seq), and Chromatin immunoprecipitation-sequencing (ChIP-seq)] and computational pipelines to identify a set of >20 000 enhancers that are active at various time points in gestation. Changes in the activity of these enhancers correlate with changes in gene expression. In addition, some of these enhancers encode risk for adverse pregnancy outcomes. We further show that integrating enhancer activity, transcription factor motif analysis, and transcription factor expression can identify distinct sets of transcription factors predicted to be more active either in early pregnancy or at term. Knockdown of selected identified transcription factors in a trophoblast stem cell culture model altered the expression of key placental marker genes. These observations provide a framework for future mechanistic studies of individual enhancer-transcription factor-target gene interactions and have the potential to inform genetic risk prediction for adverse pregnancy outcomes.

Atlas of nascent RNA transcripts reveals enhancer to gene linkages

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 chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Transcriptional responses of cancer cells to heat shock-inducing stimuli involve amplification of robust HSF1 binding

Responses of cells to stimuli are increasingly discovered to involve the binding of sequence-specific transcription factors outside of known target genes. We wanted to determine to what extent the genome-wide binding and function of a transcription factor are shaped by the cell type versus the stimulus. To do so, we induced the Heat Shock Response pathway in two different cancer cell lines with two different stimuli and related the binding of its master regulator HSF1 to nascent RNA and chromatin accessibility. Here, we show that HSF1 binding patterns retain their identity between basal conditions and under different magnitudes of activation, so that common HSF1 binding is globally associated with distinct transcription outcomes. HSF1-induced increase in DNA accessibility was modest in scale, but occurred predominantly at remote genomic sites. Apart from regulating transcription at existing elements including promoters and enhancers, HSF1 binding amplified during responses to stimuli may engage inactive chromatin.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

eRNA co-expression network uncovers TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of "bivalent" TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

Genome-wide measurement of RNA dissociation from chromatin classifies transcripts by their dynamics and reveals rapid dissociation of enhancer lncRNAs

Long non-coding RNAs (lncRNAs) are involved in gene expression regulation in cis. Although enriched in the cell chromatin fraction, to what degree this defines their regulatory potential remains unclear. Furthermore, the factors underlying lncRNA chromatin tethering, as well as the molecular basis of efficient lncRNA chromatin dissociation and its impact on enhancer activity and target gene expression, remain to be resolved. Here, we developed chrTT-seq, which combines the pulse-chase metabolic labeling of nascent RNA with chromatin fractionation and transient transcriptome sequencing to follow nascent RNA transcripts from their transcription on chromatin to release and allows the quantification of dissociation dynamics. By incorporating genomic, transcriptomic, and epigenetic metrics, as well as RNA-binding protein propensities, in machine learning models, we identify features that define transcript groups of different chromatin dissociation dynamics. Notably, lncRNAs transcribed from enhancers display reduced chromatin retention, suggesting that, in addition to splicing, their chromatin dissociation may shape enhancer activity.

Characterizing intergenic transcription at RNA polymerase II binding sites in normal and cancer tissues

Intergenic transcription in normal and cancerous tissues is pervasive but incompletely understood. To investigate this, we constructed an atlas of over 180,000 consensus RNA polymerase II (RNAPII)-bound intergenic regions from 900 RNAPII chromatin immunoprecipitation sequencing (ChIP-seq) experiments in normal and cancer samples. Through unsupervised analysis, we identified 51 RNAPII consensus clusters, many of which mapped to specific biotypes and revealed tissue-specific regulatory signatures. We developed a meta-clustering methodology to integrate our RNAPII atlas with active transcription across 28,797 RNA sequencing (RNA-seq) samples from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Encyclopedia of DNA Elements (ENCODE). This analysis revealed strong tissue- and disease-specific interconnections between RNAPII occupancy and transcriptional activity. We demonstrate that intergenic transcription at RNAPII-bound regions is a novel per-cancer and pan-cancer biomarker. This biomarker displays genomic and clinically relevant characteristics, distinguishing cancer subtypes and linking to overall survival. Our results demonstrate the effectiveness of coherent data integration to uncover intergenic transcriptional activity in normal and cancer tissues.

Hormone-induced enhancer assembly requires an optimal level of hormone receptor multivalent interactions

Transcription factors (TFs) activate enhancers to drive cell-specific gene programs in response to signals, but our understanding of enhancer assembly during signaling events is incomplete. Here, we show that androgen receptor (AR) forms condensates through multivalent interactions mediated by its N-terminal intrinsically disordered region (IDR) to orchestrate enhancer assembly in response to androgen signaling. AR IDR can be substituted by IDRs from selective proteins for AR condensation capacity and its function on enhancers. Expansion of the poly(Q) track within AR IDR results in a higher AR condensation propensity as measured by multiple methods, including live-cell single-molecule microscopy. Either weakening or strengthening AR condensation propensity impairs its heterotypic multivalent interactions with other enhancer components and diminishes its transcriptional activity. Our work reveals the requirement of an optimal level of AR condensation in mediating enhancer assembly and suggests that alteration of the fine-tuned multivalent IDR-IDR interactions might underlie AR-related human pathologies.

Serine starvation silences estrogen receptor signaling through histone hypoacetylation

Loss of estrogen receptor (ER) pathway activity promotes breast cancer progression, yet how this occurs remains poorly understood. Here, we show that serine starvation, a metabolic stress often found in breast cancer, represses estrogen receptor alpha (ERα) signaling by reprogramming glucose metabolism and epigenetics. Using isotope tracing and time-resolved metabolomic analyses, we demonstrate that serine is required to maintain glucose flux through glycolysis and the TCA cycle to support acetyl-CoA generation for histone acetylation. Consequently, limiting serine depletes histone H3 lysine 27 acetylation (H3K27ac), particularly at the promoter region of ER pathway genes including the gene encoding ERα, ESR1. Mechanistically, serine starvation impairs acetyl-CoA-dependent gene expression by inhibiting the entry of glycolytic carbon into the TCA cycle and down-regulating the mitochondrial citrate exporter SLC25A1, a critical enzyme in the production of nucleocytosolic acetyl-CoA from glucose. Consistent with this model, total H3K27ac and ERα expression are suppressed by SLC25A1 inhibition and restored by acetate, an alternate source of acetyl-CoA, in serine-free conditions. We thus uncover an unexpected role for serine in sustaining ER signaling through the regulation of acetyl-CoA metabolism.

Perturbed epigenetic transcriptional regulation in AML with IDH mutations causes increased susceptibility to NK cells

Isocitrate dehydrogenase (IDH) mutations are found in 20% of acute myeloid leukemia (AML) patients. However, only 30-40% of the patients respond to IDH inhibitors (IDHi). We aimed to identify a molecular vulnerability to tailor novel therapies for AML patients with IDH mutations. We characterized the transcriptional and epigenetic landscape with the IDH2i AG-221, using an IDH2 mutated AML cell line model and AML patient cohorts, and discovered a perturbed transcriptional regulatory network involving myeloid transcription factors that were partly restored after AG-221 treatment. In addition, hypermethylation of the HLA cluster caused a down-regulation of HLA class I genes, triggering an enhanced natural killer (NK) cell activation and an increased susceptibility to NK cell-mediated responses. Finally, analyses of DNA methylation data from IDHi-treated patients showed that non-responders still harbored hypermethylation in HLA class I genes. In conclusion, this study provides new insights suggesting that IDH mutated AML is particularly sensitive to NK cell-based personalized immunotherapy.

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.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Interactive enhancer hubs (iHUBs) mediate transcriptional reprogramming and adaptive resistance in pancreatic cancer

OBJECTIVE

Pancreatic ductal adenocarcinoma (PDAC) displays a remarkable propensity towards therapy resistance. However, molecular epigenetic and transcriptional mechanisms enabling this are poorly understood. In this study, we aimed to identify novel mechanistic approaches to overcome or prevent resistance in PDAC.

DESIGN

We used in vitro and in vivo models of resistant PDAC and integrated epigenomic, transcriptomic, nascent RNA and chromatin topology data. We identified a JunD-driven subgroup of enhancers, called interactive hubs (iHUBs), which mediate transcriptional reprogramming and chemoresistance in PDAC.

RESULTS

iHUBs display characteristics typical for active enhancers (H3K27ac enrichment) in both therapy sensitive and resistant states but exhibit increased interactions and production of enhancer RNA (eRNA) in the resistant state. Notably, deletion of individual iHUBs was sufficient to decrease transcription of target genes and sensitise resistant cells to chemotherapy. Overlapping motif analysis and transcriptional profiling identified the activator protein 1 (AP1) transcription factor JunD as a master transcription factor of these enhancers. JunD depletion decreased iHUB interaction frequency and transcription of target genes. Moreover, targeting either eRNA production or signaling pathways upstream of iHUB activation using clinically tested small molecule inhibitors decreased eRNA production and interaction frequency, and restored chemotherapy responsiveness in vitro and in vivo. Representative iHUB target genes were found to be more expressed in patients with poor response to chemotherapy compared with responsive patients.

CONCLUSION

Our findings identify an important role for a subgroup of highly connected enhancers (iHUBs) in regulating chemotherapy response and demonstrate targetability in sensitisation to chemotherapy.

macroH2A2 antagonizes epigenetic programs of stemness in glioblastoma

Self-renewal is a crucial property of glioblastoma cells that is enabled by the choreographed functions of chromatin regulators and transcription factors. Identifying targetable epigenetic mechanisms of self-renewal could therefore represent an important step toward developing effective treatments for this universally lethal cancer. Here we uncover an epigenetic axis of self-renewal mediated by the histone variant macroH2A2. With omics and functional assays deploying patient-derived in vitro and in vivo models, we show that macroH2A2 shapes chromatin accessibility at enhancer elements to antagonize transcriptional programs of self-renewal. macroH2A2 also sensitizes cells to small molecule-mediated cell death via activation of a viral mimicry response. Consistent with these results, our analyses of clinical cohorts indicate that high transcriptional levels of this histone variant are associated with better prognosis of high-grade glioma patients. Our results reveal a targetable epigenetic mechanism of self-renewal controlled by macroH2A2 and suggest additional treatment approaches for glioblastoma patients.

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.

Estrogen receptors as potential therapeutic target in endometrial cancer

Endometrial cancer (EC) is one of the most common gynecological carcinomas in both developed and developing countries. Majority of the gynecological malignancies are hormonally driven where estrogen signaling acts as an oncogenic signal. Estrogen's effects are mediated via classical nuclear estrogen receptors; estrogen receptor alpha and beta (ERα and ERβ) and a trans-membrane G protein-coupled estrogen receptor (GPR30 and GPER). ERs and GPER through ligand binding triggers multiple downstream signaling pathways causing cell cycle regulation, cell differentiation, migration, and apoptosis in various tissues including endometrium. Although the molecular aspect of estrogen function in ER-mediated signaling is now partly understood, the same is not true for GPER-mediated signaling in endometrial malignancies. Understanding the physiological roles of ERα and GPER in EC biology therefore leads to the identification of some novel therapeutic targets. Here we review the effect of estrogen signaling through ERα-and GPER in EC, major types, and some affordable treatment approaches for endometrial tumor patients which has interesting implications in understanding uterine cancer progression.

Emerging regulatory mechanisms of noncoding RNAs in topologically associating domains

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

Revealing eRNA interactions: TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced from transcriptional enhancers that are highly correlated with their activities. Using capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across individuals, we identified inter-individual variation of expression in over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers interact with each other and with promoters. Mid-to-long range interactions showed distance-dependent decay, which was modified by TF occupancy. In particular, we found a class of 'bivalent' TFs, including Cohesin, which both facilitates and insulates the interaction between enhancers and/or promoters depending on the topology. In short ranges, we observed strand specific interactions between nearby eRNAs in both convergent or divergent orientations. Our finding supports a cooperative convergent eRNA model, which is compatible with eRNA remodeling neighboring enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions depending on the distance, orientation, and the binding landscapes of TFs.

Enhancer-promoter entanglement explains their transcriptional interdependence

Enhancers not only activate target promoters to stimulate messenger RNA (mRNA) synthesis, but they themselves also undergo transcription to produce enhancer RNAs (eRNAs), the significance of which is not well understood. Transcription at the participating enhancer-promoter pair appears coordinated, but it is unclear why and how. Here, we employ cell-free transcription assays using constructs derived from the human GREB1 locus to demonstrate that transcription at an enhancer and its target promoter is interdependent. This interdependence is observable under conditions where direct enhancer-promoter contact (EPC) takes place. We demonstrate that transcription activation at a participating enhancer-promoter pair is dependent on i) the mutual availability of the enhancer and promoter, ii) the state of transcription at both the enhancer and promoter, iii) local abundance of both eRNA and mRNA, and iv) direct EPC. Our results suggest transcriptional interdependence between the enhancer and the promoter as the basis of their transcriptional concurrence and coordination throughout the genome. We propose a model where transcriptional concurrence, coordination and interdependence are possible if the participating enhancer and promoter are entangled in the form of EPC, reside in a proteinaceous bubble, and utilize shared transcriptional resources and regulatory inputs.

Estrogen-Inducible LncRNA BNAT1 Functions as a Modulator for Estrogen Receptor Signaling in Endocrine-Resistant Breast Cancer Cells

Recent advances in RNA studies have revealed that functional long noncoding RNAs (lncRNAs) contribute to the biology of cancers. In breast cancer, estrogen receptor α (ERα) is an essential transcription factor that primarily promotes the growth of luminal-type cancer, although only a small number of lncRNAs are identified as direct ERα targets and modulators for ERα signaling. In this study, we performed RNA-sequencing for ER-positive breast cancer cells and identified a novel estrogen-inducible antisense RNA in the COL18A1 promoter region, named breast cancer natural antisense transcript 1 (BNAT1). In clinicopathological study, BNAT1 may have clinical relevance as a potential diagnostic factor for prognoses of ER-positive breast cancer patients based on an in situ hybridization study for breast cancer specimens. siRNA-mediated BNAT1 silencing significantly inhibited the in vitro and in vivo growth of tamoxifen-resistant ER-positive breast cancer cells. Notably, BNAT1 silencing repressed cell cycle progression whereas it promoted apoptosis. Microarray analysis revealed that BNAT1 silencing in estrogen-sensitive breast cancer cells repressed estrogen signaling. We showed that BNAT1 knockdown decreased ERα expression and repressed ERα transactivation. RNA immunoprecipitation showed that BNAT1 physically binds to ERα protein. In summary, BNAT1 would play a critical role in the biology of ER-positive breast cancer by modulating ERα-dependent transcription regulation. We consider that BNAT1 could be a potential molecular target for diagnostic and therapeutic options targeting luminal-type and endocrine-resistant breast cancer.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Recognition RNA Structure and Dynamics &gt; RNA Structure, Dynamics and Chemistry RNA Methods &gt; RNA Analyses in Cells.

Xrn2 substrate mapping identifies torpedo loading sites and extensive premature termination of RNA pol II transcription

The exonuclease torpedo Xrn2 loads onto nascent RNA 5'-PO₄ ends and chases down pol II to promote termination downstream from polyA sites. We report that Xrn2 is recruited to preinitiation complexes and "travels" to 3' ends of genes. Mapping of 5'-PO₄ ends in nascent RNA identified Xrn2 loading sites stabilized by an active site mutant, Xrn2(D235A). Xrn2 loading sites are approximately two to 20 bases downstream from where CPSF73 cleaves at polyA sites and histone 3' ends. We propose that processing of all mRNA 3' ends comprises cleavage and limited 5'-3' trimming by CPSF73, followed by handoff to Xrn2. A similar handoff occurs at tRNA 3' ends, where cotranscriptional RNase Z cleavage generates novel Xrn2 substrates. Exonuclease-dead Xrn2 increased transcription in 3' flanking regions by inhibiting polyA site-dependent termination. Surprisingly, the mutant Xrn2 also rescued transcription in promoter-proximal regions to the same extent as in 3' flanking regions. eNET-seq revealed Xrn2-mediated degradation of sense and antisense nascent RNA within a few bases of the TSS, where 5'-PO₄ ends may be generated by decapping or endonucleolytic cleavage. These results suggest that a major fraction of pol II complexes terminates prematurely close to the start site under normal conditions by an Xrn2-mediated torpedo mechanism.

Subtype and cell type specific expression of lncRNAs provide insight into breast cancer

Long non-coding RNAs (lncRNAs) are involved in breast cancer pathogenesis through chromatin remodeling, transcriptional and post-transcriptional gene regulation. We report robust associations between lncRNA expression and breast cancer clinicopathological features in two population-based cohorts: SCAN-B and TCGA. Using co-expression analysis of lncRNAs with protein coding genes, we discovered three distinct clusters of lncRNAs. In silico cell type deconvolution coupled with single-cell RNA-seq analyses revealed that these three clusters were driven by cell type specific expression of lncRNAs. In one cluster lncRNAs were expressed by cancer cells and were mostly associated with the estrogen signaling pathways. In the two other clusters, lncRNAs were expressed either by immune cells or fibroblasts of the tumor microenvironment. To further investigate the cis-regulatory regions driving lncRNA expression in breast cancer, we identified subtype-specific transcription factor (TF) occupancy at lncRNA promoters. We also integrated lncRNA expression with DNA methylation data to identify long-range regulatory regions for lncRNA which were validated using ChiA-Pet-Pol2 loops. lncRNAs play an important role in shaping the gene regulatory landscape in breast cancer. We provide a detailed subtype and cell type-specific expression of lncRNA, which improves the understanding of underlying transcriptional regulation in breast cancer.

RNA Pol II pausing facilitates phased pluripotency transitions by buffering transcription

Promoter-proximal RNA Pol II pausing is a critical step in transcriptional control. Pol II pausing has been predominantly studied in tissue culture systems. While Pol II pausing has been shown to be required for mammalian development, the phenotypic and mechanistic details of this requirement are unknown. Here, we found that loss of Pol II pausing stalls pluripotent state transitions within the epiblast of the early mouse embryo. Using Nelfb -/- mice and a NELFB degron mouse pluripotent stem cell model, we show that embryonic stem cells (ESCs) representing the naïve state of pluripotency successfully initiate a transition program but fail to balance levels of induced and repressed genes and enhancers in the absence of NELF. We found an increase in chromatin-associated NELF during transition from the naïve to later pluripotent states. Overall, our work defines the acute and long-term molecular consequences of NELF loss and reveals a role for Pol II pausing in the pluripotency continuum as a modulator of cell state transitions.

Novel estrogen-responsive genes (ERGs) for the evaluation of estrogenic activity

Estrogen action is mediated by various genes, including estrogen-responsive genes (ERGs). ERGs have been used as reporter-genes and markers for gene expression. Gene expression profiling using a set of ERGs has been used to examine statistically reliable transcriptomic assays such as DNA microarray assays and RNA sequencing (RNA-seq). However, the quality of ERGs has not been extensively examined. Here, we obtained a set of 300 ERGs that were newly identified by six sets of RNA-seq data from estrogen-treated and control human breast cancer MCF-7 cells. The ERGs exhibited statistical stability, which was based on the coefficient of variation (CV) analysis, correlation analysis, and examination of the functional association with estrogen action using database searches. A set of the top 30 genes based on CV ranking were further evaluated quantitatively by RT-PCR and qualitatively by a functional analysis using the GO and KEGG databases and by a mechanistic analysis to classify ERα/β-dependent or ER-independent types of transcriptional regulation. The 30 ERGs were characterized according to (1) the enzymes, such as metabolic enzymes, proteases, and protein kinases, (2) the genes with specific cell functions, such as cell-signaling mediators, tumor-suppressors, and the roles in breast cancer, (3) the association with transcriptional regulation, and (4) estrogen-responsiveness. Therefore, the ERGs identified here represent various cell functions and cell signaling pathways, including estrogen signaling, and thus, may be useful to evaluate estrogenic activity.

Enhancer RNAs step forward: new insights into enhancer function

Enhancers confer precise spatiotemporal patterns of gene expression in response to developmental and environmental stimuli. Over the last decade, the transcription of enhancer RNAs (eRNAs) – nascent RNAs transcribed from active enhancers – has emerged as a key factor regulating enhancer activity. eRNAs are relatively short-lived RNA species that are transcribed at very high rates but also quickly degraded. Nevertheless, eRNAs are deeply intertwined within enhancer regulatory networks and are implicated in a number of transcriptional control mechanisms. Enhancers show changes in function and sequence over evolutionary time, raising questions about the relationship between enhancer sequences and eRNA function. Moreover, the vast majority of single nucleotide polymorphisms associated with human complex diseases map to the non-coding genome, with causal disease variants enriched within enhancers. In this Primer, we survey the diverse roles played by eRNAs in enhancer-dependent gene expression, evaluating different models for eRNA function. We also explore questions surrounding the genetic conservation of enhancers and how this relates to eRNA function and dysfunction.

Summary: This Primer evaluates the ideas that underpin developing models for eRNA function, exploring cases in which perturbed eRNA function contributes to disease.