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 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.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification.

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

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.

Unusual Base Pair between Two 2-Thiouridines and Its Implication for Nonenzymatic RNA Copying.

2-Thiouridine (s2U) is a nucleobase modification that confers enhanced efficiency and fidelity both on modern tRNA codon translation and on nonenzymatic and ribozyme-catalyzed RNA copying. We have discovered an unusual base pair between two 2-thiouridines that stabilizes an RNA duplex to a degree that is comparable to that of a native A:U base pair. High-resolution crystal structures indicate similar base-pairing geometry and stacking interactions in duplexes containing s2U:s2U compared to those with U:U pairs. Notably, the C═O···H-N hydrogen bond in the U:U pair is replaced with a C═S···H-N hydrogen bond in the s2U:s2U base pair. The thermodynamic stability of the s2U:s2U base pair suggested that this self-pairing might lead to an increased error frequency during nonenzymatic RNA copying. However, competition experiments show that s2U:s2U base-pairing induces only a low level of misincorporation during nonenzymatic RNA template copying because the correct A:s2U base pair outcompetes the slightly weaker s2U:s2U base pair. In addition, even if an s2U is incorrectly incorporated, the addition of the next base is greatly hindered. This strong stalling effect would further increase the effective fidelity of nonenzymatic RNA copying with s2U. Our findings suggest that s2U may enhance the rate and extent of nonenzymatic copying with only a minimal cost in fidelity.

Internally controlled RNA sequencing comparisons using nucleoside recoding chemistry.

Quantitative comparisons of RNA levels from different samples can lead to new biological understanding if they are able to distinguish biological variation from variable sample preparation. These challenges are pronounced in comparisons that require complex biochemical manipulations (e.g. isolating polysomes to study translation). Here, we present Transcript Regulation Identified by Labeling with Nucleoside Analogues in Cell Culture (TILAC), an internally controlled approach for quantitative comparisons of RNA content. TILAC uses two metabolic labels, 4-thiouridine (s4U) and 6-thioguanosine (s6G), to differentially label RNAs in cells, allowing experimental and control samples to be pooled prior to downstream biochemical manipulations. TILAC leverages nucleoside recoding chemistry to generate characteristic sequencing signatures for each label and uses statistical modeling to compare the abundance of RNA transcripts between samples. We verified the performance of TILAC in transcriptome-scale experiments involving RNA polymerase II inhibition and heat shock. We then applied TILAC to quantify changes in mRNA association with actively translating ribosomes during sodium arsenite stress and discovered a set of transcripts that are translationally upregulated, including MCM2 and DDX5. TILAC is broadly applicable to uncover differences between samples leading to improved biological insights.

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

Gene expression is determined by a balance between RNA synthesis and RNA degradation. To elucidate the underlying regulatory mechanisms and principles of this, simultaneous measurements of RNA synthesis and degradation are required. Here, we report the development of "Dyrec-seq," which uses 4-thiouridine and 5-bromouridine to simultaneously quantify RNA synthesis and degradation rates. Dyrec-seq enabled the quantification of RNA synthesis and degradation rates of 4702 genes in HeLa cells. Functional enrichment analysis showed that the RNA synthesis and degradation rates of genes are actually determined by the genes' biological functions. A comparison of theoretical and experimental analyses revealed that the amount of RNA is determined by the ratio of RNA synthesis to degradation rates, whereas the rapidity of responses to external stimuli is determined only by the degradation rate. This study emphasizes that not only RNA synthesis but also RNA degradation is important in shaping gene expression patterns.

Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling.

The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT-, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT- to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.

Tracking Distinct RNA Populations Using Efficient and Reversible Covalent Chemistry.

We describe a chemical method to label and purify 4-thiouridine (s(4)U)-containing RNA. We demonstrate that methanethiosulfonate (MTS) reagents form disulfide bonds with s(4)U more efficiently than the commonly used HPDP-biotin, leading to higher yields and less biased enrichment. This increase in efficiency allowed us to use s(4)U labeling to study global microRNA (miRNA) turnover in proliferating cultured human cells without perturbing global miRNA levels or the miRNA processing machinery. This improved chemistry will enhance methods that depend on tracking different populations of RNA, such as 4-thiouridine tagging to study tissue-specific transcription and dynamic transcriptome analysis (DTA) to study RNA turnover.

A non-radioactive, improved PAR-CLIP and small RNA cDNA library preparation protocol.

Crosslinking and immunoprecipitation (CLIP) methods are powerful techniques to interrogate direct protein-RNA interactions and dissect posttranscriptional gene regulatory networks. One widely used CLIP variant is photoactivatable ribonucleoside enhanced CLIP (PAR-CLIP) that involves in vivo labeling of nascent RNAs with the photoreactive nucleosides 4-thiouridine (4SU) or 6-thioguanosine (6SG), which can efficiently crosslink to interacting proteins using UVA and UVB light. Crosslinking of 4SU or 6SG to interacting amino acids changes their base-pairing properties and results in characteristic mutations in cDNA libraries prepared for high-throughput sequencing, which can be computationally exploited to remove abundant background from non-crosslinked sequences and help pinpoint RNA binding protein binding sites at nucleotide resolution on a transcriptome-wide scale. Here we present a streamlined protocol for fluorescence-based PAR-CLIP (fPAR-CLIP) that eliminates the need to use radioactivity. It is based on direct ligation of a fluorescently labeled adapter to the 3'end of crosslinked RNA on immobilized ribonucleoproteins, followed by isolation of the adapter-ligated RNA and efficient conversion into cDNA without the previously needed size fractionation on denaturing polyacrylamide gels. These improvements cut the experimentation by half to 2 days and increases sensitivity by 10-100-fold.

