Identifying transcription start sites and active enhancer elements using BruUV-seq

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Cdc73 protects Notch-induced T-cell leukemia cells from DNA damage and mitochondrial stress

ABSTRACT

Activated Notch signaling is highly prevalent in T-cell acute lymphoblastic leukemia (T-ALL), but pan-Notch inhibitors showed excessive toxicity in clinical trials. To find alternative ways to target Notch signals, we investigated cell division cycle 73 (Cdc73), which is a Notch cofactor and key component of the RNA polymerase-associated transcriptional machinery, an emerging target in T-ALL. Although we confirmed previous work that CDC73 interacts with NOTCH1, we also found that the interaction in T-ALL was context-dependent and facilitated by the transcription factor ETS1. Using mouse models, we showed that Cdc73 is important for Notch-induced T-cell development and T-ALL maintenance. Mechanistically, chromatin and nascent gene expression profiling showed that Cdc73 intersects with Ets1 and Notch at chromatin within enhancers to activate expression of known T-ALL oncogenes through its enhancer functions. Cdc73 also intersects with these factors within promoters to activate transcription of genes that are important for DNA repair and oxidative phosphorylation through its gene body functions. Consistently, Cdc73 deletion induced DNA damage and apoptosis and impaired mitochondrial function. The CDC73-induced DNA repair expression program co-opted by NOTCH1 is more highly expressed in T-ALL than in any other cancer. These data suggest that Cdc73 might induce a gene expression program that was eventually intersected and hijacked by oncogenic Notch to augment proliferation and mitigate the genotoxic and metabolic stresses of elevated Notch signaling. Our report supports studying factors such as CDC73 that intersect with Notch to derive a basic scientific understanding on how to combat Notch-dependent cancers without directly targeting the Notch complex.

Enhanced frequency of transcription pre-initiation complexes assembly after exposure to UV irradiation results in increased repair activity and reduced probabilities for mutagenesis

In addition to being essential for gene expression, transcription is crucial for the maintenance of genome integrity. Here, we undertook a systematic approach, to monitor the assembly kinetics of the pre-initiating RNA Polymerase (Pol) II at promoters at steady state and different stages during recovery from UV irradiation-stress, when pre-initiation and initiation steps have been suggested to be transiently shut down. Taking advantage of the reversible dissociation of pre-initiating Pol II after high salt treatment, we found that de novo recruitment of the available Pol II molecules at active promoters not only persists upon UV at all times tested but occurs significantly faster in the early phase of recovery (2 h) than in unexposed human fibroblasts at the majority of active genes. Our method unveiled groups of genes with significantly different pre-initiation complex (PIC) assembly dynamics after UV that present distinct rates of UV-related mutational signatures in melanoma tumours, providing functional relevance to the importance of keeping transcription initiation active during UV recovery. Our findings uncover novel mechanistic insights further detailing the multilayered transcriptional response to genotoxic stress and link PIC assembly dynamics after exposure to genotoxins with cancer mutational landscapes.

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.

Cdc73 protects Notch-induced T-cell leukemia cells from DNA damage and mitochondrial stress

Activated Notch signaling is highly prevalent in T-cell acute lymphoblastic leukemia (T-ALL) but pan-Notch inhibitors were toxic in clinical trials. To find alternative ways to target Notch signals, we investigated Cell division cycle 73 (Cdc73), which is a Notch cofactor and component of transcriptional machinery, a potential target in T-ALL. While we confirmed previous work that CDC73 interacts with NOTCH1, we also found that the interaction in T-ALL was context-dependent and facilitated by the lymphoid transcription factor ETS1. Using mouse models, we showed that Cdc73 is important for Notch-induced T-cell development and T-ALL maintenance. Mechanistically, Cdc73, Ets1, and Notch intersect chromatin at promoters and enhancers to activate oncogenes and genes that are important for DNA repair and oxidative phosphorylation. Consistently, Cdc73 deletion in T-ALL cells induced DNA damage and impaired mitochondrial function. Our data suggests that Cdc73 might promote a gene expression program that was eventually intersected by Notch to mitigate the genotoxic and metabolic stresses of elevated Notch signaling. We also provide mechanistic support for testing inhibitors of DNA repair, oxidative phosphorylation, and transcriptional machinery. Inhibiting pathways like Cdc73 that intersect with Notch at chromatin might constitute a strategy to weaken Notch signals without directly targeting the Notch complex.

Interplay Between the Histone Variant H2A.Z and the Epigenome in Pancreatic Cancer

