"RNA therapeutics" refers to a disease treatment or drug that utilizes RNA as a component. In this context, RNA may be the direct target of a small-molecule drug or RNA itself may be the drug, designed to bind to a protein, or to mimic or target another RNA. RNA has gained attention in the drug-development world, as recent clinical successes and breakthrough technologies have revolutionized the drug-like qualities of the molecule or its usefulness as a drug target. In this special issue of RNA, we gathered expert perspectives on the past, present, and future of the field, to serve as a primer and also a challenge to the broad scientific community to incorporate RNA into their experimental design and problem-solving process, and to imagine and realize the potential of RNA as a therapeutic drug or target.

Darwin's assertion that "it is mere rubbish thinking, at present, of origin of life" is no longer valid. By synthesizing origin of life (OoL) research from its inception to recent findings, with a focus on (i) proof-of-principle prebiotically plausible syntheses and (ii) molecular relics of the ancient RNA World, we present a comprehensive up-to-date description of science's understanding of the OoL and the RNA World hypothesis. Based on these observations, we solidify the consensus that RNA evolved before coded proteins and DNA genomes, such that the biosphere began with an RNA core where much of the translation apparatus and related RNA architecture arose before RNA transcription and DNA replication. This supports the conclusion that the OoL was a gradual process of chemical evolution involving a series of transitional forms between prebiotic chemistry and the last universal common ancestor (LUCA) during which RNA played a central role, and that many of the events and their relative order of occurrence along this pathway are known. The integrative nature of this synthesis also extends previous descriptions and concepts and should help inform future questions and experiments about the ancient RNA World and the OoL.

RNA sequencing has spurred a significant number of research areas in recent years. Most protocols rely on synthesizing a more stable complementary DNA (cDNA) copy of the RNA molecule during the reverse transcription reaction. The resulting cDNA pool is often wrongfully assumed to be quantitatively and molecularly similar to the original RNA input. Sadly, biases and artifacts confound the resulting cDNA mixture. These issues are often overlooked or ignored in the literature by those that rely on the reverse transcription process. In this review, we confront the reader with intra- and intersample biases and artifacts caused by the reverse transcription reaction during RNA sequencing experiments. To fight the reader's despair, we also provide solutions to most issues and inform on good RNA sequencing practices. We hope the reader can use this review to their advantage, thereby contributing to scientifically sound RNA studies.

Synthetic RNA oligonucleotides composed of canonical and modified ribonucleotides are highly effective for RNA antisense therapeutics and RNA-based genome engineering applications utilizing CRISPR-Cas9. Yet, synthesis of synthetic RNA using phosphoramidite chemistry is highly inefficient and expensive relative to DNA oligonucleotides, especially for relatively long RNA oligonucleotides. Thus, new biotechnologies are needed to significantly reduce costs, while increasing synthesis rates and yields of synthetic RNA. Here, we engineer human DNA polymerase theta (Polθ) variants and demonstrate their ability to synthesize long (95-200 nt) RNA oligonucleotides with canonical ribonucleotides and ribonucleotide analogs commonly used for stabilizing RNA for therapeutic and genome engineering applications. In contrast to natural promoter-dependent RNA polymerases, Polθ variants synthesize RNA by initiating from DNA or RNA primers, which enables the production of RNA without short abortive byproducts. Remarkably, Polθ variants show the lower capacity to misincorporate ribonucleotides compared to T7 RNA polymerase. Automation of this enzymatic RNA synthesis technology can potentially increase yields while reducing costs of synthetic RNA oligonucleotide production.

The ribosome is the universally conserved ribozyme that translates DNA coded instructions into proteins with the assistance of other RNA molecules, including transfer and messenger RNAs. Of particular interest is the segmentation phenomena, which is found in trypanosomatids and other protists. In these organisms, the large subunit ribosomal RNA is assembled from multiple smaller RNAs. This phenomenon posits several challenges to the folding and stabilization of such ribosomes to retain functionality and efficiency. In earlier studies, RNA/protein interactions were suggested to fully compensate for the fragmentation. Recently, several conserved RNA/RNA interaction regions were described in the cryo-EM structures of segmented ribosomes from trypanosomatids. These regions also seemed to aid in the folding and stabilization of such ribosomes, even before the ribosomal proteins start their association. In the present study, the existence of conserved RNA/RNA interaction regions shared between trypanosomatid and Euglena gracilis segmented ribosomes was confirmed, despite differences in segmentation patterns. Analysis of the crystallographic structures of unsegmented ribosomes from other Eukaryotes, Bacteria, and Archaea allowed us to estimate the relative age of highly conserved RNA/RNA interaction regions. These results strongly suggest that common interaction regions likely date far back into the ribosomes of the last common ancestor. Results also revealed that single hydrogen bonds are overwhelmingly facilitated by the 2'OH, a distinctive RNA feature. This supports the notion that RNA predates DNA and places some constraints on alternative nucleic acids proposals.

DEAD-box RNA helicases are implicated in most aspects of RNA biology, where these enzymes unwind short RNA duplexes in an ATP-dependent manner. During the central step of the unwinding cycle, the two domains of the helicase core form a distinct closed conformation that destabilizes the RNA duplex, which ultimately leads to duplex melting. Despite the importance of this step for the unwinding process no high-resolution structures of this state are available. Here, I used nuclear magnetic resonance spectroscopy and X-ray crystallography to determine structures of the DEAD-box helicase DbpA in the closed conformation, complexed with substrate duplexes and single-stranded unwinding product. These structures reveal that DbpA initiates duplex unwinding by interacting with up to three base-paired nucleotides and a 5' single-stranded RNA duplex overhang. These high-resolution snapshots, together with biochemical assays, rationalize the destabilization of the RNA duplex and are integrated into a conclusive model of the unwinding process.

The study of eukaryotic tRNA processing has given rise to an explosion of new information and insights in the last several years. We now have unprecedented knowledge of each step in the tRNA processing pathway, revealing unexpected twists in biochemical pathways, multiple new connections with regulatory pathways, and numerous biological effects of defects in processing steps that have profound consequences throughout eukaryotes, leading to growth phenotypes in the yeast Saccharomyces cerevisiae and to neurological and other disorders in humans. This review highlights seminal new results within the pathways that comprise the life of a tRNA, from its birth after transcription until its death by decay. We focus on new findings and revelations in each step of the pathway including the end-processing and splicing steps, many of the numerous modifications throughout the main body and anticodon loop of tRNA that are so crucial for tRNA function, the intricate tRNA trafficking pathways, and the quality control decay pathways, as well as the biogenesis and biology of tRNA-derived fragments. We also describe the many interactions of these pathways with signaling and other pathways in the cell.

Although RNA plays a vital role in gene expression, it is less used as an in situ biomarker for clinical diagnostics than DNA and protein. This is mainly due to technical challenges caused by the low expression level and easy degradation of RNA molecules. To tackle this issue, methods that are sensitive and specific are needed. Here, we present an RNA single-molecule chromogenic in situ hybridization assay based on DNA probe proximity ligation and rolling circle amplification. When the DNA probes hybridize into close proximity to the RNA molecules, they form a V-shape structure and mediate the circularization of circle probes. Thus, our method was termed vsmCISH. We successfully applied our method to assess HER2 mRNA expression status in invasive breast cancer tissue and investigated the utility of albumin mRNA ISH for differentiating primary from metastatic liver cancer. The promising results on clinical samples indicate that our method has great potential for application in diagnosing diseases using RNA biomarkers.

As we continue to find new regulatory roles for RNAs, a theme is emerging in which regulation may not be mediated through the actions of a specific RNA, as one typically thinks of a regulator and target, but rather through the collective nature of many RNAs, each contributing a small degree of the regulatory load. This mechanism has been termed "crowd-control" and may apply broadly to miRNAs and to RNAs that bind and regulate protein activity. This provides an alternative way of thinking about how RNAs can act as biological regulators and has repercussions, both for the understanding of biological systems, and for the interpretation of results in which individual members of the "crowd" can replicate the effects of the crowd when overexpressed, but are not individually significant biological regulators.

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

