4-thiouridine inhibits rRNA synthesis and causes a nucleolar stress response

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

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

Nuclear export is a limiting factor in eukaryotic mRNA metabolism

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

Translation-dependent and -independent mRNA decay occur through mutually exclusive pathways defined by ribosome density during T cell activation

mRNA translation and decay are tightly interconnected processes both in the context of mRNA quality-control pathways and for the degradation of functional mRNAs. Cotranslational mRNA degradation through codon usage, ribosome collisions, and the recruitment of specific proteins to ribosomes is an important determinant of mRNA turnover. However, the extent to which translation-dependent mRNA decay (TDD) and translation-independent mRNA decay (TID) pathways participate in the degradation of mRNAs has not been studied yet. Here we describe a comprehensive analysis of basal and signal-induced TDD and TID in mouse primary CD4+ T cells. Our results indicate that most cellular transcripts are decayed to some extent in a translation-dependent manner. Our analysis further identifies the length of untranslated regions, the density of ribosomes, and GC3 content as important determinants of TDD magnitude. Consistently, all transcripts that undergo changes in ribosome density within their coding sequence upon T cell activation display a corresponding change in their TDD level. Moreover, we reveal a dynamic modulation in the relationship between GC3 content and TDD upon T cell activation, with a reversal in the impact of GC3- and AU3-rich codons. Altogether, our data show a strong and dynamic interconnection between mRNA translation and decay in mammalian primary cells.

Correcting 4sU induced quantification bias in nucleotide conversion RNA-seq data

Nucleoside analogues like 4-thiouridine (4sU) are used to metabolically label newly synthesized RNA. Chemical conversion of 4sU before sequencing induces T-to-C mismatches in reads sequenced from labelled RNA, allowing to obtain total and labelled RNA expression profiles from a single sequencing library. Cytotoxicity due to extended periods of labelling or high 4sU concentrations has been described, but the effects of extensive 4sU labelling on expression estimates from nucleotide conversion RNA-seq have not been studied. Here, we performed nucleotide conversion RNA-seq with escalating doses of 4sU with short-term labelling (1h) and over a progressive time course (up to 2h) in different cell lines. With high concentrations or at later time points, expression estimates were biased in an RNA half-life dependent manner. We show that bias arose by a combination of reduced mappability of reads carrying multiple conversions, and a global, unspecific underrepresentation of labelled RNA emerging during library preparation and potentially global reduction of RNA synthesis. We developed a computational tool to rescue unmappable reads, which performed favourably compared to previous read mappers, and a statistical method, which could fully remove remaining bias. All methods developed here are freely available as part of our GRAND-SLAM pipeline and grandR package.

A kinetic dichotomy between mitochondrial and nuclear gene expression processes

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS subunit biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared with nuclear mRNAs, mt-mRNAs were produced 1,100-fold more, degraded 7-fold faster, and accumulated to 160-fold higher levels. Quantitative modeling and depletion of mitochondrial factors LRPPRC and FASTKD5 identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

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

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

TudS desulfidases recycle 4-thiouridine-5'-monophosphate at a catalytic [4Fe-4S] cluster

In all domains of life, transfer RNAs (tRNAs) contain post-transcriptionally sulfur-modified nucleosides such as 2- and 4-thiouridine. We have previously reported that a recombinant [4Fe-4S] cluster-containing bacterial desulfidase (TudS) from an uncultured bacterium catalyzes the desulfuration of 2- and 4-thiouracil via a [4Fe-5S] cluster intermediate. However, the in vivo function of TudS enzymes has remained unclear and direct evidence for substrate binding to the [4Fe-4S] cluster during catalysis was lacking. Here, we provide kinetic evidence that 4-thiouridine-5'-monophosphate rather than sulfurated tRNA, thiouracil, thiouridine or 4-thiouridine-5'-triphosphate is the preferred substrate of TudS. The occurrence of sulfur- and substrate-bound catalytic intermediates was uncovered from the observed switch of the S = 3/2 spin state of the catalytic [4Fe-4S] cluster to a S = 1/2 spin state upon substrate addition. We show that a putative gene product from Pseudomonas putida KT2440 acts as a TudS desulfidase in vivo and conclude that TudS-like enzymes are widespread desulfidases involved in recycling and detoxifying tRNA-derived 4-thiouridine monophosphate nucleosides for RNA synthesis.

A kinetic dichotomy between mitochondrial and nuclear gene expression drives OXPHOS biogenesis

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared to nuclear mRNAs, mt-mRNAs were produced 700-fold higher, degraded 5-fold faster, and accumulated to 170-fold higher levels. Quantitative modeling and depletion of mitochondrial factors, LRPPRC and FASTKD5, identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

grandR: a comprehensive package for nucleotide conversion RNA-seq data analysis

Metabolic labeling of RNA is a powerful technique for studying the temporal dynamics of gene expression. Nucleotide conversion approaches greatly facilitate the generation of data but introduce challenges for their analysis. Here we present grandR, a comprehensive package for quality control, differential gene expression analysis, kinetic modeling, and visualization of such data. We compare several existing methods for inference of RNA synthesis rates and half-lives using progressive labeling time courses. We demonstrate the need for recalibration of effective labeling times and introduce a Bayesian approach to study the temporal dynamics of RNA using snapshot experiments.

Cycloaddition enabled mutational profiling of 5-vinyluridine in RNA

We report the detection of 5-vinyluridine (5-VUrd) in RNA at single nucleotide resolution via mutational profiling. Maleimide cycloadducts with 5-VUrd in RNA cause a stop in primer extension during reverse transcription, and the full-length cDNA product from reverse transcription contains misincorporation across the cycloadduct site.

Studying miRNA-mRNA Interactions: An Optimized CLIP-Protocol for Endogenous Ago2-Protein

Transcriptome-wide analysis of RNA-binding partners is commonly achieved using UV crosslinking and immunoprecipitation (CLIP). Individual-nucleotide-resolution CLIP (iCLIP)enables identification of the specific position of the protein-RNA interaction. In addition to RNA-binding proteins (RBPs), microRNA (miRNA)-mRNA interactions also play a crucial role in the regulation of gene expression. Argonaute-2 (Ago2) mediates miRNA binding to a multitude of mRNA target sites, enabling the identification of miRNA-mRNA interactions by employing modified Ago2-CLIP protocols. Here, we describe an Ago2-specific CLIP protocol optimized for the use of small quantities of cell material, targeting endogenous Ago2 while avoiding possible methodological biases such as metabolic labeling or Ago2 overexpression and applying the latest advances in CLIP library preparation, the iCLIP2 protocol. In particular, we focus on the optimization of lysis conditions and improved radioactive labeling of the 5' end of the miRNA.

