RNA

5-fluorouracil (5-FU) has been used for chemotherapy for colorectal and other cancers for over 50 years. The prevailing view of its mechanism of action is inhibition of thymidine synthase leading to defects in DNA replication and repair. However, 5-FU is also incorporated into RNA causing toxicity due to defects in RNA metabolism, inhibition of pseudouridine modification, and altered ribosome function. Here, we examine the impact of 5-FU on the expression and export of small RNAs (sRNAs) into small extracellular vesicles (sEVs). Moreover, we assess the role of 5-FU in regulation of post-transcriptional sRNA modifications (PTxM) using mass spectrometry approaches. EVs are secreted by all cells and contain a variety of proteins and RNAs that can function in cell-cell communication. PTxMs on cellular and extracellular sRNAs provide yet another layer of gene regulation. We found that treatment of the colorectal cancer (CRC) cell line DLD-1 with 5-FU led to surprising differential export of miRNA snRNA, and snoRNA transcripts. Strikingly, 5-FU treatment significantly decreased the levels of pseudouridine on both cellular and secreted EV sRNAs. In contrast, 5-FU exposure led to increased levels of cellular sRNAs containing a variety of methyl-modified bases. Our results suggest that 5-FU exposure leads to altered expression, base modifications, and mislocalization of EV base-modified sRNAs.

N6-methyladenosine is one of the most abundant epitranscriptomic signatures that can affect every aspect of RNA biology, from structure and stability to intra- and intermolecular interactions. The accurate quantitative assessment of RNA stoichiometry at single-nucleotide resolution is a prerequisite to evaluate the biological significance of m6A in the context of specific RNA. We have developed a new method, termed 4-Selenothymidine 5-triphosphate reverse transcription and Ligation Assisted PCR analysis (SLAP), for quantitative and unbiased assessment of the m6A fraction on target RNA. The inclusion of thymidine triphosphate derivative during reverse transcription discourages base pair formation with m6A resulting in the reactions cessation, while maintaining normal A-T base pairing. The site-specific ligation of the resulting cDNAs with adapters, followed by amplification, generates two distinct products that reflect the modified and unmodified fraction of the analyzed RNA. These PCR products are subsequently separated by gel electrophoresis and quantified using densitometric analysis. We applied the SLAP to verify the position and assess the frequency of m6A sites present on two exemplary long non-coding RNAs. We assessed the SLAP specificity, accuracy, and sensitivity, proving the applicability of this method for the m6A analysis on less abundant transcripts. Overall, this method constitutes an extension of the birds-eye view of RNA m6A landscape provided by epitranscriptome-wide analyses by delivering quantitative assessment of modification frequency and can therefore aid the understanding of the consequences of m6A on biological processes. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/460679v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@15845aeorg.highwire.dtl.DTLVardef@46a95aorg.highwire.dtl.DTLVardef@118633dorg.highwire.dtl.DTLVardef@1b52f9_HPS_FORMAT_FIGEXP M_FIG C_FIG

Understanding the intricate roles of RNA molecules in virulence and host-pathogen interactions can provide valuable insights into combatting infections and improving human health. Although much progress has been achieved in understanding transcriptional regulation during host-pathogen interactions in diverse species, more is needed to know about the structure of pathogen RNAs. This is particularly true for fungal pathogens, including pathogenic yeasts of the Candida genus, which are the leading cause of hospital-acquired fungal infections. Deciphering the relation between RNA structure and their biology remains a significant gap. Despite advancements in transcriptional regulation studies, especially for fungal pathogens like Candida, the structural aspects of pathogenic RNAs remain understudied. Our work addresses this gap by employing genome-wide structure probing to comprehensively explore the structural landscape of mRNAs and long non-coding RNAs (lncRNAs) in the four major Candida pathogens. Specifically focusing on mRNA, we observe a robust correlation between sequence conservation and structural characteristics in orthologous transcripts, significantly when sequence identity exceeds 50%, highlighting structural feature conservation among closely related species. We investigate the impact of single nucleotide polymorphisms (SNPs) on mRNA secondary structure. SNPs within 5 untranslated regions (UTRs) tend to occur in less structured positions, suggesting structural constraints influencing transcript regulation. Furthermore, we compare the structural properties of coding regions and UTRs, noting that coding regions are generally more structured than UTRs, consistent with similar trends in other species. Additionally, we provide the first experimental characterization of lncRNA structures in Candida species. Most lncRNAs form independent subdomains, similar to human lncRNAs. Notably, we identify hairpin-like structures in lncRNAs, a feature known to be functionally significant. Comparing hairpin prevalence between lncRNAs and protein-coding genes, we find enrichment in lncRNAs across Candida species, humans, and Arabidopsis thaliana, suggesting a conserved role for these structures. In summary, our study offers valuable insights into the interplay between RNA sequence, structure, and function in Candida pathogens, with implications for gene expression regulation and potential therapeutic strategies against Candida infections.

