A transient transcriptional activation governs unpolarized-to-polarized morphogenesis during embryo implantation.

During implantation, embryos undergo an unpolarized-to-polarized transition to initiate postimplantation morphogenesis. However, the underlying molecular mechanism is unknown. Here, we identify a transient transcriptional activation governing embryonic morphogenesis and pluripotency transition during implantation. In naive pluripotent embryonic stem cells (ESCs), which represent preimplantation embryos, we find that the microprocessor component DGCR8 can recognize stem-loop structures within nascent mRNAs to sequester transcriptional coactivator FLII to suppress transcription directly. When mESCs exit from naive pluripotency, the ERK/RSK/P70S6K pathway rapidly activates, leading to FLII phosphorylation and disruption of DGCR8/FLII interaction. Phosphorylated FLII can bind to transcription factor JUN, activating cell migration-related genes to establish poised pluripotency akin to implanting embryos. Resequestration of FLII by DGCR8 drives poised ESCs into formative pluripotency. In summary, we identify a DGCR8/FLII/JUN-mediated transient transcriptional activation mechanism. Disruption of this mechanism inhibits naive-poised-formative pluripotency transition and the corresponding unpolarized-to-polarized transition during embryo implantation, which are conserved in mice and humans.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes.

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.

Butt-seq: a new method for facile profiling of transcription.

A wide range of sequencing methods has been developed to assess nascent RNA transcription and resolve the single-nucleotide position of RNA polymerase genome-wide. These techniques are often burdened with high input material requirements and lengthy protocols. We leveraged the template-switching properties of thermostable group II intron reverse transcriptase (TGIRT) and developed Butt-seq (bulk analysis of nascent transcript termini sequencing), which can produce libraries from purified nascent RNA in 6 h and from as few as 10,000 cells-an improvement of at least 10-fold over existing techniques. Butt-seq shows that inhibition of the superelongation complex (SEC) causes promoter-proximal pausing to move upstream in a fashion correlated with subnucleosomal fragments. To address transcriptional regulation in a tissue, Butt-seq was used to measure the circadian regulation of transcription from fly heads. All the results indicate that Butt-seq is a simple and powerful technique to analyze transcription at a high level of resolution.

Dynamic control of chromatin-associated m6A methylation regulates nascent RNA synthesis.

N6-methyladenosine (m6A) methylation is co-transcriptionally deposited on mRNA, but a possible role of m6A on transcription remains poorly understood. Here, we demonstrate that the METTL3/METTL14/WTAP m6A methyltransferase complex (MTC) is localized to many promoters and enhancers and deposits the m6A modification on nascent transcripts, including pre-mRNAs, promoter upstream transcripts (PROMPTs), and enhancer RNAs. PRO-seq analyses demonstrate that nascent RNAs originating from both promoters and enhancers are significantly decreased in the METTL3-depleted cells. Furthermore, genes targeted by the Integrator complex for premature termination are depleted of METTL3, suggesting a potential antagonistic relationship between METTL3 and Integrator. Consistently, we found the Integrator complex component INTS11 elevated at promoters and enhancers upon loss of MTC or nuclear m6A binders. Taken together, our findings suggest that MTC-mediated m6A modification protects nascent RNAs from Integrator-mediated termination and promotes productive transcription, thus unraveling an unexpected layer of gene regulation imposed by RNA m6A modification.

Co-transcriptional splicing regulates 3' end cleavage during mammalian erythropoiesis.

