Metabolic labeling of RNA using multiple ribonucleoside analogs enables the simultaneous evaluation of RNA synthesis and degradation rates

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Missense mutations in Myc Box I influence MYC cellular localization, mRNA partitioning and turnover to promote leukemogenesis

Somatic missense mutations in the phosphodegron domain of the MYC gene ( M YC Box I) are detected in the dominant clones of a subset of acute myeloid leukemia (AML) patients, but the mechanisms by which they contribute to AML are unknown. To unveil unique proprieties of MBI MYC mutant proteins, we systematically compared the cellular and molecular consequences of expressing similar oncogenic levels of wild type and MBI mutant MYC. We found that MBI MYC mutants can accelerate leukemia by driving unique transcriptional signatures in highly selected, myeloid progenitor subpopulations. Although these mutations increase MYC stability, they overall dampen MYC chromatin localization and lead to a cytoplasmic accumulation of the mutant proteins. This phenotype is coupled with increased translation of RNA binding proteins and nuclear export machinery, which results in altered RNA partitioning and accelerated decay of select transcripts encoding proapoptotic and proinflammatory genes. Heterozygous knockin mice harboring the germline MBI mutation Myc p.T73N exhibit cytoplasmic MYC localization, myeloid progenitors' expansion with similar transcriptional signatures to the overexpression model, and eventually develop hematological malignancies. This study uncovers that MBI MYC mutations alter MYC localization and disrupt mRNA subcellular distribution and turnover of select transcripts to accelerate tumor initiation and growth.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Simultaneous measurement of nascent transcriptome and translatome using 4-thiouridine metabolic RNA labeling and translating ribosome affinity purification

Regulation of gene expression in response to various biological processes, including extracellular stimulation and environmental adaptation requires nascent RNA synthesis and translation. Analysis of the coordinated regulation of dynamic RNA synthesis and translation is required to determine functional protein production. However, reliable methods for the simultaneous measurement of nascent RNA synthesis and translation at the gene level are limited. Here, we developed a novel method for the simultaneous assessment of nascent RNA synthesis and translation by combining 4-thiouridine (4sU) metabolic RNA labeling and translating ribosome affinity purification (TRAP) using a monoclonal antibody against evolutionarily conserved ribosomal P-stalk proteins. The P-stalk-mediated TRAP (P-TRAP) technique recovered endogenous translating ribosomes, allowing easy translatome analysis of various eukaryotes. We validated this method in mammalian cells by demonstrating that acute unfolded protein response (UPR) in the endoplasmic reticulum (ER) induces dynamic reprogramming of nascent RNA synthesis and translation. Our nascent P-TRAP (nP-TRAP) method may serve as a simple and powerful tool for analyzing the coordinated regulation of transcription and translation of individual genes in various eukaryotes.

PerturbSci-Kinetics: Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

Here we described PerturbSci-Kinetics, a novel combinatorial indexing method for capturing three-layer single-cell readout (i.e., whole transcriptomes, nascent transcriptomes, sgRNA identities) across hundreds of genetic perturbations. Through PerturbSci-Kinetics profiling of pooled CRISPR screens targeting a variety of biological processes, we were able to decipher the complexity of RNA regulations at multiple levels (e.g., synthesis, processing, degradation), and revealed key regulators involved in miRNA and mitochondrial RNA processing pathways. Our technique opens the possibility of systematically decoding the genome-wide regulatory network underlying RNA temporal dynamics at scale and cost-effectively.

Single-molecule imaging reveals translation-dependent destabilization of mRNAs

The relationship between mRNA translation and decay is incompletely understood, with conflicting reports suggesting that translation can either promote decay or stabilize mRNAs. The effect of translation on mRNA decay has mainly been studied using ensemble measurements and global transcription and translation inhibitors, which can have pleiotropic effects. We developed a single-molecule imaging approach to control the translation of a specific transcript that enabled simultaneous measurement of translation and mRNA decay. Our results demonstrate that mRNA translation reduces mRNA stability, and mathematical modeling suggests that this process is dependent on ribosome flux. Furthermore, our results indicate that miRNAs mediate efficient degradation of both translating and non-translating target mRNAs and reveal a predominant role for mRNA degradation in miRNA-mediated regulation. Simultaneous observation of translation and decay of single mRNAs provides a framework to directly study how these processes are interconnected in cells.

The genetic and biochemical determinants of mRNA degradation rates in mammals

BACKGROUND

Degradation rate is a fundamental aspect of mRNA metabolism, and the factors governing it remain poorly characterized. Understanding the genetic and biochemical determinants of mRNA half-life would enable more precise identification of variants that perturb gene expression through post-transcriptional gene regulatory mechanisms.

RESULTS

