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 (m5C), using 2D-HELS MS Seq.

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 tRNAPhe 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 tRNAPhe 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 N2, N2-dimethylguanosine (m22G), a 7-methylguanosine (m7G), a 1-methyladenosine (m1A), and a wybutosine (Y), suggesting numerous post-transcriptional regulations in tRNA. Two truncated isoforms at the 3'-CCA tail of the tRNAPhe (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 tRNAPhe (76 nt, 17% abundance). We discovered a new isoform with A-G transitions/editing at the 44 and 45 positions in the tRNAPhe 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.