DUSP11 activity on triphosphorylated transcripts promotes Argonaute association with noncanonical viral microRNAs and regulates steady-state levels of cellular noncoding RNAs.

RNA silencing is a conserved eukaryotic gene expression regulatory mechanism mediated by small RNAs. In Caenorhabditis elegans, the accumulation of a distinct class of siRNAs synthesized by an RNA-dependent RNA polymerase (RdRP) requires the PIR-1 phosphatase. However, the function of PIR-1 in RNAi has remained unclear. Since mammals lack an analogous siRNA biogenesis pathway, an RNA silencing role for the mammalian PIR-1 homolog (dual specificity phosphatase 11 [DUSP11]) was unexpected. Here, we show that the RNA triphosphatase activity of DUSP11 promotes the RNA silencing activity of viral microRNAs (miRNAs) derived from RNA polymerase III (RNAP III) transcribed precursors. Our results demonstrate that DUSP11 converts the 5' triphosphate of miRNA precursors to a 5' monophosphate, promoting loading of derivative 5p miRNAs into Argonaute proteins via a Dicer-coupled 5' monophosphate-dependent strand selection mechanism. This mechanistic insight supports a likely shared function for PIR-1 in C. elegans Furthermore, we show that DUSP11 modulates the 5' end phosphate group and/or steady-state level of several host RNAP III transcripts, including vault RNAs and Alu transcripts. This study shows that steady-state levels of select noncoding RNAs are regulated by DUSP11 and defines a previously unknown portal for small RNA-mediated silencing in mammals, revealing that DUSP11-dependent RNA silencing activities are shared among diverse metazoans.

Team effort by TRAPP forces a nucleotide fumble.

TRAPPI is a multisubunit protein complex on the Golgi that activates the small GTPase Ypt1p to facilitate the receipt of transport vesicles inbound from the endoplasmic reticulum. Cai et al. (2008) now present structural and biochemical analyses of yeast TRAPPI in a complex with Ypt1p revealing a unique mechanism by which TRAPPI catalyzes guanine nucleotide exchange.

BCDIN3D regulates tRNAHis 3' fragment processing.

5' ends are important for determining the fate of RNA molecules. BCDIN3D is an RNA phospho-methyltransferase that methylates the 5' monophosphate of specific RNAs. In order to gain new insights into the molecular function of BCDIN3D, we performed an unbiased analysis of its interacting RNAs by Thermostable Group II Intron Reverse Transcriptase coupled to next generation sequencing (TGIRT-seq). Our analyses showed that BCDIN3D interacts with full-length phospho-methylated tRNAHis and miR-4454. Interestingly, we found that miR-4454 is not synthesized from its annotated genomic locus, which is a primer-binding site for an endogenous retrovirus, but rather by Dicer cleavage of mature tRNAHis. Sequence analysis revealed that miR-4454 is identical to the 3' end of tRNAHis. Moreover, we were able to generate this 'miRNA' in vitro through incubation of mature tRNAHis with Dicer. As found previously for several pre-miRNAs, a 5'P-tRNAHis appears to be a better substrate for Dicer cleavage than a phospho-methylated tRNAHis. Moreover, tRNAHis 3'-fragment/'miR-4454' levels increase in cells depleted for BCDIN3D. Altogether, our results show that in addition to microRNAs, BCDIN3D regulates tRNAHis 3'-fragment processing without negatively affecting tRNAHis's canonical function of aminoacylation.

Simultaneous sequencing of coding and noncoding RNA reveals a human transcriptome dominated by a small number of highly expressed noncoding genes.

