Native Elongating Transcript Sequencing Research Articles

Size fractionated NET-Seq reveals a conserved architecture of transcription units around yeast genes.

Genomes from yeast to humans are subject to pervasive transcription. A single round of pervasive transcription is sufficient to alter local chromatin conformation, nucleosome dynamicsand gene expression, but is hard to distinguish from background signals. Size fractionated native elongating transcript sequencing (sfNET-Seq) was developed to precisely map nascent transcripts independent of expression levels. RNAPII-associated nascent transcripts are fractionation into different size ranges before library construction. When anchored to the transcription start sites (TSS) of annotated genes, the combined pattern of the output metagenes gives the expected reference pattern. Bioinformatic pattern matching to the reference pattern identified 9542 transcription units in Saccharomyces cerevisiae, of which 47% are coding and 53% are noncoding.In total, 3113 (33%) are unannotated noncoding transcription units. Anchoring all transcription units to the TSS or polyadenylation site (PAS) of annotated genes reveals distinctive architectures of linked pairs of divergent transcripts approximately 200nt apart. The Reb1 transcription factor is enriched 30nt downstream of the PAS only when an upstream (TSS -60nt with respect to PAS) noncoding transcription unit co-occurs with a downstream (TSS+150nt) coding transcription unit and acts to limit levels of upstream antisense transcripts. The potential for extensive transcriptional interference is evident from low abundance unannotated transcription units with variable TSS (median -240nt) initiating within a 500nt window upstream of, and transcribing over, the promoters of protein-coding genes. This study confirms a highly interleaved yeast genome with different types of transcription units altering the chromatin landscape in distinctive ways, with the potential to exert extensive regulatory control.

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

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

Rate of transcription elongation and sequence-specific pausing by RNA polymerase I directly influence rRNA processing

One of the first steps in ribosome biogenesis is transcription of the ribosomal DNA by RNA polymerase I (Pol I). Processing of the resultant rRNA begins cotranscriptionally, and perturbation of Pol I transcription elongation results in defective rRNA processing. Mechanistic insight regarding the link between transcription elongation and ribosome assembly is lacking because of limited invivo methods to assay Pol I transcription. Here, we use native elongating transcript sequencing (NET-Seq) with a strain of Saccharomyces cerevisiae containing a point mutation in Pol I, rpa190-F1205H, which results in impaired rRNA processing and ribosome assembly. We previously demonstrated that this mutation caused a mild reduction in the transcription elongation rate of Pol I invitro; however, transcription elongation by the mutant has not been characterized invivo. Here, our findings demonstrate that the mutant Pol I has an increased pause propensity during processive transcription elongation both invitro and invivo. NET-Seq reveals that rpa190-F1205H Pol I displays alternative pause site preferences invivo. Specifically, the mutant is sensitized to A/G residues in the RNA:DNA hybrid and at the last incorporated nucleotide position. Furthermore, both NET-Seq and EM analysis of Miller chromatin spreads reveal pileups of rpa190-F1205H Pol I throughout the ribosomal DNA, particularly at the 5' end of the 35S gene. This combination of invitro and invivo analyses of a Pol I mutant provides novel insights into Pol I elongation properties and indicates how these properties are crucial for efficient cotranscriptional rRNA processing and ribosome assembly.

Transcription elongation is finely tuned by dozens of regulatory factors.

Understanding the complex network that regulates transcription elongation requires the quantitative analysis of RNA polymerase II (Pol II) activity in a wide variety of regulatory environments. We performed native elongating transcript sequencing (NET-seq) in 41 strains of Saccharomyces cerevisiae lacking known elongation regulators, including RNA processing factors, transcription elongation factors, chromatin modifiers, and remodelers. We found that the opposing effects of these factors balance transcription elongation and antisense transcription. Different sets of factors tightly regulate Pol II progression across gene bodies so that Pol II density peaks at key points of RNA processing. These regulators control where Pol II pauses with each obscuring large numbers of potential pause sites that are primarily determined by DNA sequence and shape. Antisense transcription varies highly across the regulatory landscapes analyzed, but antisense transcription in itself does not affect sense transcription at the same locus. Our findings collectively show that a diverse array of factors regulate transcription elongation by precisely balancing Pol II activity.