We establish a compendium of 39 human and 27 mouse transcriptome-wide mRNA decay rate datasets. A meta-analysis of these data identified a prevalence of technical noise and measurement bias, induced partially by the underlying experimental strategy. Correcting for these biases allowed us to derive more precise, consensus measurements of half-life which exhibit enhanced consistency between species. We trained substantially improved statistical models based upon genetic and biochemical features to better predict half-life and characterize the factors molding it. Our state-of-the-art model, Saluki, is a hybrid convolutional and recurrent deep neural network which relies only upon an mRNA sequence annotated with coding frame and splice sites to predict half-life (r=0.77). The key novel principle learned by Saluki is that the spatial positioning of splice sites, codons, and RNA-binding motifs within an mRNA is strongly associated with mRNA half-life. Saluki predicts the impact of RNA sequences and genetic mutations therein on mRNA stability, in agreement with functional measurements derived from massively parallel reporter assays.

CONCLUSIONS

Our work produces a more robust ground truth for transcriptome-wide mRNA half-lives in mammalian cells. Using these revised measurements, we trained Saluki, a model that is over 50% more accurate in predicting half-life from sequence than existing models. Saluki succinctly captures many of the known determinants of mRNA half-life and can be rapidly deployed to predict the functional consequences of arbitrary mutations in the transcriptome.

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.

IL12RB1 allele bias in human TH cells is regulated by functional SNPs in its 3'UTR

Allele bias is an epigenetic mechanism wherein only the maternal- or paternal-derived allele of a gene is preferentially expressed. Allele bias is used by T cells to regulate expression of numerous genes, including those which govern their development and response to cytokines. Here we demonstrate that human T_H cell expression of the cytokine receptor gene IL12RB1 is subject to allele bias, and the extent to which this bias occurs is influenced by cells' differentiation status and two polymorphic sites in the IL12RB1 3'UTR. The single nucleotide polymorphisms (SNPs) at these sites, rs3746190 and rs404733, function to increase expression of their encoding allele. Modeling suggests this is due to a stabilizing effect of these SNPs on the predicted mRNA secondary structure. The SNP rs3746190 is also proximal to the predicted binding site of microRNA miR-1277, raising the possibility that miR-1277 cannot exert suppression in the presence of rs3746190. Functional experiments demonstrate, however, that miR-1277 suppression of IL12RB1 3'UTR expression-which itself has not been previously reported-is nevertheless independent of rs3746190. Collectively, these data demonstrate that rs3746190 and rs404733 are functional SNPs which regulate IL12RB1 allele bias in human T_H cells.

Eukaryotic mRNA Decapping Activation

The 5′-terminal cap is a fundamental determinant of eukaryotic gene expression which facilitates cap-dependent translation and protects mRNAs from exonucleolytic degradation. Enzyme-directed hydrolysis of the cap (decapping) decisively affects mRNA expression and turnover, and is a heavily regulated event. Following the identification of the decapping holoenzyme (Dcp1/2) over two decades ago, numerous studies revealed the complexity of decapping regulation across species and cell types. A conserved set of Dcp1/2-associated proteins, implicated in decapping activation and molecular scaffolding, were identified through genetic and molecular interaction studies, and yet their exact mechanisms of action are only emerging. In this review, we discuss the prevailing models on the roles and assembly of decapping co-factors, with considerations of conservation across species and comparison across physiological contexts. We next discuss the functional convergences of decapping machineries with other RNA-protein complexes in cytoplasmic P bodies and compare current views on their impact on mRNA stability and translation. Lastly, we review the current models of decapping activation and highlight important gaps in our current understanding.

Direct detection of RNA modifications and structure using single-molecule nanopore sequencing

Modifications are present on many classes of RNA, including tRNA, rRNA, and mRNA. These modifications modulate diverse biological processes such as genetic recoding and mRNA export and folding. In addition, modifications can be introduced to RNA molecules using chemical probing strategies that reveal RNA structure and dynamics. Many methods exist to detect RNA modifications by short-read sequencing; however, limitations on read length inherent to short-read-based methods dissociate modifications from their native context, preventing single-molecule modification analysis. Here, we demonstrate direct RNA nanopore sequencing to detect endogenous and exogenous RNA modifications on long RNAs at the single-molecule level. We detect endogenous 2'-O-methyl and base modifications across E. coli and S. cerevisiae ribosomal RNAs as shifts in current signal and dwell times distally through interactions with the helicase motor protein. We further use the 2'-hydroxyl reactive SHAPE reagent acetylimidazole to probe RNA structure at the single-molecule level with readout by direct nanopore sequencing.

Stephenson et al. employ direct RNA nanopore sequencing to detect endogenous and exogenous modifications on single RNA molecules. The authors demonstrate detection of endogenous 2'-O-methylation (Nm) on native ribosomal RNAs, confirming known modification patterns. They describe the development of nanoSHAPE, a method that involves exogenously labeling RNA with a small-adduct-generating chemical probe that can reveal RNA structure using long-read sequencing.

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.

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.

Dynamics of transcriptional and post-transcriptional regulation

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