Nascent RNA Sequencing Research Articles

A role for condensin in mediating transcriptional adaptation to environmental stimuli.

Nuclear organisation shapes gene regulation; however, the principles by which three-dimensional genome architecture influences gene transcription are incompletely understood. Condensin is a key architectural chromatin constituent, best known for its role in mitotic chromosome condensation. Yet at least a subset of condensin is bound to DNA throughout the cell cycle. Studies in various organisms have reported roles for condensin in transcriptional regulation, but no unifying mechanism has emerged. Here, we use rapid conditional condensin depletion in the budding yeast Saccharomyces cerevisiae to study its role in transcriptional regulation. We observe a large number of small gene expression changes, enriched at genes located close to condensin-binding sites, consistent with a possible local effect of condensin on gene expression. Furthermore, nascent RNA sequencing reveals that transcriptional down-regulation in response to environmental stimuli, in particular to heat shock, is subdued without condensin. Our results underscore the multitude by which an architectural chromosome constituent can affect gene regulation and suggest that condensin facilitates transcriptional reprogramming as part of adaptation to environmental changes.

Nascent RNA sequencing identifies a widespread sigma70-dependent pausing regulated by Gre factors in bacteria

Promoter-proximal pausing regulates eukaryotic gene expression and serves as checkpoints to assemble elongation/splicing machinery. Little is known how broadly this type of pausing regulates transcription in bacteria. We apply nascent elongating transcript sequencing combined with RNase I footprinting for genome-wide analysis of σ70-dependent transcription pauses in Escherichia coli. Retention of σ70 induces strong backtracked pauses at a 10−20-bp distance from many promoters. The pauses in the 10−15-bp register of the promoter are dictated by the canonical −10 element, 6−7 nt spacer and “YR+1Y” motif centered at the transcription start site. The promoters for the pauses in the 16−20-bp register contain an additional −10-like sequence recognized by σ70. Our in vitro analysis reveals that DNA scrunching is involved in these pauses relieved by Gre cleavage factors. The genes coding for transcription factors are enriched in these pauses, suggesting that σ70 and Gre proteins regulate transcription in response to changing environmental cues.

Fast and furious: insights of back splicing regulation during nascent RNA synthesis.

Alternative splicing of eukaryotic precursor (messenger) RNAs in the nucleus not only increases transcriptomic complexity, but also expands proteomic and functional diversity. In addition to basic types of alternative splicing, recent transcriptome-wide analyses have also suggested other new types of non-canonical splicing, such as back splicing and recursive splicing, and their widespread expression across species Increasing lines of evidence have suggested mechanisms for back splicing, including insights from analyses of nascent RNA sequencing. In this review, we discuss our current understanding of back splicing regulation, and highlight its distinct characteristics in processing during nascent RNA synthesis by taking advantage of metabolic tagging nascent RNA sequencing. Features of recursive splicing are also discussed in the perspective of nascent RNA sequencing.

Revealing nascent RNA processing dynamics with nano-COP.

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

Population-scale study of eRNA transcription reveals bipartite functional enhancer architecture

Enhancer RNAs (eRNA) are unstable non-coding RNAs, transcribed bidirectionally from active regulatory sequences, whose expression levels correlate with enhancer activity. We use capped-nascent-RNA sequencing to efficiently capture bidirectional transcription initiation across several human lymphoblastoid cell lines (Yoruba population) and detect ~75,000 eRNA transcription sites with high sensitivity and specificity. The use of nascent-RNA sequencing sidesteps the confounding effect of eRNA instability. We identify quantitative trait loci (QTLs) associated with the level and directionality of eRNA expression. High-resolution analyses of these two types of QTLs reveal distinct positions of enrichment at the central transcription factor (TF) binding regions and at the flanking eRNA initiation regions, both of which are associated with mRNA expression QTLs. These two regions—the central TF-binding footprint and the eRNA initiation cores—define a bipartite architecture of enhancers, inform enhancer function, and can be used as an indicator of the significance of non-coding regulatory variants.

