Kinetic competition during the transcription cycle results in stochastic RNA processing

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

CCAT1 lncRNA is chromatin-retained and post-transcriptionally spliced

Super-enhancers are unique gene expression regulators widely involved in cancer development. Spread over large DNA segments, they tend to be found next to oncogenes. The super-enhancer c-MYC locus forms long-range chromatin looping with nearby genes, which brings the enhancer and the genes into proximity, to promote gene activation. The colon cancer-associated transcript 1 (CCAT1) gene, which is part of the MYC locus, transcribes a lncRNA that is overexpressed in colon cancer cells through activation by MYC. Comparing different types of cancer cell lines using RNA fluorescence in situ hybridization (RNA FISH), we detected very prominent CCAT1 expression in HeLa cells, observed as several large CCAT1 nuclear foci. We found that dozens of CCAT1 transcripts accumulate on the gene locus, in addition to active transcription occurring from the gene. The accumulating transcripts are released from the chromatin during cell division. Examination of CCAT1 lncRNA expression patterns on the single-RNA level showed that unspliced CCAT1 transcripts are released from the gene into the nucleoplasm. Most of these unspliced transcripts were observed in proximity to the active gene but were not associated with nuclear speckles in which unspliced RNAs usually accumulate. At larger distances from the gene, the CCAT1 transcripts appeared spliced, implying that most CCAT1 transcripts undergo post-transcriptional splicing in the zone of the active gene. Finally, we show that unspliced CCAT1 transcripts can be detected in the cytoplasm during splicing inhibition, which suggests that there are several CCAT1 variants, spliced and unspliced, that the cell can recognize as suitable for export.

Analysis of a detailed multi-stage model of stochastic gene expression using queueing theory and model reduction

We introduce a biologically detailed, stochastic model of gene expression describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation of mRNA into protein. The processes in sub-cellular compartments are described by an arbitrary number of processing stages, thus accounting for a significantly finer molecular description of gene expression than conventional models such as the telegraph, two-stage and three-stage models of gene expression. We use two distinct tools, queueing theory and model reduction using the slow-scale linear-noise approximation, to derive exact or approximate analytic expressions for the moments or distributions of nuclear mRNA, cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their Fano factors in steady-state conditions. We use these to study the phase diagram of the stochastic model; in particular we derive parametric conditions determining three types of transitions in the properties of mRNA fluctuations: from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus to high noise in the cytoplasm, and from a monotonic increase to a monotonic decrease of the Fano factor with the number of processing stages. In contrast, protein fluctuations are always super-Poissonian and show weak dependence on the number of mRNA processing stages. Our results delineate the region of parameter space where conventional models give qualitatively incorrect results and provide insight into how the number of processing stages, e.g. the number of rate-limiting steps in initiation, splicing and mRNA degradation, shape stochastic gene expression by modulation of molecular memory.

Protocol to quantify the rate of RNA polymerase II elongation with MS2/MCP- and PP7/PCP-based live-imaging systems in early Drosophila embryos

The MS2-PP7 two-color live-imaging system provides insights into the spatiotemporal dynamics of nascent transcripts at tagged loci. Here, we present a protocol to quantitatively measure the rate of RNA polymerase II elongation for each actively transcribing nucleus in living Drosophila embryos. The elongation rate is calculated by measuring the effective distance and the time elapsed between MS2 and PP7 trajectories. We describe steps for preparing embryo samples, performing live imaging, and measuring the elongation rate. For complete details on the use and execution of this protocol, please refer to Keller et al.1.

Co-transcriptional splicing facilitates transcription of gigantic genes

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Post-transcriptional splicing can occur in a slow-moving zone around the gene

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

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.

Measuring Transcription Dynamics of Individual Genes Inside Living Cells

Transcription is a highly dynamic process, which, for many genes, occurs in stochastic bursts. Studying what regulates these stochastic bursts requires visualization and quantification of transcription dynamics in single living cells. Such measurements of bursting can be accomplished by labeling nascent transcripts of single genes fluorescently with the MS2 and PP7 RNA labeling techniques. Live-cell single-molecule microscopy of transcription in real time allows for the extraction of transcriptional bursting kinetics inside single cells. This chapter describes how to set up the MS2 or PP7 RNA labeling system of endogenous genes in both budding yeast (Saccharomyces cerevisiae) and mammalian cells (mouse embryonic stem cells). We include how to genetically engineer the cells with the MS2 and PP7 system, describe how to perform the live-microscopy experiments and discuss how to extract transcriptional bursting parameters of the genes of interest.

Stress-induced nuclear speckle reorganization is linked to activation of immediate early gene splicing