RNA modifications have a substantial impact on tRNA function, with modifications in the anticodon loop contributing to translational fidelity and modifications in the tRNA core impacting structural stability. In bacteria, tRNA modifications are crucial for responding to stress and regulating the expression of virulence factors. Although tRNA modifications are well-characterized in a few model organisms, our knowledge of tRNA modifications in human pathogens, such as Pseudomonas aeruginosa, remains limited. Here we leveraged two orthogonal approaches to build a reference landscape of tRNA modifications in E. coli, which enabled us to identify similar modifications in P. aeruginosa. Our analysis revealed a substantial degree of conservation between the two organisms, while also uncovering potential sites of tRNA modification in P. aeruginosa tRNAs that are not present in E. coli. The mutational signature at one of these sites, position 46 of tRNAGln1(UUG) is dependent on the P. aeruginosa homolog of TapT, the enzyme responsible for the 3-(3-amino-3-carboxypropyl) uridine (acp3U) modification. Identifying which modifications are present on different tRNAs will uncover the pathways impacted by the different tRNA modifying enzymes, some of which play roles in determining virulence and pathogenicity.

The parasite Trypanosoma brucei grows as bloodstream forms in mammals, and as procyclic forms in tsetse flies. Transcription is polycistronic, all mRNAs are trans spliced, and polyadenylation sites are defined by downstream splicing signals. Expression regulation therefore depends heavily on post-transcriptional mechanisms. The RNA-binding protein DRBD18 was previously implicated in the export of some mRNAs from the nucleus in procyclic forms. It copurifies the outer ring of the nuclear pore, mRNA export factors and exon-junction-complex proteins. We show that for >200 mRNAs, DRBD18 depletion caused preferential accumulation of versions with shortened 3-untranslated regions, arising from use of polyadenylation sites that were either undetectable or rarely seen in non-depleted cells. The shortened mRNAs were often, but not always, more abundant in depleted cells than the corresponding longer versions in normal cells. Their appearance was linked to the appearance of trans spliced, polyadenylated RNAs containing only downstream 3-untranslated-region-derived sequences. Experiments with one mRNA suggested that nuclear retention alone, through depletion of MEX67, did not affect mRNA length, suggesting a specific effect of DRBD18 on processing. Since DRBD18-bound mRNAs were enriched in polypyrimidine tract motifs, and it is found in both the nucleus and the cytoplasm, we suggest that DRBD18 acts in the nucleus by binding to polypyrimidine tracts in 3-UTRs. DRBD18 binding might both prevent polypyrimidine tract recognition by splicing factors, and promote export of the bound RNAs to the cytosol.

Transfer RNAs acquire a large plethora of chemical modifications. Among those, modifications of the anticodon loop play important roles in translational fidelity and tRNA stability. Four human wobble U containing tRNAs obtain 5-methoxycarbonylmethyluridine (mcm5U34) or 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U34), which play a role in decoding. This mark involves a cascade of enzymatic activities. The last step is mediated by Alkylation Repair Homolog 8 (ALKBH8). In this study, we performed a transcriptome-wide analysis of the repertoire of ALKBH8 RNA targets. Using a combination of HITS-CLIP-seq and RIP-seq analyses, we uncover ALKBH8-bound RNAs. It targets an additional wobble U-containing tRNA tRNALys(UUU). More interestingly, the spectrum of bound RNAs includes other types of non-coding RNAs, such as C/D box snoRNAs, 7SK RNA or some miRNAs, respectively.

