Multiple direct and indirect roles of the Paf1 complex in transcription elongation, splicing, and histone modifications

The polymerase-associated factor 1 (Paf1) complex (Paf1C) is a conserved protein complex with critical functions during eukaryotic transcription. Previous studies showed that Paf1C is multi-functional, controlling specific aspects of transcription ranging from RNA polymerase II (RNAPII) processivity to histone modifications. However, it is unclear how specific Paf1C subunits directly impact transcription and coupled processes. We have compared conditional depletion to steady-state deletion for each Paf1C subunit to determine the direct and indirect contributions to gene expression in Saccharomyces cerevisiae. Using nascent transcript sequencing, RNAPII profiling, and modeling of transcription elongation dynamics, we have demonstrated direct effects of Paf1C subunits on RNAPII processivity and elongation rate and indirect effects on transcript splicing and repression of antisense transcripts. Further, our results suggest that the direct transcriptional effects of Paf1C cannot be readily assigned to any particular histone modification. This work comprehensively analyzes both the immediate and the extended roles of each Paf1C subunit in transcription elongation and transcript regulation.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

In eukaryotes, genes produce a variety of distinct RNA isoforms, each with potentially unique protein products, coding potential or regulatory signals such as poly(A) tail and nucleotide modifications. Assessing the kinetics of RNA isoform metabolism, such as transcription and decay rates, is essential for unraveling gene regulation. However, it is currently impeded by lack of methods that can differentiate between individual isoforms. Here, we introduce RNAkinet, a deep convolutional and recurrent neural network, to detect nascent RNA molecules following metabolic labeling with the nucleoside analog 5-ethynyl uridine and long-read, direct RNA sequencing with nanopores. RNAkinet processes electrical signals from nanopore sequencing directly and distinguishes nascent from pre-existing RNA molecules. Our results show that RNAkinet prediction performance generalizes in various cell types and organisms and can be used to quantify RNA isoform half-lives. RNAkinet is expected to enable the identification of the kinetic parameters of RNA isoforms and to facilitate studies of RNA metabolism and the regulatory elements that influence it.

YY1-controlled regulatory connectivity and transcription are influenced by the cell cycle

Few transcription factors have been examined for their direct roles in physically connecting enhancers and promoters. Here acute degradation of Yin Yang 1 (YY1) in erythroid cells revealed its requirement for the maintenance of numerous enhancer-promoter loops, but not compartments or domains. Despite its reported ability to interact with cohesin, the formation of YY1-dependent enhancer-promoter loops does not involve stalling of cohesin-mediated loop extrusion. Integrating mitosis-to-G1-phase dynamics, we observed partial retention of YY1 on mitotic chromatin, predominantly at gene promoters, followed by rapid rebinding during mitotic exit, coinciding with enhancer-promoter loop establishment. YY1 degradation during the mitosis-to-G1-phase interval revealed a set of enhancer-promoter loops that require YY1 for establishment during G1-phase entry but not for maintenance in interphase, suggesting that cell cycle stage influences YY1's architectural function. Thus, as revealed here for YY1, chromatin architectural functions of transcription factors can vary in their interplay with CTCF and cohesin as well as by cell cycle stage.

Uncovering the functions and mechanisms of regulatory elements-associated non-coding RNAs

Over the past decade, regulatory non-coding RNAs (ncRNAs) produced by RNA Pol II have been revealed as meaningful players in various essential cellular functions. In particular, thousands of ncRNAs are produced at transcriptional regulatory elements such as enhancers and promoters, where they may exert multiple functions to regulate proper development, cellular programming, transcription or genomic stability. Here, we review the mechanisms involving these regulatory element-associated ncRNAs, and particularly enhancer RNAs (eRNAs) and PROMoter uPstream Transcripts (PROMPTs). We contextualize the mechanisms described to the processing and degradation of these short lived RNAs. We summarize recent findings explaining how ncRNAs operate locally at promoters and enhancers, or further away, either shortly after their production by RNA Pol II, or through post-transcriptional stabilization. Such discoveries lead to a converging model accounting for how ncRNAs influence cellular fate, by acting on transcription and chromatin structure, which may further involve factors participating to 3D nuclear organization.

H3.3K122A results in a neomorphic phenotype in mouse embryonic stem cells