Differential expression analysis of RNA sequencing (RNA-seq) data can identify changes in cellular RNA levels, but provides limited information about the kinetic mechanisms underlying such changes. Nucleotide recoding RNA-seq methods (NR-seq; e.g., TimeLapse-seq, SLAM-seq, etc.) address this shortcoming and are widely used approaches to identify changes in RNA synthesis and degradation kinetics. While advanced statistical models implemented in user-friendly software (e.g., DESeq2) have ensured the statistical rigor of differential expression analyses, no such tools that facilitate differential kinetic analysis with NR-seq exist. Here, we report the development of Bayesian analysis of the kinetics of RNA (bakR; https:// github.com/simonlabcode/bakR), an R package to address this need. bakR relies on Bayesian hierarchical modeling of NR-seq data to increase statistical power by sharing information across transcripts. Analyses of simulated data confirmed that bakR implementations of the hierarchical model outperform attempts to analyze differential kinetics with existing models. bakR also uncovers biological signals in real NR-seq data sets and provides improved analyses of existing data sets. This work establishes bakR as an important tool for identifying differential RNA synthesis and degradation kinetics.

The SRP9/SRP14 heterodimer is a central component of signal recognition particle (SRP) RNA (7SL) processing and Alu retrotransposition. In this study, we sought to establish the role of nuclear SRP9/SRP14 in the transcriptional regulation of 7SL and BC200 RNA. 7SL and BC200 RNA steady-state levels, rate of decay, and transcriptional activity were evaluated under SRP9/SRP14 knockdown conditions. Immunofluorescent imaging, and subcellular fractionation of MCF-7 cells, revealed a distinct nuclear localization for SRP9/SRP14. The relationship between this localization and transcriptional activity at 7SL and BC200 genes was also examined. These findings demonstrate a novel nuclear function of SRP9/SRP14 establishing that this heterodimer transcriptionally regulates 7SL and BC200 RNA expression. We describe a model in which SRP9/SRP14 cotranscriptionally regulate 7SL and BC200 RNA expression. Our model is also a plausible pathway for regulating Alu RNA transcription and is consistent with the hypothesized roles of SRP9/SRP14 transporting 7SL RNA into the nucleolus for posttranscriptional processing, and trafficking of Alu RNA for retrotransposition.

The multifunctional RNA recognition motif-containing protein Y14/RBM8A participates in mRNA metabolism and is essential for the efficient repair of DNA double-strand breaks (DSBs). Y14 contains highly charged, low-complexity sequences in both the amino- and carboxy-terminal domains. The feature of charge segregation suggests that Y14 may undergo liquid-liquid phase separation (LLPS). Recombinant Y14 formed phase-separated droplets, which were sensitive to pH and salt concentration. Domain mapping suggested that LLPS of Y14 involves multivalent electrostatic interactions and is partly determined by the net charge of its low-complexity regions. Phospho-mimicry of the carboxy-terminal arginine-serine dipeptides of Y14 suppressed phase separation. Moreover, RNA could phase separate into Y14 droplets and modulate Y14 LLPS in a concentration-dependent manner. Finally, the capacity of Y14 in LLPS and coacervation with RNA in vitro correlated with its activity in DSB repair. These results reveal a molecular rule for LLPS of Y14 in vitro and an implication for its co-condensation with RNA in genome stability.

Aptamers with fluorogenic ligands are emerging as useful tools to quantify and track RNA molecules. The RNA Mango family of aptamers have a useful combination of tight ligand binding, bright fluorescence, and small size. However, the simple structure of these aptamers, with a single base-paired stem capped by a G-quadruplex, can limit the sequence and structural modifications needed for many use-inspired designs. Here we report new structural variants of RNA Mango that have two base-paired stems attached to the quadruplex. Fluorescence saturation analysis of one of the double-stemmed constructs showed a maximum fluorescence that is ∼75% brighter than the original single-stemmed Mango I. A small number of mutations to nucleotides in the tetraloop-like linker of the second stem were subsequently analyzed. The effect of these mutations on the affinity and fluorescence suggested that the nucleobases of the second linker do not directly interact with the fluorogenic ligand (TO1-biotin), but may instead induce higher fluorescence by indirectly altering the ligand properties in the bound state. The effects of the mutations in this second tetraloop-like linker indicate the potential of this second stem for rational design and reselection experiments. Additionally, we demonstrated that a bimolecular mango designed by splitting the double-stemmed Mango can function when two RNA molecules are cotranscribed from different DNA templates in a single in vitro transcription. This bimolecular Mango has potential application in detecting RNA-RNA interactions. Together, these constructs expand the designability of the Mango aptamers to facilitate future applications of RNA imaging.

The addition of chemical groups to cellular RNA to modulate RNA fate and/or function is summarized under the term epitranscriptomic modification. More than 170 different modifications have been identified on cellular RNA, such as tRNA, rRNA and, to a lesser extent, on other RNA types. Recently, epitranscriptomic modification of viral RNA has received considerable attention as a possible additional mechanism regulating virus infection and replication. N6-methyladenosine (m6A) and C5-methylcytosine (m5C) have been most broadly studied in different RNA viruses. Various studies, however, reported varying results with regard to number and extent of the modification. Here we investigated the m5C methylome of SARS-CoV-2, and we reexamined reported m5C sites in HIV and MLV. Using a rigorous bisulfite-sequencing protocol and stringent data analysis, we found no evidence for the presence of m5C in these viruses. The data emphasize the necessity for optimizing experimental conditions and bioinformatic data analysis.

RNA-binding proteins (RBPs) are key regulators of gene expression. Small molecules targeting these RBP-RNA interactions are a rapidly emerging class of therapeutics for treating a variety of diseases. Ro-08-2750 (Ro) is a small molecule identified as a competitive inhibitor of Musashi (MSI)-RNA interactions. Here, we show that multiple Ro-dependent cellular phenotypes, specifically adrenocortical steroid production and cell viability, are Musashi-2 (MSI2)-independent. Using an unbiased proteome-wide approach, we discovered Ro broadly interacts with RBPs, many containing RRM domains. To confirm this finding, we leveraged the large-scale ENCODE data to identify a subset of RBPs whose depletion phenocopies Ro inhibition, indicating Ro is a promiscuous inhibitor of multiple RBPs. Consistent with broad disruption of ribonucleoprotein complexes, Ro treatment leads to stress granule formation. This strategy represents a generalizable framework for validating the specificity and identifying targets of RBP inhibitors in a cellular context.

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

The GIGYF proteins interact with 4EHP and RNA-associated proteins to elicit transcript-specific translational repression. However, the mechanism by which the GIGYF1/2-4EHP complex is recruited to its target transcripts remain unclear. Here, we report the crystal structures of the GYF domains from GIGYF1 and GIGYF2 in complex with proline-rich sequences from the miRISC-binding proteins TNRC6C and TNRC6A, respectively. The TNRC6 proline-rich motifs bind to a conserved array of aromatic residues on the surface of the GIGYF1/2 GYF domains, thereby bridging 4EHP to Argonaute-miRNA complexes. Our structures also reveal a phenylalanine residue conserved from yeast to human GYF domains that contributes to GIGYF2 thermostability. The molecular details we outline here are likely to be conserved between GIGYF1/2 and other RNA-binding proteins to elicit 4EHP-mediated repression in different biological contexts.

Understanding the frequency and structural context of discrete noncovalent interactions between nucleotides is of pivotal significance in establishing the rules that govern RNA structure and dynamics. Although T-shaped contacts (i.e., perpendicular stacking contacts) between aromatic amino acids and nucleobases at the nucleic acid-protein interface have recently garnered attention, the analogous contacts within the nucleic acid structures have not been discussed. In this work, we have developed an automated method for identifying and unambiguously classifying T-shaped interactions between nucleobases. Using this method, we identified a total of 3261 instances of T-shaped (perpendicular stacking) contacts between two nucleobases in an array of RNA structures from an up-to-date data set of ≤3.5 Å resolution crystal structures deposited in the Protein Data Bank.

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

Ribosomal RNA (rRNA) maturation in archaea is a complex multistep process that requires well-defined endo- and exoribonuclease activities to generate fully mature linear rRNAs. However, technical challenges prevented detailed mapping of rRNA processing steps and a systematic analysis of rRNA maturation pathways across the tree of life. In this study, we used long-read (PCR)-cDNA and direct RNA nanopore-based sequencing to study rRNA maturation in three archaeal model organisms, namely the Euryarchaea Haloferax volcanii and Pyrococcus furiosus and the Crenarchaeon Sulfolobus acidocaldarius Compared to standard short-read protocols, nanopore sequencing facilitates simultaneous readout of 5'- and 3'-positions, which is required for the classification of rRNA processing intermediates. More specifically, we (i) accurately detect and describe rRNA maturation stages by analysis of terminal read positions of cDNA reads and thereupon (ii) explore the stage-dependent installation of the KsgA-mediated dimethylations in H. volcanii using base-calling and signal characteristics of direct RNA reads. Due to the single-molecule sequencing capacity of nanopore sequencing, we could detect hitherto unknown intermediates with high confidence, revealing details about the maturation of archaea-specific circular rRNA intermediates. Taken together, our study delineates common principles and unique features of rRNA processing in euryarchaeal and crenarchaeal representatives, thereby significantly expanding our understanding of rRNA maturation pathways in archaea.

k-Junctions are elaborated forms of kink turns with an additional helix on the nonbulged strand, thus forming a three-way helical junction. Two were originally identified in the structures of Arabidopsis and Escherichia coli thiamine pyrophosphate (TPP) riboswitches, and another called DUF-3268 was tentatively identified from sequence information. In this work we show that the Arabidopsis and E. coli riboswitch k-junctions fold in response to the addition of magnesium or sodium ions, and that atomic mutations that should disrupt key hydrogen bonding interactions greatly impair folding. Using X-ray crystallography, we have determined the structure of the DUF-3268 RNA and thus confirmed that it is a k-junction. It also folds upon the addition of metal ions, though requiring a 40-fold lower concentration of either divalent or monovalent ions. The key difference between the DUF-3268 and riboswitch k-junctions is the lack of nucleotides inserted between G1b and A2b in the former. We show that this insertion is primarily responsible for the difference in folding properties. Finally, we show that the DUF-3268 can functionally substitute for the k-junction in the E. coli TPP riboswitch such that the chimera can bind the TPP ligand, although less avidly.

Pseudouridine is the most frequently naturally occurring RNA modification, found in all classes of biologically functional RNAs. Compared to uridine, pseudouridine contains an additional hydrogen bond donor group and is therefore widely regarded as a structure stabilizing modification. However, the effects of pseudouridine modifications on the structure and dynamics of RNAs have so far only been investigated in a limited number of different structural contexts. Here, we introduced pseudouridine modifications into the U-turn motif and the adjacent U:U closing base pair of the neomycin-sensing riboswitch (NSR)-an extensively characterized model system for RNA structure, ligand binding, and dynamics. We show that the effects of replacing specific uridines with pseudouridines on RNA dynamics crucially depend on the exact location of the replacement site and can range from destabilizing to locally or even globally stabilizing. By using a combination of NMR spectroscopy, MD simulations and QM calculations, we rationalize the observed effects on a structural and dynamical level. Our results will help to better understand and predict the consequences of pseudouridine modifications on the structure and function of biologically important RNAs.

RT-PCR and northern blots have long been used to study RNA isoforms usage for single genes. Recent advancements in long-read sequencing have yielded unprecedented information about the usage and abundance of these RNA isoforms. However, visualization of long-read sequencing data remains challenging due to the high information density. To alleviate these issues, we have developed NanoBlot, an open-source R-package that generates northern blot and RT-PCR-like images from long-read sequencing data. NanoBlot requires aligned, positionally sorted and indexed BAM files. Plotting is based around ggplot2 and is easily customizable. Advantages of NanoBlot include a robust system for designing probes to visualize isoforms including excluding reads based on the presence or absence of a specified region, an elegant solution to representing isoforms with continuous variations in length, and the ability to overlay multiple genes in the same plot using different colors. We present examples of nanoblots compared to actual northern blot data. In addition to traditional gel-like images, the NanoBlot package can also output other visualizations such as violin plots and 3'-RACE-like plots focused on 3'-end isoforms visualization. The use of the NanoBlot package should provide a simple answer to some of the challenges of visualizing long-read RNA-sequencing data.

The ribosome is a large ribonucleoprotein assembly that uses diverse and complex molecular interactions to maintain proper folding. In vivo assembled ribosomes have been isolated using MS2 tags installed in either the 16S or 23S ribosomal RNAs (rRNAs), to enable studies of ribosome structure and function in vitro. RNA tags in the Escherichia coli 50S subunit have commonly been inserted into an extended helix H98 in 23S rRNA, as this addition does not affect cellular growth or in vitro ribosome activity. Here, we find that E. coli 50S subunits with MS2 tags inserted in H98 are destabilized compared to wild-type (WT) 50S subunits. We identify the loss of RNA-RNA tertiary contacts that bridge helices H1, H94, and H98 as the cause of destabilization. Using cryogenic electron microscopy (cryo-EM), we show that this interaction is disrupted by the addition of the MS2 tag and can be restored through the insertion of a single adenosine in the extended H98 helix. This work establishes ways to improve MS2 tags in the 50S subunit that maintain ribosome stability and investigates a complex RNA tertiary structure that may be important for stability in various bacterial ribosomes.

The RNA interference (RNAi) pathway has evolved numerous functionalities in eukaryotes, with many on display in Kingdom Fungi. RNAi can regulate gene expression, facilitate drug resistance, or even be altogether lost to improve growth potential in some fungal pathogens. In the WHO fungal priority pathogen, Aspergillus fumigatus, the RNAi system is known to be intact and functional. To extend our limited understanding of A. fumigatus RNAi, we first investigated the genetic variation in RNAi-associated genes in a collection of 217 environmental and 83 clinical genomes, where we found that RNAi components are conserved even in clinical strains. Using endogenously expressed inverted-repeat transgenes complementary to a conditionally essential gene (pabA) or a nonessential gene (pksP), we determined that a subset of the RNAi componentry is active in inverted-repeat transgene silencing in conidia and mycelium. Analysis of mRNA-seq data from RNAi double-knockout strains linked the A. fumigatus dicer-like enzymes (DclA/B) and RNA-dependent RNA polymerases (RrpA/B) to regulation of conidial ribosome biogenesis genes; however, surprisingly few endogenous small RNAs were identified in conidia that could explain this broad change. Although RNAi was not clearly linked to growth or stress response defects in the RNAi knockouts, serial passaging of RNAi knockout strains for six generations resulted in lineages with diminished spore production over time, indicating that loss of RNAi can exert a fitness cost on the fungus. Cumulatively, A. fumigatus RNAi appears to play an active role in defense against double-stranded RNA species alongside a previously unappreciated housekeeping function in regulation of conidial ribosomal biogenesis genes.

RNA-directed transposon silencing operates in the mammalian soma and germline to safeguard genomic integrity. The piRNA pathway and the HUSH complex identify active transposons through recognition of their nascent transcripts, but mechanistic understanding of how these distinct pathways evolved is lacking. TASOR is an essential component of the HUSH complex. TASOR's DUF3715 domain adopts a pseudo-PARP structure and is required for transposon silencing in a manner independent of complex assembly. TEX15, an essential piRNA pathway factor, also contains the DUF3715 domain. Here, we show that TASOR's and TEX15's DUF3715 domain share extensive structural homology. We found that the DUF3715 domain arose in early eukaryotes and that in vertebrates it is restricted to TEX15, TASOR, and TASORB orthologs. While TASOR-like proteins are found throughout metazoa, TEX15 is vertebrate-specific. The branching of TEX15 and the TASOR-like DUF3715 domain likely occurred in early metazoan evolution. Remarkably, despite this vast evolutionary distance, the DUF3715 domain from divergent TEX15 sequences can functionally substitute the DUF3715 domain of TASOR and mediates transposon silencing. We have thus termed this domain of unknown function as the RNA-directed pseudo-PARP transposon silencing (RDTS) domain. In summary, we show an unexpected functional link between these critical transposon silencing pathways.

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

RNA-binding proteins (RBPs) are critical regulators of gene expression. An RBP typically binds to multiple mRNAs and modulates their expression. Although loss-of-function experiments on an RBP can infer how it regulates a specific target mRNA, the results are confounded by potential secondary effects due to the attenuation of all other interactions of the target RBP. For example, regarding the interaction between Trim71, an evolutionarily conserved RBP, and Ago2 mRNA, although Trim71 binds to Ago2 mRNA and overexpression of Trim71 represses Ago2 mRNA translation, it is puzzling that AGO2 protein levels are not altered in the Trim71 knockdown/knockout cells. To address this, we adapted the dTAG (degradation tag) system for determining the direct effects of the endogenous Trim71. Specifically, we knocked in the dTAG to the Trim71 locus, enabling inducible rapid Trim71 protein degradation. We observed that following the induction of Trim71 degradation, Ago2 protein levels first increased, confirming the Trim71-mediated repression, and then returned to the original levels after 24 h post-induction, revealing that the secondary effects from the Trim71 knockdown/knockout counteracted its direct effects on Ago2 mRNA. These results highlight a caveat in interpreting the results from loss-of-function studies on RBPs and provide a method to determine the primary effect(s) of RBPs on their target mRNAs.

RNA binding proteins (RBPs) perform a myriad of functions and are implicated in numerous neurological diseases. To identify the targets of RBPs in small numbers of cells, we developed TRIBE, in which the catalytic domain of the RNA editing enzyme ADAR (ADARcd) is fused to an RBP. When the RBP binds to an mRNA, ADAR catalyzes A to G modifications in the target mRNA that can be easily identified in standard RNA sequencing. In STAMP, the concept is the same except the ADARcd is replaced by the RNA editing enzyme APOBEC. Here we compared TRIBE and STAMP side-by-side in human and Drosophila cells. The goal is to learn the pros and cons of each method so that researchers can choose the method best suited to their RBP and system. In human cells, TRIBE and STAMP were performed using the RBP TDP-43. Although they both identified TDP-43 target mRNAs, combining the two methods more successfully identified high-confidence targets. In Drosophila cells, RBP-APOBEC fusions generated only low numbers of editing sites, comparable to the level of control editing. This was true for two different RBPs, Hrp48 and Thor (Drosophila EIF4E-BP), indicating that STAMP does not work well in Drosophila.

Structure probing combined with next-generation sequencing (NGS) has provided novel insights into RNA structure-function relationships. To date, such studies have focused largely on bacteria and eukaryotes, with little attention given to the third domain of life, archaea. Furthermore, functional RNAs have not been extensively studied in archaea, leaving open questions about RNA structure and function within this domain of life. With archaeal species being diverse and having many similarities to both bacteria and eukaryotes, the archaea domain has the potential to be an evolutionary bridge. In this study, we introduce a method for probing RNA structure in vivo in the archaea domain of life. We investigated the structure of ribosomal RNA (rRNA) from Methanosarcina acetivorans, a well-studied anaerobic archaeal species, grown with either methanol or acetate. After probing the RNA in vivo with dimethyl sulfate (DMS), Structure-seq2 libraries were generated, sequenced, and analyzed. We mapped the reactivity of DMS onto the secondary structure of the ribosome, which we determined independently with comparative analysis, and confirmed the accuracy of DMS probing in M. acetivorans Accessibility of the rRNA to DMS in the two carbon sources was found to be quite similar, although some differences were found. Overall, this study establishes the Structure-seq2 pipeline in the archaea domain of life and informs about ribosomal structure within M. acetivorans.

We have devised a single pot, low-cost method to add azide groups to unmodified nucleic acids without the need for enzymes or chemically modified nucleoside triphosphates. This involves reacting an azide-containing sulfinate salt with the nucleic acid, leading to replacement of C-H bonds on the nucleobase aromatic rings with C-R, where R is the azide-containing linker derived from the original sulfinate salt. With the addition of azide functional groups, the modified nucleic acid can easily be reacted with any alkyne-labeled compound of interest, including fluorescent dyes as shown in this work. This methodology enables the fluorescent labeling of a wide variety of nucleic acids, including natively folded RNAs, under mild conditions with minimal effects upon biochemical function and ribozyme catalysis. To demonstrate this, we show that a pair of labeled complementary ssDNA oligonucleotides (oligos) can hybridize to form dsDNA, even when labeled with multiple fluorophores per oligo. In addition, we also demonstrate that two different group II introns can splice when prelabeled internally with fluorophores, using our method. Broadly, this demonstrates that sulfinate modification of RNA is compatible with ribozyme function and Watson-Crick pairing, while preserving the labile backbone.

Unique chemical and physical properties are introduced by inserting selenocysteine (Sec) at specific sites within proteins. Recombinant and facile production of eukaryotic selenoproteins would benefit from a yeast expression system; however, the selenoprotein biosynthetic pathway was lost in the evolution of the kingdom Fungi as it diverged from its eukaryotic relatives. Based on our previous development of efficient selenoprotein production in bacteria, we designed a novel Sec biosynthesis pathway in Saccharomyces cerevisiae using Aeromonas salmonicida translation components. S. cerevisiae tRNASer was mutated to resemble A. salmonicida tRNASec to allow recognition by S. cerevisiae seryl-tRNA synthetase as well as A. salmonicida selenocysteine synthase (SelA) and selenophosphate synthetase (SelD). Expression of these Sec pathway components was then combined with metabolic engineering of yeast to enable the production of active methionine sulfate reductase enzyme containing genetically encoded Sec. Our report is the first demonstration that yeast is capable of selenoprotein production by site-specific incorporation of Sec.

Preadenylated single-stranded DNA ligation adaptors are essential reagents in many next generation RNA sequencing library preparation protocols. These oligonucleotides can be adenylated enzymatically or chemically. Enzymatic adenylation reactions have high yield but are not amendable to scale up. In chemical adenylation, adenosine 5'-phosphorimidazolide (ImpA) reacts with 5' phosphorylated DNA. It is easily scalable but gives poor yields, requiring labor-intensive cleanup steps. Here, we describe an improved chemical adenylation method using 95% formamide as the solvent, which results in the adenylation of oligonucleotides with >90% yield. In standard conditions, with water as the solvent, hydrolysis of the starting material to adenosine monophosphate limits the yields. To our surprise, we find that rather than increasing adenylation yields by decreasing the rate of ImpA hydrolysis, formamide does so by increasing the reaction rate between ImpA and 5'-phosphorylated DNA by ∼10-fold. The method described here enables straightforward preparation of chemically adenylated adapters with higher than 90% yield, simplifying reagent preparation for NGS.

Noncoding 6S RNAs regulate transcription by binding to the active site of bacterial RNA polymerase holoenzymes. Processing and decay of 6S-1 and 6S-2 RNA were investigated in Bacillus subtilis by northern blot and RNA-seq analyses using different RNase knockout strains, as well as by in vitro processing assays. For both 6S RNA paralogs, we identified a key-but mechanistically different-role of RNase J1. RNase J1 catalyzes 5'-end maturation of 6S-1 RNA, yet relatively inefficient and possibly via the enzyme's "sliding endonuclease" activity. 5'-end maturation has no detectable effect on 6S-1 RNA function, but rather regulates its decay: The generated 5'-monophosphate on matured 6S-1 RNA propels endonucleolytic cleavage in its apical loop region. The major 6S-2 RNA degradation pathway is initiated by endonucleolytic cleavage in the 5'-central bubble to trigger 5'-to-3'-exoribonucleolytic degradation of the downstream fragment by RNase J1. The four 3'-exonucleases of B. subtilis-RNase R, PNPase, YhaM, and particularly RNase PH-are involved in 3'-end trimming of both 6S RNAs, degradation of 6S-1 RNA fragments, and decay of abortive transcripts (so-called product RNAs, ∼14 nt in length) synthesized on 6S-1 RNA during outgrowth from stationary phase. In the case of the growth-retarded RNase Y deletion strain, we were unable to infer a specific role of RNase Y in 6S RNA decay. Yet, a participation of RNase Y in 6S RNA decay still remains possible, as evidence for such a function may have been obscured by overlapping substrate specificities of RNase Y, RNase J1, and RNase J2.

The design of new RNA sequences that retain the function of a model RNA structure is a challenge in bioinformatics because of the structural complexity of these molecules. RNA can fold into its secondary and tertiary structures by forming stem-loops and pseudoknots. A pseudoknot is a set of base pairs between a region within a stem-loop and nucleotides outside of this stem-loop; this motif is very important for numerous functional structures. It is important for any computational design algorithm to take into account these interactions to give a reliable result for any structures that include pseudoknots. In our study, we experimentally validated synthetic ribozymes designed by Enzymer, which implements algorithms allowing for the design of pseudoknots. Enzymer is a program that uses an inverse folding approach to design pseudoknotted RNAs; we used it in this study to design two types of ribozymes. The ribozymes tested were the hammerhead and the glmS, which have a self-cleaving activity that allows them to liberate the new RNA genome copy during rolling-circle replication or to control the expression of the downstream genes, respectively. We demonstrated the efficiency of Enzymer by showing that the pseudoknotted hammerhead and glmS ribozymes sequences it designed were extensively modified compared to wild-type sequences and were still active.

As one of the most prevalent RNA modifications in animals, adenosine-to-inosine (A-to-I) RNA editing facilitates the environmental adaptation of organisms by diversifying the proteome in a temporal-spatial manner. In flies and bees, the editing enzyme Adar has independently gained two different autorecoding sites that form an autofeedback loop, stabilizing the overall editing efficiency. This ensures cellular homeostasis by keeping the normal function of target genes. However, in a broader range of insects, the evolutionary dynamics and significance of this Adar autoregulatory mechanism are unclear. We retrieved the genomes of 377 arthropod species covering the five major insect orders (Hemiptera, Hymenoptera, Coleoptera, Diptera, and Lepidoptera) and aligned the Adar autorecoding sites across all genomes. We found that the two autorecoding sites underwent compensatory gains and losses during the evolution of two orders with the most sequenced species (Diptera and Hymenoptera), and that the two editing sites were mutually exclusive among them: One editable site is significantly linked to another uneditable site. This autorecoding mechanism of Adar could flexibly diversify the proteome and stabilize global editing activity. Many insects independently selected different autorecoding sites to achieve a feedback loop and regulate the global RNA editome, revealing an interesting phenomenon during evolution. Our study reveals the evolutionary force acting on accurate regulation of RNA editing activity in insects and thus deepens our understanding of the functional importance of RNA editing in environmental adaptation and evolution.

Inhibition of eukaryotic translation initiation through unscheduled RNA clamping of the DEAD-box (DDX) RNA helicases eIF4A1 and eIF4A2 has been documented for pateamine A (PatA) and rocaglates-two structurally different classes of compounds that share overlapping binding sites on eIF4A. Clamping of eIF4A to RNA causes steric blocks that interfere with ribosome binding and scanning, rationalizing the potency of these molecules since not all eIF4A molecules need to be engaged to elicit a biological effect. In addition to targeting translation, PatA and analogs have also been shown to target the eIF4A homolog, eIF4A3-a helicase necessary for exon junction complex (EJC) formation. EJCs are deposited on mRNAs upstream of exon-exon junctions and, when present downstream from premature termination codons (PTCs), participate in nonsense-mediated decay (NMD), a quality control mechanism aimed at preventing the production of dominant-negative or gain-of-function polypeptides from faulty mRNA transcripts. We find that rocaglates can also interact with eIF4A3 to induce RNA clamping. Rocaglates also inhibit EJC-dependent NMD in mammalian cells, but this does not appear to be due to induced eIF4A3-RNA clamping, but rather a secondary consequence of translation inhibition incurred by clamping eIF4A1 and eIF4A2 to mRNA.

The main goal of the origin of life (OoL) hypothesis is to reconstruct the missing link between the primordial soup and the extant biology. However, the OoL itself is just the initial part of the link representing the bootstrapping operation of Darwinian evolution. The rest of the link is the emergence of the evolution to the present day primary biological system-the ribosome-based translation apparatus. A valid hypothesis must (i) not invoke Darwinian evolution in the bootstrapping and (ii) transform the ab initio life form into the translation apparatus without violating the principle of continuity (i.e., only incremental steps without foresight). Currently, no such hypothesis exists. Here, I discuss the Quadruplex World hypothesis, which fully complies with these requirements and suggests a spontaneous emergence of the ab initio life form. The spontaneity of OoL arises from the physicochemical properties of guanine monomers in a manner of causal determinism: each step of the process (i.e., scaffolding, polymerization, and folding) is caused by the most recent past step such that in the end only the specific 3D architecture forms. The architecture (i) has a length-independent folding pattern; (ii) can play the role of the predecessor of tRNA and single-handedly conduct a primitive form of translation; and (iii) can evolve into the extant translation apparatus without any paradoxes.

Messenger RNA (mRNA) localization is an important mechanism controlling local protein synthesis. In budding yeast, asymmetric localization of transcripts such as ASH1 mRNA to the bud tip depends on the She2 RNA-binding protein. She2 assembles as a tetramer to bind RNA, but the regulation of this process as part of the mRNA locasome is still unclear. Here, we performed a phosphoproteomic analysis of She2 in vivo and identified new phosphosites, several of which are located at the dimerization or tetramerization interfaces of She2. Remarkably, phosphomimetic mutations at these residues disrupt the capacity of She2 to promote Ash1 asymmetric accumulation. A detailed analysis of one of these residues, T109, shows that a T109D mutation inhibits She2 oligomerization and its interaction with She3 and the importin-α Srp1. She2 proteins harboring the T109D mutation also display reduced expression. More importantly, this phosphomimetic mutation strongly impairs the capacity of She2 to bind RNA and disrupts ASH1 mRNA localization. These results demonstrate that the control of She2 oligomerization by phosphorylation constitutes an important regulatory step in the mRNA localization pathway.

The gene expression pathway from DNA sequence to functional protein is not as straightforward as simple depictions of the central dogma might suggest. Each step is highly regulated, with complex and only partially understood molecular mechanisms at play. Translation is one step where the "one gene-one protein" paradigm breaks down, as often a single mature eukaryotic mRNA leads to more than one protein product. One way this occurs is through translation reinitiation, in which a ribosome starts making protein from one initiation site, translates until it terminates at a stop codon, but then escapes normal recycling steps and subsequently reinitiates at a different downstream site. This process is now recognized as both important and widespread, but we are only beginning to understand the interplay of factors involved in termination, recycling, and initiation that cause reinitiation events. There appear to be several ways to subvert recycling to achieve productive reinitiation, different types of stresses or signals that trigger this process, and the mechanism may depend in part on where the event occurs in the body of an mRNA. This perspective reviews the unique characteristics and mechanisms of reinitiation events, highlights the similarities and differences between three major scenarios of reinitiation, and raises outstanding questions that are promising avenues for future research.

Assemblysomes are EDTA- and RNase-resistant ribonucleoprotein (RNP) complexes of paused ribosomes with protruding nascent polypeptide chains. They have been described in yeast and human cells for the proteasome subunit Rpt1, and the disordered amino-terminal part of the nascent chain was found to be indispensable for the accumulation of the Rpt1-RNP into assemblysomes. Motivated by this, to find other assemblysome-associated RNPs we used bioinformatics to rank subunits of Saccharomyces cerevisiae protein complexes according to their amino-terminal disorder propensity. The results revealed that gene products involved in DNA repair are enriched among the top candidates. The Sgs1 DNA helicase was chosen for experimental validation. We found that indeed nascent chains of Sgs1 form EDTA-resistant RNP condensates, assemblysomes by definition. Moreover, upon exposure to UV, SGS1 mRNA shifted from assemblysomes to polysomes, suggesting that external stimuli are regulators of assemblysome dynamics. We extended our studies to human cell lines. The BLM helicase, ortholog of yeast Sgs1, was identified upon sequencing assemblysome-associated RNAs from the MCF7 human breast cancer cell line, and mRNAs encoding DNA repair proteins were overall enriched. Using the radiation-resistant A549 cell line, we observed by transmission electron microscopy that 1,6-hexanediol, an agent known to disrupt phase-separated condensates, depletes ring ribosome structures compatible with assemblysomes from the cytoplasm of cells and makes the cells more sensitive to X-ray treatment. Taken together, these findings suggest that assemblysomes may be a component of the DNA damage response from yeast to human.

Removal of introns during pre-mRNA splicing, which is central to gene expression, initiates by base pairing of U1 snRNA with a 5' splice site (5'SS). In mammals, many introns contain weak 5'SSs that are not efficiently recognized by the canonical U1 snRNP, suggesting alternative mechanisms exist. Here, we develop a cross-linking immunoprecipitation coupled to a high-throughput sequencing method, BCLIP-seq, to identify NRDE2 (nuclear RNAi-defective 2), and CCDC174 (coiled-coil domain-containing 174) as novel RNA-binding proteins in mouse ES cells that associate with U1 snRNA and 5'SSs. Both proteins bind directly to U1 snRNA independently of canonical U1 snRNP-specific proteins, and they are required for the selection and effective processing of weak 5'SSs. Our results reveal that mammalian cells use noncanonical splicing factors bound directly to U1 snRNA to effectively select suboptimal 5'SS sequences in hundreds of genes, promoting proper splice site choice, and accurate pre-mRNA splicing.