Probing Nascent RNA with Metabolic Incorporation of Modified Nucleosides

The discovery of previously unknown functional roles of RNA in biological systems has led to increased interest in revealing novel RNA molecules as therapeutic targets and the development of tools to better understand the role of RNA in cells. RNA metabolic labeling broadens the scope of studying RNA by incorporating of unnatural nucleobases and nucleosides with bioorthogonal handles that can be utilized for chemical modification of newly synthesized cellular RNA. Such labeling of RNA provides access to applications including measurement of the rates of synthesis and decay of RNA, cellular imaging for RNA localization, and selective enrichment of nascent RNA from the total RNA pool. Several unnatural nucleosides and nucleobases have been shown to be incorporated into RNA by endogenous RNA synthesis machinery of the cells. RNA metabolic labeling can also be performed in a cell-specific manner, where only cells expressing an essential enzyme incorporate the unnatural nucleobase into their RNA. Although several discoveries have been enabled by the current RNA metabolic labeling methods, some key challenges still exist: (i) toxicity of unnatural analogues, (ii) lack of RNA-compatible conjugation chemistries, and (iii) background incorporation of modified analogues in cell-specific RNA metabolic labeling. In this Account, we showcase work done in our laboratory to overcome these challenges faced by RNA metabolic labeling.To begin, we discuss the cellular pathways that have been utilized to perform RNA metabolic labeling and study the interaction between nucleosides and nucleoside kinases. Then we discuss the use of vinyl nucleosides for metabolic labeling and demonstrate the low toxicity of 5-vinyluridine (5-VUrd) compared to other widely used nucleosides. Next, we discuss cell-specific RNA metabolic labeling with unnatural nucleobases, which requires the expression of a specific phosphoribosyl transferase (PRT) enzyme for incorporation of the nucleobase into RNA. In the course of this work, we discovered the enzyme uridine monophosphate synthase (UMPS), which is responsible for nonspecific labeling with modified uracil nucleobases. We were able to overcome this background labeling by discovering a mutant uracil PRT (UPRT) that demonstrates highly specific RNA metabolic labeling with 5-vinyluracil (5-VU). Furthermore, we discuss the optimization of inverse-electron-demand Diels-Alder (IEDDA) reactions for performing chemical modification of vinyl nucleosides to achieve covalent conjugation of RNA without transcript degradation. Finally, we highlight our latest endeavor: the development of mutually orthogonal chemical reactions for selective labeling of 5-VUrd and 2-vinyladenosine (2-VAdo), which allows for potential use of multiple vinyl nucleosides for simultaneous investigation of multiple cellular processes involving RNA. We hope that our methods and discoveries encourage scientists studying biological systems to include RNA metabolic labeling in their toolkit for studying RNA and its role in biological systems.

Control of immediate early gene expression by CPEB4-repressor complex-mediated mRNA degradation

Background

Cytoplasmic polyadenylation element-binding protein 4 (CPEB4) is known to associate with cytoplasmic polyadenylation elements (CPEs) located in the 3′ untranslated region (UTR) of specific mRNAs and assemble an activator complex promoting the translation of target mRNAs through cytoplasmic polyadenylation.

Results

Here, we find that CPEB4 is part of an alternative repressor complex that mediates mRNA degradation by associating with the evolutionarily conserved CCR4-NOT deadenylase complex. We identify human CPEB4 as an RNA-binding protein (RBP) with enhanced association to poly(A) RNA upon inhibition of class I histone deacetylases (HDACs), a condition known to cause widespread degradation of poly(A)-containing mRNA. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) analysis using endogenously tagged CPEB4 in HeLa cells reveals that CPEB4 preferentially binds to the 3′UTR of immediate early gene mRNAs, at G-containing variants of the canonical U- and A-rich CPE located in close proximity to poly(A) sites. By transcriptome-wide mRNA decay measurements, we find that the strength of CPEB4 binding correlates with short mRNA half-lives and that loss of CPEB4 expression leads to the stabilization of immediate early gene mRNAs. Akin to CPEB4, we demonstrate that CPEB1 and CPEB2 also confer mRNA instability by recruitment of the CCR4-NOT complex.

Conclusions

While CPEB4 was previously known for its ability to stimulate cytoplasmic polyadenylation, our findings establish an additional function for CPEB4 as the RNA adaptor of a repressor complex that enhances the degradation of short-lived immediate early gene mRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02760-5.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

RNA–protein interactomes as invaluable resources to study RNA viruses: Insights from SARS CoV‐2 studies

Understanding the molecular mechanisms of severe respiratory syndrome coronavirus 2 (SARS‐CoV‐2) infection is essential for the successful development of therapeutic strategies against the COVID‐19 pandemic. Numerous studies have focused on the identification of host factors and cellular pathways involved in the viral replication cycle. The speed and magnitude of hijacking the translation machinery of host mRNA, and shutting down host transcription are still not well understood. Since SARS‐CoV‐2 relies on host RNA‐binding proteins for the infection progression, several efforts have been made to define the SARS‐CoV‐2 RNA‐bound proteomes (RNA–protein interactomes). Methodologies that enable the systemic capture of protein interactors of given RNA in vivo have been adapted for the identification of the SARS‐CoV‐2 RNA interactome. The obtained proteomic data aided by genome‐wide and targeted CRISPR perturbation screens, revealed host factors with either pro‐ or anti‐viral activity and highlighted cellular processes and factors involved in host response. We focus here on the recent studies on SARS‐CoV‐2 RNA–protein interactomes, with regard to both the technological aspects of RNA interactome capture methods and the obtained results. We also summarize several related studies, which were used in the interpretation of the SARS‐CoV‐2 RNA–protein interactomes. These studies provided the selection of host factors that are potentially suitable candidates for antiviral therapy. Finally, we underscore the importance of RNA–protein interactome studies in regard to the effective development of antiviral strategies against current and future threats.

This article is categorized under:

RNA Interactions with Proteins and Other Molecules > Protein‐RNA Interactions: Functional Implications

RNA in Disease and Development > RNA in Disease

RNA Methods > RNA Analyses in Cells

Metabolic incorporation of electron-rich ribonucleosides enhances APEX-seq for profiling spatially restricted nascent transcriptome

