A general LC-MS-based RNA sequencing method for direct analysis of multiple-base modifications in RNA mixtures

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Nucleic acid degradation as barrier to gene delivery: a guide to understand and overcome nuclease activity

Gene therapy is on its way to revolutionize the treatment of both inherited and acquired diseases, by transferring nucleic acids to correct a disease-causing gene in the target cells of patients. In the fight against infectious diseases, mRNA-based therapeutics have proven to be a viable strategy in the recent Covid-19 pandemic. Although a growing number of gene therapies have been approved, the success rate is limited when compared to the large number of preclinical and clinical trials that have been/are being performed. In this review, we highlight some of the hurdles which gene therapies encounter after administration into the human body, with a focus on nucleic acid degradation by nucleases that are extremely abundant in mammalian organs, biological fluids as well as in subcellular compartments. We overview the available strategies to reduce the biodegradation of gene therapeutics after administration, including chemical modifications of the nucleic acids, encapsulation into vectors and co-administration with nuclease inhibitors and discuss which strategies are applied for clinically approved nucleic acid therapeutics. In the final part, we discuss the currently available methods and techniques to qualify and quantify the integrity of nucleic acids, with their own strengths and limitations.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Current technical advancements in plant epitranscriptomic studies

The growth and development of plants are the result of the interplay between the internal developmental programming and plant-environment interactions. Gene expression regulations in plants are made up of multi-level networks. In the past few years, many studies were carried out on co- and post-transcriptional RNA modifications, which, together with the RNA community, are collectively known as the "epitranscriptome." The epitranscriptomic machineries were identified and their functional impacts characterized in a broad range of physiological processes in diverse plant species. There is mounting evidence to suggest that the epitranscriptome provides an additional layer in the gene regulatory network for plant development and stress responses. In the present review, we summarized the epitranscriptomic modifications found so far in plants, including chemical modifications, RNA editing, and transcript isoforms. The various approaches to RNA modification detection were described, with special emphasis on the recent development and application potential of third-generation sequencing. The roles of epitranscriptomic changes in gene regulation during plant-environment interactions were discussed in case studies. This review aims to highlight the importance of epitranscriptomics in the study of gene regulatory networks in plants and to encourage multi-omics investigations using the recent technical advancements.

Small RNA modifications: regulatory molecules and potential applications

Small RNAs (also referred to as small noncoding RNAs, sncRNA) are defined as polymeric ribonucleic acid molecules that are less than 200 nucleotides in length and serve a variety of essential functions within cells. Small RNA species include microRNA (miRNA), PIWI-interacting RNA (piRNA), small interfering RNA (siRNA), tRNA-derived small RNA (tsRNA), etc. Current evidence suggest that small RNAs can also have diverse modifications to their nucleotide composition that affect their stability as well as their capacity for nuclear export, and these modifications are relevant to their capacity to drive molecular signaling processes relevant to biogenesis, cell proliferation and differentiation. In this review, we highlight the molecular characteristics and cellular functions of small RNA and their modifications, as well as current techniques for their reliable detection. We also discuss how small RNA modifications may be relevant to the clinical applications for the diagnosis and treatment of human health conditions such as cancer.

Oligonucleotide mapping via mass spectrometry to enable comprehensive primary structure characterization of an mRNA vaccine against SARS-CoV-2

Oligonucleotide mapping via liquid chromatography with UV detection coupled to tandem mass spectrometry (LC-UV-MS/MS) was recently developed to support development of Comirnaty, the world's first commercial mRNA vaccine which immunizes against the SARS-CoV-2 virus. Analogous to peptide mapping of therapeutic protein modalities, oligonucleotide mapping described here provides direct primary structure characterization of mRNA, through enzymatic digestion, accurate mass determinations, and optimized collisionally-induced fragmentation. Sample preparation for oligonucleotide mapping is a rapid, one-pot, one-enzyme digestion. The digest is analyzed via LC-MS/MS with an extended gradient and resulting data analysis employs semi-automated software. In a single method, oligonucleotide mapping readouts include a highly reproducible and completely annotated UV chromatogram with 100% maximum sequence coverage, and a microheterogeneity assessment of 5' terminus capping and 3' terminus poly(A)-tail length. Oligonucleotide mapping was pivotal to ensure the quality, safety, and efficacy of mRNA vaccines by providing: confirmation of construct identity and primary structure and assessment of product comparability following manufacturing process changes. More broadly, this technique may be used to directly interrogate the primary structure of RNA molecules in general.