Impaired Transcription Elongation by RNA Polymerase I Results in Defective Ribosomal RNA Processing

RNA polymerase I (Pol I) synthesizes the majority of the ribosomal RNA (rRNA) in a eukaryotic cell, which is then assembled into the mature ribosome with ribosomal proteins. Pol I transcribes a single gene, the 35S gene in yeast, containing three gene elements (18S, 5.8S, and 25S) and four spacer regions. Previous literature demonstrates that processing events occur both co‐ and post‐transcriptionally so that the spacer regions are removed and degraded, resulting in mature 18S, 5.8S, and 25S rRNAs. While it has been established that rRNA processing begins during transcription, the relationship between Pol I transcription elongation properties and processing events have not been well defined. The objective of this project was to determine whether transcription elongation kinetics have a regulatory role in rRNA processing, and we hypothesized that if transcription elongation was impaired, this would result in defects in processing. To test this hypothesis, we generated a yeast strain containing a point mutation in RPA190, the largest subunit of Pol I. We found that there was a modest decrease in transcription elongation rate and an increase in pausing in the mutant compared to wild‐type Pol I in vitro. To analyze elongation effects of the mutant Pol I in vivo, we used native elongating transcript sequencing, which probes for Pol I occupancy on the ribosomal DNA template at single‐nucleotide resolution. Our in vivo results corroborated the in vitro findings, indicating that this mutation introduced an increase in pausing across the template during transcription elongation. Furthermore, Northern blots and polysome profiling indicated that rRNA processing and ribosome assembly were impaired when this mutation was introduced. These three independent in vitroand in vivo experimental strategies converge to demonstrate that transcription elongation and Pol I pausing may play a regulatory role in downstream processing steps.

R-loops at microRNA encoding loci promote co-transcriptional processing of pri-miRNAs in plants

In most organisms, the maturation of nascent RNAs is coupled to transcription. Unlike in animals, the RNA polymerase II (RNAPII) transcribes microRNA genes (MIRNAs) as long and structurally variable pri-miRNAs in plants. Current evidence suggests that the miRNA biogenesis complex assembly initiates early during the transcription of pri-miRNAs in plants. However, it is unknown whether miRNA processing occurs co-transcriptionally. Here, we used native elongating transcript sequencing data and imaging techniques to demonstrate that plant miRNA biogenesis occurs coupled to transcription. We found that the entire biogenesis occurs co-transcriptionally for pri-miRNAs processed from the loop of the hairpin but requires a second nucleoplasmic step for those processed from the base. Furthermore, we found that co- and post-transcriptional miRNA processing mechanisms co-exist for most miRNAs in a dynamic balance. Notably, we discovered that R-loops, formed near the transcription start site region of MIRNAs, promote co-transcriptional pri-miRNA processing. Furthermore, our results suggest the neofunctionalization of co-transcriptionally processed miRNAs, boosting countless regulatory scenarios.

Defining the Influence of the A12.2 Subunit on Transcription Elongation and Termination by RNA Polymerase I In Vivo.

Saccharomyces cerevisiae has approximately 200 copies of the 35S rDNA gene, arranged tandemly on chromosome XII. This gene is transcribed by RNA polymerase I (Pol I) and the 35S rRNA transcript is processed to produce three of the four rRNAs required for ribosome biogenesis. An intergenic spacer (IGS) separates each copy of the 35S gene and contains the 5S rDNA gene, the origin of DNA replication, and the promoter for the adjacent 35S gene. Pol I is a 14-subunit enzyme responsible for the majority of rRNA synthesis, thereby sustaining normal cellular function and growth. The A12.2 subunit of Pol I plays a crucial role in cleavage, termination, and nucleotide addition during transcription. Deletion of this subunit causes alteration of nucleotide addition kinetics and read-through of transcription termination sites. To interrogate both of these phenomena, we performed native elongating transcript sequencing (NET-seq) with an rpa12Δ strain of S. cerevisiae and evaluated the resultant change in Pol I occupancy across the 35S gene and the IGS. Compared to wild-type (WT), we observed template sequence-specific changes in Pol I occupancy throughout the 35S gene. We also observed rpa12Δ Pol I occupancy downstream of both termination sites and throughout most of the IGS, including the 5S gene. Relative occupancy of rpa12Δ Pol I increased upstream of the promoter-proximal Reb1 binding site and dropped significantly downstream, implicating this site as a third terminator for Pol I transcription. Collectively, these high-resolution results indicate that the A12.2 subunit of Pol I plays an important role in transcription elongation and termination.