The COVID-19 pandemic persists despite the development of effective vaccines. As such, it remains crucial to identify new targets for antiviral therapies. The causative virus of COVID-19, SARS-CoV-2, is a positive-sense RNA virus with RNA structures that could serve as therapeutic targets. One such RNA with established function is the frameshift stimulatory element (FSE), which promotes programmed ribosomal frameshifting. To accelerate identification of additional functional RNA elements, we introduce a novel computational approach termed the Functional RNA Identification (FRID) pipeline. The guiding principle of our pipeline, which uses established component programs as well as customized component programs, is that functional RNA elements have conserved secondary and pseudoknot structures that facilitate function. To assess the presence and conservation of putative functional RNA elements in SARS-CoV-2, we compared over 6,000 SARS-CoV-2 genomic isolates. We identified 22 functional RNA elements from the SARS-CoV-2 genome, 14 of which have conserved pseudoknots and serve as potential targets for small molecule or antisense oligonucleotide therapeutics. The FRID pipeline is general and can be applied to identify pseudoknotted RNAs for targeted therapeutics in genomes or transcriptomes from any virus or organism.

The conserved endoribonuclease RNase E dominates the dynamic landscape of RNA metabolism and underpins control mediated by small regulatory RNAs in diverse bacterial species. We explored the enzymes hydrolytic mechanism, allosteric activation, and interplay with partner proteins in the multi-component RNA degradosome assembly. RNase E cleaves single-stranded RNA with preference to attack the phosphate located at the 5{square} nucleotide preceding uracil, and we corroborate key interactions that select that base. Unexpectedly, RNase E activity is impeded strongly when the recognised uracil is isomerised to 5-ribosyluracil (pseudouridine), from which we infer the detailed geometry of the hydrolytic attack process. Kinetics analyses support models for recognition of secondary structure in substrates by RNase E and for allosteric auto-regulation. The catalytic power of the enzyme is boosted when it is assembled into the multi-enzyme RNA degradosome, most likely as a consequence of substrate channeling. Our results rationalize the origins of substrate preferences of RNase E and illuminate its catalytic mechanism, supporting the roles of allosteric domain closure and cooperation with other components of the RNA degradosome complex.

Regulation of mRNA translation by eukaryotic initiation factors (eIFs) is crucial for cell survival. In humans, eIF3 stimulates translation of the JUN mRNA which encodes the transcription factor JUN, an oncogenic transcription factor involved in cell cycle progression, apoptosis, and cell proliferation. Previous studies revealed that eIF3 activates translation of the JUN mRNA by interacting with a stem loop in the 5' untranslated region (5' UTR) and with the 5' -7-methylguanosine cap structure. In addition to its interaction site with eIF3, the JUN 5' UTR is nearly one kilobase in length, and has a high degree of secondary structure, high GC content, and an upstream start codon (uAUG). This motivated us to explore the complexity of JUN mRNA translation regulation in human cells. Here we find that JUN translation is regulated in a sequence and structure-dependent manner in regions adjacent to the eIF3-interacting site in the JUN 5' UTR. Furthermore, we identify contributions of an additional initiation factor, eIF4A, in JUN regulation. We show that enhancing the interaction of eIF4A with JUN by using the compound Rocaglamide A (RocA) represses JUN translation. We also find that both the upstream AUG (uAUG) and the main AUG (mAUG) contribute to JUN translation and that they are conserved throughout vertebrates. Our results reveal additional layers of regulation for JUN translation and show the potential of JUN as a model transcript for understanding multiple interacting modes of translation regulation.