Pre-mRNA processing steps are tightly coordinated with transcription in many organisms. To determine how co-transcriptional splicing is integrated with transcription elongation and 3' end formation in mammalian cells, we performed long-read sequencing of individual nascent RNAs and precision run-on sequencing (PRO-seq) during mouse erythropoiesis. Splicing was not accompanied by transcriptional pausing and was detected when RNA polymerase II (Pol II) was within 75-300 nucleotides of 3' splice sites (3'SSs), often during transcription of the downstream exon. Interestingly, several hundred introns displayed abundant splicing intermediates, suggesting that splicing delays can take place between the two catalytic steps. Overall, splicing efficiencies were correlated among introns within the same transcript, and intron retention was associated with inefficient 3' end cleavage. Remarkably, a thalassemia patient-derived mutation introducing a cryptic 3'SS improved both splicing and 3' end cleavage of individual ß-globin transcripts, demonstrating functional coupling between the two co-transcriptional processes as a determinant of productive gene output.

Nascent RNA antagonizes the interaction of a set of regulatory proteins with chromatin.

A number of regulatory factors are recruited to chromatin by specialized RNAs. Whether RNA has a more general role in regulating the interaction of proteins with chromatin has not been determined. We used proteomics methods to measure the global impact of nascent RNA on chromatin in embryonic stem cells. Surprisingly, we found that nascent RNA primarily antagonized the interaction of chromatin modifiers and transcriptional regulators with chromatin. Transcriptional inhibition and RNA degradation induced recruitment of a set of transcriptional regulators, chromatin modifiers, nucleosome remodelers, and regulators of higher-order structure. RNA directly bound to factors, including BAF, NuRD, EHMT1, and INO80 and inhibited their interaction with nucleosomes. The transcriptional elongation factor P-TEFb directly bound pre-mRNA, and its recruitment to chromatin upon Pol II inhibition was regulated by the 7SK ribonucleoprotein complex. We postulate that by antagonizing the interaction of regulatory proteins with chromatin, nascent RNA links transcriptional output with chromatin composition.

Nascent RNA scaffolds contribute to chromosome territory architecture and counter chromatin compaction.

Nuclear chromosomes transcribe far more RNA than required to encode protein. Here we investigate whether non-coding RNA broadly contributes to cytological-scale chromosome territory architecture. We develop a procedure that depletes soluble proteins, chromatin, and most nuclear RNA from the nucleus but does not delocalize XIST, a known architectural RNA, from an insoluble chromosome "scaffold." RNA-seq analysis reveals that most RNA in the nuclear scaffold is repeat-rich, non-coding, and derived predominantly from introns of nascent transcripts. Insoluble, repeat-rich (C0T-1) RNA co-distributes with known scaffold proteins including scaffold attachment factor A (SAF-A), and distribution of these components inversely correlates with chromatin compaction in normal and experimentally manipulated nuclei. We further show that RNA is required for SAF-A to interact with chromatin and for enrichment of structurally embedded "scaffold attachment regions" prevalent in euchromatin. Collectively, the results indicate that long nascent transcripts contribute a dynamic structural role that promotes the open architecture of active chromosome territories.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing.

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

POINT technology illuminates the processing of polymerase-associated intact nascent transcripts.

Mammalian chromatin is the site of both RNA polymerase II (Pol II) transcription and coupled RNA processing. However, molecular details of such co-transcriptional mechanisms remain obscure, partly because of technical limitations in purifying authentic nascent transcripts. We present a new approach to characterize nascent RNA, called polymerase intact nascent transcript (POINT) technology. This three-pronged methodology maps nascent RNA 5' ends (POINT-5), establishes the kinetics of co-transcriptional splicing patterns (POINT-nano), and profiles whole transcription units (POINT-seq). In particular, we show by depletion of the nuclear exonuclease Xrn2 that this activity acts selectively on cleaved 5' P-RNA at polyadenylation sites. Furthermore, POINT-nano reveals that co-transcriptional splicing either occurs immediately after splice site transcription or is delayed until Pol II transcribes downstream sequences. Finally, we connect RNA cleavage and splicing with either premature or full-length transcript termination. We anticipate that POINT technology will afford full dissection of the complexity of co-transcriptional RNA processing.

The RNA tether model for human chromosomal translocation fragile zones.

One of the two chromosomal breakage events in recurring translocations in B cell neoplasms is often due to the recombination-activating gene complex (RAG complex) releasing DNA ends before end joining. The other break occurs in a fragile zone of 20-600 bp in a non-antigen receptor gene locus, with a more complex and intriguing set of mechanistic factors underlying such narrow fragile zones. These factors include activation-induced deaminase (AID), which acts only at regions of single-stranded DNA (ssDNA). Recent work leads to a model involving the tethering of AID to the nascent RNA as it emerges from the RNA polymerase. This mechanism may have relevance in class switch recombination (CSR) and somatic hypermutation (SHM), as well as broader relevance for other DNA enzymes.

Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates.

Transcription elongation rates influence RNA processing, but sequence-specific regulation is poorly understood. We addressed this in vivo, analyzing RNAPI in S. cerevisiae. Mapping RNAPI by Miller chromatin spreads or UV crosslinking revealed 5' enrichment and strikingly uneven local polymerase occupancy along the rDNA, indicating substantial variation in transcription speed. Two features of the nascent transcript correlated with RNAPI distribution: folding energy and GC content in the transcription bubble. In vitro experiments confirmed that strong RNA structures close to the polymerase promote forward translocation and limit backtracking, whereas high GC in the transcription bubble slows elongation. A mathematical model for RNAPI elongation confirmed the importance of nascent RNA folding in transcription. RNAPI from S. pombe was similarly sensitive to transcript folding, as were S. cerevisiae RNAPII and RNAPIII. For RNAPII, unstructured RNA, which favors slowed elongation, was associated with faster cotranscriptional splicing and proximal splice site use, indicating regulatory significance for transcript folding.

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.

Pre-mRNA splicing and its cotranscriptional connections.

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

The Escherichia coli ribB riboswitch senses flavin mononucleotide within a defined transcriptional window.

Riboswitches are metabolite-binding RNA regulators that modulate gene expression at the levels of transcription and translation. One of the hallmarks of riboswitch regulation is that they undergo structural changes upon metabolite binding. While a lot of effort has been put to characterize how the metabolite is recognized by the riboswitch, there is still relatively little information regarding how ligand sensing is performed within a transcriptional context. Here, we study the ligand-dependent cotranscriptional folding of the FMN-sensing ribB riboswitch of Escherichia coli. Using RNase H assays to study nascent ribB riboswitch transcripts, DNA probes targeting the P1 and sequestering stems indicate that FMN binding leads to the protection of these regions from RNase H cleavage, consistent with the riboswitch inhibiting translation initiation when bound to FMN. Our results show that ligand sensing is strongly affected by the position of elongating RNA polymerase, which is defining an FMN binding transcriptional window that is bordered in its 3' extremity by a transcriptional pause site. Also, using successively overlapping DNA probes targeting a subdomain of the riboswitch, our data suggest the presence of a previously unsuspected helical region involving the 3' strand of the P1 stem. Our results show that this helical region is conserved across bacterial species, thus suggesting that this predicted structure, the anti*-P1 stem, is involved in the FMN-free conformation of the ribB riboswitch. Overall, our study further demonstrates that intricate folding strategies may be used by riboswitches to perform metabolite sensing during the transcriptional process.

SAFB associates with nascent RNAs and can promote gene expression in mouse embryonic stem cells.

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Lampbrush chromosome studies in the post-genomic era.

Extraordinary extended lampbrush chromosomes with thousands of transcription loops are favorable objects in chromosome biology. Chromosomes become lampbrushy due to unusually high rate of transcription during oogenesis. However, until recently, the information on the spectrum of transcribed sequences as well as genomic context of individual chromomeres was mainly limited to tandemly repetitive elements. Here we briefly outline novel findings and future directions in lampbrush chromosome studies in the post-genomic era. We emphasize the fruitfulness of combining genome-wide approaches with microscopy imaging techniques using lampbrush chromosomes as a remarkable model object. We believe that new data on the spectrum of sequences transcribed on the lateral loops of lampbrush chromosomes and their structural organization push the boundaries in the discussion of their biological role. Also see the video abstract here: https://youtu.be/zexoHfzX9rM.

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.

Global Run-On sequencing to measure nascent transcription in Saccharomyces cerevisiae.

Global Run-On sequencing is a reliable and widely used approach for monitoring nascent transcription on a genomewide scale. The assay has been successfully used for studying global transcription in humans, plants, worms, flies, and fission yeast. Here we describe a GRO-seq protocol for studying transcription in budding yeast, Saccharomyces cerevisiae. Briefly, the technique involves permeabilization of actively growing yeast cells, allowing transcription to proceed in permeabilized cells in the presence of brominated UTP, affinity purification of bromo-UMP incorporated nascent transcripts followed by cDNA library construction, deep sequencing, and mapping against the reference genome. The approach maps the position of transcriptionally active RNA polymerase on a genomewide basis. In addition to identifying the complete set of transcriptionally active genes in a cell under a given set of conditions, the method can be used to determine elongation rate, termination defect and promoter directionality at the genomewide level. The approach is especially useful in identifying short-lived unstable transcripts that are rapidly degraded even before they leave the nucleus.

Human RNA Polymerase II Segregates from Genes and Nascent RNA and Transcribes in the Presence of DNA-Bound dCas9.

RNA polymerase II (Pol II) dysfunction is frequently implied in human disease. Understanding its functional mechanism is essential for designing innovative therapeutic strategies. To visualize its supra-molecular interactions with genes and nascent RNA, we generated a human cell line carrying ~335 consecutive copies of a recombinant ß-globin gene. Confocal microscopy showed that Pol II was not homogeneously concentrated around these identical gene copies. Moreover, Pol II signals partially overlapped with the genes and their nascent RNA, revealing extensive compartmentalization. Using a cell line carrying a single copy of the ß-globin gene, we also tested if the binding of catalytically dead CRISPR-associated system 9 (dCas9) to different gene regions affected Pol II transcriptional activity. We assessed Pol II localization and nascent RNA levels using chromatin immunoprecipitation and droplet digital reverse transcription PCR, respectively. Some enrichment of transcriptionally paused Pol II accumulated in the promoter region was detected in a strand-specific way of gRNA binding, and there was no decrease in nascent RNA levels. Pol II preserved its transcriptional activity in the presence of DNA-bound dCas9. Our findings contribute further insight into the complex mechanism of mRNA transcription in human cells.

Dynamics of transcriptional and post-transcriptional regulation.

Despite gene expression programs being notoriously complex, RNA abundance is usually assumed as a proxy for transcriptional activity. Recently developed approaches, able to disentangle transcriptional and post-transcriptional regulatory processes, have revealed a more complex scenario. It is now possible to work out how synthesis, processing and degradation kinetic rates collectively determine the abundance of each gene's RNA. It has become clear that the same transcriptional output can correspond to different combinations of the kinetic rates. This underscores the fact that markedly different modes of gene expression regulation exist, each with profound effects on a gene's ability to modulate its own expression. This review describes the development of the experimental and computational approaches, including RNA metabolic labeling and mathematical modeling, that have been disclosing the mechanisms underlying complex transcriptional programs. Current limitations and future perspectives in the field are also discussed.