Comparing the abundance of one RNA molecule to another is crucial for understanding cellular functions but most sequencing techniques can target only specific subsets of RNA. In this study, we used a new fragmented ribodepleted TGIRT sequencing method that uses a thermostable group II intron reverse transcriptase (TGIRT) to generate a portrait of the human transcriptome depicting the quantitative relationship of all classes of nonribosomal RNA longer than 60 nt. Comparison between different sequencing methods indicated that FRT is more accurate in ranking both mRNA and noncoding RNA than viral reverse transcriptase-based sequencing methods, even those that specifically target these species. Measurements of RNA abundance in different cell lines using this method correlate with biochemical estimates, confirming tRNA as the most abundant nonribosomal RNA biotype. However, the single most abundant transcript is 7SL RNA, a component of the signal recognition particle. Structured noncoding RNAs (sncRNAs) associated with the same biological process are expressed at similar levels, with the exception of RNAs with multiple functions like U1 snRNA. In general, sncRNAs forming RNPs are hundreds to thousands of times more abundant than their mRNA counterparts. Surprisingly, only 50 sncRNA genes produce half of the non-rRNA transcripts detected in two different cell lines. Together the results indicate that the human transcriptome is dominated by a small number of highly expressed sncRNAs specializing in functions related to translation and splicing.

Broad role for YBX1 in defining the small noncoding RNA composition of exosomes.

RNA is secreted from cells enclosed within extracellular vesicles (EVs). Defining the RNA composition of EVs is challenging due to their coisolation with contaminants, lack of knowledge of the mechanisms of RNA sorting into EVs, and limitations of conventional RNA-sequencing methods. Here we present our observations using thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) to characterize the RNA extracted from HEK293T cell EVs isolated by flotation gradient ultracentrifugation and from exosomes containing the tetraspanin CD63 further purified from the gradient fractions by immunoisolation. We found that EV-associated transcripts are dominated by full-length, mature transfer RNAs (tRNAs) and other small noncoding RNAs (ncRNAs) encapsulated within vesicles. A substantial proportion of the reads mapping to protein-coding genes, long ncRNAs, and antisense RNAs were due to DNA contamination on the surface of vesicles. Nevertheless, sequences mapping to spliced mRNAs were identified within HEK293T cell EVs and exosomes, among the most abundant being transcripts containing a 5' terminal oligopyrimidine (5' TOP) motif. Our results indicate that the RNA-binding protein YBX1, which is required for the sorting of selected miRNAs into exosomes, plays a role in the sorting of highly abundant small ncRNA species, including tRNAs, Y RNAs, and Vault RNAs. Finally, we obtained evidence for an EV-specific tRNA modification, perhaps indicating a role for posttranscriptional modification in the sorting of some RNA species into EVs. Our results suggest that EVs and exosomes could play a role in the purging and intercellular transfer of excess free RNAs, including full-length tRNAs and other small ncRNAs.

RNA-seq of human reference RNA samples using a thermostable group II intron reverse transcriptase.

Next-generation RNA sequencing (RNA-seq) has revolutionized our ability to analyze transcriptomes. Current RNA-seq methods are highly reproducible, but each has biases resulting from different modes of RNA sample preparation, reverse transcription, and adapter addition, leading to variability between methods. Moreover, the transcriptome cannot be profiled comprehensively because highly structured RNAs, such as tRNAs and snoRNAs, are refractory to conventional RNA-seq methods. Recently, we developed a new method for strand-specific RNA-seq using thermostable group II intron reverse transcriptases (TGIRTs). TGIRT enzymes have higher processivity and fidelity than conventional retroviral reverse transcriptases plus a novel template-switching activity that enables RNA-seq adapter addition during cDNA synthesis without using RNA ligase. Here, we obtained TGIRT-seq data sets for well-characterized human RNA reference samples and compared them to previous data sets obtained for these RNAs by the Illumina TruSeq v2 and v3 methods. We find that TGIRT-seq recapitulates the relative abundance of human transcripts and RNA spike-ins in ribo-depleted, fragmented RNA samples comparably to non-strand-specific TruSeq v2 and better than strand-specific TruSeq v3. Moreover, TGIRT-seq is more strand specific than TruSeq v3 and eliminates sampling biases from random hexamer priming, which are inherent to TruSeq. The TGIRT-seq data sets also show more uniform 5' to 3' gene coverage and identify more splice junctions, particularly near the 5' ends of mRNAs, than do the TruSeq data sets. Finally, TGIRT-seq enables the simultaneous profiling of mRNAs and lncRNAs in the same RNA-seq experiment as structured small ncRNAs, including tRNAs, which are essentially absent with TruSeq.