Small noncoding RNAs are an important class of regulatory RNAs in bacteria, often regulating responses to changes in environmental conditions. OxyS is a 110 nt, stable, trans-encoded small RNA found in Escherichia coli and is induced by an increased concentration of hydrogen peroxide. OxyS has an important regulatory role in cell stress response, affecting the expression of multiple genes. In this work, we investigated the structure of OxyS and the interaction with fhlA mRNA using nuclear magnetic resonance spectroscopy, small-angle X-ray scattering, and unbiased molecular dynamics simulations. We determined the secondary structures of isolated stem-loops and confirmed their structural integrity in OxyS. Unexpectedly, stem-loop SL4 was identified in the region that was predicted to be unstructured. Three-dimensional models of OxyS demonstrate that OxyS adopts an extended structure with four solvent-exposed stem-loops, which are available for interaction with other RNAs and proteins. Furthermore, we provide evidence of base-pairing between OxyS and fhlA mRNA.

Following transcription, tRNAs undergo a series of processing and modification events to become functional adaptors in protein synthesis. Eukaryotes have also evolved intracellular transport systems whereby nucleus-encoded tRNAs may travel out and into the nucleus. In trypanosomes, nearly all tRNAs are also imported from the cytoplasm into the mitochondrion, which lacks tRNA genes. Differential subcellular localization of the cytoplasmic splicing machinery and a nuclear enzyme responsible for queuosine modification at the anticodon "wobble" position appear to be important quality control mechanisms for tRNATyr, the only intron-containing tRNA in T. brucei Since tRNA-guanine transglycosylase (TGT), the enzyme responsible for Q formation, cannot act on an intron-containing tRNA, retrograde nuclear transport is an essential step in maturation. Unlike maturation/processing pathways, the general mechanisms of tRNA stabilization and degradation in T. brucei are poorly understood. Using a combination of cellular and molecular approaches, we show that tRNATyr has an unusually short half-life. tRNATyr, and in addition tRNAAsp, also show the presence of slow-migrating bands during electrophoresis; we term these conformers: alt-tRNATyr and alt-tRNAAsp, respectively. Although we do not know the chemical or structural nature of these conformers, alt-tRNATyr has a short half-life resembling that of tRNATyr; the same is not true for alt-tRNAAsp We also show that RRP44, which is usually an exosome subunit in other organisms, is involved in tRNA degradation of the only intron-containing tRNA in T. brucei and is partly responsible for its unusually short half-life.

