Bidirectional Direct Sequencing of Noncanonical RNA by Two-Dimensional Analysis of Mass Chromatograms

Mass Spectrometry-Based Direct Sequencing of tRNAs De Novo and Quantitative Mapping of Multiple RNA Modifications

Despite the extensive use of next-generation sequencing (NGS) of RNA, simultaneous direct sequencing and quantitative mapping of multiple RNA nucleotide modifications remains challenging. Mass spectrometry (MS)-based sequencing can directly sequence all RNA modifications without being limited to specific ones, but it requires a perfect MS ladder that few tRNAs can provide. Here, we describe an MS ladder complementation sequencing approach (MLC-Seq) that circumvents the perfect ladder requirement, allowing de novo MS sequencing of full-length heterogeneous cellular tRNAs with multiple nucleotide modifications at single-nucleotide precision. Unlike NGS-based methods, which lose RNA modification information, MLC-Seq preserves RNA sequence diversity and modification information, revealing new detailed stoichiometric tRNA modification profiles and their changes upon treatment with the dealkylating enzyme AlkB. It can also be combined with reference sequences to provide quantitative analysis of diverse tRNAs and modifications in total tRNA samples. MLC-Seq enables systematic, quantitative, and site-specific mapping of RNA modifications, revealing the truly complete informational content of tRNA.

Advances in methods for tRNA sequencing and quantification

In the past decade tRNA sequencing (tRNA-seq) has attracted considerable attention as an important tool for the development of novel approaches to quantify highly modified tRNA species and to propel tRNA research aimed at understanding the cellular physiology and disease and development of tRNA-based therapeutics. Many methods are available to quantify tRNA abundance while accounting for modifications and tRNA charging/acylation. Advances in both library preparation methods and bioinformatic workflows have enabled developments in next-generation sequencing (NGS) workflows. Other approaches forgo NGS applications in favor of hybridization-based approaches. In this review we provide a brief comparative overview of various tRNA quantification approaches, focusing on the advantages and disadvantages of these methods, which together facilitate reliable tRNA quantification.

tRNA renovatio: Rebirth through fragmentation

tRNA function is based on unique structures that enable mRNA decoding using anticodon trinucleotides. These structures interact with specific aminoacyl-tRNA synthetases and ribosomes using 3D shape and sequence signatures. Beyond translation, tRNAs serve as versatile signaling molecules interacting with other RNAs and proteins. Through evolutionary processes, tRNA fragmentation emerges as not merely random degradation but an act of recreation, generating specific shorter molecules called tRNA-derived small RNAs (tsRNAs). These tsRNAs exploit their linear sequences and newly arranged 3D structures for unexpected biological functions, epitomizing the tRNA "renovatio" (from Latin, meaning renewal, renovation, and rebirth). Emerging methods to uncover full tRNA/tsRNA sequences and modifications, combined with techniques to study RNA structures and to integrate AI-powered predictions, will enable comprehensive investigations of tRNA fragmentation products and new interaction potentials in relation to their biological functions. We anticipate that these directions will herald a new era for understanding biological complexity and advancing pharmaceutical engineering.

Trivalent rare earth metal cofactors confer rapid NP-DNA polymerase activity

A DNA polymerase with a single mutation and a divalent calcium cofactor catalyzes the synthesis of unnatural N3'→P5' phosphoramidate (NP) bonds to form NP-DNA. However, this template-directed phosphoryl transfer activity remains orders of magnitude slower than native phosphodiester synthesis. Here, we used time-resolved x-ray crystallography to show that NP-DNA synthesis proceeds with a single detectable calcium ion in the active site. Using insights from isotopic and elemental effects, we propose that one-metal-ion electrophilic substrate activation is inferior to the native two-metal-ion mechanism. We found that this deficiency in divalent activation could be ameliorated by trivalent rare earth and post-transition metal cations, substantially enhancing NP-DNA synthesis. Scandium(III), in particular, confers highly specific NP activity with kinetics enhanced by more than 100-fold over calcium(II), yielding NP-DNA strands up to 100 nucleotides in length.

Challenges and emerging trends in liquid chromatography-based analyses of mRNA pharmaceuticals