BACKGROUND

The oncogenic process is orchestrated by a complex network of chromatin remodeling elements that shape the cancer epigenome. Histone variant H2A.Z regulates DNA control elements such as promoters and enhancers in different types of cancer; however, the interplay between H2A.Z and the pancreatic cancer epigenome is unknown.

OBJECTIVE

This study analyzed the role of H2A.Z in different DNA regulatory elements.

METHODS

We performed Chromatin Immunoprecipitation Sequencing assays (ChiP-seq) with total H2A.Z and acetylated H2A.Z (acH2A.Z) antibodies and analyzed published data from ChIP-seq, RNA-seq, bromouridine labeling-UV and sequencing (BruUV-seq), Hi-C and ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) in the pancreatic cancer cell line PANC-1.

RESULTS

The results indicate that total H2A.Z facilitates the recruitment of RNA polymerase II and transcription factors at promoters and enhancers allowing the expression of pro-oncogenic genes. Interestingly, we demonstrated that H2A.Z is enriched in super-enhancers (SEs) contributing to the transcriptional activation of key genes implicated in tumor development. Importantly, we established that H2A.Z contributes to the three-dimensional (3D) genome organization of pancreatic cancer and that it is a component of the Topological Associated Domains (TADs) boundaries in PANC-1 and that total H2A.Z and acH2A.Z are associated with A and B compartments, respectively.

CONCLUSIONS

H2A.Z participates in the biology and development of pancreatic cancer by generating a pro-oncogenic transcriptome through its posttranslational modifications, interactions with different partners, and regulatory elements, contributing to the oncogenic 3D genome organization. These data allow us to understand the molecular mechanisms that promote an oncogenic transcriptome in pancreatic cancer mediated by H2A.Z.

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 > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

Transcription and genome integrity

Transcription can cause genome instability by promoting R-loop formation but also act as a mutation-suppressing machinery by sensing of DNA lesions leading to the activation of DNA damage signaling and transcription-coupled repair. Recovery of RNA synthesis following the resolution of repair of transcription-blocking lesions is critical to avoid apoptosis and several new factors involved in this process have recently been identified. Some DNA repair proteins are recruited to initiating RNA polymerases and this may expediate the recruitment of other factors that participate in the repair of transcription-blocking DNA lesions. Recent studies have shown that transcription of protein-coding genes does not always give rise to spliced transcripts, opening the possibility that cells may use the transcription machinery in a splicing-uncoupled manner for other purposes including surveillance of the transcribed genome.

A comparison of experimental assays and analytical methods for genome-wide identification of active enhancers

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. Here, we compared 13 genome-wide RNA sequencing assays in K562 cells and showed that the nuclear run-on followed by cap-selection assay (GRO/PRO-cap) has advantages in eRNA detection and active enhancer identification. We also introduced a tool, Peak Identifier for Nascent Transcript Starts (PINTS), to identify active promoters and enhancers genome-wide and pinpoint the precise location of the 5′ transcription start sites. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected eRNA TSSs available in 120 cell and tissue types that can be accessed at https://pints.yulab.org. With the knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the way for selection and characterization of their functions in a time- and labor-efficient manner in the future.

Integration of high-resolution promoter profiling assays reveals novel, cell type-specific transcription start sites across 115 human cell and tissue types

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type–specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3′ ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

Implications of Enhancer Transcription and eRNAs in Cancer

Despite extensive progress in developing anticancer therapies, therapy resistance remains a major challenge that promotes disease relapse. The changes that lead to therapy resistance can be intrinsically present or may be initiated during treatment. Genetic and epigenetic heterogeneity in tumors make it more challenging to deal with therapy resistance. Recent advances in genome-wide analyses have revealed that the deregulation of distal gene regulatory elements, such as enhancers, appears in several pathophysiological conditions, including cancer. Beyond the conventional function of enhancers in recruiting transcription factors to gene promoters, enhancer elements are also transcribed into noncoding RNAs known as enhancer RNAs (eRNA). Accumulating evidence suggests that uncontrolled enhancer activity with aberrant eRNA expression promotes oncogenesis. Interestingly, tissue-specific, transcribed eRNAs from active enhancers can serve as potential therapeutic targets or biomarkers in several cancer types. This review provides a comprehensive overview of the mechanisms of enhancer transcription and eRNAs as well as their potential roles in cancer and drug resistance.

Oncogenic seRNA functional activation: a novel mechanism of tumorigenesis