CLIP technologies are now widely used to study RNA-protein interactions and many data sets are now publicly available. An important first step in CLIP data exploration is the visual inspection and assessment of processed genomic data on selected genes or regions and performing comparisons: either across conditions within a particular project, or incorporating publicly available data. However, the output files produced by data processing pipelines or preprocessed files available to download from data repositories are often not suitable for direct comparison and usually need further processing. Furthermore, to derive biological insight it is usually necessary to visualize a CLIP signal alongside other data such as annotations, or orthogonal functional genomic data (e.g., RNA-seq). We have developed a simple, but powerful, command-line tool: clipplotr, which facilitates these visual comparative and integrative analyses with normalization and smoothing options for CLIP data and the ability to show these alongside reference annotation tracks and functional genomic data. These data can be supplied as input to clipplotr in a range of file formats, which will output a publication quality figure. It is written in R and can both run on a laptop computer independently or be integrated into computational workflows on a high-performance cluster. Releases, source code, and documentation are freely available at https://github.com/ulelab/clipplotr.

The potential for microRNAs (miRNAs) to regulate gene expression remains incompletely understood. DROSHA initiates the biogenesis of miRNAs while variants of Argonaute (AGO) and trinucleotide repeat containing six (TNRC6) family proteins form complexes with miRNAs to facilitate RNA recognition and gene regulation. Here we investigate the fate of miRNAs in the absence of these critical RNAi protein factors. Knockout of DROSHA expression reduces levels of some miRNAs annotated in miRBase but not others. The identity of miRNAs with reduced expression matches the identity of miRNAs previously identified by experimental approaches. The MirGeneDB resource offers the closest alignment with experimental results. In contrast, the loss of TNRC6 proteins had much smaller effects on miRNA levels. Knocking out AGO proteins, which directly contact the mature miRNA, decreased expression of the miRNAs most strongly associated with AGO2 as determined from enhanced crosslinking immunoprecipitation (AGO2-eCLIP). Evaluation of miRNA binding to endogenously expressed AGO proteins revealed that miRNA:AGO association was similar for AGO1, AGO2, AGO3, and AGO4. Our data emphasize the need to evaluate annotated miRNAs based on approximate cellular abundance, DROSHA-dependence, and physical association with AGO when forming hypotheses related to their function.