The small-molecule BMH-21 directly inhibits transcription elongation and DNA occupancy of RNA polymerase I in vivo and in vitro

Cancer cells are dependent upon an abundance of ribosomes to maintain rapid cell growth and proliferation. The rate-limiting step of ribosome biogenesis is ribosomal RNA (rRNA) synthesis by RNA polymerase I (Pol I). Therefore, a goal of the cancer therapeutic field is to develop and characterize Pol I inhibitors. Here, we elucidate the mechanism of Pol I inhibition by a first-in-class small-molecule BMH-21. To characterize the effects of BMH-21 on Pol I transcription, we leveraged high-resolution in vitro transcription assays and in vivo native elongating transcript sequencing (NET-seq). We find that Pol I transcription initiation, promoter escape, and elongation are all inhibited by BMH-21 in vitro. In particular, the transcription elongation phase is highly sensitive to BMH-21 treatment, as it causes a decrease in transcription elongation rate and an increase in paused Pols on the ribosomal DNA (rDNA) template. In vivo NET-seq experiments complement these findings by revealing a reduction in Pol I occupancy on the template and an increase in sequence-specific pausing upstream of G-rich rDNA sequences after BMH-21 treatment. Collectively, these data reveal the mechanism of action of BMH-21, which is a critical step forward in the development of this compound and its derivatives for clinical use.

Making sense of the natural antisense transcript puzzle.

In plants, thousands of genes are associated with antisense transcription, which often produces noncoding RNAs. Although widespread, sense-antisense pairs have been implicated in a limited variety of functions in plants and are often thought to form extensive dsRNA stretches triggering gene silencing. In this opinion, we show that evidence does not support gene silencing as a major role for antisense transcription. In fact, it is more likely that antisense transcripts play diverse functions in gene regulation. We propose a general framework for the initial functional dissection of antisense transcripts, suggesting testable hypotheses relying on an experiment-based decision tree. By moving beyond the gene silencing paradigm, we argue that a broad and diverse role for natural antisense transcription will emerge.

Coupling Between Production of Ribosomal RNA and Maturation: Just at the Beginning.

Ribosomal RNA (rRNA) production represents the most active transcription in the cell. Synthesis of the large rRNA precursors (35S/47S in yeast/human) is achieved by up to hundreds of RNA polymerase I (Pol I) enzymes simultaneously transcribing a single rRNA gene. In this review, we present recent advances in understanding the coupling between rRNA production and nascent rRNA folding. Mapping of the distribution of Pol I along ribosomal DNA at nucleotide resolution, using either native elongating transcript sequencing (NET-Seq) or crosslinking and analysis of cDNAs (CRAC), revealed frequent Pol I pausing, and CRAC results revealed a direct coupling between pausing and nascent RNA folding. High density of Pol I per gene imposes topological constraints that establish a defined pattern of polymerase distribution along the gene, with a persistent spacing between transcribing enzymes. RNA folding during transcription directly acts as an anti-pausing mechanism, implying that proper folding of the nascent rRNA favors elongation in vivo. Defects in co-transcriptional folding of rRNA are likely to induce Pol I pausing. We propose that premature termination of transcription, at defined positions, can control rRNA production in vivo.

Transcription and splicing dynamics during early Drosophila development.

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

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

Spt4 Promotes Pol I Processivity and Transcription Elongation.