The elucidation of the anti-inflammatory mechanism of EMO in rheumatoid arthritis through an integrative approach combining bioinformatics and experimental verification

Introduction: Emodin (EMO), a natural derivative of the anthraquinone family mainly extracted from rhubarb (Rheum palmatum), has previously been demonstrated to possess superior anti-inflammatory properties from a single target or pathway. In order to explore the underlying mechanism of action of EMO against rheumatoid arthritis (RA), a network pharmacology approach was employed. Methods: A gene expression profile from GSE55457 available from the Gene Expression Omnibus (GEO) database was used to identify the targets of EMO action. Further, single cell RNA sequencing data from GEO database of RA patients (GSE159117) were downloaded and analysed. To further investigate the anti-RA effect of EMO on MH7A cells, the expression of IL-6 and IL-1β were monitored. Finally, RNA-seq analyses were conducted on synovial fibroblasts from EMO-treated. Result: We screened the key targets of EMO against RA using network pharmacology methods, including HMGB1, STAT1, EGR1, NR3C1, EGFR, MAPK14, CASP3, CASP1, IL4, IL13, IKBKB and FN1, and their reliability was verified using ROC curve. Single-cell RNA sequencing data analysis showed that these core target proteins mainly played a role by modulating monocytes. The anti-RA effect of EMO was further verified with MH7A cells, which showed that EMO could block cell differentiation and reduce the expression of IL-6 and IL-1β. WB experiments confirmed that EMO could affect the expression of COX2, HMBG1 and the phosphorylation of p38. Finally, sequencing of synovial fibroblasts from rats treated with EMO showed consistent results with those predicted and verified, further proving the anti-inflammatory effect of EMO. Conclusion: Our research shows that EMO inhibits inflammatory response of rheumatoid arthritis (RA) by targeting HMGB1, STAT1, EGR1, NR3C1, EGFR, MAPK14, CASP3, CASP1, IL4, IL13, IKBKB, FN1 and Monocytes/macrophages.

Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection

tRNA-derived RNA fragments (tRFs) are a recently discovered family of small noncoding RNAs (sncRNAs). We previously reported that respiratory syncytial virus (RSV) infection induces functional tRFs, which are derived from a limited subset of parent tRNAs, in airway epithelial cells. Such induction is also observed in nasopharyngeal wash samples from RSV patients and correlates to RSV genome copies, suggesting a clinical significance of tRFs in RSV infection. This work also investigates whether the modification of parent tRNAs is changed by RSV to induce tRFs, using one of the most inducible tRFs as a model. We discovered that RSV infection changed the methylation modification of adenine at position 57 in tRNA glutamic acid, with a codon of CTC (tRNA-GluCTC), and the change is essential for its cleavage. AlkB homolog 1, a previously reported tRNA demethylase, appears to remove methyladenine from tRNA-GluCTC, prompting the subsequent production of tRFs from the 5'-end of tRNA-GluCTC, a regulator of RSV replication. This study demonstrates for the first time the importance of post-transcriptional modification of tRNAs in tRF biogenesis following RSV infection, providing critical insights for antiviral strategy development.

A novel method to purify small RNAs from human tissues for methylation analysis by LC-MS/MS

Methylation modification of small RNAs, including miRNA, piRNA, and tsRNA, is critical for small RNA biogenesis and biological function. Methylation of individual small RNA can be defined by liquid chromatography-coupled with mass spectrometry (LC-MS/MS). However, LC-MS/MS analysis requires a high purity of individual small RNA. Due to the difficulty of purifying specific small RNA from tissues or cells, the progress in characterizing small RNA methylation by LC-MS/MS is limited. Here, we report a novel method that can efficiently purify small RNA from human tissues for LC-MS/MS analysis. This method includes two steps: 1) pull down the target small RNA by incubating total small RNAs (18–24 nt) extracted from human tissues with a biotinylated antisense oligonucleotide of the target small RNA, followed by capturing the binding duplex of biotinylated antisense and small RNA via streptavidin magnetic beads, and 2) protect the target small RNA by pairing it with a single-strand DNA, which sequence is complementary to the target small RNA, to form a DNA/RNA hybrid double-strand, followed by sequential digestion with exonuclease I, nuclease S1, and DNase I, respectively. Furthermore, employing a mixture of four pairs of synthetic methylated and non-methylated small RNAs, we further refined this two-step method by optimizing the nuclease S1 treatment condition. With this method, we successfully purified miR-21-5p, miR-26-5p, piR-020485, and tsRNA from human lung and sperm tissue samples and analyzed their 2′-O-methylation modification at the 3′-end by LC-MS/MS.

