Nanopore Dwell Time Analysis Permits Sequencing and Conformational Assignment of Pseudouridine in SARS-CoV-2

Implications of RNA pseudouridylation for cancer biology and therapeutics: a narrative review

BACKGROUND

Pseudouridine (Ψ), a C5-glycoside isomer of uridine, stands as one of the most prevalent RNA modifications in all RNA types. Distinguishing from the C-N bond linking uridine to ribose, the link between Ψ and ribose is a C-C bond, endowing Ψ modified RNA distinct properties and functions in various biological processes. The conversion of uridine to Ψ is governed by pseudouridine synthases (PUSs). RNA pseudouridylation is implicated in cancer biology and therapeutics.

OBJECTIVES

In this review, we will summarize the methods for detecting Ψ, the process of Ψ generation, the impact of Ψ modification on RNA metabolism and gene expression, the roles of dysregulated Ψ and pseudouridine synthases in cancers, and the underlying mechanism.

METHODS

We conducted a comprehensive search of PubMed from its inception through February 2024. The search terms included "pseudouridine"; "pseudouridine synthase"; "PUS"; "dyskerin"; "cancer"; "tumor"; "carcinoma"; "malignancy"; "tumorigenesis"; "biomarker"; "prognosis" and "therapy". We included studies published in peer-reviewed journals that focused on Ψ detection, specific mechanisms involving Ψ and PUSs, and prognosis in cancer patients with high Ψ expression. We excluded studies lacking sufficient methodological details or appropriate controls.

RESULTS

Ψ has been recognized as a significant biomarker in cancer diagnosis and prognosis. Abnormal Ψ modifications mediated by various PUSs result in dysregulated RNA metabolism and impaired RNA function, promoting the development of various cancers. Overexpression of PUSs is common in cancer cells and predicts poor prognosis. PUSs inhibition arrests cell proliferation and enhances apoptosis in cancer cells, suggesting PUS-targeting cancer therapy may be a potential strategy in cancer treatment.

DISCUSSION

High Ψ levels in serum, urine, and saliva may suggest cancer, but do not specify the type, requiring additional lab markers and imaging for accurate diagnosis. Standardized detection methods are also crucial for reliable results. PUSs are linked to cancer, but more researches are needed to understand their mechanisms in different cancers. Anticancer treatments targeting PUSs are still under developed.

Nanopore signal deviations from pseudouridine modifications in RNA are sequence-specific: quantification requires dedicated synthetic controls

Chemical modifications to mRNA respond dynamically to environmental cues and are important modulators of gene expression. Nanopore direct RNA sequencing has been applied for assessing the presence of pseudouridine (ψ) modifications through basecalling errors and signal analysis. These approaches strongly depend on the sequence context around the modification, and the occupancies derived from these measurements are not quantitative. In this work, we combine direct RNA sequencing of synthetic RNAs bearing site-specific modifications and supervised machine learning models (ModQuant) to achieve near-analytical, site-specific ψ quantification. Our models demonstrate that the ionic current signal features important for accurate ψ classification are sequence dependent and encompass information extending beyond n + 2 and n - 2 nucleotides from the ψ site. This is contradictory to current models, which assume that accurate ψ classification can be achieved with signal information confined to the 5-nucleotide k-mer window (n + 2 and n - 2 nucleotides from the ψ site). We applied our models to quantitatively profile ψ occupancy in five mRNA sites in datasets from seven human cell lines, demonstrating conserved and variable sites. Our study motivates a wider pipeline that uses ground-truth RNA control sets with site-specific modifications for quantitative profiling of RNA modifications. The ModQuant pipeline and guide are freely available at https://github.com/wanunulab/ModQuant .

Characterization of mRNA therapeutics

Therapeutic messenger RNAs (mRNAs) have emerged as powerful tools in the treatment of complex diseases, especially for conditions that lack efficacious treatment. The successful application of this modality can be attributed to its ability to encode entire proteins. While the large nature of these molecules has supported their success as therapeutics, its extended size creates several analytical challenges. To further support therapeutic mRNA development and its deployment in clinical trials, appropriate methods to support their characterization must be developed. In this review, we describe current analytical methods that have been used in the characterization of RNA quality, identity, and integrity. Advantages and limitations from several analytical techniques ranging from gel electrophoresis to liquid chromatography-mass spectrometry and from shotgun sequencing to intact mass measurements are discussed. We comprehensively describe the application of analytical methods in the measurements of capping efficiency, poly A tail analysis, as well as their applicability in stability studies.

Studying the intersection of nucleoside modifications and SARS-CoV-2 RNA-dependent RNA transcription using an in vitro reconstituted system

There is growing recognition that viral RNA genomes possess enzymatically incorporated modified nucleosides. These small chemical changes are analogous to epigenomic modifications in DNA and have the potential to be similarly important modulators of viral transcription and evolution. However, the molecular level consequences of individual sites of modification remain to be broadly explored. Here we describe an in vitro assay to examine the impact of nucleoside modifications on the rate and fidelity of SARS-CoV-2 RNA transcription. Establishing the role of modified nucleotides in SARS-CoV-2 is of interest both for advancing fundamental knowledge of RNA modifications in viruses, and because modulating the modification-landscape of SARS-CoV-2 may represent a therapeutic strategy to interfere with viral RNA replication. Our approach can be used to assess the influence both of modifications present in a template RNA, as well nucleotide analog inhibitors. These methods provide a reproducible guide for generating active SARS-CoV-2 replication/transcription complexes capable of establishing how RNA modifications influence the pre-steady state rate constants of nucleotide addition by RNA-dependent RNA polymerases.

Adenosine modifications impede SARS-CoV-2 RNA-dependent RNA transcription