Widespread Transcriptional Readthrough Caused by Nab2 Depletion Leads to Chimeric Transcripts with Retained Introns

Nascent RNA sequencing has revealed that pre-mRNA splicing can occur shortly after introns emerge from RNA polymerase II (RNA Pol II). Differences in co-transcriptional splicing profiles suggest regulation by cis- and/or trans-acting factors. Here, we use single-molecule intron tracking (SMIT) to identify a cohort of regulators by machine learning in budding yeast. Of these, Nab2 displays reduced co-transcriptional splicing when depleted. Unexpectedly, these splicing defects are attributable to aberrant "intrusive" transcriptional readthrough from upstream genes, as revealed by long-read sequencing. Transcripts that originate from the intron-containing gene's own transcription start site (TSS) are efficiently spliced, indicating no direct role of Nab2 in splicing per se. This work highlights the coupling between transcription, splicing, and 3' end formation in the context of gene organization along chromosomes. We conclude that Nab2 is required for proper 3' end processing, which ensures gene-specific control of co-transcriptional RNA processing.

Characterization of Aminobenzylphenols as Protein Disulfide Isomerase Inhibitors in Glioblastoma Cell Lines.

Disulfide bond formation is a critical post-translational modification of newly synthesized polypeptides in the oxidizing environment of the endoplasmic reticulum and is mediated by protein disulfide isomerase (PDIA1). In this study, we report a series of α-aminobenzylphenol analogues as potent PDI inhibitors. The lead compound, AS15, is a covalent nanomolar inhibitor of PDI, and the combination of AS15 analogues with glutathione synthesis inhibitor buthionine sulfoximine (BSO) leads to synergistic cell growth inhibition. Using nascent RNA sequencing, we show that an AS15 analogue triggers the unfolded protein response in glioblastoma cells. A BODIPY-labeled analogue binds proteins including PDIA1, suggesting that the compounds are cell-permeable and reach the intended target. Taken together, these findings demonstrate an extensive biochemical characterization of a novel series of highly potent reactive small molecules that covalently bind to PDI.

Abstract 5854: New therapeutic targets of ER positive breast cancer are revealed through nascent RNA transcription time series