Cloning and Sequencing Eukaryotic Small RNAs

Small RNAs are ubiquitous regulators of gene expression that participate in nearly all aspects of physiology in a wide range of organisms. There are many different classes of eukaryotic small RNAs that play regulatory roles at every level of gene expression, including transcription, RNA stability, and translation. While eukaryotic small RNAs display diverse functions across and within classes, they are generally grouped functionally based on the machinery required for their biogenesis, the effector proteins they associate with, and their molecular characteristics. The development of techniques to clone and sequence small RNAs has been critical for their identification, yet the ligation-dependent addition of RNA adapters and the use of reverse transcriptase to generate cDNA in traditional library preparation protocols can be unsuitable to detect certain small RNA subtypes. In particular, 3' or 5' chemical modifications that are characteristic of specific types of small RNAs can impede the ligation-dependent addition of RNA adapters, while internal RNA modifications can interfere with accurate reverse transcription. The inability to clone certain small RNA subtypes with traditional protocols results in an inaccurate assessment of small RNA abundance and diversity, where some RNAs appear over-represented and others are not detected. This overview aims to guide users on how to design small RNA cloning workflows in eukaryotes to more accurately capture specific small RNAs of interest. Hence, we discuss the molecular biology underlying the identification and quantitation of small RNAs, explore the limitations of commonly used protocols, and detail the alternative approaches that can be used to enrich specific small RNA classes. © 2022 Wiley Periodicals LLC.

The Epitranscriptome in miRNAs: Crosstalk, Detection, and Function in Cancer

The epitranscriptome encompasses all post-transcriptional modifications that occur on RNAs. These modifications can alter the function and regulation of their RNA targets, which, if dysregulated, result in various diseases and cancers. As with other RNAs, miRNAs are highly modified by epitranscriptomic modifications such as m⁶A methylation, 2′-O-methylation, m⁵C methylation, m⁷G methylation, polyuridine, and A-to-I editing. miRNAs are a class of small non-coding RNAs that regulates gene expression at the post-transcriptional level. miRNAs have gathered high clinical interest due to their role in disease, development, and cancer progression. Epitranscriptomic modifications alter the targeting, regulation, and biogenesis of miRNAs, increasing the complexity of miRNA regulation. In addition, emerging studies have revealed crosstalk between these modifications. In this review, we will summarize the epitranscriptomic modifications—focusing on those relevant to miRNAs—examine the recent crosstalk between these modifications, and give a perspective on how this crosstalk expands the complexity of miRNA biology.

Rapid Determination of RNA Modifications in Consensus Motifs by Nuclease Protection with Ion-Tagged Oligonucleotide Probes and Matrix-Assisted Laser Desorption Ionization Mass Spectrometry

The reversible and substoichiometric modification of RNA has recently emerged as an additional layer of translational regulation in normal biological function and disease. Modifications are often enzymatically deposited in and removed from short (~5 nt) consensus motif sequences to carefully control the translational output of the cell. Although characterization of modification occupancy at consensus motifs can be accomplished using RNA sequencing methods, these approaches are generally time-consuming and do not directly detect post-transcriptional modifications. Here, we present a nuclease protection assay coupled with matrix-assisted laser desorption ionization mass spectrometry (MALDI-MS) to rapidly characterize modifications in consensus motifs, such as GGACU, which frequently harbor N6-methyladenosine (m⁶A). While conventional nuclease protection methods rely on long (~30 nt) oligonucleotide probes that preclude the global assessment of consensus motif modification stoichiometry, we investigated a series of ion-tagged oligonucleotide (ITO) probes and found that a benzylimidazolium-functionalized ITO (ABzIM-ITO) conferred significantly improved nuclease resistance for GGACU targets. After optimizing the conditions of the nuclease protection assay, we applied the ITO and MALDI-MS-based method for determining the stoichiometry of GG(m⁶A)CU and GGACU in RNA mixtures. Overall, the ITO-based nuclease protection and MALDI-MS method constitutes a rapid and promising approach for determining modification stoichiometries of consensus motifs.