The histone variant H3.3 acts in coordination with histone posttranslational modifications and other chromatin features to facilitate appropriate transcription. Canonical histone H3 and histone variant H3.3 are post-translationally modified with the genomic distribution of these marks denoting different features and with more recent evidence suggesting that these modifications may influence transcription. While the majority of posttranslational modifications occur on histone tails, there are defined modifications within the globular domain, such as acetylation of H3K122/H3.3K122. To understand the function of the residue H3.3K122 in transcriptional regulation, we attempted to generate H3.3K122A mouse embryonic stem (mES) cells but were unsuccessful. Through multi-omic profiling of mutant cell lines harboring two or three of four H3.3 targeted alleles, we have uncovered that H3.3K122A is neomorphic and results in lethality. This is surprising as prior studies demonstrate H3.3-null mES cells are viable and pluripotent, albeit with reduced differentiation capacity. Together, these studies have uncovered a novel dependence of a globular domain residue of H3.3 for viability and broadened our understanding of how histone variants contribute to transcription regulation and pluripotency in mES cells.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

Multiple Forms and Functions of Premature Termination by RNA Polymerase II

Eukaryotic genomes are widely transcribed by RNA polymerase II (pol II) both within genes and in intergenic regions. POL II elongation complexes comprising the polymerase, the DNA template and nascent RNA transcript must be extremely processive in order to transcribe the longest genes which are over 1 megabase long and take many hours to traverse. Dedicated termination mechanisms are required to disrupt these highly stable complexes. Transcription termination occurs not only at the 3' ends of genes once a full length transcript has been made, but also within genes and in promiscuously transcribed intergenic regions. Termination at these latter positions is termed "premature" because it is not triggered in response to a specific signal that marks the 3' end of a gene, like a polyA site. One purpose of premature termination is to remove polymerases from intergenic regions where they are "not wanted" because they may interfere with transcription of overlapping genes or the progress of replication forks. Premature termination has recently been appreciated to occur at surprisingly high rates within genes where it is speculated to serve regulatory or quality control functions. In this review I summarize current understanding of the different mechanisms of premature termination and its potential functions.

The Pseudo-Natural Product Tafbromin Selectively Targets the TAF1 Bromodomain 2

Phenotypic assays detect small-molecule bioactivity at functionally relevant cellular sites, and inherently cover a variety of targets and mechanisms of action. They can uncover new small molecule-target pairs and may give rise to novel biological insights. By means of an osteoblast differentiation assay which employs a Hedgehog (Hh) signaling agonist as stimulus and which monitors an endogenous marker for osteoblasts, we identified a pyrrolo[3,4-g]quinoline (PQ) pseudo-natural product (PNP) class of osteogenesis inhibitors. The most potent PQ, termed Tafbromin, impairs canonical Hh signaling and modulates osteoblast differentiation through binding to the bromodomain 2 of the TATA-box binding protein-associated factor 1 (TAF1). Tafbromin is the most selective TAF1 bromodomain 2 ligand and promises to be an invaluable tool for the study of biological processes mediated by TAF1(2) bromodomains.

Mechanisms of RNA Polymerase II Termination at the 3'-End of Genes

RNA polymerase II (RNAPII) is responsible for the synthesis of a diverse set of RNA molecules, including protein-coding messenger RNAs (mRNAs) and many short non-coding RNAs (ncRNAs). For this purpose, RNAPII relies on a multitude of factors that regulate the transcription cycle, from initiation and promoter-proximal pausing, through elongation and finally termination. RNAPII transcription termination at the end of genes ensures the release of RNAPII from the DNA template and its efficient recycling for further rounds of transcription. Termination of RNAPII is tightly coupled to 3'-end mRNA processing, which constitutes an important trigger for the subsequent transcription termination event. In this review, we discuss the current understanding of RNAPII termination mechanisms, focusing on 'canonical' termination at the 3'-end of genes. We also integrate the allosteric and 'torpedo' models into a unified model of termination, and describe the different termination factors that have been identified to date, paying special attention to the human factors and their mechanism of action at the molecular level. Indeed, in recent years the development of novel approaches in structural biology, biochemistry and cell biology have together led to a more detailed comprehension of the different mechanisms of RNAPII termination, and a better understanding of their importance in regulating gene expression, especially under cellular stress and pathological situations.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Transcription directionality is licensed by Integrator at active human promoters

A universal characteristic of eukaryotic transcription is that the promoter recruits RNA polymerase II (RNAPII) to produce both precursor mRNAs (pre-mRNAs) and short unstable promoter upstream transcripts (PROMPTs) toward the opposite direction. However, how the transcription machinery selects the correct direction to produce pre-mRNAs is largely unknown. Here, through multiple acute auxin-inducible degradation systems, we show that rapid depletion of an RNAPII-binding protein complex, Integrator, results in robust PROMPT accumulation throughout the genome. Interestingly, the accumulation of PROMPTs is compensated by the reduction of pre-mRNA transcripts in actively transcribed genes. Consistently, Integrator depletion alters the distribution of polymerase between the sense and antisense directions, which is marked by increased RNAPII-carboxy-terminal domain Tyr1 phosphorylation at PROMPT regions and a reduced Ser2 phosphorylation level at transcription start sites. Mechanistically, the endonuclease activity of Integrator is critical to suppress PROMPT production. Furthermore, our data indicate that the presence of U1 binding sites on nascent transcripts could counteract the cleavage activity of Integrator. In this process, the absence of robust U1 signal at most PROMPTs allows Integrator to suppress the antisense transcription and shift the transcriptional balance in favor of the sense direction.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