Current models posit that nuclear speckles (NSs) serve as reservoirs of splicing factors and facilitate posttranscriptional mRNA processing. Here, we discovered that ribotoxic stress induces a profound reorganization of NSs with enhanced recruitment of factors required for splice-site recognition, including the RNA-binding protein TIAR, U1 snRNP proteins and U2-associated factor 65, as well as serine 2 phosphorylated RNA polymerase II. NS reorganization relies on the stress-activated p38 mitogen-activated protein kinase (MAPK) pathway and coincides with splicing activation of both pre-existing and newly synthesized pre-mRNAs. In particular, ribotoxic stress causes targeted excision of retained introns from pre-mRNAs of immediate early genes (IEGs), whose transcription is induced during the stress response. Importantly, enhanced splicing of the IEGs ZFP36 and FOS is accompanied by relocalization of the corresponding nuclear mRNA foci to NSs. Our study reveals NSs as a dynamic compartment that is remodeled under stress conditions, whereby NSs appear to become sites of IEG transcription and efficient cotranscriptional splicing.

Regulation of the dynamic RNA Pol II elongation rate in Drosophila embryos

An increasing number of studies have shown the key role that RNA polymerase II (RNA Pol II) elongation plays in gene regulation. We systematically examine how various enhancers, promoters, and gene body composition influence the RNA Pol II elongation rate through a single-cell-resolution live imaging assay. By using reporter constructs containing 5' MS2 and 3' PP7 repeating stem loops, we quantify the rate of RNA Pol II elongation in live Drosophila embryos. We find that promoters and exonic gene lengths have no effect on elongation rate, while enhancers and the presence of long introns may significantly change how quickly RNA Pol II moves across a gene. Furthermore, we observe in multiple constructs that the RNA Pol II elongation rate accelerates after the transcriptional onset of nuclear cycle 14 in Drosophila embryos. Our study provides a single-cell view of various mechanisms that affect the dynamic RNA Pol II elongation rate, ultimately affecting the rate of mRNA production.

Stochastic modeling of the mRNA life process: A generalized master equation

The mRNA life cycle is a complex biochemical process, involving transcription initiation, elongation, termination, splicing, and degradation. Each of these molecular events is multistep and can create a memory. The effect of this molecular memory on gene expression is not clear, although there are many related yet scattered experimental reports. To address this important issue, we develop a general theoretical framework formulated as a master equation in the sense of queue theory, which can reduce to multiple previously studied gene models in limiting cases. This framework allows us to interpret experimental observations, extract kinetic parameters from experimental data, and identify how the mRNA kinetics vary under regulatory influences. Notably, it allows us to evaluate the influences of elongation processes on mature RNA distribution; e.g., we find that the non-exponential elongation time can induce the bimodal mRNA expression and there is an optimal elongation noise intensity such that the mature RNA noise achieves the lowest level. In a word, our framework can not only provide insight into complex mRNA life processes but also bridge a dialogue between theoretical studies and experimental data.

Following the Birth, Life, and Death of mRNAs in Single Cells

Recent advances in single-molecule imaging of mRNAs in fixed and living cells have enabled the lives of mRNAs to be studied with unprecedented spatial and temporal detail. These approaches have moved beyond simply being able to observe specific events and have begun to allow an understanding of how regulation is coupled between steps in the mRNA life cycle. Additionally, these methodologies are now being applied in multicellular systems and animals to provide more nuanced insights into the physiological regulation of RNA metabolism.

Full-length transcript alterations in human bronchial epithelial cells with U2AF1 S34F mutations

U2AF1 is one of the most recurrently mutated splicing factors in lung adenocarcinoma and has been shown to cause transcriptome-wide pre-mRNA splicing alterations; however, the full-length altered mRNA isoforms associated with the mutation are largely unknown. To better understand the impact U2AF1 has on full-length isoform fate and function, we conducted high-throughput long-read cDNA sequencing from isogenic human bronchial epithelial cells with and without a U2AF1 S34F mutation. We identified 49,366 multi-exon transcript isoforms, more than half of which did not match GENCODE or short-read-assembled isoforms. We found 198 transcript isoforms with significant expression and usage changes relative to WT, only 68% of which were assembled by short reads. Expression of isoforms from immune-related genes is largely down-regulated in mutant cells and without observed splicing changes. Finally, we reveal that isoforms likely targeted by nonsense-mediated decay are down-regulated in U2AF1 S34F cells, suggesting that isoform changes may alter the translational output of those affected genes. Altogether, our work provides a resource of full-length isoforms associated with U2AF1 S34F in lung cells.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the spatiotemporal arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find a close link between chromatin mobility and transcriptional status: active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

Integrating analog and digital modes of gene expression at Arabidopsis FLC