RNA polymerases (Pols) I, II, and III collectively synthesize most of the RNA in a eukaryotic cell. Transcription by Pols I, II, and III is regulated by hundreds of trans-acting factors. One such protein, Spt4, has been previously identified as a transcription factor that influences both Pols I and II. Spt4 forms a complex with Spt5, described as the Spt4/5 complex (or DSIF in mammalian cells). This complex has been shown previously to directly interact with Pol I and potentially affect transcription elongation. The previous literature identified defects in transcription by Pol I when SPT4 was deleted, but the necessary tools to characterize the mechanism of this effect were not available at the time. Here, we use a technique called Native Elongating Transcript Sequencing (NET-seq) to probe for the global occupancy of Pol I in wild-type (WT) and spt4△ Saccharomyces cerevisiae (yeast) cells at single nucleotide resolution in vivo. Analysis of NET-seq data reveals that Spt4 promotes Pol I processivity and enhances transcription elongation through regions of the ribosomal DNA that are particularly G-rich. These data suggest that Spt4/5 may directly affect transcription elongation by Pol I in vivo.

A machine learning-based framework for modeling transcription elongation

RNA polymerase II (Pol II) generally pauses at certain positions along gene bodies, thereby interrupting the transcription elongation process, which is often coupled with various important biological functions, such as precursor mRNA splicing and gene expression regulation. Characterizing the transcriptional elongation dynamics can thus help us understand many essential biological processes in eukaryotic cells. However, experimentally measuring Pol II elongation rates is generally time and resource consuming. We developed PEPMAN (polymerase II elongation pausing modeling through attention-based deep neural network), a deep learning-based model that accurately predicts Pol II pausing sites based on the native elongating transcript sequencing (NET-seq) data. Through fully taking advantage of the attention mechanism, PEPMAN is able to decipher important sequence features underlying Pol II pausing. More importantly, we demonstrated that the analyses of the PEPMAN-predicted results around various types of alternative splicing sites can provide useful clues into understanding the cotranscriptional splicing events. In addition, associating the PEPMAN prediction results with different epigenetic features can help reveal important factors related to the transcription elongation process. All these results demonstrated that PEPMAN can provide a useful and effective tool for modeling transcription elongation and understanding the related biological factors from available high-throughput sequencing data.

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.

Illuminating Enhancer Transcription at Nucleotide Resolution with Native Elongating Transcript Sequencing (NET-Seq).

Enhancers are transcribed by RNA polymerase II (Pol II). In order to study the regulation of enhancer transcription and its function in target gene control, methods are required that track genome transcription with high precision in vivo. Here, we provide step-by-step guidance for performing native elongating transcript sequencing (NET-Seq) in mammalian cells. NET-Seq allows quantitative measurements of transcription genome-wide, including enhancer transcription, with single-nucleotide and DNA strand resolution. The approach consists of capturing and efficiently converting the 3'-ends of the nascent RNA into a sequencing library followed by next-generation sequencing and computational data analysis. The protocol includes quality control measurements to monitor the success of the main steps. Following this protocol, a NET-Seq library is obtained within 5days.

NusG controls transcription pausing and RNA polymerase translocation throughout the Bacillus subtilis genome

Transcription is punctuated by RNA polymerase (RNAP) pausing. These pauses provide time for diverse regulatory events that can modulate gene expression. Transcription elongation factors dramatically affect RNAP pausing in vitro, but the genome-wide role of such factors on pausing has not been examined. Using native elongating transcript sequencing followed by RNase digestion (RNET-seq), we analyzed RNAP pausing in Bacillus subtilis genome-wide and identified an extensive role of NusG in pausing. This universally conserved transcription elongation factor is known as Spt5 in archaeal and eukaryotic organisms. B. subtilis NusG shifts RNAP to the posttranslocation register and induces pausing at 1,600 sites containing a consensus TTNTTT motif in the nontemplate DNA strand within the paused transcription bubble. The TTNTTT motif is necessary but not sufficient for NusG-dependent pausing. Approximately one-fourth of these pause sites were localized to untranslated regions and could participate in posttranscription initiation control of gene expression as was previously shown for tlrB and the trpEDCFBA operon. Most of the remaining pause sites were identified in protein-coding sequences. NusG-dependent pausing was confirmed for all 10 pause sites that we tested in vitro. Putative pause hairpins were identified for 225 of the 342 strongest NusG-dependent pause sites, and some of these hairpins were shown to function in vitro. NusG-dependent pausing in the ribD riboswitch provides time for cotranscriptional binding of flavin mononucleotide, which decreases the concentration required for termination upstream of the ribD coding sequence. Our phylogenetic analysis implicates NusG-dependent pausing as a widespread mechanism in bacteria.

