Metabolic labeling and recovery of nascent RNA to accurately quantify mRNA stability

A rapid inducible RNA decay system reveals fast mRNA decay in P-bodies

RNA decay is vital for regulating mRNA abundance and gene expression. Existing technologies lack the spatiotemporal precision or transcript specificity to capture the stochastic and transient decay process. We devise a general strategy to inducibly recruit protein factors to modulate target RNA metabolism. Specifically, we introduce a Rapid Inducible Decay of RNA (RIDR) technology to degrade target mRNAs within minutes. The fast and synchronous induction enables direct visualization of mRNA decay dynamics in cells. Applying RIDR to endogenous ACTB mRNA reveals rapid formation and dissolution of RNA granules in pre-existing P-bodies. Time-resolved RNA distribution measurements demonstrate rapid RNA decay inside P-bodies, which is further supported by knocking down P-body constituent proteins. Light and oxidative stress modulate P-body behavior, potentially reconciling the contradictory literature about P-body function. This study reveals compartmentalized RNA decay kinetics, establishing RIDR as a pivotal tool for exploring the spatiotemporal RNA metabolism in cells.

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.

Coordination of rhythmic RNA synthesis and degradation orchestrates 24- and 12-h RNA expression patterns in mouse fibroblasts

Circadian RNA expression is essential to ultimately regulate a plethora of downstream rhythmic biochemical, physiological, and behavioral processes. Both transcriptional and posttranscriptional mechanisms are considered important to drive rhythmic RNA expression; however, the extent to which each regulatory process contributes to the rhythmic RNA expression remains controversial. To systematically address this, we monitored RNA dynamics using metabolic RNA labeling technology during a circadian cycle in mouse fibroblasts. We find that rhythmic RNA synthesis is the primary contributor of 24-h RNA rhythms, while rhythmic degradation is more important for 12-h RNA rhythms. These rhythms were predominantly regulated by Bmal1 and/or the core clock mechanism, and the interplay between rhythmic synthesis and degradation has a significant impact in shaping rhythmic RNA expression patterns. Interestingly, core clock RNAs are regulated by multiple rhythmic processes and have the highest amplitude of synthesis and degradation, presumably critical to sustain robust rhythmicity of cell-autonomous circadian rhythms. Our study yields invaluable insights into the temporal dynamics of both 24- and 12-h RNA rhythms in mouse fibroblasts.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

Quantification of the dynamics of RNA metabolism is essential for understanding gene regulation in health and disease. Existing methods rely on metabolic labeling of nascent RNAs and physical separation or inference of labeling through PCR-generated mutations, followed by short-read sequencing. However, these methods are limited in their ability to identify transient decay intermediates or co-analyze RNA decay with cis-regulatory elements of RNA stability such as poly(A) tail length and modification status, at single molecule resolution. Here we use 5-ethynyl uridine (5EU) to label nascent RNA followed by direct RNA sequencing with nanopores. We developed RNAkinet, a deep convolutional and recurrent neural network that processes the electrical signal produced by nanopore sequencing to identify 5EU-labeled nascent RNA molecules. RNAkinet demonstrates generalizability to distinct cell types and organisms and reproducibly quantifies RNA kinetic parameters allowing the combined interrogation of RNA metabolism and cis-acting RNA regulatory elements.

Coordination of rhythmic RNA synthesis and degradation orchestrates 24-hour and 12-hour RNA expression patterns in mouse fibroblasts

Circadian RNA expression is essential to ultimately regulate a plethora of downstream rhythmic biochemical, physiological, and behavioral processes. Both transcriptional and post-transcriptional mechanisms are considered important to drive rhythmic RNA expression, however, the extent to which each regulatory process contributes to the rhythmic RNA expression remains controversial. To systematically address this, we monitored RNA dynamics using metabolic RNA labeling technology during a circadian cycle in mouse fibroblasts. We find that rhythmic RNA synthesis is the primary contributor of 24 hr RNA rhythms, while rhythmic degradation is more important for 12 hr RNA rhythms. These rhythms were predominantly regulated by Bmal1 and/or the core clock mechanism, and interplay between rhythmic synthesis and degradation has a significant impact in shaping rhythmic RNA expression patterns. Interestingly, core clock RNAs are regulated by multiple rhythmic processes and have the highest amplitude of synthesis and degradation, presumably critical to sustain robust rhythmicity of cell-autonomous circadian rhythms. Our study yields invaluable insights into the temporal dynamics of both 24 hr and 12 hr RNA rhythms in mouse fibroblasts.

A Rapid Inducible RNA Decay system reveals fast mRNA decay in P-bodies

RNA decay plays a crucial role in regulating mRNA abundance and gene expression. Modulation of RNA degradation is imperative to investigate an RNA's function. However, information regarding where and how RNA decay occurs remains scarce, partially because existing technologies fail to initiate RNA decay with the spatiotemporal precision or transcript specificity required to capture this stochastic and transient process. Here, we devised a general method that employs inducible tethering of regulatory protein factors to target RNAs and modulate their metabolism. Specifically, we established a Rapid Inducible Decay of RNA (RIDR) technology to degrade target mRNA within minutes. The fast and synchronous induction enabled direct visualization of mRNA decay dynamics in cells with spatiotemporal precision previously unattainable. When applying RIDR to endogenous ACTB mRNA, we observed rapid formation and disappearance of RNA granules, which coincided with pre-existing processing bodies (P-bodies). We measured the time-resolved RNA distribution in P-bodies and cytoplasm after induction, and compared different models of P-body function. We determined that mRNAs rapidly decayed in P-bodies upon induction. Additionally, we validated the functional role of P-bodies by knocking down specific a P-body constituent protein and RNA degradation enzyme. This study determined compartmentalized RNA decay kinetics for the first time. Together, RIDR provides a valuable and generalizable tool to study the spatial and temporal RNA metabolism in cells.