seRNA is a noncoding RNA (ncRNA) transcribed from active super-enhancer (SE), through which SE exerts biological functions and participates in various physiological and pathological processes. seRNA recruits cofactor, RNA polymerase II and mediator to constitute and stabilize chromatin loop SE and promoter region, which regulates target genes transcription. In tumorigenesis, DNA insertion, deletion, translocation, focal amplification and carcinogen factor mediate oncogenic SE generation, meanwhile, oncogenic SE transcribes into tumor-related seRNA, termed as oncogenic seRNA. Oncogenic seRNA participates in tumorigenesis through activating various signal-pathways. The recent reports showed that oncogenic seRNA implicates in a widespread range of cytopathological processes in cancer progression including cell proliferation, apoptosis, autophagy, epithelial-mesenchymal transition, extracellular matrix stiffness and angiogenesis. In this article, we comprehensively summarized seRNA’s characteristics and functions, and emphatically introduced inducible formation of oncogenic seRNA and its functional mechanisms. Lastly, some research strategies on oncogenic seRNA were introduced, and the perspectives on cancer therapy that targets oncogenic seRNA were also discussed.

Characterization of novel primary miRNA transcription units in human cells using Bru-seq nascent RNA sequencing

MicroRNAs (miRNAs) are key contributors to gene regulatory networks. Because miRNAs are processed from RNA polymerase II transcripts, insight into miRNA regulation requires a comprehensive understanding of the regulation of primary miRNA transcripts. We used Bru-seq nascent RNA sequencing and hidden Markov model segmentation to map primary miRNA transcription units (TUs) across 32 human cell lines, allowing us to describe TUs encompassing 1443 miRNAs from miRBase and 438 from MirGeneDB. We identified TUs for 61 miRNAs with an unknown CAGE TSS signal for MirGeneDB miRNAs. Many primary transcripts containing miRNA sequences failed to generate mature miRNAs, suggesting that miRNA biosynthesis is under both transcriptional and post-transcriptional control. In addition to constitutive and cell-type specific TU expression regulated by differential promoter usage, miRNA synthesis can be regulated by transcription past polyadenylation sites (transcriptional read through) and promoter divergent transcription (PROMPTs). We identified 197 miRNA TUs with novel promoters, 97 with transcriptional read-throughs and 3 miRNA TUs that resemble PROMPTs in at least one cell line. The miRNA TU annotation data resource described here reveals a greater complexity in miRNA regulation than previously known and provides a framework for identifying cell-type specific differences in miRNA transcription in cancer and cell transition states.

Continuous transcription initiation guarantees robust repair of all transcribed genes and regulatory regions

Inhibition of transcription caused by DNA damage-impaired RNA polymerase II (Pol II) elongation conceals a local increase in de novo transcription, slowly progressing from Transcription Start Sites (TSSs) to gene ends. Although associated with accelerated repair of Pol II-encountered lesions and limited mutagenesis, it is still unclear how this mechanism is maintained during genotoxic stress-recovery. Here we uncover a widespread gain in chromatin accessibility and preservation of the active H3K27ac mark after UV-irradiation. The concomitant increase in Pol II escape from promoter-proximal pause (PPP) sites of most active genes, PROMPTs and enhancer RNAs favors unrestrained initiation, as evidenced by the synthesis of nascent RNAs including start RNAs. Accordingly, drug-inhibition of PPP-release replenishes levels of pre-initiating Pol II at TSSs after UV. Our data show that such continuous engagement of Pol II molecules ensures maximal transcription-driven repair throughout expressed genes and regulatory loci. Importantly, revealing this unanticipated regulatory layer of UV-response provides physiological relevant traction to the emerging concept that Pol II initiation rate is determined by pause-release dynamics.

A transcription-driven cellular response is activated upon UV stress. Here the authors reveal mechanistic insights into the regulatory process affecting transcription and chromatin dynamics, showing how maintaining Pol II firing safeguards the integrity of cells’ transcriptome

Methods for the analysis of transcriptome dynamics