SARS-CoV-2, the causative virus of the COVID-19 pandemic, follows SARS and MERS as recent zoonotic coronaviruses causing severe respiratory illness and death in humans. The recurrent impact of zoonotic coronaviruses demands a better understanding of their fundamental molecular biochemistry. Nucleoside modifications, which modulate many steps of the RNA life cycle, have been found in SARS-CoV-2 RNA, although whether they confer a pro- or antiviral effect is unknown. Regardless, the viral RNA-dependent RNA polymerase will encounter these modifications as it transcribes through the viral genomic RNA. We investigated the functional consequences of nucleoside modification on the pre-steady state kinetics of SARS-CoV-2 RNA-dependent RNA transcription using an in vitro reconstituted transcription system with modified RNA templates. Our findings show that N 6-methyladenosine and 2'-O-methyladenosine modifications slow the rate of viral transcription at magnitudes specific to each modification, which has the potential to impact SARS-CoV-2 genome maintenance.

Multicellular, IVT-derived, unmodified human transcriptome for nanopore-direct RNA analysis

Nanopore direct RNA sequencing (DRS) enables measurements of RNA modifications. Modification-free transcripts are a practical and targeted control for DRS, providing a baseline measurement for canonical nucleotides within a matched and biologically-derived sequence context. However, these controls can be challenging to generate and carry nanopore-specific nuances that can impact analyses. We produced DRS datasets using modification-free transcripts from in vitro transcription of cDNA from six immortalized human cell lines. We characterized variation across cell lines and demonstrated how these may be interpreted. These data will serve as a versatile control and resource to the community for RNA modification analyses of human transcripts.

Multicellular, IVT-derived, unmodified human transcriptome for nanopore-direct RNA analysis

Nanopore direct RNA sequencing (DRS) enables measurements of RNA modifications. Modification-free transcripts are a practical and targeted control for DRS, providing a baseline measurement for canonical nucleotides within a matched and biologically derived sequence context. However, these controls can be challenging to generate and carry nanopore-specific nuances that can impact analysis. We produced DRS datasets using modification-free transcripts from in vitro transcription (IVT) of cDNA from six immortalized human cell lines. We characterized variation across cell lines and demonstrated how these may be interpreted. These data will serve as a versatile control and resource to the community for RNA modification analysis of human transcripts.

Determination of Acetylamantadine by γ-Cyclodextrin-Assisted α-HL Nanopore for Potential Cancer Prediagnosis

The high expression of Spermidine/spermine N1-acetyltransferase (SSAT-1) is an important indicator in early cancer diagnosis. Here, we developed a nanopore-based methodology with γ-cyclodextrin as an adaptor to detect and quantify acetylamantadine, the specific SSAT-1-catalyzed product from amantadine, to accordingly reflect the activity of SSAT-1. We employ γ-cyclodextrin and report that amantadine cannot cause any secondary signals in γ-cyclodextrin-assisted α-HL nanopore, while its acetylation product, acetylamantadine, does. This allows γ-cyclodextrin to practically detect acetylamantadine in the interference of excessive amantadine, superior to the previously reported β-cyclodextrin. The quantification of acetylamantadine was not interfered with even a 50-fold amantadine and displayed no interference in artificial urine sample analysis, which indicates the good feasibility of this nanopore-based methodology in painless cancer prediagnosis. In addition, the discrimination mechanism is also explored by 2-D nuclear magnetic resonance (NMR) and nanopore experiments with a series of adamantane derivatives with different hydrophilic and hydrophobic groups. We found that both the hydrophobic region matching effect and hydrophilic interactions play a synergistic effect in forming a host-guest complex to further generate the characteristic signals, which may provide insights for the subsequent design and study of drug-cyclodextrin complexes.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

New directions for Ψ and m1A decoding in mRNA: deciphering the stoichiometry and function

Over the past decade, advancements in epitranscriptomics have significantly enhanced our understanding of mRNA metabolism and its role in human development and diseases. This period has witnessed breakthroughs in sequencing technologies and the identification of key proteins involved in RNA modification processes. Alongside the well-studied m6A, Ψ and m1A have emerged as key epitranscriptomic markers. Initially identified through transcriptome-wide profiling, these modifications are now recognized for their broad impact on RNA metabolism and gene expression. In this Perspective, we focus on the detections and functions of Ψ and m1A modifications in mRNA and discuss previous discrepancies and future challenges. We summarize recent advances and highlight the latest sequencing technologies for stoichiometric detection and their mechanistic investigations for functional unveiling in mRNA as the new research directions.

Probing enzyme-dependent pseudouridylation using direct RNA sequencing to assess neuronal epitranscriptome plasticity

Chemical modifications in mRNAs such as pseudouridine (psi) can regulate gene expression, although our understanding of the functional impact of individual psi modifications, especially in neuronal cells, is limited. We apply nanopore direct RNA sequencing to investigate psi dynamics under cellular perturbations in SH-SY5Y cells. We assign sites to psi synthases using siRNA-based knockdown. A steady-state enzyme-substrate model reveals a strong correlation between psi synthase and mRNA substrate levels and psi modification frequencies. Next, we performed either differentiation or lead-exposure to SH-SY5Y cells and found that, upon lead exposure, not differentiation, the modification frequency is less dependent on enzyme levels suggesting translational control. Finally, we compared the plasticity of psi sites across cellular states and found that plastic sites can be condition-dependent or condition-independent; several of these sites fall within transcripts encoding proteins involved in neuronal processes. Our psi analysis and validation enable investigations into the dynamics and plasticity of RNA modifications.

Utilization of nanopore direct RNA sequencing to analyze viral RNA modifications