SLAM-seq reveals early transcriptomic response mechanisms upon glutamine deprivation in Chinese hamster ovary cells

Mammalian cells frequently encounter subtle perturbations during recombinant protein production. Identifying the genetic factors that govern the cellular stress response can facilitate targeted genetic engineering to obtain production cell lines that demonstrate a higher stress tolerance. To simulate nutrient stress, Chinese hamster ovary (CHO) cells were transferred into a glutamine(Q)-free medium and transcriptional dynamics using thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) along with standard RNA-seq of stressed and unstressed cells were investigated. The SLAM-seq method allows differentiation between actively transcribed, nascent mRNA, and total (previously present) mRNA in the sample, adding an additional, time-resolved layer to classic RNA-sequencing. The cells tackle amino acid (AA) limitation by inducing the integrated stress response (ISR) signaling pathway, reflected in Atf4 overexpression in the early hours post Q deprivation, leading to subsequent activation of its targets, Asns, Atf3, Ddit3, Eif4ebp1, Gpt2, Herpud1, Slc7a1, Slc7a11, Slc38a2, Trib3, and Vegfa. The GCN2-eIF2α-ATF4 pathway is confirmed by a significant halt in transcription of translation-related genes at 24 h post Q deprivation. The downregulation of lipid synthesis indicates the inhibition of the mTOR pathway, further confirmed by overexpression of Sesn2. Furthermore, SLAM-seq detects short-lived transcription factors, such as Egr1, that would have been missed in standard experimental designs with RNA-seq. Our results describe the successful establishment of SLAM-seq in CHO cells and therefore facilitate its future use in other scenarios where dynamic transcriptome profiling in CHO cells is essential.

Computational insights into mRNA and protein dynamics underlying synaptic plasticity rules

Recent advances in experimental techniques provide an unprecedented peek into the intricate molecular dynamics inside synapses and dendrites. The experimental insights into the molecular turnover revealed that such processes as diffusion, active transport, spine uptake, and local protein synthesis could dynamically modulate the copy numbers of plasticity-related molecules in synapses. Subsequently, theoretical models were designed to understand the interaction of these processes better and to explain how local synaptic plasticity cues can up or down-regulate the molecular copy numbers across synapses. In this review, we discuss the recent advances in experimental techniques and computational models to highlight how these complementary approaches can provide insight into molecular cross-talk across synapses, ultimately allowing us to develop biologically-inspired neural network models to understand brain function.

A selective and sensitive detection system for 4-thiouridine modification in RNA

4-Thiouridine (s4U) is a modified nucleoside, found at positions 8 and 9 in tRNA from eubacteria and archaea. Studies of the biosynthetic pathway and physiological role of s4U in tRNA are ongoing in the tRNA modification field. s4U has also recently been utilized as a biotechnological tool for analysis of RNAs. Therefore, a selective and sensitive system for the detection of s4U is essential for progress in the fields of RNA technologies and tRNA modification. Here, we report the use of biotin-coupled 2-aminoethyl-methanethiosulfonate (MTSEA biotin-XX) for labeling of s4U and demonstrate that the system is sensitive and quantitative. This technique can be used without denaturation; however, addition of a denaturation step improves the limit of detection. Thermus thermophilus tRNAs, which abundantly contain 5-methyl-2-thiouridine, were tested to investigate the selectivity of the MTSEA biotin-XX s4U detection system. The system did not react with 5-methyl-2-thiouridine in tRNAs from a T. thermophilus tRNA 4-thiouridine synthetase (thiI) gene deletion strain. Thus, the most useful advantage of the MTSEA biotin-XX s4U detection system is that MTSEA biotin-XX reacts only with s4U and not with other sulfur-containing modified nucleosides such as s2U derivatives in tRNAs. Furthermore, the MTSEA biotin-XX s4U detection system can analyze multiple samples in a short time span. The MTSEA biotin-XX s4U detection system can also be used for the analysis of s4U formation in tRNA. Finally, we demonstrate that the MTSEA biotin-XX system can be used to visualize newly transcribed tRNAs in S. cerevisiae cells.

PUMILIO competes with AUF1 to control DICER1 RNA levels and miRNA processing

DICER1 syndrome is a cancer pre-disposition disorder caused by mutations that disrupt the function of DICER1 in miRNA processing. Studying the molecular, cellular and oncogenic effects of these mutations can reveal novel mechanisms that control cell homeostasis and tumor biology. Here, we conduct the first analysis of pathogenic DICER1 syndrome allele from the DICER1 3'UTR. We find that the DICER1 syndrome allele, rs1252940486, abolishes interaction with the PUMILIO RNA binding protein with the DICER1 3'UTR, resulting in the degradation of the DICER1 mRNA by AUF1. This single mutational event leads to diminished DICER1 mRNA and protein levels, and widespread reprogramming of miRNA networks. The in-depth characterization of the rs1252940486 DICER1 allele, reveals important post-transcriptional regulatory events that control DICER1 levels.

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

Enhancers regulate 3' end processing activity to control expression of alternative 3'UTR isoforms