The yeast exoribonuclease Xrn1 and associated factors modulate RNA polymerase II processivity in 5‘ and 3‘ gene regions

mRNA levels are determined by the balance between mRNA synthesis and decay. Protein factors that mediate both processes, including the 5'-3' exonuclease Xrn1, are responsible for a cross-talk between the two processes that buffers steady-state mRNA levels. However, the roles of these proteins in transcription remain elusive and controversial. Applying native elongating transcript sequencing (NET-seq) to yeast cells, we show that Xrn1 functions mainly as a transcriptional activator and that its disruption manifests as a reduction of RNA polymerase II (Pol II) occupancy downstream of transcription start sites. By combining our sequencing data and mathematical modeling of transcription, we found that Xrn1 modulates transcription initiation and elongation of its target genes. Furthermore, Pol II occupancy markedly increased near cleavage and polyadenylation sites in xrn1Δ cells, whereas its activity decreased, a characteristic feature of backtracked Pol II. We also provide indirect evidence that Xrn1 is involved in transcription termination downstream of polyadenylation sites. We noted that two additional decay factors, Dhh1 and Lsm1, seem to function similarly to Xrn1 in transcription, perhaps as a complex, and that the decay factors Ccr4 and Rpb4 also perturb transcription in other ways. Interestingly, the decay factors could differentiate between SAGA- and TFIID-dominated promoters. These two classes of genes responded differently to XRN1 deletion in mRNA synthesis and were differentially regulated by mRNA decay pathways, raising the possibility that one distinction between these two gene classes lies in the mechanisms that balance mRNA synthesis with mRNA decay.

Downstream sequence-dependent RNA cleavage and pausing by RNA polymerase I.

The sequence of the DNA template has long been thought to influence the rate of transcription by DNA-dependent RNA polymerases, but the influence of DNA sequence on transcription elongation properties of eukaryotic RNA polymerase I (Pol I) from Saccharomyces cerevisiae has not been defined. In this study, we observe changes in dinucleotide production, transcription elongation complex stability, and Pol I pausing in vitro in response to downstream DNA. In vitro studies demonstrate that AT-rich downstream DNA enhances pausing by Pol I and inhibits Pol I nucleolytic cleavage activity. Analysis of Pol I native elongating transcript sequencing data in Saccharomyces cerevisiae suggests that these downstream sequence elements influence Pol I in vivo Native elongating transcript sequencing studies reveal that Pol I occupancy increases as downstream AT content increases and decreases as downstream GC content increases. Collectively, these data demonstrate that the downstream DNA sequence directly impacts the kinetics of transcription elongation prior to the sequence entering the active site of Pol I both in vivo and in vitro.

Analysis of Mammalian Native Elongating Transcript sequencing (mNET-seq) high-throughput data.

Mammalian Native Elongating Transcript sequencing (mNET-seq) is a recently developed technique that generates genome-wide profiles of nascent transcripts associated with RNA polymerase II (Pol II) elongation complexes. The ternary transcription complexes formed by Pol II, DNA template and nascent RNA are first isolated, without crosslinking, by immunoprecipitation with antibodies that specifically recognize the different phosphorylation states of the polymerase large subunit C-terminal domain (CTD). The coordinate of the 3' end of the RNA in the complexes is then identified by high-throughput sequencing. The main advantage of mNET-seq is that it provides global, bidirectional maps of Pol II CTD phosphorylation-specific nascent transcripts and coupled RNA processing at single nucleotide resolution. Here we describe the general pipeline to prepare and analyse high-throughput data from mNET-seq experiments.