The spatial arrangement of newly synthesized transcriptome in eukaryotic cells underlies various biological processes including cell proliferation and differentiation. In this study, we combine metabolic incorporation of electron-rich ribonucleosides (e.g., 6-thioguanosine and 4-thiouridine) with a peroxidase-mediated proximity-dependent RNA labeling technique (APEX-seq) to develop a sensitive method, termed MERR APEX-seq, for selectively profiling newly transcribed RNAs at specific subcellular locations in live cells. We demonstrate that MERR APEX-seq is 20-fold more efficient than APEX-seq and offers both high spatial specificity and high coverage in mitochondrial matrix. At the ER membrane, 91% of the transcripts captured by MERR APEX-seq encode for secretory pathway proteins, thus demonstrating the high spatial specificity of MERR APEX-seq in open subcellular compartments. Application of MERR APEX-seq to the nuclear lamina of human cells reveals a local transcriptome of 1,012 RNAs, many of which encode for nuclear proteins involved in histone modification, chromosomal structure maintenance, and RNA processing.

Transcriptome-Wide Identification of Coding and Noncoding RNA-Binding Proteins Defines the Comprehensive RNA Interactome of Leishmania mexicana

Proteomic profiling of RNA-binding proteins in Leishmania is currently limited to polyadenylated mRNA-binding proteins, leaving proteins that interact with nonadenylated RNAs, including noncoding RNAs and pre-mRNAs, unidentified. Using a combination of unbiased orthogonal organic phase separation methodology and tandem mass tag-labeling-based high resolution quantitative proteomic mass spectrometry, we robustly identified 2,417 RNA-binding proteins, including 1289 putative novel non-poly(A)-RNA-binding proteins across the two main Leishmania life cycle stages. Eight out of 20 Leishmania deubiquitinases, including the recently characterized L. mexicana DUB2 with an elaborate RNA-binding protein interactome were exclusively identified in the non-poly(A)-RNA-interactome. Additionally, an increased representation of WD40 repeat domains were observed in the Leishmania non-poly(A)-RNA-interactome, thus uncovering potential involvement of this protein domain in RNA-protein interactions in Leishmania. We also characterize the protein-bound RNAs using RNA-sequencing and show that in addition to protein coding transcripts ncRNAs are also enriched in the protein-RNA interactome. Differential gene expression analysis revealed enrichment of 142 out of 195 total L. mexicana protein kinase genes in the protein-RNA-interactome, suggesting important role of protein-RNA interactions in the regulation of the Leishmania protein kinome. Additionally, we characterize the quantitative changes in RNA-protein interactions in hundreds of Leishmania proteins following inhibition of heat shock protein 90 (Hsp90). Our results show that the Hsp90 inhibition in Leishmania causes widespread disruption of RNA-protein interactions in ribosomal proteins, proteasomal proteins and translation factors in both life cycle stages, suggesting downstream effect of the inhibition on protein synthesis and degradation pathways in Leishmania. This study defines the comprehensive RNA interactome of Leishmania and provides in-depth insight into the widespread involvement of RNA-protein interactions in Leishmania biology. IMPORTANCE Advances in proteomics and mass spectrometry have revealed the mRNA-binding proteins in many eukaryotic organisms, including the protozoan parasites Leishmania spp., the causative agents of leishmaniasis, a major infectious disease in over 90 tropical and subtropical countries. However, in addition to mRNAs, which constitute only 2 to 5% of the total transcripts, many types of non-coding RNAs participate in crucial biological processes. In Leishmania, RNA-binding proteins serve as primary gene regulators. Therefore, transcriptome-wide identification of RNA-binding proteins is necessary for deciphering the distinctive posttranscriptional mechanisms of gene regulation in Leishmania. Using a combination of highly efficient orthogonal organic phase separation method and tandem mass tag-labeling-based quantitative proteomic mass spectrometry, we provide unprecedented comprehensive molecular definition of the total RNA interactome across the two main Leishmania life cycle stages. In addition, we characterize for the first time the quantitative changes in RNA-protein interactions in Leishmania following inhibition of heat shock protein 90, shedding light into hitherto unknown large-scale downstream molecular effect of the protein inhibition in the parasite. This work provides insight into the importance of total RNA-protein interactions in Leishmania, thus significantly expanding our knowledge of the emergence of RNA-protein interactions in Leishmania biology.

Enhanced immunoprecipitation techniques for the identification of RNA-binding protein partners: IGF2BP1 interactions in mammary epithelial cells

RNA-binding proteins (RBPs) regulate the expression of large cohorts of RNA species to produce programmatic changes in cellular phenotypes. To describe the function of RBPs within a cell, it is key to identify their mRNA-binding partners. This is often done by crosslinking nucleic acids to RBPs, followed by chemical release of the nucleic acid fragments for analysis. However, this methodology is lengthy, which involves complex processing with attendant sample losses, thus large amounts of starting materials and prone to artifacts. To evaluate potential alternative technologies, we tested “exclusion-based” purification of immunoprecipitates (IFAST or SLIDE) and report here that these methods can efficiently, rapidly, and specifically isolate RBP–RNA complexes. The analysis requires less than 1% of the starting material required for techniques that include crosslinking. Depending on the antibody used, 50% to 100% starting protein can be retrieved, facilitating the assay of endogenous levels of RBPs; the isolated ribonucleoproteins are subsequently analyzed using standard techniques, to provide a comprehensive portrait of RBP complexes. Using exclusion-based techniques, we show that the mRNA-binding partners for RBP IGF2BP1 in cultured mammary epithelial cells are enriched in mRNAs important for detoxifying superoxides (specifically glutathione peroxidase [GPX]-1 and GPX-2) and mRNAs encoding mitochondrial proteins. We show that these interactions are functionally significant, as loss of function of IGF2BP1 leads to destabilization of GPX mRNAs and reduces mitochondrial membrane potential and oxygen consumption. We speculate that this underlies a consistent requirement for IGF2BP1 for the expression of clonogenic activity in vitro.

Non-radioactive In Vivo Labeling of RNA with 4-Thiouracil