Modifications on viral RNAs (vRNAs), either genomic RNAs or RNA transcripts, have complex effects on the viral life cycle and cellular responses to viral infection. The advent of Oxford Nanopore Technologies Direct RNA Sequencing provides a new strategy for studying RNA modifications. To this end, multiple computational tools have been developed, but a systemic evaluation of their performance in mapping vRNA modifications is lacking. Here, 10 computational tools were tested using the Sindbis virus (SINV) RNAs isolated from infected mammalian (BHK-21) or mosquito (C6/36) cells, with in vitro-transcribed RNAs serving as modification-free control. Three single-mode approaches were shown to be inapplicable in the viral context, and three out of seven comparative methods required cutoff adjustments to reduce false-positive predictions. Utilizing optimized cutoffs, an integrated analysis of comparative tools suggested that the intersected predictions of Tombo_com and xPore were significantly enriched compared with the background. Consequently, a pipeline integrating Tombo_com and xPore was proposed for vRNA modification detection; the performance of which was supported by N6-methyladenosine prediction in severe acute respiratory syndrome coronavirus 2 RNAs using publicly available data. When applied to SINV RNAs, this pipeline revealed more intensive modifications in subgenomic RNAs than in genomic RNAs. Modified uridines were frequently identified, exhibiting substantive overlapping between vRNAs generated in different cell lines. On the other hand, the interpretation of other modifications remained unclear, underlining the limitations of the current computational tools despite their notable potential.IMPORTANCEComputational approaches utilizing Oxford Nanopore Technologies Direct RNA Sequencing data were almost exclusively designed to map eukaryotic epitranscriptomes. Therefore, extra caution must be exercised when using these tools to detect vRNA modifications, as in most cases, vRNA modification profiles should be regarded as unknown epitranscriptomes without prior knowledge. Here, we comprehensively evaluated the performance of 10 computational tools in detecting vRNA modification sites. All tested single-mode methods failed to differentiate native and in vitro-transcribed samples. Using optimized cutoff values, seven tested comparative tools generated very different predictions. An integrated analysis showed significant enrichment of Tombo_com and xPore predictions against the background. A pipeline for vRNA modification detection was proposed accordingly and applied to Sindbis virus RNAs. In conclusion, our study underscores the need for the careful application of computational tools to analyze viral epitranscriptomics. It also offers insights into alphaviral RNA modifications, although further validation is required.

Simultaneous nanopore profiling of mRNA m6A and pseudouridine reveals translation coordination

N6-methyladenosine (m6A) and pseudouridine (Ψ) are the two most abundant modifications in mammalian messenger RNA, but the coordination of their biological functions remains poorly understood. We develop a machine learning-based nanopore direct RNA sequencing method (NanoSPA) that simultaneously analyzes m6A and Ψ in the human transcriptome. Applying NanoSPA to polysome profiling, we reveal opposing transcriptomic co-occurrence of m6A and Ψ and synergistic, hierarchical effects of m6A and Ψ on the polysome.

Why U matters: detection and functions of pseudouridine modifications in mRNAs

The uridine modifications pseudouridine (Ψ), dihydrouridine, and 5-methyluridine are present in eukaryotic mRNAs. Many uridine-modifying enzymes are associated with human disease, underscoring the importance of uncovering the functions of uridine modifications in mRNAs. These modified uridines have chemical properties distinct from those of canonical uridines, which impact RNA structure and RNA-protein interactions. Ψ, the most abundant of these uridine modifications, is present across (pre-)mRNAs. Recent work has shown that many Ψs are present at intermediate to high stoichiometries that are likely conducive to function and at locations that are poised to influence pre-/mRNA processing. Technological innovations and mechanistic investigations are unveiling the functions of uridine modifications in pre-mRNA splicing, translation, and mRNA stability, which are discussed in this review.

Unveiling the role of PUS7-mediated pseudouridylation in host protein interactions specific for the SARS-CoV-2 RNA genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, engages in complex interactions with host cell proteins throughout its life cycle. While these interactions enable the host to recognize and inhibit viral replication, they also facilitate essential viral processes such as transcription, translation, and replication. Many aspects of these virus-host interactions remain poorly understood. Here, we employed the catRAPID algorithm and utilized the RNA-protein interaction detection coupled with mass spectrometry technology to predict and validate the host proteins that specifically bind to the highly structured 5' and 3' terminal regions of the SARS-CoV-2 RNA. Among the interactions identified, we prioritized pseudouridine synthase PUS7, which binds to both ends of the viral RNA. Using nanopore direct RNA sequencing, we discovered that the viral RNA undergoes extensive post-transcriptional modifications. Modified consensus regions for PUS7 were identified at both terminal regions of the SARS-CoV-2 RNA, including one in the viral transcription regulatory sequence leader. Collectively, our findings offer insights into host protein interactions with the SARS-CoV-2 UTRs and highlight the likely significance of pseudouridine synthases and other post-transcriptional modifications in the viral life cycle. This new knowledge enhances our understanding of virus-host dynamics and could inform the development of targeted therapeutic strategies.

Advantages and challenges associated with bisulfite-assisted nanopore direct RNA sequencing for modifications

Nanopore direct RNA sequencing is a technology that allows sequencing for epitranscriptomic modifications with the possibility of a quantitative assessment. In the present work, pseudouridine (Ψ) was sequenced with the nanopore before and after the pH 7 bisulfite reaction that yields stable ribose adducts at C1' of Ψ. The adducted sites produced greater base call errors in the form of deletion signatures compared to Ψ. Sequencing studies on E. coli rRNA and tmRNA before and after the pH 7 bisulfite reaction demonstrated that using chemically-assisted nanopore sequencing has distinct advantages for minimization of false positives and false negatives in the data. The rRNA from E. coli has 19 known U/C sequence variations that give similar base call signatures as Ψ, and therefore, are false positives when inspecting base call data; however, these sites are refractory to reacting with bisulfite as is easily observed in nanopore data. The E. coli tmRNA has a low occupancy Ψ in a pyrimidine-rich sequence context that is called a U representing a false negative; partial occupancy by Ψ is revealed after the bisulfite reaction. In a final study, 5-methylcytidine (m5C) in RNA can readily be observed after the pH 5 bisulfite reaction in which the parent C deaminates to U and the modified site does not react. This locates m5C when using bisulfite-assisted nanopore direct RNA sequencing, which is otherwise challenging to observe. The advantages and challenges of the overall approach are discussed.