KAP1 negatively regulates RNA polymerase II elongation kinetics to activate signal-induced transcription

Signal-induced transcriptional programs regulate critical biological processes through the precise spatiotemporal activation of Immediate Early Genes (IEGs); however, the mechanisms of transcription induction remain poorly understood. By combining an acute depletion system with several genomics approaches to interrogate synchronized, temporal transcription, we reveal that KAP1/TRIM28 is a first responder that fulfills the temporal and heightened transcriptional demand of IEGs. Acute KAP1 loss triggers an increase in RNA polymerase II elongation kinetics during early stimulation time points. This elongation defect derails the normal progression through the transcriptional cycle during late stimulation time points, ultimately leading to decreased recruitment of the transcription apparatus for re-initiation thereby dampening IEGs transcriptional output. Collectively, KAP1 plays a counterintuitive role by negatively regulating transcription elongation to support full activation across multiple transcription cycles of genes critical for cell physiology and organismal functions.

Notch induces transcription by stimulating release of paused RNA Polymerase II

Notch proteins undergo ligand-induced proteolysis to release a nuclear effector that influences a wide range of cellular processes by regulating transcription. Despite years of study, however, how Notch induces the transcription of its target genes remains unclear. Here, we comprehensively examined the response to human Notch1 across a time course of activation using high-resolution genomic assays of chromatin accessibility and nascent RNA production. Our data reveal that Notch induces target gene transcription primarily by releasing paused RNA polymerase II (RNAPII). Moreover, in contrast to prevailing models suggesting that Notch acts by promoting chromatin accessibility, we found that open chromatin was established at Notch-responsive regulatory elements prior to Notch signal induction, through SWI/SNF-mediated remodeling. Together, these studies show that the nuclear response to Notch signaling is dictated by the pre-existing chromatin state and RNAPII distribution at the time of signal activation.

Direct and selective pharmacological disruption of the YAP-TEAD interface by IAG933 inhibits Hippo-dependent and RAS-MAPK-altered cancers

The YAP-TEAD protein-protein interaction mediates YAP oncogenic functions downstream of the Hippo pathway. To date, available YAP-TEAD pharmacologic agents bind into the lipid pocket of TEAD, targeting the interaction indirectly via allosteric changes. However, the consequences of a direct pharmacological disruption of the interface between YAP and TEADs remain largely unexplored. Here, we present IAG933 and its analogs as potent first-in-class and selective disruptors of the YAP-TEAD protein-protein interaction with suitable properties to enter clinical trials. Pharmacologic abrogation of the interaction with all four TEAD paralogs resulted in YAP eviction from chromatin and reduced Hippo-mediated transcription and induction of cell death. In vivo, deep tumor regression was observed in Hippo-driven mesothelioma xenografts at tolerated doses in animal models as well as in Hippo-altered cancer models outside mesothelioma. Importantly this also extended to larger tumor indications, such as lung, pancreatic and colorectal cancer, in combination with RTK, KRAS-mutant selective and MAPK inhibitors, leading to more efficacious and durable responses. Clinical evaluation of IAG933 is underway.

Technology to the rescue: how to uncover the role of transposable elements in preimplantation development

Transposable elements (TEs) are highly expressed in preimplantation development. Preimplantation development is the phase when the cells of the early embryo undergo the first cell fate choice and change from being totipotent to pluripotent. A range of studies have advanced our understanding of TEs in preimplantation, as well as their epigenetic regulation and functional roles. However, many questions remain about the implications of TE expression during early development. Challenges originate first due to the abundance of TEs in the genome, and second because of the limited cell numbers in preimplantation. Here we review the most recent technological advancements promising to shed light onto the role of TEs in preimplantation development. We explore novel avenues to identify genomic TE insertions and improve our understanding of the regulatory mechanisms and roles of TEs and their RNA and protein products during early development.

CDK7 kinase activity promotes RNA polymerase II promoter escape by facilitating initiation factor release