RNA molecules and their expression dynamics play essential roles in the establishment of complex cellular phenotypes and/or in the rapid cellular adaption to environmental changes. Accordingly, analyzing RNA expression remains an important step to understand the molecular basis controlling the formation of cellular phenotypes, cellular homeostasis or disease progression. Steady-state RNA levels in the cells are controlled by the sum of highly dynamic molecular processes contributing to RNA expression and can be classified in transcription, maturation and degradation. The main goal of analyzing RNA dynamics is to disentangle the individual contribution of these molecular processes to the life cycle of a given RNA under different physiological conditions. In the recent years, the use of nonradioactive nucleotide/nucleoside analogs and improved chemistry, in combination with time-dependent and high-throughput analysis, have greatly expanded our understanding of RNA metabolism across various cell types, organisms, and growth conditions.In this chapter, we describe a step-by-step protocol allowing pulse labeling of RNA with the nonradioactive nucleotide analog, 4-thiouracil , in the eukaryotic model organism Saccharomyces cerevisiae and the model archaeon Haloferax volcanii .

A Modified Cross-Linking Analysis of cDNAs (CRAC ) Protocol for Detecting RNA-Protein Interactions and Transcription at Single-Nucleotide Resolution

Detecting protein-RNA interactions in vivo is essential for deciphering many important cellular pathways. Several methods have been described for this purpose, among which cross-linking analysis of cDNA, CRAC. This method relies on a first step of UV cross-linking of living yeast cells and several subsequent steps of purification of the protein-RNA complexes, some of which under denaturing condition. Without altering the general principle of the method, we have modified and improved the protocol, with the specific aim of sequencing the nascent RNA isolated from transcription complexes and generate high-resolution and directional transcription maps.

Measuring mRNA Decay with Roadblock-qPCR

The control of mRNA stability is fundamental to gene regulation, and a deeper understanding of this post-transcriptional regulatory step can provide key insights into gene function. Measuring mRNA half-lives directly, however, is challenging. The most common strategies for evaluating mRNA stability and decay involve blocking general transcription and then measuring the decline in mRNA levels over time. The downside of these approaches, however, is that they severely impact cell function and viability, indirectly perturbing gene expression. Here, we describe Roadblock-qPCR, a simple method for measuring mRNA decay kinetics in living cells that is both economical and quick. Cells are first incubated with the nucleoside analog 4-thiouridine (4sU), which is readily incorporated into nascent mRNAs during transcription. RNA is then extracted and treated with N-ethylmaleimide (NEM), a sulfhydryl alkylating agent that selectively modifies 4sU, before proceeding to cDNA synthesis. Because the NEM-modified 4sU creates a chemical "roadblock" that interferes with reverse transcription, this treatment ultimately results in the depletion of the nascent 4sU-containing transcripts from the cDNA pool. As such, the decay rate of the non-4sU-labeled pre-existing mRNAs can be monitored by quantitative PCR (qPCR). In combination with spike-in standards, this approach can be used to efficiently and accurately measure the half-lives of endogenous mRNAs with a wide range of stabilities, while avoiding the artifacts of transcription shutoff strategies. © 2022 Wiley Periodicals LLC. Basic Protocol: Roadblock-qPCR Support Protocol: Synthesis of spike-in mRNA.

The influence of 4-thiouridine labeling on pre-mRNA splicing outcomes

Metabolic labeling is a widely used tool to investigate different aspects of pre-mRNA splicing and RNA turnover. The labeling technology takes advantage of native cellular machineries where a nucleotide analog is readily taken up and incorporated into nascent RNA. One such analog is 4-thiouridine (4sU). Previous studies demonstrated that the uptake of 4sU at elevated concentrations (>50μM) and extended exposure led to inhibition of rRNA synthesis and processing, presumably induced by changes in RNA secondary structure. Thus, it is possible that 4sU incorporation may also interfere with splicing efficiency. To test this hypothesis, we carried out splicing analyses of pre-mRNA substrates with varying levels of 4sU incorporation (0-100%). We demonstrate that increased incorporation of 4sU into pre-mRNAs decreased splicing efficiency. The overall impact of 4sU labeling on pre-mRNA splicing efficiency negatively correlates with the strength of splice site signals such as the 3' and the 5' splice sites. Introns with weaker splice sites are more affected by the presence of 4sU. We also show that transcription by T7 polymerase and pre-mRNA degradation kinetics were impacted at the highest levels of 4sU incorporation. Increased incorporation of 4sU caused elevated levels of abortive transcripts, and fully labeled pre-mRNA is more stable than its uridine-only counterpart. Cell culture experiments show that a small number of alternative splicing events were modestly, but statistically significantly influenced by metabolic labeling with 4sU at concentrations considered to be tolerable (40 μM). We conclude that at high 4sU incorporation rates small, but noticeable changes in pre-mRNA splicing can be detected when splice sites deviate from consensus. Given these potential 4sU artifacts, we suggest that appropriate controls for metabolic labeling experiments need to be included in future labeling experiments.

Interrogating the transcriptome with metabolically incorporated ribonucleosides

RNA is a central player in biological processes, but there remain major gaps in our understanding of transcriptomic processes and the underlying biochemical mechanisms regulating RNA in cells. A powerful strategy to facilitate molecular analysis of cellular RNA is the metabolic incorporation of chemical probes. In this review, we discuss current approaches for RNA metabolic labeling with modified ribonucleosides and their integration with Next-Generation Sequencing, mass spectrometry-based proteomics, and fluorescence microscopy in order to interrogate RNA behavior in its native context.

Capture of the newly transcribed RNA interactome using click chemistry

Application of synthetic nucleoside analogues to capture newly transcribed RNAs has unveiled key features of RNA metabolism. Whether this approach could be adapted to isolate the RNA-bound proteome (RNA interactome) was, however, unexplored. We have developed a new method (capture of the newly transcribed RNA interactome using click chemistry, or RICK) for the systematic identification of RNA-binding proteins based on the incorporation of 5-ethynyluridine into newly transcribed RNAs followed by UV cross-linking and click chemistry-mediated biotinylation. The RNA-protein adducts are then isolated by affinity capture using streptavidin-coated beads. Through high-throughput RNA sequencing and mass spectrometry, the RNAs and proteins can be elucidated globally. A typical RICK experimental procedure takes only 1 d, excluding the steps of cell preparation, 5-ethynyluridine labeling, validation (silver staining, western blotting, quantitative reverse-transcription PCR (qRT-PCR) or RNA sequencing (RNA-seq)) and proteomics. Major advantages of RICK are the capture of RNA-binding proteins interacting with any type of RNA and, particularly, the ability to discern between newly transcribed and steady-state RNAs through controlled labeling. Thanks to its versatility, RICK will facilitate the characterization of the total and newly transcribed RNA interactome in different cell types and conditions.