Quantitative gene regulation at the cell population level can be achieved by two fundamentally different modes of regulation at individual gene copies. A 'digital' mode involves binary ON/OFF expression states, with population-level variation arising from the proportion of gene copies in each state, while an 'analog' mode involves graded expression levels at each gene copy. At the Arabidopsis floral repressor FLOWERING LOCUS C (FLC), 'digital' Polycomb silencing is known to facilitate quantitative epigenetic memory in response to cold. However, whether FLC regulation before cold involves analog or digital modes is unknown. Using quantitative fluorescent imaging of FLC mRNA and protein, together with mathematical modeling, we find that FLC expression before cold is regulated by both analog and digital modes. We observe a temporal separation between the two modes, with analog preceding digital. The analog mode can maintain intermediate expression levels at individual FLC gene copies, before subsequent digital silencing, consistent with the copies switching OFF stochastically and heritably without cold. This switch leads to a slow reduction in FLC expression at the cell population level. These data present a new paradigm for gradual repression, elucidating how analog transcriptional and digital epigenetic memory pathways can be integrated.

Transcription factor clusters enable target search but do not contribute to target gene activation

Many transcription factors (TFs) localize in nuclear clusters of locally increased concentrations, but how TF clustering is regulated and how it influences gene expression is not well understood. Here, we use quantitative microscopy in living cells to study the regulation and function of clustering of the budding yeast TF Gal4 in its endogenous context. Our results show that Gal4 forms clusters that overlap with the GAL loci. Cluster number, density and size are regulated in different growth conditions by the Gal4-inhibitor Gal80 and Gal4 concentration. Gal4 truncation mutants reveal that Gal4 clustering is facilitated by, but does not completely depend on DNA binding and intrinsically disordered regions. Moreover, we discover that clustering acts as a double-edged sword: self-interactions aid TF recruitment to target genes, but recruited Gal4 molecules that are not DNA-bound do not contribute to, and may even inhibit, transcription activation. We propose that cells need to balance the different effects of TF clustering on target search and transcription activation to facilitate proper gene expression.

Bayesian analysis dissects kinetic modulation during non-stationary gene expression

Labelling of nascent stem loops with fluorescent proteins has fostered the visualization of transcription in living cells. Quantitative analysis of recorded fluorescence traces can shed light on kinetic transcription parameters and regulatory mechanisms. However, existing methods typically focus on steady state dynamics. Here, we combine a stochastic process transcription model with a hierarchical Bayesian method to infer global as well locally shared parameters for groups of cells and recover unobserved quantities such as initiation times and polymerase loading of the gene. We apply our approach to the cyclic response of the yeast CUP1 locus to heavy metal stress. Within the previously described slow cycle of transcriptional activity on the scale of minutes, we discover fast time-modulated bursting on the scale of seconds. Model comparison suggests that slow oscillations of transcriptional output are regulated by the amplitude of the bursts. Several polymerases may initiate during a burst.

DNA supercoiling restricts the transcriptional bursting of neighboring eukaryotic genes

DNA supercoiling has emerged as a major contributor to gene regulation in bacteria, but how DNA supercoiling impacts transcription dynamics in eukaryotes is unclear. Here, using single-molecule dual-color nascent transcription imaging in budding yeast, we show that transcriptional bursting of divergent and tandem GAL genes is coupled. Temporal coupling of neighboring genes requires rapid release of DNA supercoils by topoisomerases. When DNA supercoils accumulate, transcription of one gene inhibits transcription at its adjacent genes. Transcription inhibition of the GAL genes results from destabilized binding of the transcription factor Gal4. Moreover, wild-type yeast minimizes supercoiling-mediated inhibition by maintaining sufficient levels of topoisomerases. Overall, we discover fundamental differences in transcriptional control by DNA supercoiling between bacteria and yeast and show that rapid supercoiling release in eukaryotes ensures proper gene expression of neighboring genes.

Dynamic epistasis analysis reveals how chromatin remodeling regulates transcriptional bursting

Transcriptional bursting has been linked to the stochastic positioning of nucleosomes. However, how bursting is regulated by the remodeling of promoter nucleosomes is unknown. Here, we use single-molecule live-cell imaging of GAL10 transcription in Saccharomyces cerevisiae to measure how bursting changes upon combined perturbations of chromatin remodelers, the transcription factor Gal4 and preinitiation complex components. Using dynamic epistasis analysis, we reveal how the remodeling of different nucleosomes regulates transcriptional bursting parameters. At the nucleosome covering the Gal4 binding sites, RSC and Gal4 binding synergistically facilitate each burst. Conversely, nucleosome remodeling at the TATA box controls only the first burst upon galactose induction. At canonical TATA boxes, the nucleosomes are displaced by TBP binding to allow for transcription activation even in the absence of remodelers. Overall, our results reveal how promoter nucleosome remodeling together with Gal4 and preinitiation complex binding regulates transcriptional bursting.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the 4D arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find that alterations in chromatin mobility, not promoter-enhancer distance, is more informative about transcriptional status. Active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