High-throughput sequencing of human plasma RNA by using thermostable group II intron reverse transcriptases.

Next-generation RNA-sequencing (RNA-seq) has revolutionized transcriptome profiling, gene expression analysis, and RNA-based diagnostics. Here, we developed a new RNA-seq method that exploits thermostable group II intron reverse transcriptases (TGIRTs) and used it to profile human plasma RNAs. TGIRTs have higher thermostability, processivity, and fidelity than conventional reverse transcriptases, plus a novel template-switching activity that can efficiently attach RNA-seq adapters to target RNA sequences without RNA ligation. The new TGIRT-seq method enabled construction of RNA-seq libraries from <1 ng of plasma RNA in <5 h. TGIRT-seq of RNA in 1-mL plasma samples from a healthy individual revealed RNA fragments mapping to a diverse population of protein-coding gene and long ncRNAs, which are enriched in intron and antisense sequences, as well as nearly all known classes of small ncRNAs, some of which have never before been seen in plasma. Surprisingly, many of the small ncRNA species were present as full-length transcripts, suggesting that they are protected from plasma RNases in ribonucleoprotein (RNP) complexes and/or exosomes. This TGIRT-seq method is readily adaptable for profiling of whole-cell, exosomal, and miRNAs, and for related procedures, such as HITS-CLIP and ribosome profiling.

RUTBC2 protein, a Rab9A effector and GTPase-activating protein for Rab36.

Rab GTPases regulate vesicle budding, motility, docking, and fusion. In cells, their cycling between active, GTP-bound states and inactive, GDP-bound states is regulated by the action of opposing enzymes called guanine nucleotide exchange factors and GTPase-activating proteins (GAPs). The substrates for most RabGAPs are unknown, and the potential for cross-talk between different membrane trafficking pathways remains uncharted territory. Rab9A and its effectors regulate recycling of mannose 6-phosphate receptors from late endosomes to the trans Golgi network. We show here that RUTBC2 is a TBC domain-containing protein that binds to Rab9A specifically both in vitro and in cultured cells but is not a GAP for Rab9A. Biochemical screening of Rab protein substrates for RUTBC2 revealed highest GAP activity toward Rab34 and Rab36. In cells, membrane-associated RUTBC2 co-localizes with Rab36, and expression of wild type RUTBC2, but not the catalytically inactive, RUTBC2 R829A mutant, decreases the amount of membrane-associated Rab36 protein. These data show that RUTBC2 can act as a Rab36 GAP in cells and suggest that RUTBC2 links Rab9A function to Rab36 function in the endosomal system.

Arg-tRNA synthetase links inflammatory metabolism to RNA splicing and nuclear trafficking via SRRM2.

Cells respond to perturbations such as inflammation by sensing changes in metabolite levels. Especially prominent is arginine, which has known connections to the inflammatory response. Aminoacyl-tRNA synthetases, enzymes that catalyse the first step of protein synthesis, can also mediate cell signalling. Here we show that depletion of arginine during inflammation decreased levels of nuclear-localized arginyl-tRNA synthetase (ArgRS). Surprisingly, we found that nuclear ArgRS interacts and co-localizes with serine/arginine repetitive matrix protein 2 (SRRM2), a spliceosomal and nuclear speckle protein, and that decreased levels of nuclear ArgRS correlated with changes in condensate-like nuclear trafficking of SRRM2 and splice-site usage in certain genes. These splice-site usage changes cumulated in the synthesis of different protein isoforms that altered cellular metabolism and peptide presentation to immune cells. Our findings uncover a mechanism whereby an aminoacyl-tRNA synthetase cognate to a key amino acid that is metabolically controlled during inflammation modulates the splicing machinery.