Mutually Orthogonal Bioconjugation of Vinyl Nucleosides for RNA Metabolic Labeling

We report a strategy for the orthogonal conjugation of the vinyl nucleosides, 5-vinyluridine (5-VU) and 2-vinyladenosine (2-VA), via selective reactivity with maleimide and tris(2-carboxyethyl)phosphine (TCEP), respectively. The orthogonality was investigated using density functional theory (DFT) and confirmed by reactions with vinyl nucleosides. Further, these chemistries were used to modify RNA for fluorescent cell imaging. These reactions allow for the expanded use of RNA metabolic labeling to study nascent RNA expression within different RNA populations.

[Advances in technologies for large-scale enrichment and identification of ribonucleic acid-protein complexes]

核糖核酸(RNA)在细胞中并非单独存在,从它们产生到被降解的过程中与大量蛋白质发生相互作用,RNA结合蛋白(RNA-binding proteins, RBPs)能与RNA结合形成RNA-蛋白质复合物(RP复合物),并以这种复合物的形式发挥生理功能。RNAs或RBPs任一组分的异常与缺失都会影响RP复合物的正常生理功能,从而导致疾病的发生,如代谢异常、肌肉萎缩症、自身免疫性疾病和癌症。因此,定性定量分析RBPs及其在正常细胞和肿瘤细胞中与RNAs靶标之间的复杂相互作用网络有助于挖掘RP复合物在肿瘤发生发展中的作用,开发肿瘤生物标志物和新的治疗方式。要深入研究和理解RNAs与RBPs的相互作用网络,须依赖组学技术对RP复合物进行大规模鉴定。而作为在组学层面系统性解析RP复合物组成、含量和功能的第一步,大规模富集RP复合物极具挑战性。为了解决这一难题,研究者们发展了各种富集鉴定策略。该文针对RP复合物富集策略的最新进展进行了综述,包括紫外光交联和免疫沉淀(crosslinking and immunoprecipitation, CLIP)及其衍生技术、基于“点击化学”的富集策略和基于相分离的富集策略,比较分析了它们的技术原理、优缺点,以方便研究者们选择合适的策略来解决感兴趣的生物学问题。该文最后总结了当前的RP复合物富集方法仍然存在富集效率低和操作繁琐等亟需解决的技术挑战,为富集策略的发展提供了研究方向。

In vivo 5-ethynyluridine (EU) labelling detects reduced transcription in Purkinje cell degeneration mouse mutants, but can itself induce neurodegeneration

Fluorescent staining of newly transcribed RNA via metabolic labelling with 5-ethynyluridine (EU) and click chemistry enables visualisation of changes in transcription, such as in conditions of cellular stress. Here, we tested whether EU labelling can be used to examine transcription in vivo in mouse models of nervous system disorders. We show that injection of EU directly into the cerebellum results in reproducible labelling of newly transcribed RNA in cerebellar neurons and glia, with cell type-specific differences in relative labelling intensities, such as Purkinje cells exhibiting the highest levels. We also observed EU-labelling accumulating into cytoplasmic inclusions, indicating that EU, like other modified uridines, may introduce non-physiological properties in labelled RNAs. Additionally, we found that EU induces Purkinje cell degeneration nine days after EU injection, suggesting that EU incorporation not only results in abnormal RNA transcripts, but also eventually becomes neurotoxic in highly transcriptionally-active neurons. However, short post-injection intervals of EU labelling in both a Purkinje cell-specific DNA repair-deficient mouse model and a mouse model of spinocerebellar ataxia 1 revealed reduced transcription in Purkinje cells compared to controls. We combined EU labelling with immunohistology to correlate altered EU staining with pathological markers, such as genotoxic signalling factors. These data indicate that the EU-labelling method provided here can be used to identify changes in transcription in vivo in nervous system disease models.

3' End sequencing of pA+ and pA- RNAs

The identity and metabolism of RNAs are often governed by their 5' and 3' ends. Single gene loci produce a variety of transcript isoforms, varying primarily in their RNA 3' end status and consequently facing radically different cellular fates. Knowledge about RNA termini is therefore key to understanding the diverse RNA output from individual transcription units. In addition, the 3' end of a nascent RNA at the catalytic center of RNA polymerase provides a precise and strand-specific measure of the transcription process. Here, we describe a modified RNA 3' end sequencing method, that utilizes the in vivo metabolic labeling of RNA followed by its purification and optional in vitro polyadenylation to provide a comprehensive view of all RNA 3' ends. The strategy offers the advantages of (i) nucleotide resolution mapping of RNA 3' ends, (ii) increased sequencing depth of lowly abundant RNA and (iii) inference of RNA 3' end polyadenylation status. We have used the method to study RNA decay and transcription termination mechanisms with the potential utility to a wider range of biological questions.

Early nucleolar responses differentiate mechanisms of cell death induced by oxaliplatin and cisplatin

Recent reports provide evidence that the platinum chemotherapeutic oxaliplatin causes cell death via ribosome biogenesis stress, while cisplatin causes cell death via the DNA damage response (DDR). Underlying differences in mechanisms that might initiate disparate routes to cell death by these two broadly used platinum compounds have not yet been carefully explored. Additionally, prior studies had demonstrated that cisplatin can also inhibit ribosome biogenesis. Therefore, we sought to directly compare the initial influences of oxaliplatin and cisplatin on nucleolar processes and on the DDR. Using pulse-chase experiments, we found that at equivalent doses, oxaliplatin but not cisplatin significantly inhibited ribosomal RNA (rRNA) synthesis by Pol I, but neither compound affected rRNA processing. Inhibition of rRNA synthesis occurred as early as 90 min after oxaliplatin treatment in A549 cells, concurrent with the initial redistribution of the nucleolar protein nucleophosmin (NPM1). We observed that the nucleolar protein fibrillarin began to redistribute by 6 h after oxaliplatin treatment and formed canonical nucleolar caps by 24 h. In cisplatin-treated cells, DNA damage, as measured by γH2AX immunofluorescence, was more extensive, whereas nucleolar organization was unaffected. Taken together, our results demonstrate that oxaliplatin causes early nucleolar disruption via inhibition of rRNA synthesis accompanied by NPM1 relocalization and subsequently causes extensive nucleolar reorganization, while cisplatin causes early DNA damage without significant nucleolar disruption. These data support a model in which, at clinically relevant doses, cisplatin kills cells via the canonical DDR, and oxaliplatin kills cells via ribosome biogenesis stress, specifically via rapid inhibition of rRNA synthesis.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and β-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