Diving into the Structural Details of In Vitro Transcribed mRNA Using Liquid Chromatography–Mass Spectrometry-Based Oligonucleotide Profiling

The production process of in vitro transcribed messenger RNA (IVT-mRNA)-based vaccines has matured in recent years, partly due to the fight against infectious diseases such as COVID-19. One key to success has been the use of modified, next to canonical, nucleotides and the efficient addition of a Cap-structure and poly A tail to the 5’ and 3’ end, respectively, of this massive biomolecule. These important features affect mRNA stability and impact translation efficiency, consequently boosting the optimization and implementation of liquid chromatography–mass spectrometry (LC–MS)-based oligonucleotide profiling methods for their characterization. This article will provide an overview of these LC–MS methods at a fundamental and application level. It will be shown how LC–MS is implemented in mRNA-based vaccine analysis to determine the capping efficiency and the poly A tail length, and how it allows, via RNA mapping, (i) to determine the mRNA sequence, (ii) to screen the fidelity of the manufactured modifications, and (iii) to identify and quantify unwanted modifications resulting from manufacturing or storage, and sequence variants resulting from mutation or transcription errors.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

Sequencing the Origins of Life

One goal of origins of life research is to understand how primitive informational and catalytic biopolymers emerged and evolved. Recently, a number of sequencing techniques have been applied to analysis of replicating and evolving primitive biopolymer systems, providing a sequence-specific and high-resolution view of primitive chemical processes. Here, we review application of sequencing techniques to analysis of synthetic and primitive nucleic acids and polypeptides. This includes next-generation sequencing of primitive polymerization and evolution processes, followed by discussion of other novel biochemical techniques that could contribute to sequence analysis of primitive biopolymer driven chemical systems. Further application of sequencing to origins of life research, perhaps as a life detection technology, could provide insight into the origin and evolution of informational and catalytic biopolymers on early Earth or elsewhere.

Direct detection of RNA modifications and structure using single-molecule nanopore sequencing

Modifications are present on many classes of RNA, including tRNA, rRNA, and mRNA. These modifications modulate diverse biological processes such as genetic recoding and mRNA export and folding. In addition, modifications can be introduced to RNA molecules using chemical probing strategies that reveal RNA structure and dynamics. Many methods exist to detect RNA modifications by short-read sequencing; however, limitations on read length inherent to short-read-based methods dissociate modifications from their native context, preventing single-molecule modification analysis. Here, we demonstrate direct RNA nanopore sequencing to detect endogenous and exogenous RNA modifications on long RNAs at the single-molecule level. We detect endogenous 2'-O-methyl and base modifications across E. coli and S. cerevisiae ribosomal RNAs as shifts in current signal and dwell times distally through interactions with the helicase motor protein. We further use the 2'-hydroxyl reactive SHAPE reagent acetylimidazole to probe RNA structure at the single-molecule level with readout by direct nanopore sequencing.

Stephenson et al. employ direct RNA nanopore sequencing to detect endogenous and exogenous modifications on single RNA molecules. The authors demonstrate detection of endogenous 2'-O-methylation (Nm) on native ribosomal RNAs, confirming known modification patterns. They describe the development of nanoSHAPE, a method that involves exogenously labeling RNA with a small-adduct-generating chemical probe that can reveal RNA structure using long-read sequencing.

Advances in mRNA 5-methylcytosine modifications: Detection, effectors, biological functions, and clinical relevance

5-methylcytosine (m⁵C) post-transcriptional modifications affect the maturation, stability, and translation of the mRNA molecule. These modifications play an important role in many physiological and pathological processes, including stress response, tumorigenesis, tumor cell migration, embryogenesis, and viral replication. Recently, there has been a better understanding of the biological implications of m⁵C modification owing to the rapid development and optimization of detection technologies, including liquid chromatography-tandem mass spectrometry (LC-MS/MS) and RNA-BisSeq. Further, predictive models (such as PEA-m⁵C, m⁵C-PseDNC, and DeepMRMP) for the identification of potential m⁵C modification sites have also emerged. In this review, we summarize the current experimental detection methods and predictive models for mRNA m⁵C modifications, focusing on their advantages and limitations. We systematically surveyed the latest research on the effectors related to mRNA m⁵C modifications and their biological functions in multiple species. Finally, we discuss the physiological effects and pathological significance of m⁵C modifications in multiple diseases, as well as their therapeutic potential, thereby providing new perspectives for disease treatment and prognosis.