The serine and arginine-rich splicing factor SRSF1 is an evolutionarily conserved, essential pre-mRNA splicing factor. Through a global protein-RNA interaction survey we discovered SRSF1 binding sites 25-50nt upstream from hundreds of pre-miRNAs. Using primary miRNA-10b as a model we demonstrate that SRSF1 directly regulates microRNA biogenesis both in vitro and in vivo. Selective 2 hydroxyl acylation analyzed by primer extension (SHAPE) defined a structured RNA element located upstream of the precursor miRNA-10b stem loop. Our data support a model where SRSF1 promotes initial steps of microRNA biogenesis by relieving the repressive effects of cis-regulatory elements within the leader sequence.

In eukaryotes, the nucleolus is the site of ribosome biosynthesis, an essential process in all cells. While human ribosome assembly is largely evolutionarily conserved, many of the regulatory details underlying its control and function have not yet been well-defined. The nucleolar protein RSL24D1 was originally identified as a factor important for ribosome biogenesis, and as an interactor with the PeBoW complex (PES1, BOP1, WDR12) in high-throughput affinity purifications. The PeBoW complex has been shown to be required for pre-28S rRNA processing. In this study, we show that RSL24D1 depletion impairs both pre-ribosomal RNA (pre-rRNA) transcription and mature 28S rRNA production, leading to decreased protein synthesis and p53 stabilization in mammalian cells. Surprisingly, each of the PeBoW complex members is also required for pre-rRNA transcription. We also demonstrate that RSL24D1 is physically complexed with RNA polymerase I, revealing a connection between large ribosomal subunit biogenesis and rDNA transcription. These results uncover the dual role of RSL24D1 and the PeBoW complex in multiple steps of ribosome biogenesis, and provide evidence implicating large subunit biogenesis factors in pre-rRNA transcription control.

The potential for miRNAs to regulate gene expression remains controversial. DROSHA initiates the biogenesis of miRNAs while Argonaute (AGO) and TNRC6 proteins form complexes with miRNAs that recognize RNA. Here we investigate the fate of miRNAs in the absence of critical RNAi protein factors. Knockout of DROSHA expression reduced levels of some miRNAs, but not others. Knocking out AGO proteins, which directly contact the mature miRNA, decreased expression of miRNAs. Quantitative analysis indicates compensation to maintain the overall pool of AGO after knockout of AGO variants. Evaluation of miRNA binding to AGO proteins revealed that association between AGO and miRNAs was similar for AGO1 - 4. Contrary to the assumptions underlying many peer-reviewed reports, not all annotated miRNAs have equal potential as biological regulators. Cellular abundance, DROSHA dependence, and physical association with AGO must be considered when forming hypotheses related to their function. Our data prioritize sixty miRNAs - under two percent of the overall annotated miRNA repertoire - as being most likely to function as robust gene regulators. Our approach will facilitate identifying biologically active miRNAs.

Alternative pre-mRNA splicing allows one gene to encode multiple spliced messenger RNAs and, in turn, multiple proteins from a single gene transcript. This process is tightly regulated by cis elements within the pre-mRNA and trans-acting RNA binding proteins that recognize and bind to these elements, thus influencing the spliceosome assembly at adjacent splice sites. Thus, chemical modifications in either the cis-elements or trans factors or both can significantly alter splicing patterns and, thereby, the cellular proteome. Recent studies highlight that many RNA binding proteins (RBPs) are modified at multiple lysine side chains via acetylation, which neutralizes the formal positive charge and disrupts RPBs ability to participate in RNA recognition, binding and protein-protein interactions. This suggests that lysine acetylation of RPBs may be a novel mode of eukaryotic gene regulation during pre-mRNA processing. To test this, we used the well-characterized polypyrimidine tract binding protein (which is acetylated at several lysine side chains) as a model system to investigate the role of reversible RBP acetylation in regulating alternative-pre mRNA splicing. Using multiple sequence analysis, structure-based electrostatic modeling of RNA-protein interactions, and multi-site glutamine (acetyllysine mimic) and arginine (deacetyllysine mimic) mutants, we show for the first time that for a subset of PTBP1-regulated exons, acetylation at RNA-interacting lysine side chains significantly alters PTBP1 splicing activity.