Direct Nanopore Sequencing for the 17 RNA Modification Types in 36 Locations in the E. coli Ribosome Enables Monitoring of Stress-Dependent Changes

The bacterium Escherichia coli possesses 16S and 23S rRNA strands that have 36 chemical modification sites with 17 different structures. Nanopore direct RNA sequencing using a protein nanopore sensor and helicase brake, which is also a sensor, was applied to the rRNAs. Nanopore current levels, base calling profile, and helicase dwell times for the modifications relative to unmodified synthetic rRNA controls found signatures for nearly all modifications. Signatures for clustered modifications were determined by selective sequencing of writer knockout E. coli and sequencing of synthetic RNAs utilizing some custom-synthesized nucleotide triphosphates for their preparation. The knowledge of each modification's signature, apart from 5-methylcytidine, was used to determine how metabolic and cold-shock stress impact rRNA modifications. Metabolic stress resulted in either no change or a decrease, and one site increased in modification occupancy, while cold-shock stress led to either no change or a decrease. The double modification m4Cm1402 resides in 16S rRNA, and it decreased with both stressors. Using the helicase dwell time, it was determined that the N4 methyl group is lost during both stressors, and the 2'-OMe group remained. In the ribosome, this modification stabilizes binding to the mRNA codon at the P-site resulting in increased translational fidelity that is lost during stress. The E. coli genome has seven rRNA operons (rrn), and the earlier studies aligned the nanopore reads to a single operon (rrnA). Here, the reads were aligned to all seven operons to identify operon-specific changes in the 11 pseudouridines. This study demonstrates that direct sequencing for >16 different RNA modifications in a strand is achievable.

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA

Nucleophilic addition of bisulfite to pyrimidine bases has been known for a half century, and the reaction has been in use for at least a quarter of a century for identifying 5-methylcytidine in DNA. This account focuses on the chemistry of bisulfite with pseudouridine, an isomer of the RNA nucleoside uridine in which the uracil base is connected to C1' of ribose via C5 instead of N1. Pseudouridine, Ψ, is the most common nucleotide modification found in cellular RNA overall, in part due to its abundance in rRNAs and tRNAs. It has a stabilizing influence on RNA structure because N1 is now available for additional hydrogen bonding and because the heterocycle is slightly better at π stacking. The isomerization of U to Ψ in RNA strands is catalyzed by 13 different enzymes in humans and 11 in E. coli; some of these enzymes are implicated in disease states which is testament to the biological importance of pseudouridine in cells. Recently, pseudouridine came into the limelight as the key modification that, after N1 methylation, enables mRNA vaccines to be delivered efficiently into human tissue with minimal generation of a deleterious immunogenic response. Here we describe the bisulfite reaction with pseudouridine which gives rise to a chemical sequencing method to map the modified base in the epitranscriptome. Unlike the reaction with cytidine, the addition of bisulfite to Ψ leads irreversibly to form an adduct that is bypassed during cDNA synthesis by reverse transcriptases yielding a characteristic deletion signature. Although there were hints to the structure of the bisulfite adduct(s) 30 to 50 years ago, it took modern spectroscopic and computational methods to solve the mystery. Raman spectroscopy along with extensive NMR, ECD, and computational work led to the assignment of the major product as the (R) diastereomer of an oxygen adduct at C1' of a ring-opened pseudouridine. Mechanistically, this arose from a succession of conjugate addition, E2 elimination, and a [2,3] sigmatropic rearrangement, all of which are stereodefined reactions. A minor reaction with excess bisulfite led to the (S) isomer of a S-adducted SO3- group. Understanding structure and mechanism aided the design of a Ψ-specific sequencing reaction and guided attempts to improve the utility and specificity of the method. Separately, we have been investigating the use of nanopore direct RNA sequencing, a single-molecule method that directly analyzes RNA strands isolated from cells after end-ligation of adaptor sequences. By combining the electrical current and base-calling data from the nanopore with dwell-time analysis from the helicase employed to deliver RNA to the nanopore, we were able to map Ψ sites in nearly all sequence contexts. This analysis was employed to find Ψ residues in the SARS-CoV-2 vRNA, to analyze the sequence context effects of mRNA vaccine synthesis via in vitro transcription, and to evaluate the impact of stress on chemical modifications in the E. coli ribosome. Most recently, we found that bisulfite treatment of RNA leading to Ψ adducts could modulate the nanopore signal to help in mapping modifications of low occupancy.

Identification and validation of key biomarkers based on RNA methylation genes in sepsis

Background

RNA methylation is closely involved in immune regulation, but its role in sepsis remains unknown. Here, we aim to investigate the role of RNA methylation-associated genes (RMGs) in classifying and diagnosing of sepsis.

Methods

Five types of RMGs (m1A, m5C, m6Am, m7G and Ψ) were used to identify sepsis subgroups based on gene expression profile data obtained from the GEO database (GSE57065, GSE65682, and GSE95233). Unsupervised clustering analysis was used to identify distinct RNA modification subtypes. The CIBERSORT, WGCNA, GO and KEGG analysis were performed to explore immune infiltration pattern and biological function of each cluster. RF, SVM, XGB, and GLM algorithm were applied to identify the diagnostic RMGs in sepsis. Finally, the expression levels of the five key RMGs were verified by collecting PBMCs from septic patients using qRT-PCR, and their diagnostic efficacy for sepsis was verified in combination with clinical data using ROC analysis.

