Time-Resolved Small RNA Sequencing Unravels the Molecular Principles of MicroRNA Homeostasis.

Argonaute-bound microRNAs silence mRNA expression in a dynamic and regulated manner to control organismal development, physiology, and disease. We employed metabolic small RNA sequencing for a comprehensive view on intracellular microRNA kinetics in Drosophila. Based on absolute rate of biogenesis and decay, microRNAs rank among the fastest produced and longest-lived cellular transcripts, disposing up to 105 copies per cell at steady-state. Mature microRNAs are produced within minutes, revealing tight intracellular coupling of biogenesis that is selectively disrupted by pre-miRNA-uridylation. Control over Argonaute protein homeostasis generates a kinetic bottleneck that cooperates with non-coding RNA surveillance to ensure faithful microRNA loading. Finally, regulated small RNA decay enables the selective rapid turnover of Ago1-bound microRNAs, but not of Ago2-bound small interfering RNAs (siRNAs), reflecting key differences in the robustness of small RNA silencing pathways. Time-resolved small RNA sequencing opens new experimental avenues to deconvolute the timescales, molecular features, and regulation of small RNA silencing pathways in living cells.

Approaching the Golden Fleece a Molecule at a Time: Biophysical Insights into Argonaute-Instructed Nucleic Acid Interactions.

Argonaute proteins act at the core of nucleic acid-guided interference pathways that regulate gene expression and defend organisms against foreign genetic elements in all domains of life. Here, we review recent biophysical studies on how Argonaute proteins instruct oligonucleotides in the process of target finding, binding, cleavage, and release, as measured at high spatiotemporal resolution by single-molecule approaches. In the context of previous structural, biochemical, and computational studies, a model emerges for how Argonaute proteins manipulate the thermodynamic rules for nucleic acid hybridization to convey efficiency and specificity to RNA- and DNA-guided regulatory processes.

Uridylation of RNA Hairpins by Tailor Confines the Emergence of MicroRNAs in Drosophila.

Uridylation of RNA species represents an emerging theme in post-transcriptional gene regulation. In the microRNA pathway, such modifications regulate small RNA biogenesis and stability in plants, worms, and mammals. Here, we report Tailor, an uridylyltransferase that is required for the majority of 3' end modifications of microRNAs in Drosophila and predominantly targets precursor hairpins. Uridylation modulates the characteristic two-nucleotide 3' overhang of microRNA hairpins, which regulates processing by Dicer-1 and destabilizes RNA hairpins. Tailor preferentially uridylates mirtron hairpins, thereby impeding the production of non-canonical microRNAs. Mirtron selectivity is explained by primary sequence specificity of Tailor, selecting substrates ending with a 3' guanosine. In contrast to mirtrons, conserved Drosophila precursor microRNAs are significantly depleted in 3' guanosine, thereby escaping regulatory uridylation. Our data support the hypothesis that evolutionary adaptation to Tailor-directed uridylation shapes the nucleotide composition of precursor microRNA 3' ends. Hence, hairpin uridylation may serve as a barrier for the de novo creation of microRNAs in Drosophila.

SLAM-ITseq: sequencing cell type-specific transcriptomes without cell sorting.

Cell type-specific transcriptome analysis is an essential tool for understanding biological processes in which diverse types of cells are involved. Although cell isolation methods such as fluorescence-activated cell sorting (FACS) in combination with transcriptome analysis have widely been used so far, their time-consuming and harsh procedures limit their applications. Here, we report a novel in vivo metabolic RNA sequencing method, SLAM-ITseq, which metabolically labels RNA with 4-thiouracil in a specific cell type in vivo followed by detection through an RNA-seq-based method that specifically distinguishes the thiolated uridine by base conversion. This method has successfully identified the cell type-specific transcriptome in three different tissues: endothelial cells in brain, epithelial cells in intestine and adipocytes in white adipose tissue. As this method does not require isolation of cells or RNA prior to the transcriptomic analysis, SLAM-ITseq provides an easy yet accurate snapshot of the transcriptional state in vivo.

Thiol-linked alkylation of RNA to assess expression dynamics.