The uprise of RNA biology in neuroendocrine neoplasms: altered splicing and RNA species unveil translational opportunities

Neuroendocrine neoplasms (NENs) comprise a highly heterogeneous group of tumors arising from the diffuse neuroendocrine system. NENs mainly originate in gastrointestinal, pancreatic, and pulmonary tissues, and despite being rare, show rising incidence. The molecular mechanisms underlying NEN development are still poorly understood, although recent studies are unveiling their genomic, epigenomic and transcriptomic landscapes. RNA was originally considered as an intermediary between DNA and protein. Today, compelling evidence underscores the regulatory relevance of RNA processing, while new RNA molecules emerge with key functional roles in core cell processes. Indeed, correct functioning of the interrelated complementary processes comprising RNA biology, its processing, transport, and surveillance, is essential to ensure adequate cell homeostasis, and its misfunction is related to cancer at multiple levels. This review is focused on the dysregulation of RNA biology in NENs. In particular, we survey alterations in the splicing process and available information implicating the main RNA species and processes in NENs pathology, including their role as biomarkers, and their functionality and targetability. Understanding how NENs precisely (mis)behave requires a profound knowledge at every layer of their heterogeneity, to help improve NEN management. RNA biology provides a wide spectrum of previously unexplored processes and molecules that open new avenues for NEN detection, classification and treatment. The current molecular biology era is rapidly evolving to facilitate a detailed comprehension of cancer biology and is enabling the arrival of personalized, predictive and precision medicine to rare tumors like NENs.

Inferring transcriptional bursting kinetics from single-cell snapshot data using a generalized telegraph model

Gene expression has inherent stochasticity resulting from transcription's burst manners. Single-cell snapshot data can be exploited to rigorously infer transcriptional burst kinetics, using mathematical models as blueprints. The classical telegraph model (CTM) has been widely used to explain transcriptional bursting with Markovian assumptions. However, growing evidence suggests that the gene-state dwell times are generally non-exponential, as gene-state switching is a multi-step process in organisms. Therefore, interpretable non-Markovian mathematical models and efficient statistical inference methods are urgently required in investigating transcriptional burst kinetics. We develop an interpretable and tractable model, the generalized telegraph model (GTM), to characterize transcriptional bursting that allows arbitrary dwell-time distributions, rather than exponential distributions, to be incorporated into the ON and OFF switching process. Based on the GTM, we propose an inference method for transcriptional bursting kinetics using an approximate Bayesian computation framework. This method demonstrates an efficient and scalable estimation of burst frequency and burst size on synthetic data. Further, the application of inference to genome-wide data from mouse embryonic fibroblasts reveals that GTM would estimate lower burst frequency and higher burst size than those estimated by CTM. In conclusion, the GTM and the corresponding inference method are effective tools to infer dynamic transcriptional bursting from static single-cell snapshot data.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

Position-dependent effects of RNA-binding proteins in the context of co-transcriptional splicing

Alternative splicing is an important step in eukaryotic mRNA pre-processing which increases the complexity of gene expression programs, but is frequently altered in disease. Previous work on the regulation of alternative splicing has demonstrated that splicing is controlled by RNA-binding proteins (RBPs) and by epigenetic DNA/histone modifications which affect splicing by changing the speed of polymerase-mediated pre-mRNA transcription. The interplay of these different layers of splicing regulation is poorly understood. In this paper, we derived mathematical models describing how splicing decisions in a three-exon gene are made by combinatorial spliceosome binding to splice sites during ongoing transcription. We additionally take into account the effect of a regulatory RBP and find that the RBP binding position within the sequence is a key determinant of how RNA polymerase velocity affects splicing. Based on these results, we explain paradoxical observations in the experimental literature and further derive rules explaining why the same RBP can act as inhibitor or activator of cassette exon inclusion depending on its binding position. Finally, we derive a stochastic description of co-transcriptional splicing regulation at the single-cell level and show that splicing outcomes show little noise and follow a binomial distribution despite complex regulation by a multitude of factors. Taken together, our simulations demonstrate the robustness of splicing outcomes and reveal that quantitative insights into kinetic competition of co-transcriptional events are required to fully understand this important mechanism of gene expression diversity.

Extracting, quantifying, and comparing dynamical and biomechanical properties of living matter through single particle tracking