It is estimated that nearly 50% of mammalian transcripts contain at least one upstream open reading frame (uORF), which are typically one to two orders of magnitude smaller than the downstream main ORF. Most uORFs are thought to be inhibitory as they sequester the scanning ribosome, but in some cases allow for translation reinitiation. However, termination in the 5' UTR at the end of uORFs resembles premature termination that is normally sensed by the nonsense-mediated mRNA decay (NMD) pathway. Translation reinitiation has been proposed as a method for mRNAs to prevent NMD. Here, we test how uORF length influences translation reinitiation and mRNA stability in HeLa cells. Using custom 5' UTRs and uORF sequences, we show that reinitiation can occur on heterologous mRNA sequences, favors small uORFs, and is supported when initiation occurs with more initiation factors. After determining reporter mRNA half-lives in HeLa cells and mining available mRNA half-life data sets for cumulative predicted uORF length, we conclude that translation reinitiation after uORFs is not a robust method for mRNAs to prevent NMD. Together, these data suggest that the decision of whether NMD ensues after translating uORFs occurs before reinitiation in mammalian cells.

The level of each RNA species depends on the balance between its rates of production and decay. Although previous studies have measured RNA decay across the genome in tissue culture and single-celled organisms, few experiments have been performed in intact complex tissues and organs. It is therefore unclear whether the determinants of RNA decay found in cultured cells are preserved in an intact tissue, and whether they differ between neighboring cell types and are regulated during development. To address these questions, we measured RNA synthesis and decay rates genome wide via metabolic labeling of whole cultured Drosophila larval brains using 4-thiouridine. Our analysis revealed that decay rates span a range of more than 100-fold, and that RNA stability is linked to gene function, with mRNAs encoding transcription factors being much less stable than mRNAs involved in core metabolic functions. Surprisingly, among transcription factor mRNAs there was a clear demarcation between more widely used transcription factors and those that are expressed only transiently during development. mRNAs encoding transient transcription factors are among the least stable in the brain. These mRNAs are characterized by epigenetic silencing in most cell types, as shown by their enrichment with the histone modification H3K27me3. Our data suggest the presence of an mRNA destabilizing mechanism targeted to these transiently expressed transcription factors to allow their levels to be regulated rapidly with high precision. Our study also demonstrates a general method for measuring mRNA transcription and decay rates in intact organs or tissues, offering insights into the role of mRNA stability in the regulation of complex developmental programs.

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