Direct Immunodetection of Global A-to-I RNA Editing Activity with a Chemiluminescent Bioassay

Adenosine-to-inosine (A-to-I) editing is a conserved eukaryotic RNA modification that contributes to development, immune response, and overall cellular function. Here, we utilize Endonuclease V (EndoV), which binds specifically to inosine in RNA, to develop an EndoV-linked immunosorbency assay (EndoVLISA) as a rapid, plate-based chemiluminescent method for measuring global A-to-I editing signatures in cellular RNA. We first optimize and validate our assay with chemically synthesized oligonucleotides. We then demonstrate rapid detection of inosine content in treated cell lines, demonstrating equivalent performance against current standard RNA-seq approaches. Lastly, we deploy our EndoVLISA for profiling differential A-to-I RNA editing signatures in normal and diseased human tissue, illustrating the utility of our platform as a diagnostic bioassay. Together, the EndoVLISA method is cost-effective, straightforward, and utilizes common laboratory equipment, offering a highly accessible new approach for studying A-to-I editing. Moreover, the multi-well plate format makes this the first assay amenable for direct high-throughput quantification of A-to-I editing for applications in disease detection and drug development.

A General LC-MS-Based Method for Direct and De Novo Sequencing of RNA Mixtures Containing both Canonical and Modified Nucleotides

Mass spectrometry (MS)-based sequencing has advantages in direct sequencing of RNA, compared to cDNA-based RNA sequencing methods, as it is completely independent of enzymes and base complementarity errors in sample preparation. In addition, it allows for sequencing of different RNA modifications in a single study, rather than just one specific modification type per study. However, many technical challenges remain in de novo MS sequencing of RNA, making it difficult to MS sequence mixed RNAs or to differentiate isomeric modifications such as pseudouridine (Ψ) from uridine (U). Our recent study incorporates a two-dimensional hydrophobic end labeling strategy into MS-based sequencing (2D-HELS MS Seq) to systematically address the current challenges in MS sequencing of RNA, making it possible to directly and de novo sequence purified single RNA and mixed RNA containing both canonical and modified nucleotides. Here, we describe the method to sequence representative single-RNA and mixed-RNA oligonucleotides, each with a different sequence and/or containing modified nucleotides such as Ψ and 5-methylcytosine (m⁵C), using 2D-HELS MS Seq.

A Compositional Model to Predict the Aggregated Isotope Distribution for Average DNA and RNA Oligonucleotides

Structural modifications of DNA and RNA molecules play a pivotal role in epigenetic and posttranscriptional regulation. To characterise these modifications, more and more MS and MS/MS- based tools for the analysis of nucleic acids are being developed. To identify an oligonucleotide in a mass spectrum, it is useful to compare the obtained isotope pattern of the molecule of interest to the one that is theoretically expected based on its elemental composition. However, this is not straightforward when the identity of the molecule under investigation is unknown. Here, we present a modelling approach for the prediction of the aggregated isotope distribution of an average DNA or RNA molecule when a particular (monoisotopic) mass is available. For this purpose, a theoretical database of all possible DNA/RNA oligonucleotides up to a mass of 25 kDa is created, and the aggregated isotope distribution for the entire database of oligonucleotides is generated using the BRAIN algorithm. Since this isotope information is compositional in nature, the modelling method is based on the additive log-ratio analysis of Aitchison. As a result, a univariate weighted polynomial regression model of order 10 is fitted to predict the first 20 isotope peaks for DNA and RNA molecules. The performance of the prediction model is assessed by using a mean squared error approach and a modified Pearson's χ2 goodness-of-fit measure on experimental data. Our analysis has indicated that the variability in spectral accuracy contributed more to the errors than the approximation of the theoretical isotope distribution by our proposed average DNA/RNA model. The prediction model is implemented as an online tool. An R function can be downloaded to incorporate the method in custom analysis workflows to process mass spectral data.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