Multi-UTR genes are widely transcribed and express their alternative 3'UTR isoforms in a cell type-specific manner. As transcriptional enhancers regulate mRNA expression, we investigated if they also regulate 3'UTR isoform expression. Endogenous enhancer deletion of the multi-UTR gene PTEN did not impair transcript production but prevented 3'UTR isoform switching which was recapitulated by silencing of an enhancer-bound transcription factor. In reporter assays, enhancers increase transcript production when paired with single-UTR gene promoters. However, when combined with multi-UTR gene promoters, they change 3'UTR isoform expression by increasing 3' end processing activity of polyadenylation sites. Processing activity of polyadenylation sites is affected by transcription factors, including NF-κB and MYC, transcription elongation factors, chromatin remodelers, and histone acetyltransferases. As endogenous cell type-specific enhancers are associated with genes that increase their short 3'UTRs in a cell type-specific manner, our data suggest that transcriptional enhancers integrate cellular signals to regulate cell type-and condition-specific 3'UTR isoform expression.

Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment

Background

Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample.

Methods

Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data.

Results

Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method.

Conclusions

In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04672-4.

Synthesis of 4-thiouridines with prodrug functionalization for RNA metabolic labeling

Metabolic labeling has emerged as a powerful tool to endow RNA with reactive handles allowing for subsequent chemical derivatization and processing. Recently, thiolated nucleosides, such as 4-thiouridine (4sU), have attracted great interest in metabolic labeling-based RNA sequencing approaches (TUC-seq, SLAM-seq, TimeLapse-seq) to study cellular RNA expression and decay dynamics. For these and other applications (e.g. PAR-CLIP), thus far only the naked nucleoside 4sU has been applied. Here we examined the concept of derivatizing 4sU into a 5′-monophosphate prodrug that would allow for cell permeation and potentially improve labeling efficiency by bypassing the rate-limiting first step of 5′ phosphorylation of the nucleoside into the ultimately bioactive 4sU triphosphate (4sUTP). To this end, we developed robust synthetic routes towards diverse 4sU monophosphate prodrugs. Using metabolic labeling assays, we found that most of the newly introduced 4sU prodrugs were well tolerated by the cells. One derivative, the bis(4-acetyloxybenzyl) 5′-monophosphate of 4sU, was also efficiently incorporated into nascent RNA.

Synthetic access to 4-thiouridine (4sU) derivatives with monophosphate prodrug patterns creates additional possibilities for metabolic labeling of RNA for different applications.

Addition-Elimination Mechanism-Activated Nucleotide Transition Sequencing for RNA Dynamics Profiling

Dynamic information of intracellular transcripts is essential to understand their functional roles. Routine RNA-sequencing (RNA-seq) methods only measure RNA species at a steady state and do not provide RNA dynamic information. Here, we develop addition-elimination mechanism-activated nucleotide transition sequencing (AENT-seq) for transcriptome-wide profiling of RNA dynamics. In AENT-seq, nascent transcripts are metabolically labeled with 4-thiouridine (4sU). The total RNA is treated with N₂H₄·H₂O under aqueous conditions. N₂H₄·H₂O is demonstrated to convert 4sU to 4-hydrazino cytosine (C*) based on an addition-elimination chemistry. C* is regarded as cytosine (C) during the DNA extension process. This 4sU-to-C transition marks nascent transcripts, so it enables sequencing analysis of RNA dynamics. We apply our AENT-seq to investigate transcript dynamic information of several genes involved in cancer progression and metastasis. This method uses a simple chemical reaction in aqueous solutions and will be rapidly disseminated with extensive applications.

Transcriptional Regulation of MECP2E1-E2 Isoforms and BDNF by Metformin and Simvastatin through Analyzing Nascent RNA Synthesis in a Human Brain Cell Line

Methyl CpG binding protein 2 (MeCP2) is the main DNA methyl-binding protein in the brain that binds to 5-methylcytosine and 5-hydroxymethyl cytosine. MECP2 gene mutations are the main origin of Rett Syndrome (RTT), a neurodevelopmental disorder in young females. The disease has no existing cure, however, metabolic drugs such as metformin and statins have recently emerged as potential therapeutic candidates. In addition, induced MECP2-BDNF homeostasis regulation has been suggested as a therapy avenue. Here, we analyzed nascent RNA synthesis versus steady state total cellular RNA to study the transcriptional effects of metformin (an anti-diabetic drug) on MECP2 isoforms (E1 and E2) and BNDF in a human brain cell line. Additionally, we investigated the impact of simvastatin (a cholesterol lowering drug) on transcriptional regulation of MECP2E1/E2-BDNF. Metformin was capable of post-transcriptionally inducing BDNF and/or MECP2E1, while transcriptionally inhibiting MECP2E2. In contrast simvastatin significantly inhibited BDNF transcription without significantly impacting MECP2E2 transcripts. Further analysis of ribosomal RNA transcripts confirmed that the drug neither individually nor in combination affected these fundamentally important transcripts. Experimental analysis was completed in conditions of the presence or absence of serum starvation that showed minimal impact for serum deprival, although significant inhibition of steady state MECP2E1 by simvastatin was only detected in non-serum starved cells. Taken together, our results suggest that metformin controls MECP2E1/E2-BDNF transcriptionally and/or post-transcriptionally, and that simvastatin is a potent transcriptional inhibitor of BDNF. The transcriptional effect of these drugs on MECP2E1/E2-BDNF were not additive under these tested conditions, however, either drug may have potential application for related disorders.

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.

Endogenous p53 expression in human and mouse is not regulated by its 3'UTR