Cyclin-dependent kinase 7 (CDK7), part of the general transcription factor TFIIH, promotes gene transcription by phosphorylating the C-terminal domain of RNA polymerase II (RNA Pol II). Here, we combine rapid CDK7 kinase inhibition with multi-omics analysis to unravel the direct functions of CDK7 in human cells. CDK7 inhibition causes RNA Pol II retention at promoters, leading to decreased RNA Pol II initiation and immediate global downregulation of transcript synthesis. Elongation, termination, and recruitment of co-transcriptional factors are not directly affected. Although RNA Pol II, initiation factors, and Mediator accumulate at promoters, RNA Pol II complexes can also proceed into gene bodies without promoter-proximal pausing while retaining initiation factors and Mediator. Further downstream, RNA Pol II phosphorylation increases and initiation factors and Mediator are released, allowing recruitment of elongation factors and an increase in RNA Pol II elongation velocity. Collectively, CDK7 kinase activity promotes the release of initiation factors and Mediator from RNA Pol II, facilitating RNA Pol II escape from the promoter.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

FACT maintains chromatin architecture and thereby stimulates RNA polymerase II pausing during transcription in vivo

Facilitates chromatin transcription (FACT) is a histone chaperone that supports transcription through chromatin in vitro, but its functional roles in vivo remain unclear. Here, we analyze the in vivo functions of FACT with the use of multi-omics analysis after rapid FACT depletion from human cells. We show that FACT depletion destabilizes chromatin and leads to transcriptional defects, including defective promoter-proximal pausing and elongation, and increased premature termination of RNA polymerase II. Unexpectedly, our analysis revealed that promoter-proximal pausing depends not only on the negative elongation factor (NELF) but also on the +1 nucleosome, which is maintained by FACT.

Global control of RNA polymerase II

RNA polymerase II (Pol II) is the multi-protein complex responsible for transcribing all protein-coding messenger RNA (mRNA). Most research on gene regulation is focused on the mechanisms controlling which genes are transcribed when, or on the mechanics of transcription. How global Pol II activity is determined receives comparatively less attention. Here, we follow the life of a Pol II molecule from 'assembly of the complex' to nuclear import, enzymatic activity, and degradation. We focus on how Pol II spends its time in the nucleus, and on the two-way relationship between Pol II abundance and activity in the context of homeostasis and global transcriptional changes.

Enhancing Transcriptome Mapping with Rapid PRO-seq Profiling of Nascent RNA

Precision nuclear run-on (PRO) sequencing (PRO-seq) is a powerful technique for mapping polymerase active sites with nucleotide resolution and measuring newly synthesized transcripts at both promoters and enhancer elements. The current PRO-seq protocol is time-intensive, technically challenging, and requires a large amount of starting material. To overcome these limitations, we developed rapid PRO-seq (rPRO-seq) which utilizes pre-adenylated single-stranded DNAs (AppDNA), a dimer blocking oligonucleotide (DBO), on-bead 5' RNA end repair, and column-based purification. These modifications enabled efficient transcriptome mapping within a single day (∼12 hours) increasing ligation efficiency, abolished adapter dimers, and reduced sample loss and RNA degradation. We demonstrate the reproducibility of rPRO-seq in measuring polymerases at promoters, gene bodies, and enhancers as compared to original PRO-seq protocols. Additionally, rPRO-seq is scalable, allowing for transcriptome mapping with as little as 25,000 cells. We apply rPRO-seq to study the role of Integrator in mouse hematopoietic stem and progenitor cell (mHSPC) homeostasis, identifying Ints11 as an essential component of transcriptional regulation and RNA processing in mHSPC homeostasis. Overall, rPRO-seq represents a significant advance in the field of nascent transcript analyses and will be a valuable tool for generating patient-specific genome-wide transcription profiles with minimal sample requirements.

Gene-specific RNA homeostasis revealed by perturbation of coactivator complexes

Transcript buffering entails the reciprocal modulation of mRNA synthesis and degradation rates to maintain stable RNA levels under varying cellular conditions. Current research supports a global, non-sequence-specific connection between mRNA synthesis and degradation, but the underlying mechanisms are still unclear. In this study, we investigated changes in RNA metabolism following acute depletion of TIP60/KAT5, the acetyltransferase subunit of the NuA4 transcriptional coactivator complex, in mouse embryonic stem cells. By combining RNA sequencing of nuclear, cytoplasmic, and newly synthesised transcript fractions with biophysical modelling, we demonstrate that TIP60 predominantly enhances transcription of numerous genes, while a smaller set of genes undergoes TIP60-dependent transcriptional repression. Surprisingly, transcription changes caused by TIP60 depletion were offset by corresponding changes in RNA nuclear export and cytoplasmic stability, indicating gene-specific buffering mechanisms. Similarly, disruption of the unrelated ATAC coactivator complex also resulted in gene-specific transcript buffering. These findings reveal that transcript buffering functions at a gene-specific level and suggest that cells dynamically adjust RNA splicing, export, and degradation in response to individual RNA synthesis alterations, thereby sustaining cellular homeostasis.