A panoply of new tools for tracking single particles and molecules has led to an explosion of experimental data, leading to novel insights into physical properties of living matter governing cellular development and function, health and disease. In this Perspective, we present tools to investigate the dynamics and mechanics of living systems from the molecular to cellular scale via single-particle techniques. In particular, we focus on methods to measure, interpret, and analyse complex data sets that are associated with forces, materials properties, transport, and emergent organisation phenomena within biological and soft-matter systems. Current approaches, challenges, and existing solutions in the associated fields are outlined in order to support the growing community of researchers at the interface of physics and the life sciences. Each section focuses not only on the general physical principles and the potential for understanding living matter, but also on details of practical data extraction and analysis, discussing limitations, interpretation, and comparison across different experimental realisations and theoretical frameworks. Particularly relevant results are introduced as examples. While this Perspective describes living matter from a physical perspective, highlighting experimental and theoretical physics techniques relevant for such systems, it is also meant to serve as a solid starting point for researchers in the life sciences interested in the implementation of biophysical methods.

Genetically encodable tagging and sensing systems for fluorescent RNA imaging

Live cell imaging of RNAs is crucial to interrogate their fundamental roles in various biological processes. The highly spatiotemporal dynamic nature of RNA abundance and localization has presented great challenges for RNA imaging. Genetically encodable tagging and sensing (GETS) systems that can be continuously produced in living systems have afforded promising tools for imaging and sensing RNA dynamics in live cells. Here we review the recent advances of GETS systems that have been developed for RNA tagging and sensing in live cells. We first describe the various GETS systems using MS2-bacteriophage-MS2 coat protein, pumilio homology domain and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9/13 for RNA labeling and tracking. The progresses of GETS systems for fluorogenic labeling and/or sensing RNAs by engineering light-up RNA aptamers, CRISPR-Cas9 systems and RNA aptamer stabilized fluorogenic proteins are then elaborated. The challenges and future perspectives in this field are finally discussed. With the continuing development, GETS systems will afford powerful tools to elucidate RNA biology in living systems.

Visualization, Quantification, and Modeling of Endogenous RNA Polymerase II Phosphorylation at a Single-copy Gene in Living Cells

In eukaryotic cells, RNA Polymerase II (RNAP2) is the enzyme in charge of transcribing mRNA from DNA. RNAP2 possesses an extended carboxy-terminal domain (CTD) that gets dynamically phosphorylated as RNAP2 progresses through the transcription cycle, therefore regulating each step of transcription from recruitment to termination. Although RNAP2 residue-specific phosphorylation has been characterized in fixed cells by immunoprecipitation-based assays, or in live cells by using tandem gene arrays, these assays can mask heterogeneity and limit temporal and spatial resolution. Our protocol employs multi-colored complementary fluorescent antibody-based (Fab) probes to specifically detect the CTD of the RNAP2 (CTD-RNAP2), and its phosphorylated form at the serine 5 residue (Ser5ph-RNAP2) at a single-copy HIV-1 reporter gene. Together with high-resolution fluorescence microscopy, single-molecule tracking analysis, and rigorous computational modeling, our system allows us to visualize, quantify, and predict endogenous RNAP2 phosphorylation dynamics and mRNA synthesis at a single-copy gene, in living cells, and throughout the transcription cycle. Graphical abstract: Schematic of the steps for visualizing, quantifying, and predicting RNAP2 phosphorylation at a single-copy gene.

Multi-color RNA imaging with CRISPR-Cas13b systems in living cells

Visualizing RNA dynamics is important for understanding RNA function. Catalytically dead (d) CRISPR-Cas13 systems have been established to image and track RNAs in living cells, but efficient dCas13 for RNA imaging is still limited. Here, we analyzed metagenomic and bacterial genomic databases to comprehensively screen Cas13 homologies for their RNA labeling capabilities in living mammalian cells. Among eight previously unreported dCas13 proteins that can be used for RNA labeling, dHgm4Cas13b and dMisCas13b displayed comparable, if not higher, efficiencies to the best-known ones when targeting endogenous MUC4 and NEAT1_2 by single guide (g) RNAs. Further examination of the labeling robustness of different dCas13 systems using the GCN4 repeats revealed that a minimum of 12 GCN4 repeats was required for dHgm4Cas13b and dMisCas13b imaging at the single RNA molecule level, while >24 GCN4 repeats were required for reported dLwaCas13a, dRfxCas13d and dPguCas13b. Importantly, by silencing pre-crRNA processing activity of dMisCas13b (ddMisCas13b) and further incorporating RNA aptamers including PP7, MS2, Pepper or BoxB to individual gRNAs, a CRISPRpalette system was developed to successfully achieve multi-color RNA visualization in living cells.

Transcriptional regulation and chromatin architecture maintenance are decoupled functions at the Sox2 locus