Results

Sepsis was divided into three subtypes (cluster 1 to 3). Cluster 1 highly expressed NSUN7 and TRMT6, with the characteristic of neutrophil activation and upregulation of MAPK signaling pathways. Cluster 2 highly expressed NSUN3, and was featured by the regulation of mRNA stability and amino acid metabolism. NSUN5 and NSUN6 were upregulated in cluster 3 which was involved in ribonucleoprotein complex biogenesis and carbohydrate metabolism pathways. In addition, we identified that five RMGs (NSUN7, NOP2, PUS1, PUS3 and FTO) could function as biomarkers for clinic diagnose of sepsis. For validation, we determined that the relative expressions of NSUN7, NOP2, PUS1 and PUS3 were upregulated, while FTO was downregulated in septic patients. The area under the ROC curve (AUC) of NSUN7, NOP2, PUS1, PUS3 and FTO was 0.828, 0.707, 0.846, 0.834 and 0.976, respectively.

Conclusions

Our study uncovered that dysregulation of RNA methylation genes (m1A, m5C, m6Am, m7G and Ψ) was closely involved in the pathogenesis of sepsis, providing new insights into the classification of sepsis endotypes. We also revealed that five hub RMGs could function as novel diagnostic biomarkers and potential targets for treatment.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair

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

No evidence for epitranscriptomic m5C modification of SARS-CoV-2, HIV and MLV viral RNA

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

Enzymatic synthesis of RNA standards for mapping and quantifying RNA modifications in sequencing analysis

Synthesis of RNA standards that contain an internal site-specific modification is important for mapping and quantification of the modified nucleotide in sequencing analysis. While RNA containing a site-specific modification can be readily synthesized by solid-state coupling for less than 100-mer nucleotides, longer RNA must be synthesized by enzymatic ligation in the presence of a DNA splint. However, long RNAs have structural heterogeneity, and those generated by in vitro transcription have 3'-end sequence heterogeneity, which together substantially reduce the yield of ligation. Here we describe a method of 3-part splint ligation that joins an in vitro transcribed left-arm RNA, an in vitro transcribed right-arm RNA, and a chemically synthesized modification-containing middle RNA, with an efficiency higher than previously reported. We report that the improved efficiency is largely attributed to the inclusion of a pair of DNA disruptors proximal to the ligation sites, and to a lesser extent to the homogeneous processing of the 3'-end of the left-arm RNA. The yields of the ligated long RNA are sufficiently high to afford purification to homogeneity for practical RNA research. We also verify the sequence accuracy at each ligation junction by nanopore sequencing.

Assessing Readability of an 8-Letter Expanded Deoxyribonucleic Acid Alphabet with Nanopores

Chemists have now synthesized new kinds of DNA that add nucleotides to the four standard nucleotides (guanine, adenine, cytosine, and thymine) found in standard Terran DNA. Such "artificially expanded genetic information systems" are today used in molecular diagnostics; to support directed evolution to create medically useful receptors, ligands, and catalysts; and to explore issues related to the early evolution of life. Further applications are limited by the inability to directly sequence DNA containing nonstandard nucleotides. Nanopore sequencing is well-suited for this purpose, as it does not require enzymatic synthesis, amplification, or nucleotide modification. Here, we take the first steps to realize nanopore sequencing of an 8-letter "hachimoji" expanded DNA alphabet by assessing its nanopore signal range using the MspA (Mycobacterium smegmatis porin A) nanopore. We find that hachimoji DNA exhibits a broader signal range in nanopore sequencing than standard DNA alone and that hachimoji single-base substitutions are distinguishable with high confidence. Because nanopore sequencing relies on a molecular motor to control the motion of DNA, we then assessed the compatibility of the Hel308 motor enzyme with nonstandard nucleotides by tracking the translocation of single Hel308 molecules along hachimoji DNA, monitoring the enzyme kinetics and premature enzyme dissociation from the DNA. We find that Hel308 is compatible with hachimoji DNA but dissociates more frequently when walking over C-glycoside nucleosides, compared to N-glycosides. C-glycocide nucleosides passing a particular site within Hel308 induce a higher likelihood of dissociation. This highlights the need to optimize nanopore sequencing motors to handle different glycosidic bonds. It may also inform designs of future alternative DNA systems that can be sequenced with existing motors and pores.

Magnipore: Prediction of differential single nucleotide changes in the Oxford Nanopore Technologies sequencing signal of SARS-CoV-2 samples

Oxford Nanopore Technologies (ONT) allows direct sequencing of ribonucleic acids (RNA) and, in addition, detection of possible RNA modifications due to deviations from the expected ONT signal. The software available so far for this purpose can only detect a small number of modifications. Alternatively, two samples can be compared for different RNA modifications. We present Magnipore, a novel tool to search for significant signal shifts between samples of Oxford Nanopore data from similar or related species. Magnipore classifies them into mutations and potential modifications. We use Magnipore to compare SARS-CoV-2 samples. Included were representatives of the early 2020s Pango lineages (n=6), samples from Pango lineages B.1.1.7 (n=2, Alpha), B.1.617.2 (n=1, Delta), and B.1.529 (n=7, Omicron). Magnipore utilizes position-wise Gaussian distribution models and a comprehensible significance threshold to find differential signals. In the case of Alpha and Delta, Magnipore identifies 55 detected mutations and 15 sites that hint at differential modifications. We predicted potential virus-variant and variant-group-specific differential modifications. Magnipore contributes to advancing RNA modification analysis in the context of viruses and virus variants.

Detection of nucleotide modifications in bacteria and bacteriophages: Strengths and limitations of current technologies and software