Lipid encapsulated messenger RNA (LNP mRNA) has garnered a significant amount of interest from the pharmaceutical industry and general public alike. This attention has been catalyzed by the clinical success of LNP mRNA for SARS-CoV-2 vaccination as well as future promises that might be fulfilled by the biotechnology pipeline, such as the in vivo delivery of a CRISPR/Cas9 complex that can edit patient cells to reduce levels of low-density lipoprotein. LNP mRNAs are comprised of various chemically diverse molecules brought together in a sophisticated intermolecular complex. This can make it challenging to achieve thorough analytical characterization. Nevertheless, liquid chromatography is becoming an increasingly relied upon technique for LNP mRNA analyses. Although there have been significant advances in all types of LNP mRNA analyses, this review focuses on recent developments and the possibilities of applying anion exchange (AEX) and ion pairing reversed phase (IP-RP) liquid chromatography for intact mRNAs as well as techniques for oligo mapping analysis, 5' endcap testing and lipid compositional assays.

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.

Locating chemical modifications in RNA sequences through ribonucleases and LC-MS based analysis

Knowledge of the structural information is essential for understanding the functional details of modified RNA. Cellular non-coding RNA such as rRNA, tRNA and even viral RNAs contain a number of post-transcriptional modifications with varied degree of diversity and density. In this chapter, we discuss the use of a combination of biochemical and analytical tools such as ribonucleases and liquid chromatography coupled with mass spectrometry approaches for characterization of modified RNA. We present the protocols and alternate strategies for obtaining confident modified sequence information to facilitate the understanding of function.

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.

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.

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.

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.

Broadly applicable oligonucleotide mass spectrometry for the analysis of RNA writers and erasers in vitro

RNAs are post-transcriptionally modified by dedicated writer or eraser enzymes that add or remove specific modifications, respectively. Mass spectrometry (MS) of RNA is a useful tool to study the modification state of an oligonucleotide (ON) in a sensitive manner. Here, we developed an ion-pairing reagent free chromatography for positive ion detection of ONs by low- and high-resolution MS, which does not interfere with other types of small compound analyses done on the same instrument. We apply ON-MS to determine the ONs from an RNase T1 digest of in vitro transcribed tRNA, which are purified after ribozyme-fusion transcription by automated size exclusion chromatography. The thus produced tRNAValAAC is substrate of the human tRNA ADAT2/3 enzyme and we confirm the deamination of adenosine to inosine and the formation of tRNAValIACin vitro by ON-MS. Furthermore, low resolution ON-MS is used to monitor the demethylation of ONs containing 1-methyladenosine by bacterial AlkB in vitro. The power of high-resolution ON-MS is demonstrated by the detection and mapping of modified ONs from native total tRNA digested with RNase T1. Overall, we present an oligonucleotide MS method which is broadly applicable to monitor in vitro RNA (de-)modification processes and native RNA.

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.

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

A complete understanding of the structural and functional potential of RNA requires understanding of chemical modifications and non-canonical bases; this in turn requires advances in current sequencing methods to be able to sequence not only canonical ribonucleotides, but at the same time directly sequence these non-standard moieties. Here, we present the first direct and modification type-independent RNA sequencing method via introduction of a 2-dimensional hydrophobic end-labeling strategy into traditional mass spectrometry-based sequencing (2D HELS MS Seq) to allow de novo sequencing of RNA mixtures and enhance sample usage efficiency. Our method can directly read out the complete sequence, while identifying, locating, and quantifying base modifications accurately in both single and mixed RNA samples containing multiple different modifications at single-base resolution. Our method can also quantify stoichiometry/percentage of modified RNA versus its canonical counterpart RNA, simulating a real biological sample where modifications exist but may not be 100% at a particular site in the RNA. This method is a critical step towards fully sequencing real complex cellular RNA samples of any type and containing any modification type and can also be used in the quality control of modified therapeutic RNAs.

Synthesis of phosphoramidate-linked DNA by a modified DNA polymerase