Gene expression profiling by high-throughput sequencing reveals qualitative and quantitative changes in RNA species at steady state but obscures the intracellular dynamics of RNA transcription, processing and decay. We developed thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM seq), an orthogonal-chemistry-based RNA sequencing technology that detects 4-thiouridine (s4U) incorporation in RNA species at single-nucleotide resolution. In combination with well-established metabolic RNA labeling protocols and coupled to standard, low-input, high-throughput RNA sequencing methods, SLAM seq enabled rapid access to RNA-polymerase-II-dependent gene expression dynamics in the context of total RNA. We validated the method in mouse embryonic stem cells by showing that the RNA-polymerase-II-dependent transcriptional output scaled with Oct4/Sox2/Nanog-defined enhancer activity, and we provide quantitative and mechanistic evidence for transcript-specific RNA turnover mediated by post-transcriptional gene regulatory pathways initiated by microRNAs and N6-methyladenosine. SLAM seq facilitates the dissection of fundamental mechanisms that control gene expression in an accessible, cost-effective and scalable manner.

Quantification of experimentally induced nucleotide conversions in high-throughput sequencing datasets.

BACKGROUND: Methods to read out naturally occurring or experimentally introduced nucleic acid modifications are emerging as powerful tools to study dynamic cellular processes. The recovery, quantification and interpretation of such events in high-throughput sequencing datasets demands specialized bioinformatics approaches. RESULTS: Here, we present Digital Unmasking of Nucleotide conversions in K-mers (DUNK), a data analysis pipeline enabling the quantification of nucleotide conversions in high-throughput sequencing datasets. We demonstrate using experimentally generated and simulated datasets that DUNK allows constant mapping rates irrespective of nucleotide-conversion rates, promotes the recovery of multimapping reads and employs Single Nucleotide Polymorphism (SNP) masking to uncouple true SNPs from nucleotide conversions to facilitate a robust and sensitive quantification of nucleotide-conversions. As a first application, we implement this strategy as SLAM-DUNK for the analysis of SLAMseq profiles, in which 4-thiouridine-labeled transcripts are detected based on T > C conversions. SLAM-DUNK provides both raw counts of nucleotide-conversion containing reads as well as a base-content and read coverage normalized approach for estimating the fractions of labeled transcripts as readout. CONCLUSION: Beyond providing a readily accessible tool for analyzing SLAMseq and related time-resolved RNA sequencing methods (TimeLapse-seq, TUC-seq), DUNK establishes a broadly applicable strategy for quantifying nucleotide conversions.

Comparison of EJC-enhanced and EJC-independent NMD in human cells reveals two partially redundant degradation pathways.

Nonsense-mediated mRNA decay (NMD) is a eukaryotic post-transcriptional gene regulation mechanism that eliminates mRNAs with the termination codon (TC) located in an unfavorable environment for efficient translation termination. The best-studied NMD-targeted mRNAs contain premature termination codons (PTCs); however, NMD regulates even many physiological mRNAs. An exon-junction complex (EJC) located downstream from a TC acts as an NMD-enhancing signal, but is not generally required for NMD. Here, we compared these "EJC-enhanced" and "EJC-independent" modes of NMD with regard to their requirement for seven known NMD factors in human cells using two well-characterized NMD reporter genes (immunoglobulin µ and ß-Globin) with or without an intron downstream from the PTC. We show that both NMD modes depend on UPF1 and SMG1, but detected transcript-specific differences with respect to the requirement for UPF2 and UPF3b, consistent with previously reported UPF2- and UPF3-independent branches of NMD. In addition and contrary to expectation, a higher sensitivity of EJC-independent NMD to reduced UPF2 and UPF3b concentrations was observed. Our data further revealed a redundancy of the endo- and exonucleolytic mRNA degradation pathways in both modes of NMD. Moreover, the relative contributions of both decay pathways differed between the reporters, with PTC-containing immunoglobulin µ transcripts being preferentially subjected to SMG6-mediated endonucleolytic cleavage, whereas ß-Globin transcripts were predominantly degraded by the SMG5/SMG7-dependent pathway. Overall, the surprising heterogeneity observed with only two NMD reporter pairs suggests the existence of several mechanistically distinct branches of NMD in human cells.

SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis.

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

Transcriptome-Wide Profiling of RNA Stability.

Gene expression is controlled at multiple levels, including RNA transcription and turnover. But determining the relative contributions of RNA biogenesis and decay to the steady-state abundance of cellular transcripts remains challenging because conventional transcriptomics approaches do not provide the temporal resolution to derive the kinetic parameters underlying steady-state gene expression.Here, we describe a protocol that combines metabolic RNA labeling by 4-thiouridine with chemical nucleoside conversion and whole-transcriptome sequencing followed by bioinformatics analysis to determine RNA stability in cultured cells at a genomic scale. Time-resolved transcriptomics by thiol (SH)-linked alkylation for the metabolic sequencing of RNA (SLAMseq) provides accurate information on transcript half-lives across annotated features in the genome, including by-products of transcription, such as introns. We provide a step-by-step instruction for time-resolved transcriptomics, which enhances traditional RNA sequencing protocols to acquire the temporal resolution required to directly measure the cellular kinetics of RNA turnover under physiological conditions.

Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets.

Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.