RNA and DNA modifications occur in eukaryotes and prokaryotes, as well as in their viruses, and serve a wide range of functions, from gene regulation to nucleic acid protection. Although the first nucleotide modification was discovered almost 100 years ago, new and unusual modifications are still being described. Nucleotide modifications have also received more attention lately because of their increased significance, but also because new sequencing approaches have eased their detection. Chiefly, third generation sequencing platforms PacBio and Nanopore offer direct detection of modified bases by measuring deviations of the signals. These unusual modifications are especially prevalent in bacteriophage genomes, the viruses of bacteria, where they mostly appear to protect DNA against degradation from host nucleases. In this Opinion article, we highlight and discuss current approaches to detect nucleotide modifications, including hardwares and softwares, and look onward to future applications, especially for studying unusual, rare, or complex genome modifications in bacteriophages. The ability to distinguish between several types of nucleotide modifications may even shed new light on metagenomic studies.

Nanopore sequencing for N1-methylpseudouridine in RNA reveals sequence-dependent discrimination of the modified nucleotide triphosphate during transcription

Direct RNA sequencing with a commercial nanopore platform was used to sequence RNA containing uridine (U), pseudouridine (Ψ) or N1-methylpseudouridine (m1Ψ) in >100 different 5-nucleotide contexts. The base calling data for Ψ or m1Ψ were similar but different from U allowing their detection. Understanding the nanopore signatures for Ψ and m1Ψ enabled a running start T7 RNA polymerase assay to study the selection of UTP versus ΨTP or m1ΨTP competing mixtures in all possible adjacent sequence contexts. A significant sequence context dependency was observed for T7 RNA polymerase with insertion yields for ΨTP versus UTP spanning a range of 20-65%, and m1ΨTP versus UTP producing variable yields that differ by 15-70%. Experiments with SP6 RNA polymerase, as well as chemically-modified triphosphates and DNA templates provide insight to explain the observations. The SP6 polymerase introduced m1ΨTP when competed with UTP with a smaller window of yields (15-30%) across all sequence contexts studied. These results may aid in future efforts that employ RNA polymerases to make therapeutic mRNAs with sub-stoichiometric amounts of m1Ψ.

Semi-quantitative detection of pseudouridine modifications and type I/II hypermodifications in human mRNAs using direct long-read sequencing

Here, we develop and apply a semi-quantitative method for the high-confidence identification of pseudouridylated sites on mammalian mRNAs via direct long-read nanopore sequencing. A comparative analysis of a modification-free transcriptome reveals that the depth of coverage and specific k-mer sequences are critical parameters for accurate basecalling. By adjusting these parameters for high-confidence U-to-C basecalling errors, we identify many known sites of pseudouridylation and uncover previously unreported uridine-modified sites, many of which fall in k-mers that are known targets of pseudouridine synthases. Identified sites are validated using 1000-mer synthetic RNA controls bearing a single pseudouridine in the center position, demonstrating systematic under-calling using our approach. We identify mRNAs with up to 7 unique modification sites. Our workflow allows direct detection of low-, medium-, and high-occupancy pseudouridine modifications on native RNA molecules from nanopore sequencing data and multiple modifications on the same strand.

Effects of natural RNA modifications on the activity of SARS-CoV-2 RNA-dependent RNA polymerase

RNA-dependent RNA polymerase (RdRp) plays a key role in the replication of RNA viruses, including SARS-CoV-2. Processive RNA synthesis by RdRp is crucial for successful genome replication and expression, especially in the case of very long coronaviral genomes. Here, we analysed the activity of SARS-CoV-2 RdRp (the nsp12-nsp7-nsp8 complex) on synthetic primer-templates of various structures, including substrates with mismatched primers or template RNA modifications. It has been shown that RdRp cannot efficiently extend RNA primers containing mismatches and has no intrinsic RNA cleavage activity to remove the primer 3'-end, thus necessitating the action of exoribonuclease for proofreading. Similar to DNA-dependent RNA polymerases, RdRp can perform processive pyrophosphorolysis of the nascent RNA product but this reaction is also blocked in the presence of mismatches. Furthermore, we have demonstrated that several natural post-transcriptional modifications in the RNA template, which do not prevent complementary interactions (N6-methyladenosine, 5-methylcytosine, inosine and pseudouridine), do not change RdRp processivity. At the same time, certain modifications of RNA bases and ribose residues strongly block RNA synthesis, either prior to nucleotide incorporation (3-methyluridine and 1-methylguanosine) or immediately after it (2'-O-methylation). The results demonstrate that the activity of SARS-CoV-2 RdRp can be strongly inhibited by common modifications of the RNA template suggesting a way to design novel antiviral compounds.

Chemical Probe-Based Nanopore Sequencing to Selectively Assess the RNA Modifications

Nanopore direct RNA sequencing (dRNA-Seq) reads reveal RNA modifications through consistent error profiles specific to a modified nucleobase. However, a null data set is required to identify actual RNA modification-associated errors for distinguishing it from confounding highly intrinsic sequencing errors. Here, we reveal that inosine creates a signature mismatch error in dRNA-Seq reads and obviates the need for a null data set by harnessing the selective reactivity of acrylonitrile for validating the presence of actual inosine modifications. Selective reactivity of acrylonitrile toward inosine altered multiple dRNA-Seq parameters like signal intensity and trace value. We also deduced the stoichiometry of inosine modification through deviation in signal intensity and trace value using this chemical biology approach. Furthermore, we devised Nano ICE-Seq, a protocol to overcome the low coverage issue associated with direct RNA sequencing. Taken together, our chemical probe-based approach may facilitate the knockout-free detection of disease-associated RNA modifications in clinical scenarios.

Nanopore Detection Assisted DNA Information Processing