Life on Earth depends on polymerases. These enzymes copy genetic information to produce the DNA and RNA strands at the core of the central dogma. Polymerases act by forming phosphodiester linkages to produce polynucleotide strands. While synthetic chemistry can generate a broad range of alternative genetic materials with unnatural linkages, polymerases have so far been limited to forming O-P bonds. Here, we show that, in fact, unnatural N-P bonds can also be formed by a modified DNA polymerase. This template-directed activity generates complementary strands linked by phosphoramidate (NP) esters, an alternative backbone linkage only known to exist in the laboratory. The emergence of NP-DNA polymerase activity implies the biochemical plausibility of alternative central dogmas for cellular life.

All known polymerases copy genetic material by catalyzing phosphodiester bond formation. This highly conserved activity proceeds by a common mechanism, such that incorporated nucleoside analogs terminate chain elongation if the resulting primer strand lacks a terminal hydroxyl group. Even conservatively substituted 3′-amino nucleotides generally act as chain terminators, and no enzymatic pathway for their polymerization has yet been found. Although 3′-amino nucleotides can be chemically coupled to yield stable oligonucleotides containing N3′→P5′ phosphoramidate (NP) bonds, no such internucleotide linkages are known to occur in nature. Here, we report that 3′-amino terminated primers are, in fact, slowly extended by the DNA polymerase from B. stearothermophilus in a template-directed manner. When its cofactor is Ca²⁺ rather than Mg²⁺, the reaction is fivefold faster, permitting multiple turnover NP bond formation to yield NP-DNA strands from the corresponding 3′-amino-2′,3′-dideoxynucleoside 5′-triphosphates. A single active site mutation further enhances the rate of NP-DNA synthesis by an additional 21-fold. We show that DNA-dependent NP-DNA polymerase activity depends on conserved active site residues and propose a likely mechanism for this activity based on a series of crystal structures of bound complexes. Our results significantly broaden the catalytic scope of polymerase activity and suggest the feasibility of a genetic transition between native nucleic acids and NP-DNA.

Analysis of oligonucleotides by ion-pair reversed-phase liquid chromatography coupled with positive mode electrospray ionization mass spectrometry

Oligonucleotides are usually analyzed by ion-pair reversed-phase liquid chromatography (IP-RPLC) coupled with negative mode electrospray ionization mass spectrometry (ESI-MS) due to their highly negative charged phosphodiester backbones. Herein, the signal suppression effect of triethylamine (TEA) adducts caused the ion-pair reagent TEA/hexafluoroisopropanol (HFIP) is greatly alleviated after improving the in-source energy in positive mode ESI-MS. This strategy is applied for different RNA sequencing through analyzing their formic acid hydrolysates via IP-RPLC MS. Comparing with negative ion mode, we demonstrate that IP-RPLC MS analysis in positive ion mode is more suitable for RNA sequencing with fewer contaminant interferences. Finally, simultaneous online separation and detection of oligonucleotides and protein digests are achieved in positive ion mode IP-RPLC MS analysis with little interference to each other.

Sequencing Mirror-Image DNA Chemically

The development of mirror-image biology systems faces a crucial barrier of lacking an L-DNA sequencing technique. Here, we developed a practical method for sequencing mirror-image DNA by adopting the Maxam-Gilbert sequencing approach, through which specific nucleobases in an end-labeled L-DNA are cleaved by achiral chemicals. This technique may facilitate the therapeutic application of nuclease-resistant L-aptamer drugs, and bring the vision of building an alternative, mirror-image self-replicating system closer to reality.

Enhancing the Mass Spectrometry Sensitivity for Oligonucleotide Detection by Organic Vapor Assisted Electrospray

There are two challenges in oligonucleotide detection by liquid chromatography coupled with mass spectrometry (LC-MS), the serious ion suppression effects caused by ion-pair reagents and the low detection sensitivity in positive mode MS. In this study, highly concentrated alcohol vapors were introduced into an enclosed electrospray ionization chamber, and oligonucleotides could be well detected in negative mode MS even with 100 mM triethylammonium acetate (TEAA) as an ion-pair reagent. The MS signal intensity was improved 600-fold (for standard oligonucleotide dT15) by the isopropanol vapor assisted electrospray, and effective ion-pair LC separation was feasibly coupled with high-sensitive MS detection. Then, oligonucleotides were successfully detected in positive mode MS with few adducts by propanoic acid vapor assisted electrospray. The signal intensity was enhanced more than 10-fold on average compared with adding acids into the electrospray solution. Finally, oligonucleotides and peptides or histones were simultaneously detected in MS with little interference with each other. Our strategy provides a useful alternative for investigating the biological functions of oligonucleotides.