The TP53 gene encodes the tumor suppressor p53 which is functionally inactivated in many human cancers. Numerous studies suggested that 3′UTR-mediated p53 expression regulation plays a role in tumorigenesis and could be exploited for therapeutic purposes. However, these studies did not investigate post-transcriptional regulation of the native TP53 gene. Here, we used CRISPR/Cas9 to delete the human and mouse TP53/Trp53 3′UTRs while preserving endogenous mRNA processing. This revealed that the endogenous 3′UTR is not involved in regulating p53 mRNA or protein expression neither in steady state nor after genotoxic stress. Using reporter assays, we confirmed the previously observed repressive effects of the isolated 3′UTR. However, addition of the TP53 coding region to the reporter had a dominant negative impact on expression as its repressive effect was stronger and abrogated the contribution of the 3′UTR. Our data highlight the importance of genetic models in the validation of post-transcriptional gene regulatory effects.

A Trans-Omics Comparison Reveals Common Gene Expression Strategies in Four Model Organisms and Exposes Similarities and Differences between Them

The ultimate goal of gene expression regulation is on the protein level. However, because the amounts of mRNAs and proteins are controlled by their synthesis and degradation rates, the cellular amount of a given protein can be attained by following different strategies. By studying omics data for six expression variables (mRNA and protein amounts, plus their synthesis and decay rates), we previously demonstrated the existence of common expression strategies (CESs) for functionally related genes in the yeast Saccharomyces cerevisiae. Here we extend that study to two other eukaryotes: the yeast Schizosaccharomyces pombe and cultured human HeLa cells. We also use genomic data from the model prokaryote Escherichia coli as an external reference. We show that six-variable profiles (6VPs) can be constructed for every gene and that these 6VPs are similar for genes with similar functions in all the studied organisms. The differences in 6VPs between organisms can be used to establish their phylogenetic relationships. The analysis of the correlations among the six variables supports the hypothesis that most gene expression control occurs in actively growing organisms at the transcription rate level, and that translation plays a minor role. We propose that living organisms use CESs for the genes acting on the same physiological pathways, especially for those belonging to stable macromolecular complexes, but CESs have been modeled by evolution to adapt to the specific life circumstances of each organism.

Comparative proteomics identifies Schlafen 5 (SLFN5) as a herpes simplex virus restriction factor that suppresses viral transcription

Intrinsic antiviral host factors confer cellular defence by limiting virus replication and are often counteracted by viral countermeasures. We reasoned that host factors that inhibit viral gene expression could be identified by determining proteins bound to viral DNA (vDNA) in the absence of key viral antagonists. Herpes simplex virus 1 (HSV-1) expresses E3 ubiquitin-protein ligase ICP0 (ICP0), which functions as an E3 ubiquitin ligase required to promote infection. Cellular substrates of ICP0 have been discovered as host barriers to infection but the mechanisms for inhibition of viral gene expression are not fully understood. To identify restriction factors antagonized by ICP0, we compared proteomes associated with vDNA during HSV-1 infection with wild-type virus and a mutant lacking functional ICP0 (ΔICP0). We identified the cellular protein Schlafen family member 5 (SLFN5) as an ICP0 target that binds vDNA during HSV-1 ΔICP0 infection. We demonstrated that ICP0 mediates ubiquitination of SLFN5, which leads to its proteasomal degradation. In the absence of ICP0, SLFN5 binds vDNA to repress HSV-1 transcription by limiting accessibility of RNA polymerase II to viral promoters. These results highlight how comparative proteomics of proteins associated with viral genomes can identify host restriction factors and reveal that viral countermeasures can overcome SLFN antiviral activity.

Direct RNA sequencing reveals m6A modifications on adenovirus RNA are necessary for efficient splicing

Adenovirus is a nuclear replicating DNA virus reliant on host RNA processing machinery. Processing and metabolism of cellular RNAs can be regulated by METTL3, which catalyzes the addition of N6-methyladenosine (m6A) to mRNAs. While m6A-modified adenoviral RNAs have been previously detected, the location and function of this mark within the infectious cycle is unknown. Since the complex adenovirus transcriptome includes overlapping spliced units that would impede accurate m6A mapping using short-read sequencing, here we profile m6A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m6A detection. Although both early and late viral transcripts contain m6A, depletion of m6A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m6A discovery within individual transcripts at nucleotide resolution, and highlight the role of m6A in regulating splicing of a viral pathogen.

Adenovirus-mediated ubiquitination alters protein-RNA binding and aids viral RNA processing

Viruses promote infection by hijacking the ubiquitin machinery of the host to counteract or redirect cellular processes. Adenovirus encodes two early proteins, E1B55K and E4orf6, that together co-opt a cellular ubiquitin ligase complex to overcome host defences and promote virus production. Adenovirus mutants lacking E1B55K or E4orf6 display defects in viral RNA processing and protein production, but previously identified substrates of the redirected ligase do not explain these phenotypes. Here, we used a quantitative proteomics approach to identify substrates of E1B55K/E4orf6-mediated ubiquitination that facilitate RNA processing. While all currently known cellular substrates of E1B55K and E4orf6 are degraded by the proteasome, we uncovered RNA-binding proteins as high-confidence substrates that are not decreased in overall abundance. We focused on two RNA-binding proteins, RALY and hnRNP-C, which we confirm are ubiquitinated without degradation. Knockdown of RALY and hnRNP-C increased levels of viral RNA splicing, protein abundance and progeny production during infection with E1B55K-deleted virus. Furthermore, infection with E1B55K-deleted virus resulted in an increased interaction of hnRNP-C with viral RNA and attenuation of viral RNA processing. These data suggest that viral-mediated ubiquitination of RALY and hnRNP-C relieves a restriction on viral RNA processing and reveal an unexpected role for non-degradative ubiquitination in the manipulation of cellular processes during virus infection.

Activation and inhibition of nonsense-mediated mRNA decay control the abundance of alternative polyadenylation products