How distal regulatory elements control gene transcription and chromatin topology is not clearly defined, yet these processes are closely linked in lineage specification during development. Through allele-specific genome editing and chromatin interaction analyses of the Sox2 locus in mouse embryonic stem cells, we found a striking disconnection between transcriptional control and chromatin architecture. We traced nearly all Sox2 transcriptional activation to a small number of key transcription factor binding sites, whose deletions have no effect on promoter-enhancer interaction frequencies or topological domain organization. Local chromatin architecture maintenance, including at the topologically associating domain (TAD) boundary downstream from the Sox2 enhancer, is widely distributed over multiple transcription factor-bound regions and maintained in a CTCF-independent manner. Furthermore, partial disruption of promoter-enhancer interactions by ectopic chromatin loop formation has no effect on Sox2 transcription. These findings indicate that many transcription factors are involved in modulating chromatin architecture independently of CTCF.

Illuminating RNA biology through imaging

RNA processing plays a central role in accurately transmitting genetic information into functional RNA and protein regulators. To fully appreciate the RNA life-cycle, tools to observe RNA with high spatial and temporal resolution are critical. Here we review recent advances in RNA imaging and highlight how they will propel the field of RNA biology. We discuss current trends in RNA imaging and their potential to elucidate unanswered questions in RNA biology.

Pre-mRNA splicing factor U2AF2 recognizes distinct conformations of nucleotide variants at the center of the pre-mRNA splice site signal

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies.

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

MYC amplifies gene expression through global changes in transcription factor dynamics

The MYC oncogene has been studied for decades, yet there is still intense debate over how this transcription factor controls gene expression. Here, we seek to answer these questions with an in vivo readout of discrete events of gene expression in single cells. We engineered an optogenetic variant of MYC (Pi-MYC) and combined this tool with single-molecule RNA and protein imaging techniques to investigate the role of MYC in modulating transcriptional bursting and transcription factor binding dynamics in human cells. We find that the immediate consequence of MYC overexpression is an increase in the duration rather than in the frequency of bursts, a functional role that is different from the majority of human transcription factors. We further propose that the mechanism by which MYC exerts global effects on the active period of genes is by altering the binding dynamics of transcription factors involved in RNA polymerase II complex assembly and productive elongation.

Next generation single-molecule techniques: Imaging, labeling, and manipulation in vitro and in cellulo

Owing to their unique abilities to manipulate, label, and image individual molecules in vitro and in cellulo, single-molecule techniques provide previously unattainable access to elementary biological processes. In imaging, single-molecule fluorescence resonance energy transfer (smFRET) and protein-induced fluorescence enhancement in vitro can report on conformational changes and molecular interactions, single-molecule pull-down (SiMPull) can capture and analyze the composition and function of native protein complexes, and single-molecule tracking (SMT) in live cells reveals cellular structures and dynamics. In labeling, the abilities to specifically label genomic loci, mRNA, and nascent polypeptides in cells have uncovered chromosome organization and dynamics, transcription and translation dynamics, and gene expression regulation. In manipulation, optical tweezers, integration of single-molecule fluorescence with force measurements, and single-molecule force probes in live cells have transformed our mechanistic understanding of diverse biological processes, ranging from protein folding, nucleic acids-protein interactions to cell surface receptor function.

Single-Molecule RNA Imaging Using Mango II Arrays

In recent years, fluorogenic RNA aptamers, such as Spinach, Broccoli, Corn, Mango, Coral, and Pepper have gathered traction as an efficient alternative labeling strategy for background-free imaging of cellular RNAs. However, their application has been somewhat limited by relatively inefficient folding and fluorescent stability. With the recent advent of novel RNA-Mango variants which are improved in both fluorescence intensity and folding stability in tandem arrays, it is now possible to image RNAs with single-molecule sensitivity. Here we discuss the protocol for imaging Mango II tagged RNAs in both fixed and live cells.

Single-Molecule Imaging of mRNA Interactions with Stress Granules

Single-molecule imaging in living cells enables the investigation of molecular dynamics and interactions underlying the physiology of a cell. We recently developed a method to visualize translation events at single-mRNA resolution in living cells. Here we describe how we apply this method to visualize mRNA interactions with stress granules in the context of translational activity during cell stress.

Understanding lncRNA-protein assemblies with imaging and single-molecule approaches

Long non-coding RNAs (lncRNAs) associate with RNA-binding proteins (RBPs) to form lncRNA-protein complexes that act in a wide range of biological processes. Understanding the molecular mechanism of how a lncRNA-protein complex is assembled and regulated is key for their cellular functions. In this mini-review, we outline molecular methods used to identify lncRNA-protein interactions from large-scale to individual levels using bulk cells as well as those recently developed imaging and single-molecule approaches that are capable of visualizing RNA-protein assemblies in single cells and in real-time. Focusing on the latter group of approaches, we discuss their applications and limitations, which nevertheless have enabled quantification and comprehensive dissection of RNA-protein interactions possible.