Sulfur Availability Impacts Accumulation of the 2-Thiouridine tRNA Modification in Bacillus subtilis.

Posttranscriptional modifications to tRNA are critical elements for the folding and functionality of these adaptor molecules. Sulfur modifications in tRNA are installed by specialized enzymes that act on cognate tRNA substrates at specific locations. Most studied organisms contain a general cysteine desulfurase to mobilize sulfur for the synthesis of S-tRNA and other thio-cofactors. Bacillus subtilis and other Gram-positive bacteria encode multiple cysteine desulfurases that partner with specific sulfur acceptors in the biosynthesis of thio-cofactors. This metabolic layout suggests an alternate mode of regulation in these biosynthetic pathways. In this study, tRNA modifications were exploited as a readout for the functionality of pathways involving cysteine desulfurases. These analyses showed that the relative abundance of 2-thiouridine-modified tRNA (s2U) responds to sulfur availability in the growth medium in a dose-dependent manner. This study found that low sulfur concentrations lead to decreased levels of the s2U cysteine desulfurase YrvO and thiouridylase MnmA, without altering the levels of other cysteine desulfurases, SufS, NifS, and NifZ. Analysis of pathway metabolites that depend on the activity of cysteine desulfurases indicates that sulfur nutrient availability specifically impacts s2U accumulation while having no effect on the levels of other S-modified tRNA or activity levels of Fe-S enzymes. Collectively, these results support a model in which s2U tRNA serves as a marker for sulfur availability in B. subtilis. IMPORTANCE The 2-thiouridine (s2U) tRNA modification is found ubiquitously across all domains of life. YrvO and MnmA, the enzymes involved in this modification, are essential in B. subtilis, confirming the well-established role of s2U in maintaining translational efficiency and, consequently, cellular viability. Herein, we show that in the model Gram-positive organism Bacillus subtilis, the levels of s2U are responsive to sulfur availability. Downregulation of the s2U biosynthetic components leads to lower s2U levels, which may serve as a signal for the slowing of the translational apparatus during cellular nutrient insufficiency. Our findings provide the basis for the identification of a potential bacterial mode of regulation during S-metabolite depletion that may use s2U as a marker of suboptimal metabolic status.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus-application of 4-thiouridine photocrosslinking.

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

An ancient type of MnmA protein is an iron-sulfur cluster-dependent sulfurtransferase for tRNA anticodons.

Transfer RNA (tRNA) is an adaptor molecule indispensable for assigning amino acids to codons on mRNA during protein synthesis. 2-thiouridine (s2U) derivatives in the anticodons (position 34) of tRNAs for glutamate, glutamine, and lysine are post-transcriptional modifications essential for precise and efficient codon recognition in all organisms. s2U34 is introduced either by (i) bacterial MnmA/eukaryote mitochondrial Mtu1 or (ii) eukaryote cytosolic Ncs6/archaeal NcsA, and the latter enzymes possess iron-sulfur (Fe-S) cluster. Here, we report the identification of novel-type MnmA homologs containing three conserved Cys residues, which could support Fe-S cluster binding and catalysis, in a broad range of bacteria, including thermophiles, Cyanobacteria, Mycobacteria, Actinomyces, Clostridium, and Helicobacter Using EPR spectroscopy, we revealed that Thermus thermophilus MnmA (TtMnmA) contains an oxygen-sensitive [4Fe-4S]-type cluster. Efficient in vitro formation of s2U34 in tRNALys and tRNAGln by holo-TtMnmA occurred only under anaerobic conditions. Mutational analysis of TtMnmA suggested that the Fe-S cluster is coordinated by the three conserved Cys residues (Cys105, Cys108, and Cys200), and is essential for its activity. Evolutionary scenarios for the sulfurtransferases, including the Fe-S cluster containing Ncs6/NcsA s2U thiouridylases and several distantly related sulfurtransferases, are proposed.

In Situ Visualizing Nascent RNA by Exploring DNA-Templated Oxidative Amination of 4-Thiouridine.

Tracking and mapping the nascent RNA molecules in cells is essential for deciphering embryonic development and neuronal differentiation. Here, we utilized 4-thiouridine (s4U) as a metabolic tag to label nascent RNA and developed a fluorescence imaging method based on the DNA-templated oxidative amination (DTOA) reaction of s4U. The DTOA reaction occurred between amine-modified DNA and s4U-containing RNA with high sequence specificity and chemical selectivity. Target nascent mRNAs in HeLa cells, including those encoding green fluorescent proteins (GFPs) and endogenous BAG-1, were thus lit up selectively by DTOA-based fluorescence in situ hybridization (DTOA FISH). We believe the DTOA conjugation chemistry shown in this study could be generally applied to investigate the spatial distribution of nascent transcription dynamics in cellular processes.