Alternative polyadenylation (APA) produces transcript 3' untranslated regions (3'UTRs) with distinct sequences, lengths, stabilities and functions. We show here that APA products include a class of cryptic nonsense-mediated mRNA decay (NMD) substrates with extended 3'UTRs that gene- or transcript-level analyses of NMD often fail to detect. Transcriptome-wide, the core NMD factor UPF1 preferentially recognizes long 3'UTR products of APA, leading to their systematic downregulation. Counteracting this mechanism, the multifunctional RNA-binding protein PTBP1 regulates the balance of short and long 3'UTR isoforms by inhibiting NMD, in addition to its previously described modulation of co-transcriptional polyadenylation (polyA) site choice. Further, we find that many transcripts with altered APA isoform abundance across multiple tumor types are controlled by NMD. Together, our findings reveal a widespread role for NMD in shaping the outcomes of APA.

Massively parallel and time-resolved RNA sequencing in single cells with scNT-seq

Single-cell RNA sequencing offers snapshots of whole transcriptomes but obscures the temporal RNA dynamics. Here we present single-cell metabolically labeled new RNA tagging sequencing (scNT-seq), a method for massively parallel analysis of newly transcribed and pre-existing mRNAs from the same cell. This droplet microfluidics-based method enables high-throughput chemical conversion on barcoded beads, efficiently marking newly transcribed mRNAs with T-to-C substitutions. Using scNT-seq, we jointly profiled new and old transcriptomes in ~55,000 single cells. These data revealed time-resolved transcription factor activities and cell-state trajectories at the single-cell level in response to neuronal activation. We further determined rates of RNA biogenesis and decay to uncover RNA regulatory strategies during stepwise conversion between pluripotent and rare totipotent two-cell embryo (2C)-like stem cell states. Finally, integrating scNT-seq with genetic perturbation identifies DNA methylcytosine dioxygenase as an epigenetic barrier into the 2C-like cell state. Time-resolved single-cell transcriptomic analysis thus opens new lines of inquiry regarding cell-type-specific RNA regulatory mechanisms.

YTHDF2 destabilizes m6A-modified neural-specific RNAs to restrain differentiation in induced pluripotent stem cells

N⁶-methyladenosine (m⁶A) is an abundant post-transcriptional modification that can impact RNA fate via interactions with m⁶A-specific RNA binding proteins. Despite accumulating evidence that m⁶A plays an important role in modulating pluripotency, the influence of m⁶A reader proteins in pluripotency is less clear. Here, we report that YTHDF2, an m⁶A reader associated with mRNA degradation, is highly expressed in induced pluripotent stem cells (iPSCs) and down-regulated during neural differentiation. Through RNA sequencing, we identified a group of m⁶A-modified transcripts associated with neural development that are directly regulated by YTDHF2. Depletion of YTHDF2 in iPSCs leads to stabilization of these transcripts, loss of pluripotency, and induction of neural-specific gene expression. Collectively, our results suggest YTHDF2 functions to restrain expression of neural-specific mRNAs in iPSCs and facilitate their rapid and coordinated up-regulation during neural induction. These effects are both achieved by destabilization of the targeted transcripts.

Global Profiling of Cellular Substrates of Human Dcp2

Decapping is the first committed step in 5'-to-3' RNA decay, and in the cytoplasm of human cells, multiple decapping enzymes regulate the stabilities of distinct subsets of cellular transcripts. However, the complete set of RNAs regulated by any individual decapping enzyme remains incompletely mapped, and no consensus sequence or property is currently known to unambiguously predict decapping enzyme substrates. Dcp2 was the first-identified and best-studied eukaryotic decapping enzyme, but it has been shown to regulate the stability of <400 transcripts in mammalian cells to date. Here, we globally profile changes in the stability of the human transcriptome in Dcp2 knockout cells via TimeLapse-seq. We find that P-body enrichment is the strongest correlate of Dcp2-dependent decay and that modification with m⁶A exhibits an additive effect with P-body enrichment for Dcp2 targeting. These results are consistent with a model in which P-bodies represent sites where translationally repressed transcripts are sorted for decay by soluble cytoplasmic decay complexes through additional molecular marks.

Immunoglobulin expression and the humoral immune response is regulated by the non-canonical poly(A) polymerase TENT5C

TENT5C is a non-canonical cytoplasmic poly(A) polymerase highly expressed by activated B cells to suppress their proliferation. Here we measure the global distribution of poly(A) tail lengths in responsive B cells using a Nanopore direct RNA-sequencing approach, showing that TENT5C polyadenylates immunoglobulin mRNAs regulating their half-life and consequently steady-state levels. TENT5C is upregulated in differentiating plasma cells by innate signaling. Compared with wild-type, Tent5c^−/− mice produce fewer antibodies and have diminished T-cell-independent immune response despite having more CD138^high plasma cells as a consequence of accelerated differentiation. B cells from Tent5c^−/− mice also have impaired capacity of the secretory pathway, with reduced ER volume and unfolded protein response. Importantly, these functions of TENT5C are dependent on its enzymatic activity as catalytic mutation knock-in mice display the same defect as Tent5c^−/−. These findings define the role of the TENT5C enzyme in the humoral immune response.

Regulating polyadenylation is important for mRNA stability, which can in turn affect B cell maturation and humoral immune responses. Here the authors use Nanopore poly(A) sequencing to explore the importance of the cytoplasmic poly(A) polymerase TENT5C, particularly in the production of immunoglobulins.

Thioguanosine Conversion Enables mRNA-Lifetime Evaluation by RNA Sequencing Using Double Metabolic Labeling (TUC-seq DUAL)