The search for RNA-binding proteins: a technical and interdisciplinary challenge

RNA-binding proteins are customarily regarded as important facilitators of gene expression. In recent years, RNA–protein interactions have also emerged as a pervasive force in the regulation of homeostasis. The compendium of proteins with provable RNA-binding function has swelled from the hundreds to the thousands astride the partnership of mass spectrometry-based proteomics and RNA sequencing. At the foundation of these advances is the adaptation of RNA-centric capture methods that can extract bound protein that has been cross-linked in its native environment. These methods reveal snapshots in time displaying an extensive network of regulation and a wealth of data that can be used for both the discovery of RNA-binding function and the molecular interfaces at which these interactions occur. This review will focus on the impact of these developments on our broader perception of post-transcriptional regulation, and how the technical features of current capture methods, as applied in mammalian systems, create a challenging medium for interpretation by systems biologists and target validation by experimental researchers.

Roadblock-qPCR: A simple and inexpensive strategy for targeted measurements of mRNA stability

The stability of mRNAs is fundamental to determining expression level and dynamics. Nonetheless, current approaches for measuring transcript half-lives (e.g. transcription shutoff) are generally toxic or technically complex. Here we describe an alternative strategy for targeted measurements of endogenous mRNA stability that is simple, inexpensive, and non-toxic. Cells are first metabolically labelled with the nucleoside analog 4-thiouridine (4sU). Extracted mRNA can then be treated with the thiol-reactive compound N-ethylmaleimide. This compound modifies 4sU nucleotides and sterically interferes with reverse transcription of 4sU-containing transcripts, disrupting their conversion into cDNA. The decay rate of non-4sU-containing pre-existing mRNA can then be monitored by quantitative PCR (qPCR). Importantly, this approach avoids the biochemical isolation of 4sU-labelled transcripts and/or RNA-seq analysis required for other metabolic-labelling strategies. In summary, our method combines the simplicity of "transcription shutoff" strategies with the accuracy of metabolic-labelling strategies for measurements of mRNA stability across a wide range of half-lives.

Global analysis of RNA-binding protein dynamics by comparative and enhanced RNA interactome capture

Interactions between RNA-binding proteins (RBPs) and RNAs are critical to cell biology. However, methods to comprehensively and quantitatively assess these interactions within cells were lacking. RNA interactome capture (RIC) uses in vivo UV crosslinking, oligo(dT) capture, and proteomics to identify RNA-binding proteomes. Recent advances have empowered RIC to quantify RBP responses to biological cues such as metabolic imbalance or virus infection. Enhanced RIC exploits the stronger binding of locked nucleic acid (LNA)-containing oligo(dT) probes to poly(A) tails to maximize RNA capture selectivity and efficiency, profoundly improving signal-to-noise ratios. The subsequent analytical use of SILAC and TMT proteomic approaches, together with high-sensitivity sample preparation and tailored statistical data analysis, substantially improves RIC's quantitative accuracy and reproducibility. This optimized approach is an extension of the original RIC protocol. It takes 3 d plus 2 weeks for proteomics and data analysis and will enable the study of RBP dynamics under different physiological and pathological conditions.

Clip for studying protein-RNA interactions that regulate virus replication

Viral and cellular RNA-binding proteins regulate numerous key steps in the replication of diverse virus genera. Viruses efficiently co-opt the host cell machinery for purposes such as transcription, splicing and subcellular localization of viral genomes. Though viral RNAs often need to resemble cellular RNAs to effectively utilize the cellular machinery, they still retain unique sequence and structural features for recognition by viral proteins for processes such as RNA polymerization, RNA export and selective packaging into virus particles. While beneficial for virus replication, distinct features of viral nucleic acids can also be recognized as foreign by several host defense proteins. Development of the crosslinking immunoprecipitation coupled with sequencing (CLIP) approach has allowed the study of viral and cellular RNA binding proteins that regulate critical aspects of viral replication in unprecedented detail. By combining immunoprecipitation of covalently crosslinked protein-RNA complexes with high-throughput sequencing, CLIP provides a global account of RNA sequences bound by RNA-binding proteins of interest in physiological settings and at near-nucleotide resolution. Here, we describe the step-by-step application of the CLIP methodology within the context of two cellular splicing regulatory proteins, hnRNP A1 and hnRNP H1 that regulate HIV-1 splicing. In principle, this versatile protocol can be applied to many other viral and cellular RNA-binding proteins.

Cell- and Polymerase-Selective Metabolic Labeling of Cellular RNA with 2'-Azidocytidine

Metabolic labeling of cellular RNA is a powerful approach to investigate RNA biology. In addition to revealing whole transcriptome dynamics, targeted labeling strategies can be used to study individual RNA subpopulations within complex systems. Here, we describe a strategy for cell- and polymerase-selective RNA labeling with 2'-azidocytidine (2'-AzCyd), a modified nucleoside amenable to bioorthogonal labeling with SPAAC chemistry. In contrast to 2'-OH-containing pyrimidine ribonucleosides, which rely upon uridine-cytidine kinase 2 (UCK2) for activation, 2'-AzCyd is phosphorylated by deoxycytidine kinase (dCK), and we find that expression of dCK mediates cell-selective 2'-AzCyd labeling. Further, 2'-AzCyd is primarily incorporated into rRNA and displays low cytotoxicity and high labeling efficiency. We apply our system to analyze the turnover of rRNA during ribophagy induced by oxidative stress or mTOR inhibition to show that 28S and 18S rRNAs undergo accelerated degradation. Taken together, our work provides a general approach for studying dynamic RNA behavior with cell and polymerase specificity and reveals fundamental insights into nucleotide and nucleic acid metabolism.

RNA Conformation Capture by Proximity Ligation

RNA proximity ligation is a set of molecular biology techniques used to analyze the conformations and spatial proximity of RNA molecules within cells. A typical experiment starts with cross-linking of a biological sample using UV light or psoralen, followed by partial fragmentation of RNA, RNA-RNA ligation, library preparation, and high-throughput sequencing. In the past decade, proximity ligation has been used to study structures of individual RNAs, networks of interactions between small RNAs and their targets, and whole RNA-RNA interactomes, in models ranging from bacteria to animal tissues and whole animals. Here, we provide an overview of the field, highlight the main findings, review the recent experimental and computational developments, and provide troubleshooting advice for new users. In the final section, we draw parallels between DNA and RNA proximity ligation and speculate on possible future research directions.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the “fifth” ribonucleotide in 1951. Since then, the ever‐increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal–Hreidarsson syndrome, Bowen–Conradi syndrome, or Williams–Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications.