KAP1 negatively regulates RNA polymerase II elongation kinetics to activate signal-induced transcription

Signal-induced transcriptional programs regulate critical biological processes through the precise spatiotemporal activation of Immediate Early Genes (IEGs); however, the mechanisms of transcription induction remain poorly understood. By combining an acute depletion system with high resolution genomics approaches to interrogate synchronized, temporal transcription, we reveal that KAP1/TRIM28 is a first responder that fulfills the temporal and heightened transcriptional demand of IEGs. Unexpectedly, acute KAP1 loss triggers an increase in RNA polymerase II elongation kinetics during early stimulation time points. This elongation defect derails the normal progression through the transcriptional cycle during late stimulation time points, ultimately leading to decreased recruitment of the transcription apparatus for re-initiation thereby dampening IEGs transcriptional output. Collectively, KAP1 plays a counterintuitive role by negatively regulating transcription elongation to support full activation across multiple transcription cycles of genes critical for cell physiology and organismal functions.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

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.

Nuclear export is a limiting factor in eukaryotic mRNA metabolism

The eukaryotic mRNA life cycle includes transcription, nuclear mRNA export and degradation. To quantify all these processes simultaneously, we perform thiol-linked alkylation after metabolic labeling of RNA with 4-thiouridine (4sU), followed by sequencing of RNA (SLAM-seq) in the nuclear and cytosolic compartments of human cancer cells. We develop a model that reliably quantifies mRNA-specific synthesis, nuclear export, and nuclear and cytosolic degradation rates on a genome-wide scale. We find that nuclear degradation of polyadenylated mRNA is negligible and nuclear mRNA export is slow, while cytosolic mRNA degradation is comparatively fast. Consequently, an mRNA molecule generally spends most of its life in the nucleus. We also observe large differences in the nuclear export rates of different 3'UTR transcript isoforms. Furthermore, we identify genes whose expression is abruptly induced upon metabolic labeling. These transcripts are exported substantially faster than average mRNAs, suggesting the existence of alternative export pathways. Our results highlight nuclear mRNA export as a limiting factor in mRNA metabolism and gene regulation.

Multiple direct and indirect roles of Paf1C in elongation, splicing, and histone post-translational modifications

Paf1C is a highly conserved protein complex with critical functions during eukaryotic transcription. Previous studies have shown that Paf1C is multi-functional, controlling specific aspects of transcription, ranging from RNAPII processivity to histone modifications. However, it is unclear how specific Paf1C subunits directly impact transcription and coupled processes. We have compared conditional depletion to steady-state deletion for each Paf1C subunit to determine the direct and indirect contributions to gene expression in Saccharomyces cerevisiae. Using nascent transcript sequencing, RNAPII profiling, and modeling of transcription elongation dynamics, we have demonstrated direct effects of Paf1C subunits on RNAPII processivity and elongation rate and indirect effects on transcript splicing and repression of antisense transcripts. Further, our results suggest that the direct transcriptional effects of Paf1C cannot be readily assigned to any particular histone modification. This work comprehensively analyzes both the immediate and extended roles of each Paf1C subunit in transcription elongation and transcript regulation.

RBM22 regulates RNA polymerase II 5' pausing, elongation rate, and termination by coordinating 7SK-P-TEFb complex and SPT5

BACKGROUND

Splicing factors are vital for the regulation of RNA splicing, but some have also been implicated in regulating transcription. The underlying molecular mechanisms of their involvement in transcriptional processes remain poorly understood.

RESULTS

Here, we describe a direct role of splicing factor RBM22 in coordinating multiple steps of RNA Polymerase II (RNAPII) transcription in human cells. The RBM22 protein widely occupies the RNAPII-transcribed gene locus in the nucleus. Loss of RBM22 promotes RNAPII pause release, reduces elongation velocity, and provokes transcriptional readthrough genome-wide, coupled with production of transcripts containing sequences from downstream of the gene. RBM22 preferentially binds to the hyperphosphorylated, transcriptionally engaged RNAPII and coordinates its dynamics by regulating the homeostasis of the 7SK-P-TEFb complex and the association between RNAPII and SPT5 at the chromatin level.

CONCLUSIONS

Our results uncover the multifaceted role of RBM22 in orchestrating the transcriptional program of RNAPII and provide evidence implicating a splicing factor in both RNAPII elongation kinetics and termination control.

Widespread impact of nucleosome remodelers on transcription at cis-regulatory elements