The nature of the modification at position 37 of tRNAPhe correlates with acquired taxol resistance

Acquired drug resistance is a major obstacle in cancer therapy. Recent studies revealed that reprogramming of tRNA modifications modulates cancer survival in response to chemotherapy. However, dynamic changes in tRNA modification were not elucidated. In this study, comparative analysis of the human cancer cell lines and their taxol resistant strains based on tRNA mapping was performed by using UHPLC-MS/MS. It was observed for the first time in all three cell lines that 4-demethylwyosine (imG-14) substitutes for hydroxywybutosine (OHyW) due to tRNA-wybutosine synthesizing enzyme-2 (TYW2) downregulation and becomes the predominant modification at the 37th position of tRNAphe in the taxol-resistant strains. Further analysis indicated that the increase in imG-14 levels is caused by downregulation of TYW2. The time courses of the increase in imG-14 and downregulation of TYW2 are consistent with each other as well as consistent with the time course of the development of taxol-resistance. Knockdown of TYW2 in HeLa cells caused both an accumulation of imG-14 and reduction in taxol potency. Taken together, low expression of TYW2 enzyme promotes the cancer survival and resistance to taxol therapy, implying a novel mechanism for taxol resistance. Reduction of imG-14 deposition offers an underlying rationale to overcome taxol resistance in cancer chemotherapy.

Regulation and Function of RNA Pseudouridylation in Human Cells

Recent advances in pseudouridine detection reveal a complex pseudouridine landscape that includes messenger RNA and diverse classes of noncoding RNA in human cells. The known molecular functions of pseudouridine, which include stabilizing RNA conformations and destabilizing interactions with varied RNA-binding proteins, suggest that RNA pseudouridylation could have widespread effects on RNA metabolism and gene expression. Here, we emphasize how much remains to be learned about the RNA targets of human pseudouridine synthases, their basis for recognizing distinct RNA sequences, and the mechanisms responsible for regulated RNA pseudouridylation. We also examine the roles of noncoding RNA pseudouridylation in splicing and translation and point out the potential effects of mRNA pseudouridylation on protein production, including in the context of therapeutic mRNAs.

2D-HELS MS Seq: A General LC-MS-Based Method for Direct and de novo Sequencing of RNA Mixtures with Different Nucleotide Modifications

Mass spectrometry (MS)-based sequencing approaches have been shown to be useful in direct sequencing RNA without the need for a complementary DNA (cDNA) intermediate. However, such approaches are rarely applied as a de novo RNA sequencing method, but used mainly as a tool that can assist in quality assurance for confirming known sequences of purified single-stranded RNA samples. Recently, we developed a direct RNA sequencing method by integrating a 2-dimensional mass-retention time hydrophobic end-labeling strategy into MS-based sequencing (2D-HELS MS Seq). This method is capable of accurately sequencing single RNA sequences as well as mixtures containing up to 12 distinct RNA sequences. In addition to the four canonical ribonucleotides (A, C, G, and U), the method has the capacity to sequence RNA oligonucleotides containing modified nucleotides. This is possible because the modified nucleobase either has an intrinsically unique mass that can help in its identification and its location in the RNA sequence, or can be converted into a product with a unique mass. In this study, we have used RNA, incorporating two representative modified nucleotides (pseudouridine (Ψ) and 5-methylcytosine (m⁵C)), to illustrate the application of the method for the de novo sequencing of a single RNA oligonucleotide as well as a mixture of RNA oligonucleotides, each with a different sequence and/or modified nucleotides. The procedures and protocols described here to sequence these model RNAs will be applicable to other short RNA samples (<35 nt) when using a standard high-resolution LC-MS system, and can also be used for sequence verification of modified therapeutic RNA oligonucleotides. In the future, with the development of more robust algorithms and with better instruments, this method could allow sequencing of more complex biological samples.

Direct Sequencing of tRNA by 2D-HELS-AA MS Seq Reveals Its Different Isoforms and Dynamic Base Modifications