The transcriptome is the complete set of transcripts in a cell or tissue and includes ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), and regulatory noncoding RNA. At steady-state, the transcriptome results from a compensatory variation of the transcription and decay rate to maintain the RNA concentration constant. RNA transcription constitutes the first stage in gene expression, and thus is a major and primary mode of gene expression control. Nevertheless, regulation of RNA decay is also a key factor in gene expression control, involving either selective RNA stabilization or enhanced degradation. Transcriptome analysis allows the identification of gene expression alterations, providing new insights regarding the pathways and mechanisms involved in physiological and pathological processes. Upon perturbation of cell homeostasis, rapid changes in gene expression are required to adapt to new conditions. Thus, to better understand the regulatory mechanisms associated with gene expression alterations, it is vital to acknowledge the relative contribution of RNA synthesis and decay to the transcriptome. To the toxicology field, the study of gene expression regulation mechanisms can help identify the early and mechanistic relevant cellular events associated with a particular response. This review aims to provide a critical comparison of the available methods used to analyze the contribution of RNA transcription and decay to gene expression dynamics. Notwithstanding, an integration of the data obtained is necessary to understand the entire repercussions of gene transcription changes at a system-level. Thus, a brief overview of the methods available for the integration and analysis of the data obtained from transcriptome analysis will also be provided.

DNA Methylation Controls Metastasis-Suppressive 14q32-Encoded miRNAs

Expression of 14q32-encoded miRNAs is a favorable prognostic factor in patients with metastatic cancer. In this study, we used genomic inhibition of DNA methylation through disruption of DNA methyltransferases DNMT1 and DNMT3B and pharmacologic inhibition with 5-Aza-2'-deoxycytidine (5-Aza-dC, decitabine) to demonstrate that DNA methylation predominantly regulates expression of metastasis-suppressive miRNAs in the 14q32 cluster. DNA demethylation facilitated CCCTC-binding factor (CTCF) recruitment to the maternally expressed gene 3 differentially methylated region (MEG3-DMR), which acts as a cis-regulatory element for 14q32 miRNA expression. 5-Aza-dC activated demethylation of the MEG3-DMR and expression of 14q32 miRNAs, which suppressed adhesion, invasion, and migration (AIM) properties of metastatic tumor cells. Cancer cells with MEG3-DMR hypomethylation exhibited constitutive expression of 14q32 miRNAs and resistance to 5-Aza-dC-induced suppression of AIM. Expression of methylation-dependent 14q32 miRNAs suppressed metastatic colonization in preclinical models of lung and liver metastasis and correlated with improved clinical outcomes in patients with metastatic cancer. These findings implicate epigenetic modification via DNA methylation in the regulation of metastatic propensity through miRNA networks and identify a previously unrecognized action of decitabine on the activation of metastasis-suppressive miRNAs. SIGNIFICANCE: This study investigates epigenetic regulation of metastasis-suppressive miRNAs and the effect on metastasis.

Enhancer-derived RNA: A Primer

Enhancer-derived RNAs (eRNAs) are a group of RNAs transcribed by RNA polymerase II from the domain of transcription enhancers, a major type of cis-regulatory elements in the genome. The correlation between eRNA production and enhancer activity has stimulated studies on the potential role of eRNAs in transcriptional regulation. Additionally, eRNA has also served as a marker for global identification of enhancers. Here I review the brief history and fascinating properties of eRNAs.

Human prostate luminal cell differentiation requires NOTCH3 induction by p38-MAPK and MYC

Many pathways dysregulated in prostate cancer are also involved in epithelial differentiation. To better understand prostate tumor initiation, we sought to investigate specific genes and mechanisms required for normal basal to luminal cell differentiation. Utilizing human prostate basal epithelial cells and an in vitro differentiation model, we tested the hypothesis that regulation of NOTCH3 by the p38 MAPK family (hereafter p38-MAPK), via MYC, is required for luminal differentiation. Inhibition (SB202190 and BIRB796) or knockdown of p38α (also known as MAPK14) and/or p38δ (also known as MAPK13) prevented proper differentiation. Additionally, treatment with a γ-secretase inhibitor (RO4929097) or knockdown of NOTCH1 and/or NOTCH3 greatly impaired differentiation and caused luminal cell death. Constitutive p38-MAPK activation through MKK6(CA) increased NOTCH3 (but not NOTCH1) mRNA and protein levels, which was diminished upon MYC inhibition (10058-F4 and JQ1) or knockdown. Furthermore, we validated two NOTCH3 enhancer elements through a combination of enhancer (e)RNA detection (BruUV-seq) and luciferase reporter assays. Finally, we found that the NOTCH3 mRNA half-life increased during differentiation or upon acute p38-MAPK activation. These results reveal a new connection between p38-MAPK, MYC and NOTCH signaling, demonstrate two mechanisms of NOTCH3 regulation and provide evidence for NOTCH3 involvement in prostate luminal cell differentiation.