RUTBC1 protein, a Rab9A effector that activates GTP hydrolysis by Rab32 and Rab33B proteins.

Rab GTPases regulate all steps of membrane trafficking. Their interconversion between active, GTP-bound states and inactive, GDP-bound states is regulated by guanine nucleotide exchange factors and GTPase-activating proteins. The substrates for most Rab GTPase-activating proteins (GAPs) are unknown. Rab9A and its effectors regulate transport of mannose 6-phosphate receptors from late endosomes to the trans-Golgi network. We show here that RUTBC1 is a Tre2/Bub2/Cdc16 domain-containing protein that binds to Rab9A-GTP both in vitro and in cultured cells, but is not a GTPase-activating protein for Rab9A. Biochemical screening of RUTBC1 Rab protein substrates revealed highest in vitro GTP hydrolysis-activating activity with Rab32 and Rab33B. Catalysis required Arg-803 of RUTBC1, and RUTBC1 could activate a catalytically inhibited Rab33B mutant (Q92A), in support of a dual finger mechanism for RUTBC1 action. Rab9A binding did not influence GAP activity of bead-bound RUTBC1 protein. In cells and cell extracts, RUTBC1 influenced the ability of Rab32 to bind its effector protein, Varp, consistent with a physiological role for RUTBC1 in regulating Rab32. In contrast, binding of Rab33B to its effector protein, Atg16L1, was not influenced by RUTBC1 in cells or extracts. The identification of a protein that binds Rab9A and inactivates Rab32 supports a model in which Rab9A and Rab32 act in adjacent pathways at the boundary between late endosomes and the biogenesis of lysosome-related organelles.

High-grade ovarian cancer associated H/ACA snoRNAs promote cancer cell proliferation and survival.

Small nucleolar RNAs (snoRNAs) are an omnipresent class of non-coding RNAs involved in the modification and processing of ribosomal RNA (rRNA). As snoRNAs are required for ribosome production, the increase of which is a hallmark of cancer development, their expression would be expected to increase in proliferating cancer cells. However, assessing the nature and extent of snoRNAs' contribution to cancer biology has been largely limited by difficulties in detecting highly structured RNA. In this study, we used a dedicated midsize non-coding RNA (mncRNA) sensitive sequencing technique to accurately survey the snoRNA abundance in independently verified high-grade serous ovarian carcinoma (HGSC) and serous borderline tumour (SBT) tissues. The results identified SNORA81, SNORA19 and SNORA56 as an H/ACA snoRNA signature capable of discriminating between independent sets of HGSC, SBT and normal tissues. The expression of the signature SNORA81 correlates with the level of ribosomal RNA (rRNA) modification and its knockdown inhibits 28S rRNA pseudouridylation and accumulation leading to reduced cell proliferation and migration. Together our data indicate that specific subsets of H/ACA snoRNAs may promote tumour aggressiveness by inducing rRNA modification and synthesis.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs.

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

Identification of protein-protected mRNA fragments and structured excised intron RNAs in human plasma by TGIRT-seq peak calling.

Human plasma contains > 40,000 different coding and non-coding RNAs that are potential biomarkers for human diseases. Here, we used thermostable group II intron reverse transcriptase sequencing (TGIRT-seq) combined with peak calling to simultaneously profile all RNA biotypes in apheresis-prepared human plasma pooled from healthy individuals. Extending previous TGIRT-seq analysis, we found that human plasma contains largely fragmented mRNAs from > 19,000 protein-coding genes, abundant full-length, mature tRNAs and other structured small non-coding RNAs, and less abundant tRNA fragments and mature and pre-miRNAs. Many of the mRNA fragments identified by peak calling correspond to annotated protein-binding sites and/or have stable predicted secondary structures that could afford protection from plasma nucleases. Peak calling also identified novel repeat RNAs, miRNA-sized RNAs, and putatively structured intron RNAs of potential biological, evolutionary, and biomarker significance, including a family of full-length excised intron RNAs, subsets of which correspond to mirtron pre-miRNAs or agotrons.