Nucleosome remodeling complexes and other regulatory factors work in concert to build a chromatin environment that directs the expression of a distinct set of genes in each cell using cis-regulatory elements (CREs), such as promoters and enhancers, that drive transcription of both mRNAs and CRE-associated non-coding RNAs (ncRNAs). Two classes of CRE-associated ncRNAs include upstream antisense RNAs (uaRNAs), which are transcribed divergently from a shared mRNA promoter, and enhancer RNAs (eRNAs), which are transcribed bidirectionally from active enhancers. The complicated network of CRE regulation by nucleosome remodelers remains only partially explored, with a focus on a select, limited number of remodelers. We endeavored to elucidate a remodeler-based regulatory network governing CRE-associated transcription (mRNA, eRNA, and uaRNA) in murine embryonic stem (ES) cells to test the hypothesis that many SNF2-family nucleosome remodelers collaborate to regulate the coding and non-coding transcriptome via alteration of underlying nucleosome architecture. Using depletion followed by transient transcriptome sequencing (TT-seq), we identified thousands of misregulated mRNAs and CRE-associated ncRNAs across the remodelers examined, identifying novel contributions by understudied remodelers in the regulation of coding and noncoding transcription. Our findings suggest that mRNA and eRNA transcription are coordinately co-regulated, while mRNA and uaRNAs sharing a common promoter are independently regulated. Subsequent mechanistic studies suggest that while remodelers SRCAP and CHD8 modulate transcription through classical mechanisms such as transcription factors and histone variants, a broad set of remodelers including SMARCAL1 indirectly contribute to transcriptional regulation through maintenance of genomic stability and proper Integrator complex localization. This study systematically examines the contribution of SNF2-remodelers to the CRE-associated transcriptome, identifying at least two classes for remodeler action.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

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.

Development and Comparison of 4-Thiouridine to Cytidine Base Conversion Reaction

To study transcriptome dynamics without harming cells, it is essential to convert chemical bases. 4-Thiouridine (4sU) is a biocompatible uridine analogue that can be converted into a cytidine analogue. Although several reactions can convert 4sU into a cytidine analogue, few studies have compared the features of these reactions. In this study, we performed three reported base conversion reactions, including osmium tetroxide, iodoacetamide, and sodium periodate treatment, as well as a new reaction using 2,4-dinitrofluorobenzene. We compared the reaction time, conversion efficacy, and effects on reverse transcription. These reactions successfully converted 4sU into a cytidine analogue quantitatively using trinucleotides. However, the conversion efficacy and effect on reverse transcription vary depending on the reaction with the RNA transcript. OsO4 treatment followed by NH4Cl treatment showed the best base-conversion efficiency. Nevertheless, each reaction has its own advantages and disadvantages as a tool for studying the transcriptome. Therefore, it is crucial to select the appropriate reaction for the target of interest.

Genome-wide kinetic profiling of pre-mRNA 3' end cleavage

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

hGRAD: A versatile "one-fits-all" system to acutely deplete RNA binding proteins from condensates

Nuclear RNA binding proteins (RBPs) are difficult to study because they often belong to large protein families and form extensive networks of auto- and crossregulation. They are highly abundant and many localize to condensates with a slow turnover, requiring long depletion times or knockouts that cannot distinguish between direct and indirect or compensatory effects. Here, we developed a system that is optimized for the rapid degradation of nuclear RBPs, called hGRAD. It comes as a "one-fits-all" plasmid, and integration into any cell line with endogenously GFP-tagged proteins allows for an inducible, rapid, and complete knockdown. We show that the nuclear RBPs SRSF3, SRSF5, SRRM2, and NONO are completely cleared from nuclear speckles and paraspeckles within 2 h. hGRAD works in various cell types, is more efficient than previous methods, and does not require the expression of exogenous ubiquitin ligases. Combining SRSF5 hGRAD degradation with Nascent-seq uncovered transient transcript changes, compensatory mechanisms, and an effect of SRSF5 on transcript stability.

Interplay between the transcription preinitiation complex and the +1 nucleosome

Eukaryotic transcription starts with the assembly of a preinitiation complex (PIC) on core promoters. Flanking this region is the +1 nucleosome, the first nucleosome downstream of the core promoter. While this nucleosome is rich in epigenetic marks and plays a key role in transcription regulation, how the +1 nucleosome interacts with the transcription machinery has been a long-standing question. Here, we summarize recent structural and functional studies of the +1 nucleosome in complex with the PIC. We specifically focus on how differently organized promoter-nucleosome templates affect the assembly of the PIC and PIC-Mediator on chromatin and result in distinct transcription initiation.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

KAS-seq profiling captures transcription dynamics during oocyte maturation