De novo sequencing of highly modified therapeutic oligonucleotides by hydrophobic tag sequencing coupled with LC-MS

Correct sequences are prerequisite for quality control of therapeutic oligonucleotides. However, there is no definitive method available for determining sequences of highly modified therapeutic RNAs, and thereby, most of the oligonucleotides have been used clinically without direct sequence determination. In this study, we developed a novel sequencing method called 'hydrophobic tag sequencing'. Highly modified oligonucleotides are sequenced by partially digesting oligonucleotides conjugated with a 5'-hydrophobic tag, followed by liquid chromatography-mass spectrometry analysis. 5'-Hydrophobic tag-printed fragments (5'-tag degradates) can be separated in order of their molecular masses from tag-free oligonucleotides by reversed-phase liquid chromatography. As models for the sequencing, the anti-VEGF aptamer (Macugen) and the highly modified 38-mer RNA sequences were analyzed under blind conditions. Most nucleotides were identified from the molecular weight of hydrophobic 5'-tag degradates calculated from monoisotopic mass in simple full mass data. When monoisotopic mass could not be assigned, the nucleotide was estimated using the molecular weight of the most abundant mass. The sequences of Macugen and 38-mer RNA perfectly matched the theoretical sequences. The hydrophobic tag sequencing worked well to obtain simple full mass data, resulting in accurate and clear sequencing. The present study provides for the first time a de novo sequencing technology for highly modified RNAs and contributes to quality control of therapeutic oligonucleotides. Copyright © 2016 John Wiley & Sons, Ltd.

Going global: the new era of mapping modifications in RNA

The post-transcriptional modification of RNA by the addition of one or more chemical groups has been known for over 50 years. These chemical modifications, once thought to be static, are now being discovered to play key regulatory roles in gene expression. The advent of massive parallel sequencing of RNA (RNA-seq) now allows us to probe the complexity of cellular RNA and how chemically altering RNA structure expands the RNA vocabulary. Here we present an overview of the various strategies and technologies that are available to profile RNA chemical modifications at the cellular level. These strategies can be characterized as targeted and untargeted approaches: targeted strategies are developed for one single chemical modification while untargeted strategies are more broadly applicable to a range of such chemical changes. Key for all of these approaches is the ability to locate modifications within the RNA sequence. While most of these methods are built upon an RNA-Seq pipeline, alternative approaches based on mass spectrometry or conventional DNA sequencing retain value in the overall analysis process. We also look forward toward future opportunities and technologies that may expand the types of modifications that can be globally profiled. Given the ever increasing recognition that these RNA chemical modifications play important biological roles, a variety of methods, preferably orthogonal approaches, will be required to globally identify, validate and quantify RNA chemical modifications found in the transcriptome. WIREs RNA 2017, 8:e1367. doi: 10.1002/wrna.1367 For further resources related to this article, please visit the WIREs website.

Pinpointing RNA-Protein Cross-Links with Site-Specific Stable Isotope-Labeled Oligonucleotides

High affinity RNA-protein interactions are critical to cellular function, but directly identifying the determinants of binding within these complexes is often difficult. Here, we introduce a stable isotope mass labeling technique to assign specific interacting nucleotides in an oligonucleotide-protein complex by photo-cross-linking. The method relies on generating site-specific oxygen-18-labeled phosphodiester linkages in oligonucleotides, such that covalent peptide-oligonucleotide cross-link sites arising from ultraviolet irradiation can be assigned to specific sequence positions in both RNA and protein simultaneously by mass spectrometry. Using Lin28A and a let-7 pre-element RNA, we demonstrate that mass labeling permits unambiguous identification of the cross-linked sequence positions in the RNA-protein complex.