Transcriptional and post-transcriptional regulation of the ionizing radiation response by ATM and p53

In response to ionizing radiation (IR), cells activate a DNA damage response (DDR) pathway to re-program gene expression. Previous studies using total cellular RNA analyses have shown that the stress kinase ATM and the transcription factor p53 are integral components required for induction of IR-induced gene expression. These studies did not distinguish between changes in RNA synthesis and RNA turnover and did not address the role of enhancer elements in DDR-mediated transcriptional regulation. To determine the contribution of synthesis and degradation of RNA and monitor the activity of enhancer elements following exposure to IR, we used the recently developed Bru-seq, BruChase-seq and BruUV-seq techniques. Our results show that ATM and p53 regulate both RNA synthesis and stability as well as enhancer element activity following exposure to IR. Importantly, many genes in the p53-signaling pathway were coordinately up-regulated by both increased synthesis and RNA stability while down-regulated genes were suppressed either by reduced synthesis or stability. Our study is the first of its kind that independently assessed the effects of ionizing radiation on transcription and post-transcriptional regulation in normal human cells.

Evolution of Brain Active Gene Promoters in Human Lineage Towards the Increased Plasticity of Gene Regulation

Adaptability to a variety of environmental conditions is a prominent feature of Homo sapiens. We hypothesize that this feature can be explained by evolutionary changes in gene promoters active in the brain prefrontal cortex leading to a more flexible gene regulation network. The genotype-dependent range of gene expression can be broader in humans than in other higher primates. Thus, we searched for specific signatures of evolutionary changes in promoter architectures of multiple hominid genes, including the genes active in human cortical neurons that may indicate an increase of variability of gene expression rather than just changes in the level of expression, such as downregulation or upregulation of the genes. We performed a whole-genome search for genetic-based alterations that may impact gene regulation "flexibility" in a process of hominids evolution, such as (i) CpG dinucleotide content, (ii) predicted nucleosome-DNA dissociation constant, and (iii) predicted affinities for TATA-binding protein (TBP) in gene promoters. We tested all putative promoter regions across the human genome and especially gene promoters in active chromatin state in neurons of prefrontal cortex, the brain region critical for abstract thinking and social and behavioral adaptation. Our data imply that the origin of modern man has been associated with an increase of flexibility of promoter-driven gene regulation in brain. In contrast, after splitting from the ancestral lineages of H. sapiens, the evolution of ape species is characterized by reduced flexibility of gene promoter functioning, underlying reduced variability of the gene expression.

Gene length as a biological timer to establish temporal transcriptional regulation

Transcriptional timing is inherently influenced by gene length, thus providing a mechanism for temporal regulation of gene expression. While gene size has been shown to be important for the expression timing of specific genes during early development, whether it plays a role in the timing of other global gene expression programs has not been extensively explored. Here, we investigate the role of gene length during the early transcriptional response of human fibroblasts to serum stimulation. Using the nascent sequencing techniques Bru-seq and BruUV-seq, we identified immediate genome-wide transcriptional changes following serum stimulation that were linked to rapid activation of enhancer elements. We identified 873 significantly induced and 209 significantly repressed genes. Variations in gene size allowed for a large group of genes to be simultaneously activated but produce full-length RNAs at different times. The median length of the group of serum-induced genes was significantly larger than the median length of all expressed genes, housekeeping genes, and serum-repressed genes. These gene length relationships were also observed in corresponding mouse orthologs, suggesting that relative gene size is evolutionarily conserved. The sizes of transcription factor and microRNA genes immediately induced after serum stimulation varied dramatically, setting up a cascade mechanism for temporal expression arising from a single activation event. The retention and expansion of large intronic sequences during evolution have likely played important roles in fine-tuning the temporal expression of target genes in various cellular response programs.

Enhancers as non-coding RNA transcription units: recent insights and future perspectives

Networks of regulatory enhancers dictate distinct cell identities and cellular responses to diverse signals by instructing precise spatiotemporal patterns of gene expression. However, 35 years after their discovery, enhancer functions and mechanisms remain incompletely understood. Intriguingly, recent evidence suggests that many, if not all, functional enhancers are themselves transcription units, generating non-coding enhancer RNAs. This observation provides a fundamental insight into the inter-regulation between enhancers and promoters, which can both act as transcription units; it also raises crucial questions regarding the potential biological roles of the enhancer transcription process and non-coding enhancer RNAs. Here, we review research progress in this field and discuss several important, unresolved questions regarding the roles and mechanisms of enhancers in gene regulation.