The RNA editing enzyme adenosine deaminase acting on RNA 1 (ADAR1) is an essential regulator of the innate immune response to both cellular and viral double-stranded RNA (dsRNA). Adenosine-to-inosine (A-to-I) editing by ADAR1 modifies the sequence and structure of endogenous dsRNA and masks it from the cytoplasmic dsRNA sensor melanoma differentiation-associated protein 5 (MDA5), preventing innate immune activation. Loss-of-function mutations in ADAR are associated with rare autoinflammatory disorders including Aicardi-Goutières syndrome (AGS), defined by a constitutive systemic up-regulation of type I interferon (IFN). The murine Adar gene encodes two protein isoforms with distinct functions: ADAR1p110 is constitutively expressed and localizes to the nucleus, whereas ADAR1p150 is primarily cytoplasmic and is inducible by IFN. Recent studies have demonstrated the critical requirement for ADAR1p150 to suppress innate immune activation by self dsRNAs. However, detailed in vivo characterization of the role of ADAR1p150 during development and in adult mice is lacking. We identified a new ADAR1p150-specific knockout mouse mutant based on a single nucleotide deletion that resulted in the loss of the ADAR1p150 protein without affecting ADAR1p110 expression. The Adar1p150 -/- died embryonically at E11.5-E12.5 accompanied by cell death in the fetal liver and an activated IFN response. Somatic loss of ADAR1p150 in adults was lethal and caused rapid hematopoietic failure, demonstrating an ongoing requirement for ADAR1p150 in vivo. The generation and characterization of this mouse model demonstrates the essential role of ADAR1p150 in vivo and provides a new tool for dissecting the functional differences between ADAR1 isoforms and their physiological contributions.