Computational approaches are often used to predict regulatory RNAs in bacteria, but their success is limited to RNAs that are highly conserved across phyla, in sequence and structure. The ANTAR regulatory system consists of a family of RNAs (the ANTAR-target RNAs) that selectively recruit ANTAR proteins. This protein-RNA complex together regulates genes at the level of translation or transcriptional elongation. Despite the widespread distribution of ANTAR proteins in bacteria, their targets RNAs havent been identified in certain bacterial phyla such as actinobacteria. Here, by using a computational search model that is tuned to actinobacterial genomes, we comprehensively identify ANTAR-target RNAs in actinobacteria. These RNA motifs lie in select transcripts, often overlapping with the ribosome binding site or start codon, to regulate translation. Transcripts harboring ANTAR-target RNAs majorly encode proteins involved in the transport and metabolism of cellular metabolites like sugars, amino acids and ions; or encode transcription factors that in turn regulate diverse genes. In this report, we substantially diversify and expand the family of ANTAR RNAs across bacteria.

U6 small nuclear (sn)RNA is the shortest and most conserved snRNA in the spliceosome and forms a substantial portion of its active site. Unlike the other four spliceosomal snRNAs, which are synthesized by RNA polymerase (RNAP) II, U6 is made by RNAP III. To determine if some aspect of U6 function is incompatible with synthesis by RNAP II, we created a U6 snRNA gene with RNAP II promoter and terminator sequences. This "U6-II" gene is functional as the sole source of U6 snRNA in yeast, but its transcript is much less stable than U6 snRNA made by RNAP III. Addition of the U4 snRNA Sm protein binding site to U6-II increased its stability and led to formation of U6-II*Sm complexes. We conclude that synthesis of U6 snRNA by RNAP III is not required for its function and that U6 snRNPs containing the Sm complex can form in vivo. The ability to synthesize U6 snRNA with RNAP II relaxes sequence restraints imposed by intragenic RNAP III promoter and terminator elements and allows facile control of U6 levels via regulators of RNAP II transcription.

RNA metabolic labeling using 4-thiouridine (s4U) captures the dynamics of RNA synthesis and decay. The power of this approach is dependent on appropriate quantification of labeled and unlabeled sequencing reads, which can be compromised by the apparent loss of s4U-labeled reads in a process we refer to as dropout. Here we show that s4U-containing transcripts can be selectively lost when RNA samples are handled under sub-optimal conditions, but that this loss can be minimized using an optimized protocol. We demonstrate a second cause of dropout in nucleotide recoding and RNA sequencing (NR-seq) experiments that is computational and downstream of library preparation. NR-seq experiments involve chemically converting s4U from a uridine analog to a cytidine analog and using the apparent T-to-C mutations to identify the populations of newly synthesized RNA. We show that high levels of T-to-C mutations can prevent read alignment with some computational pipelines, but that this bias can be overcome using improved alignment pipelines. Importantly, kinetic parameter estimates are affected by dropout independent of the NR chemistry employed, and all chemistries are practically indistinguishable in bulk, short-read RNA-seq experiments. Dropout is an avoidable problem that can be identified by including unlabeled controls, and mitigated through improved sample handing and read alignment that together improve the robustness and reproducibility of NR-seq experiments.

Ribonuclease P (RNase P), a ribozyme conserved across all domains of life, is involved in the tRNA 5' maturation. RNase P recognizes precursor tRNA based on its structure, not the tRNA sequence. This feature had been exploited to engineer RNase P to selectively target and cleave any RNA as a gene inactivation strategy. Of these, a strategy called M1GS involves tethering an Escherichia coli M1 RNA to a short stretch of guide sequence (GS), which is complementary to the RNA targeted for cleavage. Despite its simplicity and versatility, M1GS tool appears to be underutilized compared to other gene inactivation strategies. Perhaps one of the reasons is that employment of the M1GS strategy requires prior knowledge about the requirements of the M1GS target sites. To facilitate its broader use, we have developed a Python script-based user-friendly bioinformatic tool built based on the requirements of M1GS to predict its target sites for any given RNA (using either DNA or RNA sequence as an input). In this study, we first demonstrate the utility of the bioinformatic tool in predicting M1GS target sites for human 28S rRNA and then we show that the customized M1GS-mediated downregulation of 28S rRNA in a human cancer cell line. We further validate the bioinformatic tool by predicting M1GS target sites for two previously targeted RNAs. Lastly, we discuss the utility of M1GS ribozyme-mediated rRNA downregulation as a potential anticancer modality in cancers where rRNAs are upregulated.

