General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

In-Gel Cyanoethylation for Pseudouridines Mass Spectrometry Detection of Bacterial Regulatory RNA

Regulatory RNAs, as well as many RNA families, contain chemically modified nucleotides, including pseudouridines (ψ). To map nucleotide modifications, approaches based on enzymatic digestion of RNA followed by nano liquid chromatography-tandem mass spectrometry (nanoLC-MS/MS) analysis were implemented several years ago. However, detection of ψ by mass spectrometry (MS) is challenging as ψ exhibits the same mass as uridine. Thus, a chemical labeling strategy using acrylonitrile was developed to detect this mass-silent modification. Acrylonitrile reacts specifically to ψ to form 1-cyanoethylpseudouridine (Ceψ), resulting in a mass shift of ψ detectable by MS. Here, a protocol detailing the steps from the purification of RNA by polyacrylamide gel electrophoresis, including in-gel labeling of ψ, to MS data interpretation to map ψ and other modifications is proposed. To demonstrate its efficiency, the protocol was applied to bacterial regulatory RNAs from E. coli: 6S RNA and transfer-messenger RNA (tmRNA, also known as 10Sa RNA). Moreover, ribonuclease P (RNase P) was also mapped using this approach. This method enabled the detection of several ψ at single nucleotide resolution.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

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.

Non-covalent dyes in microscale thermophoresis for studying RNA ligand interactions and modifications

Microscale Thermophoresis (MST) is a powerful biophysical technique that measures the mobility of biomolecules in response to a temperature gradient, making it useful for investigating the interactions between biological molecules. This study presents a novel methodology for studying RNA-containing samples using non-covalent nucleic acid-sensitive dyes in MST. This "mix-and-measure" protocol uses non-covalent dyes, such as those from the Syto or Sybr series, which lead to the statistical binding of one fluorophore per RNA oligo showing key advantages over traditional covalent labelling approaches. This new approach has been successfully used to study the binding of ligands to RNA molecules (e.g., SAM- and PreQ1 riboswitches) and the identification of modifications (e.g., m6A) in short RNA oligos which can be written by the RNA methyltransferase METTL3/14.

Data Analysis Pipeline for Detection and Quantification of Pseudouridine (ψ) in RNA by HydraPsiSeq

Pseudouridine, a modified RNA residue formed by the isomerization of its parental U nucleotide, is prevalent in a majority of cellular RNAs; its presence was reported in tRNA, rRNA, and sn/snoRNA as well as in mRNA/lncRNA. Multiple analytical deep sequencing-based approaches have been proposed for pseudouridine detection and quantification, among which the most popular relies on the use of soluble carbodiimide (termed CMCT). Recently, we developed an alternative protocol for pseudouridine mapping and quantification. The principle is based on protection of pseudouridine against random RNA cleavage by hydrazine/aniline treatment (HydraPsiSeq protocol). This "negative" detection mode requires higher sequencing depth and provides a precise quantification of the pseudouridine content. All "wet-lab" technical details of the HydraPsiSeq protocol have been described in recent publications. Here, we describe all bioinformatics analysis steps required for data processing from raw reads to the pseudouridylation profile of known or unknown RNA.

Oligonucleotides Featuring a Covalently Mercurated 6-Phenylcarbazole Residue as High-Affinity Hybridization Probes for Thiopyrimidine-Containing Sequences

Short oligonucleotides incorporating either 1-mercuri-6-phenylcarbazole, 8-mercuri-6-phenylcarbazole, or 1,8-dimercuri-6-phenylcarbazole C-nucleoside in the middle of the chain have been synthesized and studied for their potential as hybridization probes for sequences containing thiopyrimidine nucleobases. All of these oligonucleotides formed very stable duplexes with complementary sequences pairing the organometallic moiety with either 2- or 4-thiothymine. The isomeric monomercurated oligonucleotides were also able to discriminate between 2- and 4-thiothymine based on the different melting temperatures of the respective duplexes. DFT-optimized structures of the most stable mononuclear Hg^II -mediated base pairs featured a coordinated covalent bond between Hg^II and either S2 or S4 and a hydrogen bond between the carbazole nitrogen and N3. The dinuclear Hg^II -mediated base pairs, in turn, were geometrically very similar to the one previously reported to form between 1,8-dimercuri-6-phenylcarbazole and thymine and had one Hg^II ion coordinated to a thio and the other one to an oxo substituent.

Pseudouridylation of Epstein-Barr virus noncoding RNA EBER2 facilitates lytic replication

Epstein-Barr virus (EBV) expresses two highly abundant noncoding RNAs called EBV-encoded RNA 1 (EBER1) and EBER2, which are preserved in all clinical isolates of EBV, thus underscoring their essential function in the viral life cycle. Recent epitranscriptomics studies have uncovered a vast array of distinct RNA modifications within cellular as well as viral noncoding RNAs that are instrumental in executing their function. Here we show that EBER2 is marked by pseudouridylation, and by using HydraPsiSeq the modification site was mapped to a single nucleotide within the 3' region of EBER2. The writer enzyme was identified to be the snoRNA-dependent pseudouridine synthase Dyskerin, which is the catalytic subunit of H/ACA small nucleolar ribonucleoprotein complexes, and is guided to EBER2 by SNORA22. Similar to other noncoding RNAs for which pseudouridylation has a positive effect on RNA stability, loss of EBER2 pseudouridylation results in a decrease in RNA levels. Furthermore, pseudouridylation of EBER2 is required for the prolific accumulation of progeny viral genomes, suggesting that this single modification in EBER2 is important for efficient viral lytic replication. Taken together, our findings add to the list of RNA modifications that are essential for noncoding RNAs to implement their physiological roles.

Identified Gefitinib Metabolism-Related lncRNAs can be Applied to Predict Prognosis, Tumor Microenvironment, and Drug Sensitivity in Non-Small Cell Lung Cancer

Gefitinib has shown promising efficacy in the treatment of patients with locally advanced or metastatic EGFR-mutated non-small cell lung cancer (NSCLC). Molecular biomarkers for gefitinib metabolism-related lncRNAs have not yet been elucidated. Here, we downloaded relevant genes and matched them to relevant lncRNAs. We then used univariate, LASSO, and multivariate regression to screen for significant genes to construct prognostic models. We investigated TME and drug sensitivity by risk score data. All lncRNAs with differential expression were selected for GO/KEGG analysis. Imvigor210 cohort was used to validate the value of the prognostic model. Finally, we performed a stemness indices difference analysis. lncRNA-constructed prognostic models were significant in the high-risk and low-risk subgroups. Immune pathways were identified in both groups at low risk. The higher the risk score the greater the value of exclusion, MDSC, and CAF. PRRophetic algorithm screened a total of 58 compounds. In conclusion, the prognostic model we constructed can accurately predict OS in NSCLC patients. Two groups of low-risk immune pathways are beneficial to patients. Gefitinib metabolism was again validated to be related to cytochrome P450 and lipid metabolism. Finally, drugs that might be used to treat NSCLC patients were screened.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.