Abstract Mechanistic understanding of Estrogen Receptor (ER) signaling will lead to new therapeutic targets for ER positive breast cancer. To this end, Arpeggio performed nascent RNA sequencing following estradiol (E2) treatment in MCF7 breast cancer cells every 15 minutes for 6 hours. To date, our study represents the densest time series of estrogen signaling ever reported. Arpeggio noted three transcriptional waves following E2 treatment: peaks at 30 minutes, 2.5 hours and 5 hours, respectively. We noted that the early response was enriched for canonical ER signaling, however, mid and late response were enriched for metabolic processes. Because our study was statistically powered with 24 time points, we could use delay-coordinate embedding to look for causal gene interactions. We implicated NME2 as the causal driver of the 5-hour transcriptional wave. Indeed, NME2 and Estrogen Receptor (ER) have been shown to co-localize where NME2 blocks ER effect. Our study represents a new methodology for target discovery in hormone-driven cancers. Citation Format: Joseph Azofeifa, Joel Basken, Maria Lai, Ryan Langendorf, Laura Norris, Tim Read, Adam Robbins-Pianka. New therapeutic targets of ER positive breast cancer are revealed through nascent RNA transcription time series [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5854.

RAI1 Regulates Activity-Dependent Nascent Transcription and Synaptic Scaling

Long-lasting forms of synaptic plasticity such as synaptic scaling are critically dependent on transcription. Activity-dependent transcriptional dynamics in neurons, however, remain incompletely characterized because most previous efforts relied on measurement of steady-state mRNAs. Here, we use nascent RNA sequencing to profile transcriptional dynamics of primary neuron cultures undergoing network activity shifts. We find pervasive transcriptional changes, in which ∼45% of expressed genes respond to network activity shifts. We further link retinoic acid-induced 1 (RAI1), the Smith-Magenis syndrome gene, to the transcriptional program driven by reduced network activity. Remarkable agreement among nascent transcriptomes, dynamic chromatin occupancy of RAI1, and electrophysiological properties of Rai1-deficient neurons demonstrates the essential roles of RAI1 in suppressing synaptic upscaling in the naive network, while promoting upscaling triggered by activity silencing. These results highlight the utility of bona fide transcription profiling to discover mechanisms of activity-dependent chromatin remodeling that underlie normal and pathological synaptic plasticity.

Acrylonitrile-Mediated Nascent RNA Sequencing for Transcriptome-Wide Profiling of Cellular RNA Dynamics.

RNA sequencing has greatly facilitated gene expression studies but is weak in studying temporal RNA dynamics; this issue can be addressed by analyzing nascent RNAs. A famous method for nascent RNA analysis is metabolic labeling with noncanonic nucleoside followed by affinity purification, however, purification processes can always introduce biases into data analysis. Here, a chemical method for nascent RNA sequencing that avoids affinity purification based on acrylonitrile‐mediated uridine‐to‐cytidine (U‐to‐C) conversion (AMUC‐seq) via 4‐thiouridine (s4U) cyanoethylation is presented. This method converts s4U base‐pairing with guanine through the nucleophilic addition of s4U to acrylonitrile. The high reaction efficiency permits AMUC‐seq directly and efficiently to recover nascent RNA information from total RNAs. AMUC‐seq is validated by being used to detect mRNA half‐lives and investigating the direct gene targets of a G‐quadruplex stabilizer, which can be regarded as potential anticancer drug, in human cells. Thousands of direct gene targets of this drug are verified (these genes are significantly enriched in cancer such as SRC and HRAS). AMUC‐seq also confirms G‐quadruplex stabilization that impacts RNA polyadenylation. These results show AMUC‐seq is qualified for the study of temporal RNA dynamics, and it can be a promising strategy to study the therapeutic mechanism of transcription‐modulating drugs.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores.

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Abstract LB-B15: New causal drivers of estrogen signaling revealed by dense time series of nascent RNA transcription

Abstract Mechanistic understanding of Estrogen Receptor (ER) signaling will lead to new therapeutic targets for ER positive breast cancer. To this end, Arpeggio performed nascent RNA sequencing following estradiol (E2) treatment in MCF7 breast cancer cells every 15 minutes for 6 hours. To date, our study represents the densest time series of estrogen signaling ever reported. Arpeggio noted three transcriptional waves following E2 treatment: peaks at 30 minutes, 2.5 hours and 5 hours respectively. We noted that the early response was enriched for canonical ER signaling, however mid and late response were enriched for metabolic processes. Because our study was statistically powered with 24 time points, we could use delay-coordinate embedding to look for causal gene interactions. We implicated NME2 as the causal driver of the 5-hour transcriptional wave. Indeed, NME2 and Estrogen Receptor (ER) have been shown to co-localize where NME2 blocks ER effect. Our study represents a new methodology for target discovery in hormone driven cancers. Citation Format: Ryan Langendorf, Joel Basken, Maria Lai, Joseph Azofeifa. New causal drivers of estrogen signaling revealed by dense time series of nascent RNA transcription [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference on Molecular Targets and Cancer Therapeutics; 2019 Oct 26-30; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2019;18(12 Suppl):Abstract nr LB-B15. doi:10.1158/1535-7163.TARG-19-LB-B15

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.

DNA damage and transcriptional regulation in iPSC-derived neurons from Ataxia Telangiectasia patients

Ataxia Telangiectasia (A-T) is neurodegenerative syndrome caused by inherited mutations inactivating the ATM kinase, a master regulator of the DNA damage response (DDR). What makes neurons vulnerable to ATM loss remains unclear. In this study we assessed on human iPSC-derived neurons whether the abnormal accumulation of DNA-Topoisomerase 1 adducts (Top1ccs) found in A-T impairs transcription elongation, thus favoring neurodegeneration. Furthermore, whether neuronal activity-induced immediate early genes (IEGs), a process involving the formation of DNA breaks, is affected by ATM deficiency. We found that Top1cc trapping by CPT induces an ATM-dependent DDR as well as an ATM-independent induction of IEGs and repression especially of long genes. As revealed by nascent RNA sequencing, transcriptional elongation and recovery were found to proceed with the same rate, irrespective of gene length and ATM status. Neuronal activity induced by glutamate receptors stimulation, or membrane depolarization with KCl, triggered a DDR and expression of IEGs, the latter independent of ATM. In unperturbed A-T neurons a set of genes (FN1, DCN, RASGRF1, FZD1, EOMES, SHH, NR2E1) implicated in the development, maintenance and physiology of central nervous system was specifically downregulated, underscoring their potential involvement in the neurodegenerative process in A-T patients.

Inhibition of protein disulfide isomerase in glioblastoma causes marked downregulation of DNA repair and DNA damage response genes.

Aberrant overexpression of endoplasmic reticulum (ER)-resident oxidoreductase protein disulfide isomerase (PDI) plays an important role in cancer progression. In this study, we demonstrate that PDI promotes glioblastoma (GBM) cell growth and describe a class of allosteric PDI inhibitors that are selective for PDI over other PDI family members.Methods: We performed a phenotypic screening triage campaign of over 20,000 diverse compounds to identify PDI inhibitors cytotoxic to cancer cells. From this screen, BAP2 emerged as a lead compound, and we assessed BAP2-PDI interactions with gel filtration, thiol-competition assays, and site-directed mutagenesis studies. To assess selectivity, we compared BAP2 activity across several PDI family members in the PDI reductase assay. Finally, we performed in vivo studies with a mouse xenograft model of GBM combining BAP2 and the standard of care (temozolomide and radiation), and identified affected gene pathways with nascent RNA sequencing (Bru-seq).Results: BAP2 and related analogs are novel PDI inhibitors that selectively inhibit PDIA1 and PDIp. Though BAP2 contains a weak Michael acceptor, interaction with PDI relies on Histidine 256 in the b' domain of PDI, suggesting allosteric binding. Furthermore, both in vitro and in vivo, BAP2 reduces cell and tumor growth. BAP2 alters the transcription of genes involved in the unfolded protein response, ER stress, apoptosis and DNA repair response.Conclusion: These results indicate that BAP2 has anti-tumor activity and the suppressive effect on DNA repair gene expression warrants combination with DNA damaging agents to treat GBM.

Dynamic Imaging of Nascent RNA Reveals General Principles of Transcription Dynamics And Stochastic Splice Site Selection

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex and nascent RNA sequencing further reveals widespread splice site selection within introns. We propose and validate a unified theoretical model to explain the general feature of transcription and pervasive stochastic splice site selection.

Single-molecule nascent RNA sequencing identifies regulatory domain architecture at promoters and enhancers.

Eukaryotic RNA Polymerase II has been found at both promoters and distal enhancers, suggesting additional functions beyond mRNA production. To understand this role, we sequenced nascent RNAs at single-molecule resolution to unravel the interplay between Pol2 initiation, capping, and pausing genome-wide. Our analyses reveal two pause classes that are associated with different RNA capping profiles. More proximal pausing is associated with less complete capping, less elongation, and a more enhancer-like complement of transcription factors than later pausing. Unexpectedly, Transcription Start Sites (TSSes) are predominantly found in constellations composed of multiple divergent pairs. TSS clusters are intimately associated with precise arrays of nucleosomes, and correspond with boundaries of transcription factor binding and chromatin modification at promoters and enhancers. TSS architecture is remarkably similar after the dramatic transcriptional changes caused by heat‬ stress. Together, our results suggest that promoter- and enhancer-associated Pol2 is a regulatory nexus for integrating information across TSS ensembles.

Numerous recursive sites contribute to accuracy of splicing in long introns in flies.

Recursive splicing, a process by which a single intron is removed from pre-mRNA transcripts in multiple distinct segments, has been observed in a small subset of Drosophila melanogaster introns. However, detection of recursive splicing requires observation of splicing intermediates that are inherently unstable, making it difficult to study. Here we developed new computational approaches to identify recursively spliced introns and applied them, in combination with existing methods, to nascent RNA sequencing data from Drosophila S2 cells. These approaches identified hundreds of novel sites of recursive splicing, expanding the catalog of recursively spliced fly introns by 4-fold. A subset of recursive sites were validated by RT-PCR and sequencing. Recursive sites occur in most very long (> 40 kb) fly introns, including many genes involved in morphogenesis and development, and tend to occur near the midpoints of introns. Suggesting a possible function for recursive splicing, we observe that fly introns with recursive sites are spliced more accurately than comparably sized non-recursive introns.

Nascent RNA sequencing analysis provides insights into enhancer-mediated gene regulation

BackgroundEnhancers are distal cis-regulatory elements that control gene expression. Despite an increasing appreciation of the importance of enhancers in cellular function and disease, our knowledge of enhancer-regulated transcription is very limited. Nascent RNA sequencing technologies, such as global nuclear run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), not only provide a direct and reliable measurement of enhancer activity, but also allow for quantifying transcription of enhancers and target genes simultaneously, making these technologies extremely useful for exploring enhancer-mediated regulation.ResultsNascent RNA sequencing analysis (NRSA) provides a comprehensive view of enhancer-mediated gene regulation. NRSA not only outperforms existing methods for enhancer identification, but also enables annotation and quantification of active enhancers, and prediction of their target genes. Furthermore, NRSA identifies functionally important enhancers by integrating 1) nascent transcriptional changes in enhancers and their target genes and 2) binding profiles from regulator(s) of interest. Applied to wildtype and histone deacetylase 3 (Hdac3) knockout mouse livers, NRSA showed that HDAC3 regulates RNA polymerase recruitment through both proximal (promoter) and distal (enhancer) regulatory elements. Integrating ChIP-seq with PRO-seq data, NRSA prioritized enhancers based on their potential contribution to mediating HDAC3 regulation.ConclusionsNRSA will greatly facilitate the usage of nascent RNA sequencing techniques and accelerate the study of enhancer-mediated regulation.