Nuclear envelope mechanobiology: linking the nuclear structure and function

The nucleus, central to cellular activity, relies on both direct mechanical input as well as its molecular transducers to sense external stimuli and respond by regulating intra-nuclear chromatin organization that determines cell function and fate. In mesenchymal stem cells of musculoskeletal tissues, changes in nuclear structures are emerging as a key modulator of their differentiation and proliferation programs. In this review we will first introduce the structural elements of the nucleoskeleton and discuss the current literature on how nuclear structure and signaling are altered in relation to environmental and tissue level mechanical cues. We will focus on state-of-the-art techniques to apply mechanical force and methods to measure nuclear mechanics in conjunction with DNA, RNA, and protein visualization in living cells. Ultimately, combining real-time nuclear deformations and chromatin dynamics can be a powerful tool to study mechanisms of how forces affect the dynamics of genome function.

Imaging Organization of RNA Processing within the Nucleus

Within the nucleus, messenger RNA is generated and processed in a highly organized and regulated manner. Messenger RNA processing begins during transcription initiation and continues until the RNA is translated and degraded. Processes such as 5' capping, alternative splicing, and 3' end processing have been studied extensively with biochemical methods and more recently with single-molecule imaging approaches. In this review, we highlight how imaging has helped understand the highly dynamic process of RNA processing. We conclude with open questions and new technological developments that may further our understanding of RNA processing.

SRSF protein kinase 1 modulates RAN translation and suppresses CGG repeat toxicity

Transcribed CGG repeat expansions cause neurodegeneration in Fragile X-associated tremor/ataxia syndrome (FXTAS). CGG repeat RNAs sequester RNA-binding proteins (RBPs) into nuclear foci and undergo repeat-associated non-AUG (RAN) translation into toxic peptides. To identify proteins involved in these processes, we employed a CGG repeat RNA-tagging system to capture repeat-associated RBPs by mass spectrometry in mammalian cells. We identified several SR (serine/arginine-rich) proteins that interact selectively with CGG repeats basally and under cellular stress. These proteins modify toxicity in a Drosophila model of FXTAS. Pharmacologic inhibition of serine/arginine protein kinases (SRPKs), which alter SRSF protein phosphorylation, localization, and activity, directly inhibits RAN translation of CGG and GGGGCC repeats (associated with C9orf72 ALS/FTD) and triggers repeat RNA retention in the nucleus. Lowering SRPK expression suppressed toxicity in both FXTAS and C9orf72 ALS/FTD model flies, and SRPK inhibitors suppressed CGG repeat toxicity in rodent neurons. Together, these findings demonstrate roles for CGG repeat RNA binding proteins in RAN translation and repeat toxicity and support further evaluation of SRPK inhibitors in modulating RAN translation associated with repeat expansion disorders.

In-vivo Single-Molecule Imaging in Yeast: Applications and Challenges

Single-molecule imaging has gained momentum to quantify the dynamics of biomolecules in live cells, as it provides direct real-time measurements of various cellular activities under their physiological environment. Yeast, a simple and widely used eukaryote, serves as a good model system to quantify single-molecule dynamics of various cellular processes because of its low genomic and cellular complexities, as well as its facile ability to be genetically manipulated. In the past decade, significant developments have been made regarding the intracellular labeling of biomolecules (proteins, mRNA, fatty acids), the microscopy setups to visualize single-molecules and capture their fast dynamics, and the data analysis pipelines to interpret such dynamics. In this review, we summarize the current state of knowledge for the single-molecule imaging in live yeast cells to provide a ready reference for beginners. We provide a comprehensive table to demonstrate how various labs tailored the imaging regimes and data analysis pipelines to estimate various biophysical parameters for a variety of biological processes. Lastly, we present current challenges and future directions for developing better tools and resources for single-molecule imaging in live yeast cells.

RNA transport from transcription to localized translation: a single molecule perspective

Transport of mRNAs is an important step of gene expression, which brings the genetic message from the DNA in the nucleus to a precise cytoplasmic location in a regulated fashion. Perturbation of this process can lead to pathologies such as developmental and neurological disorders. In this review, we discuss recent advances in the field of mRNA transport made using single molecule fluorescent imaging approaches. We present an overview of these approaches in fixed and live cells and their input in understanding the key steps of mRNA journey: transport across the nucleoplasm, export through the nuclear pores and delivery to its final cytoplasmic location. This review puts a particular emphasis on the coupling of mRNA transport with translation, such as localization-dependent translational regulation and translation-dependent mRNA localization. We also highlight the recently discovered translation factories, and how cellular and viral RNAs can hijack membrane transport systems to travel in the cytoplasm.