In fully grown oocytes, the genome is considered to be globally transcriptionally silenced. However, this conclusion is primarily derived from the results obtained through immunofluorescence staining or inferred from the highly condensed state of chromosomes, lacking more direct evidence. Here, by using a kethoxal-assisted single-stranded DNA sequencing (KAS-seq) approach, we investigated the landscape of single-strand DNA (ssDNA) throughout the genome and provided a readout of the activity and dynamics of transcription during oocyte meiotic maturation. In non-surrounded nucleolus (NSN) oocytes, we observed a robust KAS-seq signal, indicating the high transcriptional activity. In surrounded nucleolus (SN) oocytes, the presence of ssDNA still persists although the KAS-seq signal was relatively weak, suggesting the presence of transcription. Accompanying with the meiotic resumption, the transcriptional activity gradually decreased, and global repression was detected in matured oocytes. Moreover, we preformed the integrative genomics analysis to dissect the transcriptional dynamics during mouse oocyte maturation. In sum, the present study delineates the detailed transcriptional activity during mammalian oocyte maturation.

Internal and external normalization of nascent RNA sequencing run-on experiments

In experiments with significant perturbations to transcription, nascent RNA sequencing protocols are dependent on external spike-ins for reliable normalization. Unlike in RNA-seq, these spike-ins are not standardized and, in many cases, depend on a run-on reaction that is assumed to have constant efficiency across samples. To assess the validity of this assumption, we analyze a large number of published nascent RNA spike-ins to quantify their variability across existing normalization methods. Furthermore, we develop a new biologically-informed Bayesian model to estimate the error in spike-in based normalization estimates, which we term Virtual Spike-In (VSI). We apply this method both to published external spike-ins as well as using reads at the [Formula: see text] end of long genes, building on prior work from Mahat (Mol Cell 62(1):63-78, 2016. https://doi.org/10.1016/j.molcel.2016.02.025 ) and Vihervaara (Nat Commun 8(1):255, 2017. https://doi.org/10.1038/s41467-017-00151-0 ). We find that spike-ins in existing nascent RNA experiments are typically under sequenced, with high variability between samples. Furthermore, we show that these high variability estimates can have significant downstream effects on analysis, complicating biological interpretations of results.

PSIP1/LEDGF reduces R-loops at transcription sites to maintain genome integrity

R-loops that accumulate at transcription sites pose a persistent threat to genome integrity. PSIP1 is a chromatin protein associated with transcriptional elongation complex, possesses histone chaperone activity, and is implicated in recruiting RNA processing and DNA repair factors to transcription sites. Here, we show that PSIP1 interacts with R-loops and other proteins involved in R-loop homeostasis, including PARP1. Genome-wide mapping of PSIP1, R-loops and γ-H2AX in PSIP1-depleted human and mouse cell lines revealed an accumulation of R-loops and DNA damage at gene promoters in the absence of PSIP1. R-loop accumulation causes local transcriptional arrest and transcription-replication conflict, leading to DNA damage. PSIP1 depletion increases 53BP1 foci and reduces RAD51 foci, suggesting altered DNA repair choice. Furthermore, PSIP1 depletion increases the sensitivity of cancer cells to PARP1 inhibitors and DNA-damaging agents that induce R-loop-induced DNA damage. These findings provide insights into the mechanism through which PSIP1 maintains genome integrity at the site of transcription.

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

PRO-seq: Precise Mapping of Engaged RNA Pol II at Single-Nucleotide Resolution

Gene regulation is dependent on the production of mRNAs and a repertoire of non-coding RNAs by RNA polymerase II (RNAPII). Precision run-on sequencing (PRO-seq) maps the position of engaged RNAPII complexes at single-nucleotide resolution and can reveal direct targets of regulation, locations of enhancers, and transcription mechanisms that are difficult or impossible to measure by analysis of total cellular RNA. Briefly, this method first involves permeabilizing cells with mild detergents to remove intracellular NTPs and halt transcription. Transcription is then resumed in the presence of biotin-NTPs and sarkosyl to allow transcriptional incorporation of a single biotinylated NTP by RNAPII. The biotin moiety is then bound to streptavidin beads to stringently enrich for nascent RNAs. Sequencing libraries are then generated such that the first base read corresponds to the 3' end of the nascent transcript. Here, we describe our current protocol for generating PRO-seq libraries from metazoan cells, including adaptations of previously published protocols to incorporate unique molecular identifiers, reduce ligation bias, and improve library yields. Additional commentary describes quality control and processing of PRO-seq data and references for more advanced downstream analysis such as gene and enhancer identification. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Cell permeabilization for PRO-seq Basic Protocol 2: Construction of PRO-seq libraries Support Protocol: Adenylation of 3' adapter.

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

High-throughput detection of nascent RNA is critical for studies of transcription and much more challenging than that of mRNA. Recently, several massively parallel nascent RNA sequencing methods were established in eukaryotic cells. Here, we systematically compared 3 classes of methods on the same pure or crude nuclei preparations: GRO-seq for sequence nuclear run-on RNAs, pNET-seq for sequence RNA polymerase II-associated RNAs, and CB RNA-seq for sequence chromatin-bound (CB) RNAs in Arabidopsis (Arabidopsis thaliana). To improve the resolution of CB RNAs, 3'CB RNA-seq was established to sequence the 3' ends of CB RNAs. In addition, we modified pNET-seq to establish the Chromatin Native Elongation Transcript sequencing (ChrNET) method using chromatin as the starting material for RNA immunoprecipitation. Reproducibility, sensitivity and accuracy in detecting nascent transcripts, experimental procedures, and costs were analyzed, which revealed the strengths and weaknesses of each method. We found that pNET and GRO methods best detected active RNA polymerase II. CB RNA-seq is a simple and cost-effective alternative for nascent RNA studies, due to its high correlation with pNET-seq and GRO-seq. Compared with pNET, ChrNET has higher specificity for nascent RNA capture and lower sequencing cost. 3'CB is sensitive to transcription-coupled splicing. Using these methods, we identified 1,404 unknown transcripts, 4,482 unannotated splicing events, and 60 potential recursive splicing events. This comprehensive comparison of different nascent/chromatin RNA sequencing methods highlights the strengths of each method and serves as a guide for researchers aiming to select a method that best meets their study goals.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

Quantification of the dynamics of RNA metabolism is essential for understanding gene regulation in health and disease. Existing methods rely on metabolic labeling of nascent RNAs and physical separation or inference of labeling through PCR-generated mutations, followed by short-read sequencing. However, these methods are limited in their ability to identify transient decay intermediates or co-analyze RNA decay with cis-regulatory elements of RNA stability such as poly(A) tail length and modification status, at single molecule resolution. Here we use 5-ethynyl uridine (5EU) to label nascent RNA followed by direct RNA sequencing with nanopores. We developed RNAkinet, a deep convolutional and recurrent neural network that processes the electrical signal produced by nanopore sequencing to identify 5EU-labeled nascent RNA molecules. RNAkinet demonstrates generalizability to distinct cell types and organisms and reproducibly quantifies RNA kinetic parameters allowing the combined interrogation of RNA metabolism and cis-acting RNA regulatory elements.

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.

What is an enhancer?

Tight control of the transcription process is essential for the correct spatial and temporal gene expression pattern during development and in homeostasis. Enhancers are at the core of correct transcriptional activation. The original definition of an enhancer is straightforward: a DNA sequence that activates transcription independent of orientation and direction. Dissection of numerous enhancer loci has shown that many enhancer-like elements might not conform to the original definition, suggesting that enhancers and enhancer-like elements might use multiple different mechanisms to contribute to transcriptional activation. Here, we review methodologies to identify enhancers and enhancer-like elements and discuss pitfalls and consequences for our understanding of transcriptional regulation.

Genome-wide Mapping of 5'-monophosphorylated Ends of Mammalian Nascent RNA Transcripts

In eukaryotic cells, RNA biogenesis generally requires processing of the nascent transcript as it is being synthesized by RNA polymerase. These processing events include endonucleolytic cleavage, exonucleolytic trimming, and splicing of the growing nascent transcript. Endonucleolytic cleavage events that generate an exposed 5'-monophosphorylated (5'-PO4) end on the growing nascent transcript occur in the maturation of rRNAs, tRNAs, and mRNAs. These 5'-PO4 ends can be a target of further processing or be subjected to 5'-3' exonucleolytic digestion that may result in termination of transcription. Here, we describe how to identify 5'-PO4 ends of intermediates in nascent RNA metabolism. We capture these species via metabolic labeling with bromouridine followed by immunoprecipitation and specific ligation of 5'-PO4 RNA ends with the 3'-hydroxyl group of a 5' adaptor (5'-PO4 Bru-Seq) using RNA ligase I. These ligation events are localized at single nucleotide resolution via highthroughput sequencing, which identifies the position of 5'-PO4 groups precisely. This protocol successfully detects the 5'monophosphorylated ends of RNA processing intermediates during production of mature ribosomal, transfer, and micro RNAs. When combined with inhibition of the nuclear 5'-3' exonuclease Xrn2, 5'-PO4 Bru-Seq maps the 5' splice sites of debranched introns and mRNA and tRNA 3' end processing sites cleaved by CPSF73 and RNaseZ, respectively. Key features • Metabolic labeling for brief periods with bromouridine focuses the analysis of 5'-PO4 RNA ends on the population of nascent transcripts that are actively transcribed. • Detects 5'-PO4 RNA ends on nascent transcripts produced by all RNA polymerases. • Detects 5'-PO4 RNA ends at single nucleotide resolution.