A functional role for the GCC185 golgin in mannose 6-phosphate receptor recycling.

Mannose 6-phosphate receptors (MPRs) deliver newly synthesized lysosomal enzymes to endosomes and then recycle to the Golgi. MPR recycling requires Rab9 GTPase; Rab9 recruits the cytosolic adaptor TIP47 and enhances its ability to bind to MPR cytoplasmic domains during transport vesicle formation. Rab9-bearing vesicles then fuse with the trans-Golgi network (TGN) in living cells, but nothing is known about how these vesicles identify and dock with their target. We show here that GCC185, a member of the Golgin family of putative tethering proteins, is a Rab9 effector that is required for MPR recycling from endosomes to the TGN in living cells, and in vitro. GCC185 does not rely on Rab9 for its TGN localization; depletion of GCC185 slightly alters the Golgi ribbon but does not interfere with Golgi function. Loss of GCC185 triggers enhanced degradation of mannose 6-phosphate receptors and enhanced secretion of hexosaminidase. These data assign a specific pathway to an interesting, TGN-localized protein and suggest that GCC185 may participate in the docking of late endosome-derived, Rab9-bearing transport vesicles at the TGN.

Distinct mechanisms of microRNA sorting into cancer cell-derived extracellular vesicle subtypes.

Extracellular vesicles (EVs) encompass a variety of vesicles secreted into the extracellular space. EVs have been implicated in promoting tumor metastasis, but the molecular composition of tumor-derived EV sub-types and the mechanisms by which molecules are sorted into EVs remain mostly unknown. We report the separation of two small EV sub-populations from a metastatic breast cancer cell line, with biochemical features consistent with different sub-cellular origins. These EV sub-types use different mechanisms of miRNA sorting (selective and non-selective), suggesting that sorting occurs via fundamentally distinct processes, possibly dependent on EV origin. Using biochemical and genetic tools, we identified the Lupus La protein as mediating sorting of selectively packaged miRNAs. We found that two motifs embedded in miR-122 are responsible for high-affinity binding to Lupus La and sorting into vesicles formed in a cell-free reaction. Thus, tumor cells can simultaneously deploy multiple EV species using distinct sorting mechanisms that may enable diverse functions in normal and cancer biology.

Measuring Rab GTPase-activating protein (GAP) activity in live cells and extracts.

Mammalian cells encode a diverse set of Rab GTPases and their corresponding regulators. In vitro biochemical screening has proven invaluable in assigning particular Rabs as substrates for their cognate GTPase-activating proteins. However, in vitro activity does not always reflect substrate specificity in cells. This method describes a functional test of GAP activity in cells or extracts that takes into account the presence of other factors or conditions that might change observed in vitro specificity.

Mutant enzymes challenge all assumptions.

Enzymes called Rab GTPases that carry so-called "activating" mutations may never become activated at all.

Multiple Rab GTPase binding sites in GCC185 suggest a model for vesicle tethering at the trans-Golgi.

GCC185, a trans-Golgi network-localized protein predicted to assume a long, coiled-coil structure, is required for Rab9-dependent recycling of mannose 6-phosphate receptors (MPRs) to the Golgi and for microtubule nucleation at the Golgi via CLASP proteins. GCC185 localizes to the Golgi by cooperative interaction with Rab6 and Arl1 GTPases at adjacent sites near its C terminus. We show here by yeast two-hybrid and direct biochemical tests that GCC185 contains at least four additional binding sites for as many as 14 different Rab GTPases across its entire length. A central coiled-coil domain contains a specific Rab9 binding site, and functional assays indicate that this domain is important for MPR recycling to the Golgi complex. N-Terminal coiled-coils are also required for GCC185 function as determined by plasmid rescue after GCC185 depletion by using small interfering RNA in cultured cells. Golgi-Rab binding sites may permit GCC185 to contribute to stacking and lateral interactions of Golgi cisternae as well as help it function as a vesicle tether.