The deoxyribonucleotide (DNA) molecule is a stable carrier for large amounts of genetic information and provides an ideal storage medium for next-generation information processing technologies. Technologies that process DNA information, representing a cross-disciplinary integration of biology and computer techniques, have become attractive substitutes for technologies that process electronic information alone. The detailed applications of DNA technologies can be divided into three components: storage, computing, and self-assembly. The quality of DNA information processing relies on the accuracy of DNA reading. Nanopore detection allows researchers to accurately sequence nucleotides and is thus widely used to read DNA. In this paper, we introduce the principles and development history of nanopore detection and conduct a systematic review of recent developments and specific applications in DNA information processing involving nanopore detection and nanopore-based storage. We also discuss the potential of artificial intelligence in nanopore detection and DNA information processing. This work not only provides new avenues for future nanopore detection development, but also offers a foundation for the construction of more advanced DNA information processing technologies.

Identification of nucleoside monophosphates and their epigenetic modifications using an engineered nanopore

RNA modifications play critical roles in the regulation of various biological processes and are associated with many human diseases. Direct identification of RNA modifications by sequencing remains challenging, however. Nanopore sequencing is promising, but the current strategy is complicated by sequence decoding. Sequential nanopore identification of enzymatically cleaved nucleoside monophosphates may simultaneously provide accurate sequence and modification information. Here we show a phenylboronic acid-modified hetero-octameric Mycobacterium smegmatis porin A nanopore, with which direct distinguishing between monophosphates of canonical nucleosides, 5-methylcytidine, N⁶-methyladenosine, N⁷-methylguanosine, N¹-methyladenosine, inosine, pseudouridine and dihydrouridine was achieved. A custom machine learning algorithm, which reports an accuracy of 0.996, was also applied to the quantitative analysis of modifications in microRNA and natural transfer RNA. It is generally suitable for sensing of a variety of other nucleoside or nucleotide derivatives and may bring new insights to epigenetic RNA sequencing.

Analyzing viral epitranscriptomes using nanopore direct RNA sequencing

RNA modifications are a common occurrence across all domains of life. Several chemical modifications, including N6-methyladenosine, have also been found in viral transcripts and viral RNA genomes. Some of the modifications increase the viral replication efficiency while also helping the virus to evade the host immune system. Nonetheless, there are numerous examples in which the host's RNA modification enzymes function as antiviral factors. Although established methods like MeRIP-seq and miCLIP can provide a transcriptome- wide overview of how viral RNA is modified, it is difficult to distinguish between the complex overlapping viral transcript isoforms using the short read-based techniques. Nanopore direct RNA sequencing (DRS) provides both long reads and direct signal readings, which may carry information about the modifications. Here, we describe a refined protocol for analyzing the RNA modifications in viral transcriptomes using nanopore technology.

Role of Epitranscriptomic and Epigenetic Modifications during the Lytic and Latent Phases of Herpesvirus Infections

Herpesviruses are double-stranded DNA viruses occurring at a high prevalence in the human population and are responsible for a wide array of clinical manifestations and diseases, from mild to severe. These viruses are classified in three subfamilies (Alpha-, Beta- and Gammaherpesvirinae), with eight members currently known to infect humans. Importantly, all herpesviruses can establish lifelong latent infections with symptomatic or asymptomatic lytic reactivations. Accumulating evidence suggest that chemical modifications of viral RNA and DNA during the lytic and latent phases of the infections caused by these viruses, are likely to play relevant roles in key aspects of the life cycle of these viruses by modulating and regulating their replication, establishment of latency and evasion of the host antiviral response. Here, we review and discuss current evidence regarding epitranscriptomic and epigenetic modifications of herpesviruses and how these can influence their life cycles. While epitranscriptomic modifications such as m⁶A are the most studied to date and relate to positive effects over the replication of herpesviruses, epigenetic modifications of the viral genome are generally associated with defense mechanisms of the host cells to suppress viral gene transcription. However, herpesviruses can modulate these modifications to their own benefit to persist in the host, undergo latency and sporadically reactivate.

Pseudouridine and N1-Methylpseudouridine Display pH-Independent Reaction Rates with Bisulfite Yielding Ribose Adducts

In RNA, pseudouridine (Ψ) and 5-methylcytidine (m5C) are located by their differential reactions with NaHSO3 at pH 5. The pyrimidines were allowed to react with NaHSO3, NaN3, NaCN, or NaSCN at pH 5 to find that NaHSO3 was unique in achieving quantitative yields. Pseudouridine reaction selectivity with NaHSO3 was found at pH 7 supported by the reaction rate constants. The Ψ derivative N1-methylpseudouridine found in mRNA vaccines reacts similarly with bisulfite to yield ribose adducts.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

The plant epitranscriptome: revisiting pseudouridine and 2'-O-methyl RNA modifications

There is growing evidence that post-transcriptional RNA modifications are highly dynamic and can be used to improve crop production. Although more than 172 unique types of RNA modifications have been identified throughout the kingdom of life, we are yet to leverage upon the understanding to optimize RNA modifications in crops to improve productivity. The contributions of internal mRNA modifications such as N6-methyladenosine (m⁶ A) and 5-methylcytosine (m⁵ C) methylations to embryonic development, root development, leaf morphogenesis, flowering, fruit ripening and stress response are sufficiently known, but the roles of the two most abundant RNA modifications, pseudouridine (Ψ) and 2'-O-methylation (Nm), in the cell remain unclear due to insufficient advances in high-throughput technologies in plant development. Therefore, in this review, we discuss the latest methods and insights gained in mapping internal Ψ and Nm and their unique properties in plants and other organisms. In addition, we discuss the limitations that remain in high-throughput technologies for qualitative and quantitative mapping of these RNA modifications and highlight future challenges in regulating the plant epitranscriptome.

Nanopore-Based Detection of Viral RNA Modifications

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