Scaffold attachment factor B (SAFB) is a conserved RNA-binding protein that is essential for early mammalian development. However, the functions of SAFB in mouse embryonic stem cells (ESCs) have not been characterized. Using RNA immunoprecipitation followed by RNA-seq (RIP-seq), we examined the RNAs associated with SAFB in wild-type and SAFB/SAFB2 double-knockout ESCs. SAFB predominantly associated with introns of protein-coding genes through purine-rich motifs. The transcript most enriched in SAFB association was the lncRNA Malat1, which also contains a purine-rich region in its 5' end. Knockout of SAFB/SAFB2 led to differential expression of approximately 1000 genes associated with multiple biological processes, including apoptosis, cell division, and cell migration. Knockout of SAFB/SAFB2 also led to splicing changes in a set of genes that were largely distinct from those that exhibited changes in expression level. The spliced and nascent transcripts of many genes whose expression levels were positively regulated by SAFB also associated with high levels of SAFB, implying that SAFB binding promotes their expression. Reintroduction of SAFB into double-knockout cells restored gene expression toward wild-type levels, an effect again observable at the level of spliced and nascent transcripts. Proteomics analysis revealed a significant enrichment of nuclear speckle-associated and RS domain-containing proteins among SAFB interactors. Neither Xist nor Polycomb functions were dramatically altered in SAFB/2 knockout ESCs. Our findings suggest that among other potential functions in ESCs, SAFB promotes the expression of certain genes through its ability to bind nascent RNA.

Systematic evolution of ligands through exponential enrichment (SELEX) is widely used to identify functional nucleic acids, such as aptamers and ribozymes. Ideally, selective pressure drives the enrichment of sequences that display the function of interest (binding, catalysis, etc.). However, amplification biases from reverse transcription can overwhelm this enrichment and leave some functional sequences at a disadvantage, with cumulative effects across multiple rounds of selection. Libraries that are designed to include structural scaffolds can improve selection outcomes by sampling sequence space more strategically, but they are also susceptible to such amplification biases, particularly during reverse transcription. Therefore, we tested five reverse transcriptases (RTs)-ImProm-II, Marathon RT (MaRT), TGIRT-III, SuperScript IV (SSIV), and BST 3.0 DNA polymerase (BST)-to determine which enzymes introduced the least bias. We directly compared cDNA yield and processivity for these enzymes on RNA templates with varying degrees of structure under various reaction conditions. In these analyses, BST exhibited excellent processivity, generated large quantities of the full-length cDNA product, displayed little bias among templates with varying structure and sequence, and performed well on long, highly structured viral RNAs. Additionally, six RNA libraries containing either strong, moderate, or no incorporated structural elements were pooled and competed head-to-head in six rounds of an amplification-only selection without external selective pressure using either SSIV, ImProm-II, or BST during reverse transcription. High-throughput sequencing established that BST maintained the most neutral enrichment values, indicating low interlibrary bias over the course of six rounds, relative to SSIV and ImProm-II, and it introduced minimal mutational bias.

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

LncRNAs comprise a heterogeneous class of RNA-encoding genes typified by low expression, nuclear enrichment, high tissue-specificity, and functional diversity, but the vast majority remain uncharacterized. Here, we assembled the mouse liver noncoding transcriptome from >2000 bulk RNA-seq samples and discovered 48,261 liver-expressed lncRNAs, a majority novel. Using these lncRNAs as a single-cell transcriptomic reference set, we elucidated lncRNA dysregulation in mouse models of high fat diet-induced nonalcoholic steatohepatitis and carbon tetrachloride-induced liver fibrosis. Trajectory inference analysis revealed lncRNA zonation patterns across the liver lobule in each major liver cell population. Perturbations in lncRNA expression and zonation were common in several disease-associated liver cell types, including nonalcoholic steatohepatitis-associated macrophages, a hallmark of fatty liver disease progression, and collagen-producing myofibroblasts, a central feature of liver fibrosis. Single-cell-based gene regulatory network analysis using bigSCale2 linked individual lncRNAs to specific biological pathways, and network-essential regulatory lncRNAs with disease-associated functions were identified by their high network centrality metrics. For a subset of these lncRNAs, promoter sequences of the network-defined lncRNA target genes were significantly enriched for lncRNA triplex formation, providing independent mechanistic support for the lncRNA-target gene linkages predicted by the gene regulatory networks. These findings elucidate liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. A subset of the liver disease-associated regulatory lncRNAs identified have human orthologs and are promising candidates for biomarkers and therapeutic targets.

Rae1 is a well-conserved endoribonuclease among Gram-positive bacteria, cyanobacteria, and the chloroplasts of higher plants. We have previously shown that Rae1 cleaves the Bacillus subtilis yrzI operon mRNA in a translation-dependent manner within a short open reading frame (ORF) called S1025, encoding a 17-amino acid (aa) peptide of unknown function. Here, we map a new Rae1 cleavage site in the bmrBCD operon mRNA encoding a multidrug transporter, within an unannotated 26-aa cryptic ORF that we have named bmrX Expression of the bmrCD portion of the mRNA is ensured by an antibiotic-dependent ribosome attenuation mechanism within the upstream ORF bmrB Cleavage by Rae1 within bmrX suppresses bmrCD expression that escapes attenuation control in the absence of antibiotics. Similar to S1025, Rae1 cleavage within bmrX is both translation- and reading frame-dependent. Consistent with this, we show that translation-dependent cleavage by Rae1 promotes ribosome rescue by the tmRNA.

Riboswitches function as important translational regulators in bacteria. Comprehensive mutational analysis of transcriptional riboswitches has been used to probe the energetic intricacies of interplay between the aptamer and expression platform, but translational riboswitches have been inaccessible to massively parallel techniques. The guanidine-II (gdm-II) riboswitch is an exclusively translational class. We have integrated RelE cleavage with next-generation sequencing to quantify ligand-dependent changes in translation initiation for all single and double mutations of the Pseudomonas aeruginosa gdm-II riboswitch, a total of more than 23,000 variants. This extensive mutational analysis is consistent with the prominent features of the bioinformatic consensus. These data indicate, unexpectedly, that direct sequestration of the Shine-Dalgarno sequence is dispensable for riboswitch function. Additionally, this comprehensive data set reveals important positions not identified in previous computational and crystallographic studies. Mutations in the variable linker region stabilize alternate conformations. The double mutant data reveal the functional importance of the previously modeled P0b helix formed by the 5' and 3' tails that serves as the basis for translational control. Additional mutations to GU wobble base pairs in both P1 and P2 reveal how the apparent cooperativity of the system involves an intricate network of communication between the two binding sites. This comprehensive examination of a translational riboswitch's expression platform illuminates how the riboswitch is precisely tuned and tunable with regard to ligand sensitivity, the amplitude of expression between ON and OFF states, and the cooperativity of ligand binding.

The interactions of iron regulatory proteins (IRPs) with mRNAs containing an iron-responsive element (IRE) maintain cellular iron homeostasis and coordinate it with metabolism and possibly cellular behavior. The mRNA encoding transferrin receptor-1 (TFRC, TfR1), which is a major means of iron importation, has five IREs within its 3' UTR, and IRP interactions help maintain cytosolic iron through the protection of the TfR1 mRNA from degradation. An IRE within the 3' UTR of an mRNA splice variant encoding human cell division cycle 14A (CDC14A) has the potential to coordinate the cellular iron status with cellular behavior through a similar IRP-mediated mechanism. However, the stability of the CDC14A splice variant was reported earlier to be unaffected by the cellular iron status, which suggested that the IRE is not functional. We labeled newly synthesized mRNA in HEK293 cells with 5-ethynyl uridine and found that the stability of the CDC14A variant is responsive to iron deprivation, but there are two major differences from the regulation of TfR1 mRNA stability. First, the decay of the CDC14A mRNA does not utilize the Roquin-mediated reaction that acts on the TfR1 mRNA, indicating that there is flexibility in the degradative machinery antagonized by the IRE-IRP interactions. Second, the stabilization of the CDC14A mRNA is delayed relative to the TfR1 mRNA and does not occur until IRP binding activity has been induced. The result is consistent with a hierarchy of IRP interactions in which the maintenance of cellular iron through the stabilization of the TfR1 mRNA is initially prioritized.