Mutant U2AF1-induced alternative splicing of H2afy (macroH2A1) regulates B-lymphopoiesis in mice

Somatic mutations in spliceosome genes are found in ∼50% of patients with myelodysplastic syndromes (MDS), a myeloid malignancy associated with low blood counts. Expression of the mutant splicing factor U2AF1(S34F) alters hematopoiesis and mRNA splicing in mice. Our understanding of the functionally relevant alternatively spliced target genes that cause hematopoietic phenotypes in vivo remains incomplete. Here, we demonstrate that reduced expression of H2afy1.1, an alternatively spliced isoform of the histone H2A variant gene H2afy, is responsible for reduced B cells in U2AF1(S34F) mice. Deletion of H2afy or expression of U2AF1(S34F) reduces expression of Ebf1 (early B cell factor 1), a key transcription factor for B cell development, and mechanistically, H2AFY is enriched at the EBF1 promoter. Induced expression of H2AFY1.1 in U2AF1(S34F) cells rescues reduced EBF1 expression and B cells numbers in vivo. Collectively, our data implicate alternative splicing of H2AFY as a contributor to lymphopenia induced by U2AF1(S34F) in mice and MDS.

Probing steps in DNA transcription using single-molecule methods

Transcriptional regulation is one of the key steps in determining gene expression. Diverse single-molecule techniques have been applied to characterize the stepwise progression of transcription, yielding complementary results. These techniques include, but are not limited to, fluorescence-based microscopy with single or multiple colors, force measuring and manipulating microscopy using magnetic field or light, and atomic force microscopy. Here, we summarize and evaluate these current methodologies in studying and resolving individual steps in the transcription reaction, which encompasses RNA polymerase binding, initiation, elongation, mRNA production, and termination. We also describe the advantages and disadvantages of each method for studying transcription.

Quantitative imaging of RNA polymerase II activity in plants reveals the single-cell basis of tissue-wide transcriptional dynamics

The responses of plants to their environment are often dependent on the spatiotemporal dynamics of transcriptional regulation. While live-imaging tools have been used extensively to quantitatively capture rapid transcriptional dynamics in living animal cells, the lack of implementation of these technologies in plants has limited concomitant quantitative studies in this kingdom. Here, we applied the PP7 and MS2 RNA-labelling technologies for the quantitative imaging of RNA polymerase II activity dynamics in single cells of living plants as they respond to experimental treatments. Using this technology, we counted nascent RNA transcripts in real time in Nicotiana benthamiana (tobacco) and Arabidopsis thaliana. Examination of heat shock reporters revealed that plant tissues respond to external signals by modulating the proportion of cells that switch from an undetectable basal state to a high-transcription state, instead of modulating the rate of transcription across all cells in a graded fashion. This switch-like behaviour, combined with cell-to-cell variability in transcription rate, results in mRNA production variability spanning three orders of magnitude. We determined that cellular heterogeneity stems mainly from stochasticity intrinsic to individual alleles instead of variability in cellular composition. Together, our results demonstrate that it is now possible to quantitatively study the dynamics of transcriptional programs in single cells of living plants.

Intron removal by the spliceosome: A solo job or a team effort?

Wan et al. (2021) establish a powerful new platform to measure the dynamics of transcription and splicing of endogenous genes in single cells in real time. Combining real-time measurements with multiple deep-sequencing tools reveals an unexpectedly high amount of spliceosome activity, prompting a reconsideration of current models of how introns are removed from pre-mRNA.

Live-cell imaging reveals the spatiotemporal organization of endogenous RNA polymerase II phosphorylation at a single gene

The carboxyl-terminal domain of RNA polymerase II (RNAP2) is phosphorylated during transcription in eukaryotic cells. While residue-specific phosphorylation has been mapped with exquisite spatial resolution along the 1D genome in a population of fixed cells using immunoprecipitation-based assays, the timing, kinetics, and spatial organization of phosphorylation along a single-copy gene have not yet been measured in living cells. Here, we achieve this by combining multi-color, single-molecule microscopy with fluorescent antibody-based probes that specifically bind to different phosphorylated forms of endogenous RNAP2 in living cells. Applying this methodology to a single-copy HIV-1 reporter gene provides live-cell evidence for heterogeneity in the distribution of RNAP2 along the length of the gene as well as Serine 5 phosphorylated RNAP2 clusters that remain separated in both space and time from nascent mRNA synthesis. Computational models determine that 5 to 40 RNAP2 cluster around the promoter during a typical transcriptional burst, with most phosphorylated at Serine 5 within 6 seconds of arrival and roughly half escaping the promoter in ~1.5 minutes. Taken together, our data provide live-cell support for the notion of efficient transcription clusters that transiently form around promoters and contain high concentrations of RNAP2 phosphorylated at Serine 5.

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, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.