Temporal information about cellular RNA populations is essential to understand the functional roles of RNA. We have developed the hydrazine/NH₄ Cl/OsO₄ -based conversion of 6-thioguanosine (6sG) into A', where A' constitutes a 6-hydrazino purine derivative. A' retains the Watson-Crick base-pair mode and is efficiently decoded as adenosine in primer extension assays and in RNA sequencing. Because 6sG is applicable to metabolic labeling of freshly synthesized RNA and because the conversion chemistry is fully compatible with the conversion of the frequently used metabolic label 4-thiouridine (4sU) into C, the combination of both modified nucleosides in dual-labeling setups enables high accuracy measurements of RNA decay. This approach, termed TUC-seq DUAL, uses the two modified nucleosides in subsequent pulses and their simultaneous detection, enabling mRNA-lifetime evaluation with unprecedented precision.

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.

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.

Zika virus noncoding sfRNAs sequester multiple host-derived RNA-binding proteins and modulate mRNA decay and splicing during infection

Insect-borne flaviviruses produce a 300-500-base long noncoding RNA, termed subgenomic flavivirus RNA (sfRNA), by stalling the cellular 5'-3'-exoribonuclease 1 (XRN1) via structures located in their 3' UTRs. In this study, we demonstrate that sfRNA production by Zika virus represses XRN1 analogous to what we have previously shown for other flaviviruses. Using protein-RNA reconstitution and a stringent RNA pulldown assay with human choriocarcinoma (JAR) cells, we demonstrate that the sfRNAs from both dengue type 2 and Zika viruses interact with a common set of 21 RNA-binding proteins that contribute to the regulation of post-transcriptional processes in the cell, including splicing, RNA stability, and translation. We found that four of these sfRNA-interacting host proteins, DEAD-box helicase 6 (DDX6) and enhancer of mRNA decapping 3 (EDC3) (two RNA decay factors), phosphorylated adaptor for RNA export (a regulator of the biogenesis of the splicing machinery), and apolipoprotein B mRNA-editing enzyme catalytic subunit 3C (APOBEC3C, a nucleic acid-editing deaminase), inherently restrict Zika virus infection. Furthermore, we demonstrate that the regulations of cellular mRNA decay and RNA splicing are compromised by Zika virus infection as well as by sfRNA alone. Collectively, these results reveal the large extent to which Zika virus-derived sfRNAs interact with cellular RNA-binding proteins and highlight the potential for widespread dysregulation of post-transcriptional control that likely limits the effective response of these cells to viral infection.

On the optimal design of metabolic RNA labeling experiments

Massively parallel RNA sequencing (RNA-seq) in combination with metabolic labeling has become the de facto standard approach to study alterations in RNA transcription, processing or decay. Regardless of advances in the experimental protocols and techniques, every experimentalist needs to specify the key aspects of experimental design: For example, which protocol should be used (biochemical separation vs. nucleotide conversion) and what is the optimal labeling time? In this work, we provide approximate answers to these questions using the asymptotic theory of optimal design. Specifically, we investigate, how the variance of degradation rate estimates depends on the time and derive the optimal time for any given degradation rate. Subsequently, we show that an increase in sample numbers should be preferred over an increase in sequencing depth. Lastly, we provide some guidance on use cases when laborious biochemical separation outcompetes recent nucleotide conversion based methods (such as SLAMseq) and show, how inefficient conversion influences the precision of estimates. Code and documentation can be found at https://github.com/dieterich-lab/DesignMetabolicRNAlabeling.

Massively parallel RNA sequencing (RNA-seq) in combination with metabolic labeling has become the de facto standard approach to study alterations in RNA transcription, processing or decay. In our manuscript, we address several key aspects of experimental design: 1) The optimal labeling time, 2) the number of replicate samples over sequencing depth and 3) the choice of experimental protocol. We provide approximate answers to these questions using asymptotic theory of optimal design.

hnRNP L-dependent protection of normal mRNAs from NMD subverts quality control in B cell lymphoma

The human nonsense-mediated mRNA decay pathway (NMD) performs quality control and regulatory functions within complex post-transcriptional regulatory networks. In addition to degradation-promoting factors, efficient and accurate detection of NMD substrates involves proteins that safeguard normal mRNAs. Here, we identify hnRNP L as a factor that protects mRNAs with NMD-inducing features including long 3'UTRs. Using biochemical and transcriptome-wide approaches, we provide evidence that the susceptibility of a given transcript to NMD can be modulated by its 3'UTR length and ability to recruit hnRNP L. Integrating these findings with the previously defined role of polypyrimidine tract binding protein 1 in NMD evasion enables enhanced prediction of transcript susceptibility to NMD. Unexpectedly, this system is subverted in B cell lymphomas harboring translocations that produce BCL2:IGH fusion mRNAs. CRISPR/Cas9 deletion of hnRNP L binding sites near the BCL2 stop codon reduces expression of the fusion mRNAs and induces apoptosis. Together, our data indicate that protection by hnRNP L overrides the presence of multiple 3'UTR introns, allowing these aberrant mRNAs to evade NMD and promoting BCL2 overexpression and neoplasia.

Global analysis of RNA metabolism using bio-orthogonal labeling coupled with next-generation RNA sequencing

Many open questions in RNA biology relate to the kinetics of gene expression and the impact of RNA binding regulatory factors on processing or decay rates of particular transcripts. Steady state measurements of RNA abundance obtained from RNA-seq approaches are not able to separate the effects of transcription from those of RNA decay in the overall abundance of any given transcript, instead only giving information on the (presumed steady-state) abundances of transcripts. Through the combination of metabolic labeling and high-throughput sequencing, several groups have been able to measure both transcription rates and decay rates of the entire transcriptome of an organism in a single experiment. This review focuses on the methodology used to specifically measure RNA decay at a global level. By comparing and contrasting approaches and describing the experimental protocols in a modular manner, we intend to provide both experienced and new researchers to the field the ability to combine aspects of various protocols to fit the unique needs of biological questions not addressed by current methods.