The small self-cleaving hairpin ribozyme has served as a model for RNA structure and function and has been engineered for biotechnology applications. Hairpin ribozymes were thought to be rare with only four known examples, which limited the interpretation of their biological importance and the starting sequences for engineering efforts. Recently, a bioinformatics approach identified hundreds of different RNA sequences in metatranscriptomic data that matched a novel permutation of the hairpin ribozyme. However, the self-cleavage activity of most of these sequences has not been experimentally demonstrated. Here, a high-throughput sequencing-based approach was used to evaluate the co-transcriptional self-cleavage activity of 855 different hairpin ribozymes in parallel. The results showed that nearly all sequences are very efficient self-cleaving ribozymes, and even rare nucleotides at highly conserved positions did not prevent observable ribozyme activity. The distribution of activity observed suggests that the metatranscriptomic sequences could contain random mutations from efficient wild-type ribozymes. The results further validate the bioinformatics approach that was used for ribozyme discovery and opens further questions about the biological roles of these ribozymes in the diverse environments where they were discovered.

Cap-independent translation initiation in eukaryotes involves initiation factor (eIF) binding to a transcripts 5 untranslated region (UTR). Internal-ribosome-entry-site (IRES)-like cap-independent translation initiation does not require a free 5 end for eIF binding, as eIFs recruit the ribosome to or near the start codon. For viral mRNA, recruitment usually utilizes RNA structure, such as a pseudoknot. However, for cellular mRNA cap-independent translation, no consensus RNA structures or sequences have yet been identified for eIF binding. Fibroblast-growth factor 9 (FGF-9) is a member of a subset of mRNA that are cap-independently upregulated in breast and colorectal cancer cells using this IRES-like method. Death-associated factor 5 (DAP5), an eIF4GI homolog, binds directly to the FGF-9 5 UTR to initiate translation. However, the DAP5 binding site within the FGF-9 5 UTR is unknown. Moreover, DAP5 binds to other, dissimilar 5 UTRs, some of which need a free 5 end to stimulate cap-independent translation. We propose that a particular RNA structure involving tertiary folding, rather than a conserved sequence or secondary structure, acts as a DAP5 binding site. Using SHAPE-seq, we modeled the FGF-9 5 UTR RNAs complex secondary and tertiary structure in vitro. Further, DAP5 footprinting and toeprinting experiments show DAP5s preference for one face of this structure. DAP5 binding appears to stabilize a higher-energy RNA fold that frees the 5 end to solvent and brings the start codon close to the recruited ribosome. Our findings offer a fresh perspective in the hunt for cap-independent translational enhancers. Structural, rather than sequence-specific, eIF binding sites may act as attractive chemotherapeutic targets or as dosage tools for mRNA-based therapies.

The outbreak of the COVID-19 pandemic has led to intensive studies of both the structure and replication mechanism of SARS-CoV-2. In spite of some secondary structure experiments being carried out, the 3D structure of the key function regions of the viral RNA has not yet been well understood. At the beginning of COVID-19 breakout, RNA-Puzzles community attempted to envisage the three-dimensional structure of 5'- and 3'-Un-Translated Regions (UTRs) of the SARS-CoV-2 genome. Here, we report the results of this prediction challenge, presenting the methodologies developed by six participating groups and discussing 100 RNA 3D models (60 models of 5'-UTR and 40 of 3'-UTR) predicted through applying both human experts and automated server approaches. We describe the original protocol for the reference-free comparative analysis of RNA 3D structures designed especially for this challenge. We elaborate on the deduced consensus structure and the reliability of the predicted structural motifs. All the computationally simulated models, as well as the development and the testing of computational tools dedicated to 3D structure analysis, are available for further study.