Post-transcriptional modifications are intrinsic to RNA structure and function. However, methods to sequence RNA typically require a cDNA intermediate and are either not able to sequence these modifications or are tailored to sequence one specific nucleotide modification only. Interestingly, some of these modifications occur with <100% frequency at their particular sites, and site-specific quantification of their stoichiometries is another challenge. Here, we report a direct method for sequencing tRNA^Phe without cDNA by integrating a two-dimensional hydrophobic RNA end-labeling strategy with an anchor-based algorithm in mass spectrometry-based sequencing (2D-HELS-AA MS Seq). The entire tRNA^Phe was sequenced and the identity, location, and stoichiometry of all eleven different RNA modifications was determined, five of which were not 100% modified, including a 2'-O-methylated G (Gm) in the wobble anticodon position as well as an N², N²-dimethylguanosine (m²₂G), a 7-methylguanosine (m⁷G), a 1-methyladenosine (m¹A), and a wybutosine (Y), suggesting numerous post-transcriptional regulations in tRNA. Two truncated isoforms at the 3'-CCA tail of the tRNA^Phe (75 nt with a 3'-CC tail (80% abundance) and 74 nt with a 3'-C tail (3% abundance)) were identified in addition to the full-length 3'-CCA-tailed tRNA^Phe (76 nt, 17% abundance). We discovered a new isoform with A-G transitions/editing at the 44 and 45 positions in the tRNA^Phe variable loop, and discuss possible mechanisms related to the emergence and functions of the isoforms with these base transitions or editing. Our method revealed new isoforms, base modifications, and RNA editing as well as their stoichiometries in the tRNA that cannot be determined by current cDNA-based methods, opening new opportunities in the field of epitranscriptomics.

Naturally occurring modified ribonucleosides

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the “fifth” ribonucleotide in 1951. Since then, the ever‐increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal–Hreidarsson syndrome, Bowen–Conradi syndrome, or Williams–Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications.

This article is categorized under:

RNA Processing > RNA Editing and Modification

Ribonucleic Acid Sequence Characterization by Negative Electron Transfer Dissociation Mass Spectrometry

Modified oligonucleotides represent a promising avenue for drug development, with small interfering RNAs (siRNA) and microRNAs gaining traction in the therapeutic market. Mass spectrometry (MS)-based analysis offers many benefits for characterizing modified nucleic acids. Negative electron transfer dissociation (NETD) has proven valuable in sequencing oligonucleotide anions, particularly because it can retain modifications while generating sequence-informative fragments. We show that NETD can be successfully implemented on a widely available quadrupole-Orbitrap-linear ion trap mass spectrometer that uses a front-end glow discharge source to generate radical fluoranthene reagent cations. We characterize both unmodified and modified ribonucleic acids and present the first application of activated-ion negative electron transfer dissociation (AI-NETD) to nucleic acids. AI-NETD achieved 100% sequence coverage for both a 6-mer (5'-rGmUrArCmUrG-3') with 2'-O-methyl modifications and a 21-mer (5'-rCrArUrCrCrUrCrUrArGrArGrGrArUrArGrArArUrG-3'), the luciferase antisense siRNA. Both NETD and AI-NETD afforded complete sequence coverage of these molecules while maintaining a relatively low degree of undesired base-loss products and internal products relative to collision-based methods.

Rapid, Unbiased PRRSV Strain Detection Using MinION Direct RNA Sequencing and Bioinformatics Tools

Prompt detection and effective control of porcine reproductive and respiratory syndrome virus (PRRSV) during outbreaks is important given its immense adverse impact on the swine industry. However, the diagnostic process can be challenging due to the high genetic diversity and high mutation rate of PRRSV. A diagnostic method that can provide more detailed genetic information about pathogens is urgently needed. In this study, we evaluated the ability of Oxford Nanopore MinION direct RNA sequencing to generate a PRRSV whole genome sequence and detect and discriminate virus at the strain-level. A nearly full length PRRSV genome was successfully generated from raw sequence reads, achieving an accuracy of 96% after consensus genome generation. Direct RNA sequencing reliably detected the PRRSV strain present with an accuracy of 99.9% using as few as 5 raw sequencing reads and successfully differentiated multiple co-infecting strains present in a sample. In addition, PRRSV strain information was obtained from clinical samples containing 10⁴ to 10⁶ viral copies or more within 6 hours of sequencing. Overall, direct viral RNA sequencing followed by bioinformatic analysis proves to be a promising approach for identification of the viral strain or strains involved in clinical infections, allowing for more precise prevention and control strategies during PRRSV outbreaks.