Current methods for detecting unlabeled antisense oligonucleotide (ASO) drugs rely on immunohistochemistry (IHC) and/or conjugated molecules, which lack sufficient sensitivity, specificity, and resolution to fully investigate their biodistribution. Our aim was to demonstrate the qualitative and quantitative distribution of unlabeled bepirovirsen, a clinical stage ASO, in livers and kidneys of dosed mice using novel staining and imaging technologies at subcellular resolution. ASOs were detected in formalin-fixed paraffin-embedded (FFPE) and frozen tissues using an automated chromogenic in situ hybridization (ISH) assay: miRNAscope. This was then combined with immunohistochemical detection of cell lineage markers. ASO distribution in hepatocytes versus nonparenchymal cell lineages was quantified using HALO AI image analysis. To complement this, hyperspectral coherent anti-Stokes Raman scattering (HS-CARS) imaging microscopy was used to specifically detect the unique cellular Raman spectral signatures following ASO treatment. Bepirovirsen was localized primarily in nonparenchymal liver cells and proximal renal tubules. Codetection of ASO with distinct cell lineage markers of liver and kidney populations aided target cell identity facilitating quantification. Positive liver signal was quantified using HALO AI, with 12.9% of the ASO localized to the hepatocytes and 87.1% in nonparenchymal cells. HS-CARS imaging specifically detected ASO fingerprints based on the unique vibrational signatures following unlabeled ASO treatment in a totally nonperturbative manner at subcellular resolution. Together, these novel detection and imaging modalities represent a significant increase in our ability to detect unlabeled ASOs in tissues, demonstrating improved levels of specificity and resolution. These methods help us understand their underlying mechanisms of action and ultimately improve the therapeutic potential of these important drugs for treating globally significant human diseases.

Fluorescence in situ hybridization (FISH) is a widely used tool for quantifying gene expression and determining the location of RNA molecules in cells. We present an improved method for FISH probe production that yields high-purity probes with a wide range of fluorophores using standard laboratory equipment at low cost. The method modifies an earlier protocol that uses terminal deoxynucleotidyl transferase to add fluorescently labeled nucleotides to synthetic deoxyoligonucleotides. In our protocol, amino-11-ddUTP is joined to an oligonucleotide pool prior to its conjugation to a fluorescent dye, thereby generating pools of probes ready for a variety of modifications. This order of reaction steps allows for high labeling efficiencies regardless of the GC content or terminal base of the oligonucleotides. The degree of labeling (DOL) for spectrally distinct fluorophores (Quasar, ATTO, and Alexa dyes) was mostly >90%, comparable with commercial probes. The ease and low cost of production allowed the generation of probe sets targeting a wide variety of RNA molecules. Using these probes, FISH assays in C2C12 cells showed the expected subcellular localization of mRNAs and pre-mRNAs for Polr2a (RNA polymerase II subunit 2a) and Gapdh, and of the long noncoding RNAs Malat1 and Neat1 Developing FISH probe sets for several transcripts containing retained introns, we found that retained introns in the Gabbr1 and Noc2l transcripts are present in subnuclear foci separate from their sites of synthesis and partially coincident with nuclear speckles. This labeling protocol should have many applications in RNA biology.

N6-methyladenosine (m6A) in mRNA regulates almost every stage in the mRNA life cycle, and the development of methodologies for the high-throughput detection of methylated sites in mRNA using m6A-specific methylated RNA immunoprecipitation with next-generation sequencing (MeRIPSeq) or m6A individual-nucleotide-resolution cross-linking and immunoprecipitation (miCLIP) have revolutionized the m6A research field. Both of these methods are based on immunoprecipitation of fragmented mRNA. However, it is well documented that antibodies often have nonspecific activities, thus verification of identified m6A sites using an antibody-independent method would be highly desirable. We mapped and quantified the m6A site in the chicken β-actin zipcode based on the data from chicken embryo MeRIPSeq results and our RNA-Epimodification Detection and Base-Recognition (RedBaron) antibody-independent assay. We also demonstrated that methylation of this site in the β-actin zipcode enhances ZBP1 binding in vitro, while methylation of a nearby adenosine abolishes binding. This suggests that m6A may play a role in regulating localized translation of β-actin mRNA, and the ability of m6A to enhance or inhibit a reader protein's RNA binding highlights the importance of m6A detection at nucleotide resolution.

Fission yeast phosphate homeostasis gene pho1 is actively repressed during growth in phosphate-rich medium by transcription in cis of a long noncoding (lnc) RNA from the 5' flanking prt(nc-pho1) gene. Pho1 expression is: (i) derepressed by genetic maneuvers that favor precocious lncRNA 3'-processing and termination, in response to DSR and PAS signals in prt; and (ii) hyperrepressed in genetic backgrounds that dampen 3'-processing/termination efficiency. Governors of 3'-processing/termination include the RNA polymerase CTD code, the CPF (cleavage and polyadenylation factor) complex, termination factors Seb1 and Rhn1, and the inositol pyrophosphate signaling molecule 1,5-IP8 Here, we present genetic and biochemical evidence that fission yeast Duf89, a metal-dependent phosphatase/pyrophosphatase, is an antagonist of precocious 3'-processing/termination. We show that derepression of pho1 in duf89Δ cells correlates with squelching the production of full-length prt lncRNA and is erased or attenuated by: (i) DSR/PAS mutations in prt; (ii) loss-of-function mutations in components of the 3'-processing and termination machinery; (iii) elimination of the CTD Thr4-PO4 mark; (iv) interdicting CTD prolyl isomerization by Pin1; (v) inactivating the Asp1 kinase that synthesizes IP8; and (vi) loss of the putative IP8 sensor Spx1. The findings that duf89Δ is synthetically lethal with pho1-derepressive mutations CTD-S7A and aps1Δ-and that this lethality is rescued by CTD-T4A, CPF/Rhn1/Pin1 mutations, and spx1Δ-implicate Duf89 more broadly as a collaborator in cotranscriptional regulation of essential fission yeast genes. The duf89-D252A mutation, which abolishes Duf89 phosphohydrolase activity, phenocopied duf89 +, signifying that duf89Δ phenotypes are a consequence of Duf89 protein absence, not absence of Duf89 catalysis.

Under certain circumstances, any of the three termination codons can be read through by a near-cognate tRNA; i.e., a tRNA whose two out of three anticodon nucleotides base pair with those of the stop codon. Unless programed to synthetize C-terminally extended protein variants with expanded physiological roles, readthrough represents an undesirable translational error. On the other side of a coin, a significant number of human genetic diseases is associated with the introduction of nonsense mutations (premature termination codons [PTCs]) into coding sequences, where stopping is not desirable. Here, the tRNA's ability to induce readthrough opens up the intriguing possibility of mitigating the deleterious effects of PTCs on human health. In yeast, the UGA and UAR stop codons were described to be read through by four readthrough-inducing rti-tRNAs-tRNATrp and tRNACys, and tRNATyr and tRNAGln, respectively. The readthrough-inducing potential of tRNATrp and tRNATyr was also observed in human cell lines. Here, we investigated the readthrough-inducing potential of human tRNACys in the HEK293T cell line. The tRNACys family consists of two isoacceptors, one with ACA and the other with GCA anticodons. We selected nine representative tRNACys isodecoders (differing in primary sequence and expression level) and tested them using dual luciferase reporter assays. We found that at least two tRNACys can significantly elevate UGA readthrough when overexpressed. This indicates a mechanistically conserved nature of rti-tRNAs between yeast and human, supporting the idea that they could be used in the PTC-associated RNA therapies.

Initiation of translation on many viral mRNAs occurs by noncanonical mechanisms that involve 5' end-independent binding of ribosomes to an internal ribosome entry site (IRES). The ∼190-nt-long intergenic region (IGR) IRES of dicistroviruses such as cricket paralysis virus (CrPV) initiates translation without Met-tRNAi Met or initiation factors. Advances in metagenomics have revealed numerous dicistrovirus-like genomes with shorter, structurally distinct IGRs, such as nedicistrovirus (NediV) and Antarctic picorna-like virus 1 (APLV1). Like canonical IGR IRESs, the ∼165-nt-long NediV-like IGRs comprise three domains, but they lack key canonical motifs, including L1.1a/L1.1b loops (which bind to the L1 stalk of the ribosomal 60S subunit) and the apex of stem-loop V (SLV) (which binds to the head of the 40S subunit). Domain 2 consists of a compact, highly conserved pseudoknot (PKIII) that contains a UACUA loop motif and a protruding CrPV-like stem--loop SLIV. In vitro reconstitution experiments showed that NediV-like IRESs initiate translation from a non-AUG codon and form elongation-competent 80S ribosomal complexes in the absence of initiation factors and Met-tRNAi Met Unlike canonical IGR IRESs, NediV-like IRESs bind directly to the peptidyl (P) site of ribosomes leaving the aminoacyl (A) site accessible for decoding. The related structures of NediV-like IRESs and their common mechanism of action indicate that they exemplify a distinct class of IGR IRES.