Expanding the Nucleoside Recoding Toolkit: Revealing RNA Population Dynamics with 6-Thioguanosine

RNA-sequencing (RNA-seq) measures RNA abundance in a biological sample but does not provide temporal information about the sequenced RNAs. Metabolic labeling can be used to distinguish newly made RNAs from pre-existing RNAs. Mutations induced from chemical recoding of the hydrogen bonding pattern of the metabolic label can reveal which RNAs are new in the context of a sequencing experiment. These nucleotide recoding strategies have been developed for a single uridine analogue, 4-thiouridine (s4U), limiting the scope of these experiments. Here we report the first use of nucleoside recoding with a guanosine analogue, 6-thioguanosine (s6G). Using TimeLapse sequencing (TimeLapse-seq), s6G can be recoded under RNA-friendly oxidative nucleophilic-aromatic substitution conditions to produce adenine analogues (substituted 2-aminoadenosines). We demonstrate the first use of s6G recoding experiments to reveal transcriptome-wide RNA population dynamics.

Gaining insight into transcriptome-wide RNA population dynamics through the chemistry of 4-thiouridine

Cellular RNA levels are the result of a juggling act between RNA transcription, processing, and degradation. By tuning one or more of these parameters, cells can rapidly alter the available pool of transcripts in response to stimuli. While RNA sequencing (RNA-seq) is a vital method to quantify RNA levels genome-wide, it is unable to capture the dynamics of different RNA populations at steady-state or distinguish between different mechanisms that induce changes to the steady-state (i.e., altered rate of transcription vs. degradation). The dynamics of different RNA populations can be studied by targeted incorporation of noncanonical nucleosides. 4-Thiouridine (s⁴ U) is a commonly used and versatile RNA metabolic label that allows the study of many properties of RNA metabolism from synthesis to degradation. Numerous experimental strategies have been developed that leverage the power of s⁴ U to label newly transcribed RNA in whole cells, followed by enrichment with activated disulfides or chemistry to induce C mutations at sites of s⁴ U during sequencing. This review presents existing methods to study RNA population dynamics genome-wide using s⁴ U metabolic labeling, as well as a discussion of considerations and challenges when designing s⁴ U metabolic labeling experiments. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Turnover and Surveillance > Regulation of RNA Stability.

Contributions of regulated transcription and mRNA decay to the dynamics of gene expression

Organisms have acquired sophisticated regulatory networks that control gene expression in response to cellular perturbations. Understanding of the mechanisms underlying the coordinated changes in gene expression in response to external and internal stimuli is a fundamental issue in biology. Recent advances in high-throughput technologies have enabled the measurement of diverse biological information, including gene expression levels, kinetics of gene expression, and interactions among gene expression regulatory molecules. By coupling these technologies with quantitative modeling, we can now uncover the biological roles and mechanisms of gene regulation at the system level. This review consists of two parts. First, we focus on the methods using uridine analogs that measure synthesis and decay rates of RNAs, which demonstrate how cells dynamically change the regulation of gene expression in response to both internal and external cues. Second, we discuss the underlying mechanisms of these changes in kinetics, including the functions of transcription factors and RNA-binding proteins. Overall, this review will help to clarify a system-level view of gene expression programs in cells. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Regulation of RNA Stability RNA Methods > RNA Analyses in vitro and In Silico.

Solid phase chemistry to covalently and reversibly capture thiolated RNA

Here, we describe an approach to enrich newly transcribed RNAs from primary mouse neurons using 4-thiouridine (s4U) metabolic labeling and solid phase chemistry. This one-step enrichment procedure captures s4U-RNA by using highly efficient methane thiosulfonate (MTS) chemistry in an immobilized format. Like solution-based methods, this solid-phase enrichment can distinguish mature RNAs (mRNA) with differential stability, and can be used to reveal transient RNAs such as enhancer RNAs (eRNAs) and primary microRNAs (pri-miRNAs) from short metabolic labeling. Most importantly, the efficiency of this solid-phase chemistry made possible the first large scale measurements of RNA polymerase II (RNAPII) elongation rates in mouse cortical neurons. Thus, our approach provides the means to study regulation of RNA metabolism in specific tissue contexts as a means to better understand gene expression in vivo.

A length-dependent evolutionarily conserved pathway controls nuclear export of circular RNAs

Here, Huang et al. investigated how circRNA localization or nuclear export is controlled and, using RNAi screening, found that depletion of the Drosophila DExH/D-box helicase Hel25E results in nuclear accumulation of long (>800-nt), but not short, circRNAs. Their findings suggest that the lengths of mature circRNAs are measured to dictate the mode of nuclear export.

Circular RNAs (circRNAs) are generated from many protein-coding genes. Most accumulate in the cytoplasm, but how circRNA localization or nuclear export is controlled remains unclear. Using RNAi screening, we found that depletion of the Drosophila DExH/D-box helicase Hel25E results in nuclear accumulation of long (>800-nucleotide), but not short, circRNAs. The human homologs of Hel25E similarly regulate circRNA localization, as depletion of UAP56 (DDX39B) or URH49 (DDX39A) causes long and short circRNAs, respectively, to become enriched in the nucleus. These data suggest that the lengths of mature circRNAs are measured to dictate the mode of nuclear export.

ICE1 promotes the link between splicing and nonsense-mediated mRNA decay