This article is categorized under:

RNA Processing > RNA Editing and Modification

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

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

Metabolic Labeling of RNAs Uncovers Hidden Features and Dynamics of the Arabidopsis Transcriptome

Transcriptome analysis by RNA sequencing (RNA-seq) has become an indispensable research tool in modern plant biology. Virtually all RNA-seq studies provide a snapshot of the steady state transcriptome, which contains valuable information about RNA populations at a given time but lacks information about the dynamics of RNA synthesis and degradation. Only a few specialized sequencing techniques, such as global run-on sequencing, have been used to provide information about RNA synthesis rates in plants. Here, we demonstrate that RNA labeling with the modified, nontoxic uridine analog 5-ethynyl uridine (5-EU) in Arabidopsis (Arabidopsis thaliana) seedlings provides insight into plant transcriptome dynamics. Pulse labeling with 5-EU revealed nascent and unstable RNAs, RNA processing intermediates generated by splicing, and chloroplast RNAs. Pulse-chase experiments with 5-EU allowed us to determine RNA stabilities without the need for chemical transcription inhibitors such as actinomycin and cordycepin. Inhibitor-free, genome-wide analysis of polyadenylated RNA stability via 5-EU pulse-chase experiments revealed RNAs with shorter half-lives than those reported after chemical inhibition of transcription. In summary, our results indicate that the Arabidopsis nascent transcriptome contains unstable RNAs and RNA processing intermediates and suggest that polyadenylated RNAs have low stability in plants. Our technique lays the foundation for easy, affordable, nascent transcriptome analysis and inhibitor-free analysis of RNA stability in plants.

Organic phase separation opens up new opportunities to interrogate the RNA-binding proteome

Protein-RNA interactions regulate all aspects of RNA metabolism and are crucial to the function of catalytic ribonucleoproteins. Until recently, the available technologies to capture RNA-bound proteins have been biased toward poly(A) RNA-binding proteins (RBPs) or involve molecular labeling, limiting their application. With the advent of organic-aqueous phase separation-based methods, we now have technologies that efficiently enrich the complete suite of RBPs and enable quantification of RBP dynamics. These flexible approaches to study RBPs and their bound RNA open up new research avenues for systems-level interrogation of protein-RNA interactions.

A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells

The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.

An optimized chemical-genetic method for cell-specific metabolic labeling of RNA

Tissues and organs are composed of diverse cell types, which poses a major challenge for cell-specific gene expression profiling. Current metabolic labeling methods rely on the inability of mammalian cells to incorporate exogenous pyrimidine analogs, which are then co-opted by ectopically-expressed enzymes. We demonstrate that mammalian cells can incorporate uracil analogs and characterize the enzymatic pathways responsible for high background incorporation. To overcome these limitations, we developed a novel small-molecule/enzyme pair consisting of uridine-cytidine kinase 2 (UCK2) and 2’-azidouridine (2’AzUd). We demonstrate that 2’AzUd is only incorporated in UCK2-expressing cells and characterize selectivity mechanisms using molecular dynamics and X-ray crystallography. Furthermore, this pair can be used to purify and track RNA from specific cellular populations, making it ideal for high-resolution cell-specific RNA labeling. Overall, these results reveal novel aspects of mammalian salvage pathways and serve as a new benchmark for designing, characterizing and evaluating cell-specific biomolecule labeling methodologies.

Using TTchem-seq for profiling nascent transcription and measuring transcript elongation

The dynamics of transcription can be studied genome wide by high-throughput sequencing of nascent and newly synthesized RNA. 4-thiouridine (4SU) labeling in vivo enables the specific capture of such new transcripts, with 4SU residues being tagged by biotin linkers and captured using streptavidin beads before library production and high-throughput sequencing. To achieve high-resolution profiles of transcribed regions, an RNA fragmentation step before biotin tagging was introduced, in an approach known as transient transcriptome sequencing (TT-seq). We recently introduced a chemical approach for RNA fragmentation that we refer to as TTchem-seq. We describe how TTchem-seq can be used in combination with transient inhibition of early elongation using the reversible CDK9 inhibitor, 5,6-dichlorobenzimidazole 1-β-D-ribofuranoside (DRB), to measure RNA polymerase II (RNAPII) elongation rates in vivo, a technique we call DRB/TTchem-seq. Here, we provide detailed protocols for carrying out TTchem-seq and DRB/TTchem-seq, including computational analysis. Experiments and data analysis can be performed over a period of 10-13 d and require molecular biology and bioinformatics skills.

Thiouridine-to-Cytidine Conversion Sequencing (TUC-Seq) to Measure mRNA Transcription and Degradation Rates

The study of RNA dynamics, specifically RNA transcription and decay rates, has gained increasing attention in recent years because various mechanisms have been discovered that affect mRNA half-life, thereby ultimately controlling protein output. Therefore, there is a need for methods enabling minimally invasive, simple and high-throughput determination of RNA stability that can be applied to determine RNA transcription and decay rates in cells and organisms. We have recently developed a protocol which we named TUC-seq to directly distinguish newly synthesized transcripts from the preexisting pool of transcripts by metabolic labeling of nascent RNAs with 4-thiouridine (4sU) followed by osmium tetroxide-mediated conversion of 4sU to cytidine (C) and direct sequencing. In contrast to other related methods (SLAM-seq, TimeLapse-seq), TUC-seq converts 4sU to a native C instead of an alkylated or otherwise modified nucleoside derivative. TUC-seq can be applied to any cell type that is amenable to 4sU labeling. By employing different labeling strategies (pulse or pulse-chase labeling), it is suitable for a broad field of applications and provides a fast and highly efficient means to determine mRNA transcription and decay rates.

Alterations of signaling pathways in response to chemical perturbations used to measure mRNA decay rates in yeast