Robust synthesis of 2'-azido modified RNA from 2'-amino precursors by diazotransfer reaction.

Azides are versatile bioorthogonal reporter moieties that are commonly used for site-specific labeling and functionalization of RNA to probe its biology. The preparation of azido modified nucleic acids by solid-phase synthesis is problematic due to the inherent reactivity of P(III) species with azides according to the Staudinger reaction. Various strategies have been developed to bypass this limitation and are often time-consuming, low-yielding and labor-intensive. In particular, the synthesis of RNA with internal 2'-azido modifications is restricted to a single approach that employs P(V) chemistry instead of the widely used P(III) phosphoramidite chemistry. To fill this methodological gap, we present a novel convenient path toward 2'-azido RNA from readily accessible 2'-amino RNA through treatment with the diazotizing reagent fluorosulfuryl azide (FSO2N3). A diazotransfer reaction was established for oligoribonucleotides of different lengths and secondary structures. The robustness of the approach was further demonstrated for RNAs containing multiple 2'-azido moieties and for RNAs containing other sensitive modifications such as thiouridine or methylated nucleobases with a positive charge. The synthetic ease of generating 2'-azido RNA will pave the way for biotechnological applications, in particular for siRNA technologies and for referencing the growing number of RNA metabolic labeling approaches that rely on 2'-azido nucleosides.

mRNA regions where 80S ribosomes pause during translation elongation in vivo interact with protein uS19, a component of the decoding site.

In eukaryotic ribosomes, the conserved protein uS19, formerly known as S15, extends with its C-terminal tail to the decoding site. The cross-linking of uS19 to the A site codon has been detected using synthetic mRNAs bearing 4-thiouridine (s4U) residues. Here, we showed that the A-site tRNA prevents this cross-linking and that the P site codon does not contact uS19. Next, we focused on determining uS19-mRNA interactions in vivo by applying the photoactivatable-ribonucleoside enhancing cross-linking and immunoprecipitation method to a stable HEK293 cell line producing FLAG-tagged uS19 and grown in a medium containing s4U. We found that when translation was stopped by cycloheximide, uS19 was efficiently cross-linked to mRNA regions with a high frequency of Glu, Lys and, more rarely, Arg codons. The results indicate that the complexes, in which the A site codon is not involved in the formation of the mRNA-tRNA duplex, are present among the cycloheximide-arrested 80S complexes, which implies pausing of elongating ribosomes at the above mRNA regions. Thus, our findings demonstrate that the human ribosomal protein uS19 interacts with mRNAs during translation elongation and highlight the regions of mRNAs where ribosome pausing occurs, bringing new structural and functional insights into eukaryotic translation in vivo.

The RNA degradosome promotes tRNA quality control through clearance of hypomodified tRNA.

The factors and mechanisms that govern tRNA stability in bacteria are not well understood. Here, we investigated the influence of posttranscriptional modification of bacterial tRNAs (tRNA modification) on tRNA stability. We focused on ThiI-generated 4-thiouridine (s4U), a modification found in bacterial and archaeal tRNAs. Comprehensive quantification of Vibrio cholerae tRNAs revealed that the abundance of some tRNAs is decreased in a ΔthiI strain in a stationary phase-specific manner. Multiple mechanisms, including rapid degradation of a subset of hypomodified tRNAs, account for the reduced abundance of tRNAs in the absence of thiI Through transposon insertion sequencing, we identified additional tRNA modifications that promote tRNA stability and bacterial viability. Genetic analysis of suppressor mutants as well as biochemical analyses revealed that rapid degradation of hypomodified tRNA is mediated by the RNA degradosome. Elongation factor Tu seems to compete with the RNA degradosome, protecting aminoacyl tRNAs from decay. Together, our observations describe a previously unrecognized bacterial tRNA quality control system in which hypomodification sensitizes tRNAs to decay mediated by the RNA degradosome.

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 N2H4·H2O under aqueous conditions. N2H4·H2O 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.

4-Thiouridine-Enhanced Peroxidase-Generated Biotinylation of RNA.

Peroxidase-generated proximity labeling is in widespread use to study subcellular proteomes and the protein interaction networks in living cells, but the development of subcellular RNA labeling is limited. APEX-seq has emerged as a new method to study subcellular RNA in living cells, but the labeling of RNA still has room to improve. In this work, we describe 4-thiouridine (s4 U)-enhanced peroxidase-generated biotinylation of RNA with high efficiency. The incorporation of s4 U could introduce additional sites for RNA labeling, enhanced biotinylation was observed on monomer, model oligo RNA and total RNA. Through the s4 U metabolic approach, the in vivo RNA biotinylation efficiency by peroxidase catalysis was also dramatically increased, which will benefit RNA isolation and study for the spatial transcriptome.

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.