Building a Science Gateway For Processing and Modeling Sequencing Data Via Apache Airavata.

The amount of DNA sequencing data has been exponentially growing during the past decade due to advances in sequencing technology. Processing and modeling large amounts of sequencing data can be computationally intractable for desktop computing platforms. High performance computing (HPC) resources offer advantages in terms of computing power, and can be a general solution to these problems. Using HPCs directly for computational needs requires skilled users who know their way around HPCs and acquiring such skills take time. Science gateways acts as the middle layer between users and HPCs, providing users with the resources to accomplish compute-intensive tasks without requiring specialized expertise. We developed a web-based computing platform for genome biologists by customizing the PHP Gateway for Airavata (PGA) framework that accesses publicly accessible HPC resources via Apache Airavata. This web computing platform takes advantage of the Extreme Science and Engineering Discovery Environment (XSEDE) which provides the resources for gateway development, including access to CPU, GPU, and storage resources. We used this platform to develop a gateway for the dREG algorithm, an online computing tool for finding functional regions in mammalian genomes using nascent RNA sequencing data. The dREG gateway provides its users a free, powerful and user-friendly GPU computing resource based on XSEDE, circumventing the need of specialized knowledge about installation, configuration, and execution on an HPC for biologists. The dREG gateway is available at: https://dREG.dnasequence.org/.