Determining mRNA Stability by Metabolic RNA Labeling and Chemical Nucleoside Conversion.

The varying rates at which mRNAs decay are tightly coordinated with transcriptional changes to shape gene expression during development and disease. But currently available RNA sequencing approaches lack the temporal information to determine the relative contribution of RNA biogenesis, processing and turnover to the establishment of steady-state gene expression profiles.Here, we describe a protocol that combines metabolic RNA labeling with chemical nucleoside conversion by thiol-linked alkylation of 4-thiouridine to determine RNA stability in cultured cells (SLAMseq). When coupled to cost-effective mRNA 3' end sequencing approaches, SLAMseq determines the half-life of polyadenylated transcripts in a global and transcript-specific manner using untargeted or targeted cDNA library preparation protocols.We provide a step-by-step instruction for time-resolved mRNA 3' end sequencing, which augments traditional RNA-seq approaches to acquire the temporal resolution necessary to study the molecular principles that control gene expression.

Publisher Correction: Sequencing cell-type-specific transcriptomes with SLAM-ITseq.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sequencing cell-type-specific transcriptomes with SLAM-ITseq.

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

Small RNAs Are Trafficked from the Epididymis to Developing Mammalian Sperm.

The biogenesis of the RNA payload of mature sperm is of great interest, because RNAs delivered to the zygote at fertilization can affect early development. Here, we tested the hypothesis that small RNAs are trafficked to mammalian sperm during the process of post-testicular maturation in the epididymis. By characterizing small RNA dynamics during germ cell maturation in mice, we confirm and extend prior observations that sperm undergo a dramatic switch in the RNA payload from piRNAs to tRNA fragments (tRFs) upon exiting the testis and entering the epididymis. Small RNA delivery to sperm could be recapitulated in vitro by incubating testicular spermatozoa with caput epididymosomes. Finally, tissue-specific metabolic labeling of RNAs in intact mice definitively shows that mature sperm carry RNAs that were originally synthesized in the epididymal epithelium. These data demonstrate that soma-germline RNA transfer occurs in male mammals, most likely via vesicular transport from the epididymis to maturing sperm.

SLAM-seq defines direct gene-regulatory functions of the BRD4-MYC axis.

Defining direct targets of transcription factors and regulatory pathways is key to understanding their roles in physiology and disease. We combined SLAM-seq [thiol(SH)-linked alkylation for the metabolic sequencing of RNA], a method for direct quantification of newly synthesized messenger RNAs (mRNAs), with pharmacological and chemical-genetic perturbation in order to define regulatory functions of two transcriptional hubs in cancer, BRD4 and MYC, and to interrogate direct responses to BET bromodomain inhibitors (BETis). We found that BRD4 acts as general coactivator of RNA polymerase II-dependent transcription, which is broadly repressed upon high-dose BETi treatment. At doses triggering selective effects in leukemia, BETis deregulate a small set of hypersensitive targets including MYC. In contrast to BRD4, MYC primarily acts as a selective transcriptional activator controlling metabolic processes such as ribosome biogenesis and de novo purine synthesis. Our study establishes a simple and scalable strategy to identify direct transcriptional targets of any gene or pathway.

A strand-specific switch in noncoding transcription switches the function of a Polycomb/Trithorax response element.

Polycomb/Trithorax response elements (PRE/TREs) can switch their function reversibly between silencing and activation by mechanisms that are poorly understood. Here we show that a switch in forward and reverse noncoding transcription from the Drosophila melanogaster vestigial (vg) PRE/TRE switches the status of the element between silencing (induced by the forward strand) and activation (induced by the reverse strand). In vitro, both noncoding RNAs inhibit PRC2 histone methyltransferase activity, but, in vivo, only the reverse strand binds PRC2. Overexpression of the reverse strand evicts PRC2 from chromatin and inhibits its enzymatic activity. We propose that the interaction of RNAs with PRC2 is differentially regulated in vivo, allowing regulated inhibition of local PRC2 activity. Genome-wide analysis shows that strand switching of noncoding RNAs occurs at several hundred Polycomb-binding sites in fly and vertebrate genomes. This work identifies a previously unreported and potentially widespread class of PRE/TREs that switch function by switching the direction of noncoding RNA transcription.