Assessing variations in mRNA stability typically involves inhibiting transcription either globally or in a gene-specific manner. Alternatively, mRNA pulse-labeling strategies offer a means to calculate mRNA stability without inhibiting transcription. However, key stress-responsive cell signaling pathways, which affect mRNA stability, may themselves be perturbed by the approaches used to measure mRNA stability, leading to artifactual results. Here, we have focused on common strategies to measure mRNA half-lives in yeast and determined that commonly used transcription inhibitors thiolutin and 1,10 phenanthroline inhibit TORC1 signaling, PKC signaling, and partially activate HOG signaling. Additionally, 4-thiouracil (4tU), a uracil analog used in mRNA pulse-labeling approaches, modestly induces P-bodies, mRNA-protein granules implicated in storage and decay of nontranslating mRNA. Thiolutin also induces P-bodies, whereas phenanthroline has no effect. Doxycycline, which controls "Tet On/Tet Off" regulatable promoters, shows no impact on the above signaling pathways or P-bodies. In summary, our data argues that broad-acting transcriptional inhibitors are problematic for determining mRNA half-life, particularly if studying the impacts of the TORC1, HOG, or PKC pathway on mRNA stability. Regulatable promoter systems are a preferred approach for individual mRNA half-life studies, with 4tU labeling representing a good approach to global mRNA half-life analysis, despite modestly inducing P-bodies.

LYAR potentiates rRNA synthesis by recruiting BRD2/4 and the MYST-type acetyltransferase KAT7 to rDNA

Activation of ribosomal RNA (rRNA) synthesis is pivotal during cell growth and proliferation, but its aberrant upregulation may promote tumorigenesis. Here, we demonstrate that the candidate oncoprotein, LYAR, enhances ribosomal DNA (rDNA) transcription. Our data reveal that LYAR binds the histone-associated protein BRD2 without involvement of acetyl-lysine-binding bromodomains and recruits BRD2 to the rDNA promoter and transcribed regions via association with upstream binding factor. We show that BRD2 is required for the recruitment of the MYST-type acetyltransferase KAT7 to rDNA loci, resulting in enhanced local acetylation of histone H4. In addition, LYAR binds a complex of BRD4 and KAT7, which is then recruited to rDNA independently of the BRD2-KAT7 complex to accelerate the local acetylation of both H4 and H3. BRD2 also helps recruit BRD4 to rDNA. By contrast, LYAR has no effect on rDNA methylation or the binding of RNA polymerase I subunits to rDNA. These data suggest that LYAR promotes the association of the BRD2-KAT7 and BRD4-KAT7 complexes with transcription-competent rDNA loci but not to transcriptionally silent rDNA loci, thereby increasing rRNA synthesis by altering the local acetylation status of histone H3 and H4.

Methods for the analysis of transcriptome dynamics

The transcriptome is the complete set of transcripts in a cell or tissue and includes ribosomal RNA (rRNA), messenger RNA (mRNA), transfer RNA (tRNA), and regulatory noncoding RNA. At steady-state, the transcriptome results from a compensatory variation of the transcription and decay rate to maintain the RNA concentration constant. RNA transcription constitutes the first stage in gene expression, and thus is a major and primary mode of gene expression control. Nevertheless, regulation of RNA decay is also a key factor in gene expression control, involving either selective RNA stabilization or enhanced degradation. Transcriptome analysis allows the identification of gene expression alterations, providing new insights regarding the pathways and mechanisms involved in physiological and pathological processes. Upon perturbation of cell homeostasis, rapid changes in gene expression are required to adapt to new conditions. Thus, to better understand the regulatory mechanisms associated with gene expression alterations, it is vital to acknowledge the relative contribution of RNA synthesis and decay to the transcriptome. To the toxicology field, the study of gene expression regulation mechanisms can help identify the early and mechanistic relevant cellular events associated with a particular response. This review aims to provide a critical comparison of the available methods used to analyze the contribution of RNA transcription and decay to gene expression dynamics. Notwithstanding, an integration of the data obtained is necessary to understand the entire repercussions of gene transcription changes at a system-level. Thus, a brief overview of the methods available for the integration and analysis of the data obtained from transcriptome analysis will also be provided.

Chemo-enzymatic treatment of RNA to facilitate analyses

Labeling RNA is a recurring problem to make RNA compatible with state-of-the-art methodology and comes in many flavors. Considering only cellular applications, the spectrum still ranges from site-specific labeling of individual transcripts, for example, for live-cell imaging of mRNA trafficking, to metabolic labeling in combination with next generation sequencing to capture dynamic aspects of RNA metabolism on a transcriptome-wide scale. Combining the specificity of RNA-modifying enzymes with non-natural substrates has emerged as a valuable strategy to modify RNA site- or sequence-specifically with functional groups suitable for subsequent bioorthogonal reactions and thus label RNA with reporter moieties such as affinity or fluorescent tags. In this review article, we will cover chemo-enzymatic approaches (a) for in vitro labeling of RNA for application in cells, (b) for treatment of total RNA, and (c) for metabolic labeling of RNA. This article is categorized under: RNA Processing < RNA Editing and Modification RNA Methods < RNA Analyses in vitro and In Silico RNA Methods < RNA Analyses in Cells.

Sequencing cell-type-specific transcriptomes with SLAM-ITseq

Analysis of cell-type-specific transcriptomes is vital for understanding the biology of tissues and organs in the context of multicellular organisms. In this Protocol Extension, we combine a previously developed cell-type-specific metabolic RNA labeling method (thiouracil (TU) tagging) and a pipeline to detect the labeled transcripts by a novel RNA sequencing (RNA-seq) method, SLAMseq (thiol (SH)-linked alkylation for the metabolic sequencing of RNA). By injecting a uracil analog, 4-thiouracil, into transgenic mice that express cell-type-specific uracil phosphoribosyltransferase (UPRT), an enzyme required for 4-thiouracil incorporation into newly synthesized RNA, only cells expressing UPRT synthesize thiol-containing RNA. Total RNA isolated from a tissue of interest is then sequenced with SLAMseq, which introduces thymine to cytosine (T>C) conversions at the sites of the incorporated 4-thiouracil. The resulting sequencing reads are then mapped with the T>C-aware alignment software, SLAM-DUNK, which allows mapping of reads containing T>C mismatches. The number of T>C conversions per transcript is further analyzed to identify which transcripts are synthesized in the UPRT-expressing cells. Thus, our method, SLAM-ITseq (SLAMseq in tissue), enables cell-specific transcriptomics without laborious FACS-based cell sorting or biochemical isolation of the labeled transcripts used in TU tagging. In the murine tissues we assessed previously, this method identified ~5,000 genes that are expressed in a cell type of interest from the total RNA pool from the tissue. Any laboratory with access to a high-throughput sequencer and high-power computing can adapt this protocol with ease, and the entire pipeline can be completed in <5 d.