The nonsense-mediated mRNA decay (NMD) pathway detects aberrant transcripts containing premature termination codons (PTCs) and regulates expression of 5–10% of non-aberrant human mRNAs. To date, most proteins involved in NMD have been identified by genetic screens in model organisms; however, the increased complexity of gene expression regulation in human cells suggests that additional proteins may participate in the human NMD pathway. To identify proteins required for NMD, we performed a genome-wide RNAi screen against >21,000 genes. Canonical members of the NMD pathway were highly enriched as top hits in the siRNA screen, along with numerous candidate NMD factors, including the conserved ICE1/KIAA0947 protein. RNAseq studies reveal that depletion of ICE1 globally enhances accumulation and stability of NMD-target mRNAs. Further, our data suggest that ICE1 uses a putative MIF4G domain to interact with exon junction complex (EJC) proteins and promotes the association of the NMD protein UPF3B with the EJC.

The DNA in our cells contains the hereditary information that makes each of us unique. Molecules called mRNAs are copies of this information and are used as templates for making proteins. When a strand of incorrectly copied mRNA, or one including errors from the original DNA template, is recognized, our cells destroy the mRNA to prevent it from producing a damaged protein. Organisms from yeast to humans have evolved a mechanism for finding and destroying faulty mRNAs, called mRNA surveillance. Animals are particularly reliant on mRNA surveillance, as their proteins are often made from cutting and pasting together mRNA from different portions of DNA, in a process known as splicing. Despite being a vital process, we still lack a good understanding of how mRNA surveillance works.

Now, Baird et al. used human kidney cells that produced an error-containing mRNA that could be tracked. To investigate how efficient RNA surveillance is under different conditions, the levels of individual proteins were reduced one at a time. By tracking the amount of faulty mRNA, it was possible to find out if a single protein plays a role in human mRNA surveillance. If the number of faulty mRNAs is high when a protein is reduced, it suggests that this protein may be involved in mRNA surveillance.

Baird et al. screened more than 21,000 proteins, the majority of proteins made in human cells. Many of the proteins that stood out as important in mRNA surveillance were the ones already known to be relevant in yeast and worm cells. But the experiments also identified new proteins that appear to play a role specifically in human RNA surveillance. One of the proteins, ICE1, is essential for the relationship between mRNA splicing and mRNA surveillance. Without ICE1, the mRNA surveillance machinery can no longer find and destroy faulty mRNAs.

Nearly one-third of genetic diseases are caused by mutations that result in faulty mRNAs, which can be detected by mRNA surveillance pathways. Depending on the disease, destroying these error-containing mRNAs can either improve or worsen disease symptoms. A better understanding of the factors that control human RNA surveillance could one day help to develop treatments that affect mRNA surveillance to improve disease outcomes.

TimeLapse-seq: adding a temporal dimension to RNA sequencing through nucleoside recoding

RNA sequencing (RNA-seq) offers a snapshot of cellular RNA populations, but not temporal information about the sequenced RNA. Here we report TimeLapse-seq, which uses oxidative-nucleophilic-aromatic substitution to convert 4-thiouridine into cytidine analogs, yielding apparent U-to-C mutations that mark new transcripts upon sequencing. TimeLapse-seq is a single-molecule approach that is adaptable to many applications and reveals RNA dynamics and induced differential expression concealed in traditional RNA-seq.

Uncovering cell type-specific complexities of gene expression and RNA metabolism by TU-tagging and EC-tagging

Cell type-specific transcription is a key determinant of cell fate and function. An ongoing challenge in biology is to develop robust and stringent biochemical methods to explore gene expression with cell type specificity. This challenge has become even greater as researchers attempt to apply high-throughput RNA analysis methods under in vivo conditions. TU-tagging and EC-tagging are in vivo biosynthetic RNA tagging techniques that allow spatial and temporal specificity in RNA purification. Spatial specificity is achieved through targeted expression of pyrimidine salvage enzymes (uracil phosphoribosyltransferase and cytosine deaminase) and temporal specificity is achieved by controlling exposure to bioorthogonal substrates of these enzymes (4-thiouracil and 5-ethynylcytosine). Tagged RNAs can be purified from total RNA extracted from an animal or tissue and used in transcriptome profiling analyses. In addition to identifying cell type-specific mRNA profiles, these techniques are applicable to noncoding RNAs and can be used to measure RNA transcription and decay. Potential applications of TU-tagging and EC-tagging also include fluorescent RNA imaging and selective definition of RNA-protein interactions. TU-tagging and EC-tagging hold great promise for supporting research at the intersection of RNA biology and developmental biology. This article is categorized under: Technologies > Analysis of the Transcriptome.

Osmium-Mediated Transformation of 4-Thiouridine to Cytidine as Key To Study RNA Dynamics by Sequencing

To understand the functional roles of RNA in the cell, it is essential to elucidate the dynamics of their production, processing and decay. A recent method for assessing mRNA dynamics is metabolic labeling with 4-thiouridine (4sU), followed by thio-selective attachment of affinity tags. Detection of labeled transcripts by affinity purification and hybridization to microarrays or by deep sequencing then reveals RNA expression levels. Here, we present a novel sequencing method (TUC-seq) that eliminates affinity purification and allows for direct assessment of 4sU-labeled RNA. It employs an OsO4 -mediated transformation to convert 4sU into cytosine. We exemplify the utility of the new method for verification of endogenous 4sU in tRNAs and for the detection of pulse-labeled mRNA of seven selected genes in mammalian cells to determine the relative abundance of the new transcripts. The results prove TUC-seq as a straight-forward and highly versatile method for studies of cellular RNA dynamics.