Chemistry solutions to facilitate nanopore detection and analysis

Characteristic ionic current modulations will be produced in a single molecule manner during the communication of individual molecules with a nanopore. Hence, the information regarding the length, composition, and structure of a molecule can be extracted from deciphering the electrical message. However, until now, achieving a satisfactory resolution for observation and quantification of a target analyte in a complex system remains a nontrivial task. In this review, we summarize the progress and especially the recent advance in utilizing chemistry solutions to facilitate nanopore detection and analysis. The discussed chemistry solutions are classified into several major categories, including covalent/non-covalent chemistry, redox chemistry, displacement chemistry, back titration chemistry, chelation chemistry, hydrolysis-chemistry, and click chemistry. Considering the significant success of using chemical reaction-assisted nanopore sensing strategies to improve sensor sensitivity & selectivity and to study various topics, other non-chemistry based methodologies can undoubtedly be employed by nanopore sensors to explore new applications in the interdisciplinary area of chemistry, biology, materials, and nanotechnology.

Epitranscriptomics of SARS-CoV-2 Infection

Recent studies on the epitranscriptomic code of SARS-CoV-2 infection have discovered various RNA modifications, such as N6-methyladenosine (m6A), pseudouridine (Ψ), and 2′-O-methylation (Nm). The effects of RNA methylation on SARS-CoV-2 replication and the enzymes involved in this mechanism are emerging. In this review, we summarize the advances in this emerging field and discuss the role of various players such as readers, writers, and erasers in m6A RNA methylation, the role of pseudouridine synthase one and seven in epitranscriptomic modification Ψ, an isomer of uridine, and role of nsp16/nsp10 heterodimer in 2′-O-methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We also discuss RNA expression levels of various enzymes involved in RNA modifications in blood cells of SARS-CoV-2 infected individuals and their impact on host mRNA modification. In conclusion, these observations will facilitate the development of novel strategies and therapeutics for targeting RNA modification of SARS-CoV-2 RNA to control SARS-CoV-2 infection.

Techniques to detect epitranscriptomic marks

Similar to epigenetic DNA modification, RNA can be methylated and altered for stability and processing. RNA modifications, i.e., epitranscriptomes involve three functions, that is, writing, erasing, and reading of marks. Methods for measurement and position detection are useful for the assessment of cellular function and human disease biomarkers. Since the first detection of pyrimidine 5-methylcytosine hundred years ago, numerous techniques have been developed to study the modifications of nucleotides, including RNAs. Recent studies focused on high throughput and direct measurements to investigate the precise function of epitranscriptomes, including the characterization of SARS-CoV-2. The current work presents an overview of the development of detection techniques for epitranscriptomic marks and updates recent progress on the related field.

Promising Assays for Examining a Putative Role of Ribosomal Heterogeneity in COVID-19 Susceptibility and Severity

The heterogeneity of ribosomes, characterized by structural variations, arises from differences in types, numbers, and/or post-translational modifications of participating ribosomal proteins (RPs), ribosomal RNAs (rRNAs) sequence variants plus post-transcriptional modifications, and additional molecules essential for forming a translational machinery. The ribosomal heterogeneity within an individual organism or a single cell leads to preferential translations of selected messenger RNA (mRNA) transcripts over others, especially in response to environmental cues. The role of ribosomal heterogeneity in SARS-CoV-2 coronavirus infection, propagation, related symptoms, or vaccine responses is not known, and a technique to examine these has not yet been developed. Tools to detect ribosomal heterogeneity or to profile translating mRNAs independently cannot identify unique or specialized ribosome(s) along with corresponding mRNA substrate(s). Concurrent characterizations of RPs and/or rRNAs with mRNA substrate from a single ribosome would be critical to decipher the putative role of ribosomal heterogeneity in the COVID-19 disease, caused by the SARS-CoV-2, which hijacks the host ribosome to preferentially translate its RNA genome. Such a protocol should be able to provide a high-throughput screening of clinical samples in a large population that would reach a statistical power for determining the impact of a specialized ribosome to specific characteristics of the disease. These characteristics may include host susceptibility, viral infectivity and transmissibility, severity of symptoms, antiviral treatment responses, and vaccine immunogenicity including its side effect and efficacy. In this study, several state-of-the-art techniques, in particular, chemical probing of ribosomal components or rRNA structures, proximity ligation to generate rRNA-mRNA chimeras for sequencing, nanopore gating of individual ribosomes, nanopore RNA sequencing and/or structural analyses, single-ribosome mass spectrometry, and microfluidic droplets for separating ribosomes or indexing rRNAs/mRNAs, are discussed. The key elements for further improvement and proper integration of the above techniques to potentially arrive at a high-throughput protocol for examining individual ribosomes and their mRNA substrates in a clinical setting are also presented.

Guide RNA acrobatics: the one-for-two shuffle

RNA modifications are crucial for the proper function of the RNAs. The sites of pseudouridines are often specified by dual hairpin guide RNAs, with one or both hairpins identifying a target uridine. In this issue of Genes & Development, Jády and colleagues (pp. 70-83) identify a novel mechanism by which a single guide RNA hairpin can specify two uridines adjacent to each other or separated by 1 nt; i.e., one for two or guide RNA acrobatics.

Interferon inducible pseudouridine modification in human mRNA by quantitative nanopore profiling

Pseudouridine (Ψ) is an abundant mRNA modification in mammalian transcriptome, but its functions have remained elusive due to the difficulty of transcriptome-wide mapping. We develop a nanopore native RNA sequencing method for quantitative Ψ prediction (NanoPsu) that utilizes native content training, machine learning modeling, and single-read linkage analysis. Biologically, we find interferon inducible Ψ modifications in interferon-stimulated gene transcripts which are consistent with a role of Ψ in enabling efficacy of mRNA vaccines.