RNA cytosine methylation analysis by bisulfite sequencing

Chemical Modifications to RNA: A New Layer of Gene Expression Regulation

The first chemical modification to RNA was discovered nearly 60 years ago; to date, more than 100 chemically distinct modifications have been identified in cellular RNA. With the recent development of novel chemical and/or biochemical methods, dynamic modifications to RNA have been identified in the transcriptome, including N⁶-methyladenosine (m⁶A), inosine (I), 5-methylcytosine (m⁵C), pseudouridine (Ψ), 5-hydroxymethylcytosine (hm⁵C), and N¹-methyladenosine (m¹A). Collectively, the multitude of RNA modifications are termed epitranscriptome, leading to the emerging field of epitranscriptomics. In this review, we primarily focus on recently reported chemical modifications to mRNA; we discuss their chemical properties, biological functions, and mechanisms with an emphasis on their high-throughput detection methods. We also envision that future tools, particularly novel chemical biology methods, could further facilitate and enable studies in the field of epitranscriptomics.

Distinct 5-methylcytosine profiles in poly(A) RNA from mouse embryonic stem cells and brain

BACKGROUND

Recent work has identified and mapped a range of posttranscriptional modifications in mRNA, including methylation of the N6 and N1 positions in adenine, pseudouridylation, and methylation of carbon 5 in cytosine (m5C). However, knowledge about the prevalence and transcriptome-wide distribution of m5C is still extremely limited; thus, studies in different cell types, tissues, and organisms are needed to gain insight into possible functions of this modification and implications for other regulatory processes.

RESULTS

We have carried out an unbiased global analysis of m5C in total and nuclear poly(A) RNA of mouse embryonic stem cells and murine brain. We show that there are intriguing differences in these samples and cell compartments with respect to the degree of methylation, functional classification of methylated transcripts, and position bias within the transcript. Specifically, we observe a pronounced accumulation of m5C sites in the vicinity of the translational start codon, depletion in coding sequences, and mixed patterns of enrichment in the 3' UTR. Degree and pattern of methylation distinguish transcripts modified in both embryonic stem cells and brain from those methylated in either one of the samples. We also analyze potential correlations between m5C and micro RNA target sites, binding sites of RNA binding proteins, and N6-methyladenosine.

CONCLUSION

Our study presents the first comprehensive picture of cytosine methylation in the epitranscriptome of pluripotent and differentiated stages in the mouse. These data provide an invaluable resource for future studies of function and biological significance of m5C in mRNA in mammals.

Post-transcriptional gene regulation by mRNA modifications

The recent discovery of reversible mRNA methylation has opened a new realm of post-transcriptional gene regulation in eukaryotes. The identification and functional characterization of proteins that specifically recognize RNA N6-methyladenosine (m6A) unveiled it as a modification that cells utilize to accelerate mRNA metabolism and translation. N6-adenosine methylation directs mRNAs to distinct fates by grouping them for differential processing, translation and decay in processes such as cell differentiation, embryonic development and stress responses. Other mRNA modifications, including N1-methyladenosine (m1A), 5-methylcytosine (m5C) and pseudouridine, together with m6A form the epitranscriptome and collectively code a new layer of information that controls protein synthesis.

Transcriptome-Wide Detection of 5-Methylcytosine by Bisulfite Sequencing

While low-throughput RNA bisulfite sequencing is the method of choice to assess the methylation status of specific cytosines in candidate RNAs, the combination of bisulfite treatment of RNA with today's high-throughput sequencing techniques opens the door to methylation studies at nucleotide resolution on a transcriptome-wide scale. Below we describe a protocol for the transcriptome-wide analysis of total or fractionated poly(A)RNA in cells and tissues. Although the nature of the bisulfite sequencing protocol makes it comparably easy to translate from a low to a high-throughput approach, several critical points require attention before starting such a project. We describe a step-by-step protocol for planning and performing the experiment and analyzing the data.

Analysis of High-Throughput RNA Bisulfite Sequencing Data

Methylation of the 5-cytosine (m⁵C) is a common but not well-understood RNA modification, which can be detected by sequencing of bisulfite-treated transcripts (RNA-BSseq). In this Chapter, we discuss computational RNA-BSseq data analysis methods for transcriptome-wide identification and quantification of m⁵C.

Detection of 5-Methylcytosine in Specific Poly(A) RNAs by Bisulfite Sequencing

RNA bisulfite sequencing (RNA-BS-seq) represents a method for the detection of methylated cytosines in RNA. Developed originally for the analysis of DNA methylation, a modified version of this method can be used for the analysis of methylated cytosine in RNA. Treatment of nucleic acids with HSO₃-ions under acidic conditions results in deamination of cytosine (C) to uracil, while 5-methylcytosine (m5C) or 5-hydroxymethylcytosine (hm5C) exhibit low reactivity in this reaction and remain unchanged. Subsequent PCR amplification and sequencing of specific targets allows for the assessment of the methylation status of single Cs in their native sequence context at nucleotide resolution. Here, we describe the application of this method for the analysis of cytosine methylation in low abundance poly(A)RNA using a combination of commercially available kits and standard lab methods to ensure reproducible results. Furthermore, useful information on optimizing the method, suitable controls for almost all steps, and general troubleshooting guides are provided.

Capture and sequencing of NAD-capped RNA sequences with NAD captureSeq

Here we describe a protocol for NAD captureSeq that allows for the identification of nicotinamide-adenine dinucleotide (NAD)-capped RNA sequences in total RNA samples from different organisms. NAD-capped RNA is first chemo-enzymatically biotinylated with high efficiency, permitting selective capture on streptavidin beads. Then, a highly efficient library preparation protocol tailored to immobilized, 5'-modified RNA is applied, with adaptor ligation to the RNA's 3' terminus and reverse transcription (RT) performed on-bead. Then, cDNA is released into solution, tailed, ligated to a second adaptor and PCR-amplified. After next-generation sequencing (NGS) of the DNA library, enriched sequences are identified by comparison with a control sample in which the first step of chemo-enzymatic biotinylation is omitted. Because the downstream protocol does not necessarily rely on NAD-modified but on 'clickable' or biotin-modified RNA, it can be applied to other RNA modifications or RNA-biomolecule interactions. The central part of this protocol can be completed in ∼7 d, excluding preparatory steps, sequencing and bioinformatic analysis.

The emerging biology of RNA post-transcriptional modifications

RNA modifications have long been known to be central in the proper function of tRNA and rRNA. While chemical modifications in mRNA were discovered decades ago, their function has remained largely mysterious until recently. Using enrichment strategies coupled to next generation sequencing, multiple modifications have now been mapped on a transcriptome-wide scale in a variety of contexts. We now know that RNA modifications influence cell biology by many different mechanisms - by influencing RNA structure, by tuning interactions within the ribosome, and by recruiting specific binding proteins that intersect with other signaling pathways. They are also dynamic, changing in distribution or level in response to stresses such as heat shock and nutrient deprivation. Here, we provide an overview of recent themes that have emerged from the substantial progress that has been made in our understanding of chemical modifications across many major RNA classes in eukaryotes.

Detection of nucleic acid modifications by chemical reagents

Nucleic acids, especially RNA, naturally contain a diversity of chemically modified nucleosides. To understand the biological role of these modified nucleosides, nucleic acid scientists need tools to specifically label, detect and enrich modified nucleic acids. These tools comprise a diverse set of chemical reagents which have been established in the early years of nucleic acid research. Recent developments in high-throughput sequencing and mass spectrometry utilize these chemical labeling strategies to efficiently detect and localize modifications in nucleic acids. As a consequence the transcriptome-wide distribution of modified nucleosides, especially 5-methylcytosine and pseudouridine, in all domains of life could be analyzed. With the help of these techniques and the gained knowledge, it becomes possible to understand the functions of modifications and even study their connections to human health and disease. Here, the differential chemical reactivity of modified nucleosides and their canonical counterpart is reviewed and discussed.

Small RNA Modifications: Integral to Function and Disease

Small RNAs have the potential to store a secondary layer of labile biological information in the form of modified nucleotides. Emerging evidence has shown that small RNAs including microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs) and tRNA-derived small RNAs (tsRNAs) harbor a diversity of RNA modifications. These findings highlight the importance of RNA modifications in the modulation of basic properties such as RNA stability and other complex physiological processes involved in stress responses, metabolism, immunity, and epigenetic inheritance of environmentally acquired traits, among others. High-resolution, high-throughput methods for detecting, mapping and screening these small RNA modifications now provide opportunities to uncover their diagnostic potential as sensitive disease markers.

CoverageAnalyzer (CAn): A Tool for Inspection of Modification Signatures in RNA Sequencing Profiles

Combination of reverse transcription (RT) and deep sequencing has emerged as a powerful instrument for the detection of RNA modifications, a field that has seen a recent surge in activity because of its importance in gene regulation. Recent studies yielded high-resolution RT signatures of modified ribonucleotides relying on both sequence-dependent mismatch patterns and reverse transcription arrests. Common alignment viewers lack specialized functionality, such as filtering, tailored visualization, image export and differential analysis. Consequently, the community will profit from a platform seamlessly connecting detailed visual inspection of RT signatures and automated screening for modification candidates. CoverageAnalyzer (CAn) was developed in response to the demand for a powerful inspection tool. It is freely available for all three main operating systems. With SAM file format as standard input, CAn is an intuitive and user-friendly tool that is generally applicable to the large community of biomedical users, starting from simple visualization of RNA sequencing (RNA-Seq) data, up to sophisticated modification analysis with significance-based modification candidate calling.

The RNA methyltransferase Dnmt2 methylates DNA in the structural context of a tRNA

The amino acid sequence of Dnmt2 is very similar to the catalytic domains of bacterial and eukaryotic DNA-(cytosine 5)-methyltransferases, but it efficiently catalyzes tRNA methylation, while its DNA methyltransferase activity is the subject of controversial reports with rates varying between zero and very weak. By using composite nucleic acid molecules as substrates, we surprisingly found that DNA fragments, when presented as covalent DNA-RNA hybrids in the structural context of a tRNA, can be more efficiently methylated than the corresponding natural tRNA substrate. Furthermore, by stepwise development of tRNA^Asp, we showed that this natural Dnmt2 substrate could be engineered to employ RNAs that act like guide RNAs in vitro. The 5’-half of tRNA^Asp was able to efficiently guide methylation toward a single stranded tRNA fragment as would result from tRNA cleavage by tRNA specific nucleases. In a more artificial setting, a composite system of guide RNAs could ultimately be engineered to enable the enzyme to perform cytidine methylation on single stranded DNA in vitro.

Next-generation sequencing technologies for detection of modified nucleotides in RNAs

Our ability to map and quantify RNA modifications at a genome-wide scale have revolutionized our understanding of the pervasiveness and dynamic regulation of diverse RNA modifications. Recent efforts in the field have demonstrated the presence of modified residues in almost any type of cellular RNA. Next-generation sequencing (NGS) technologies are the primary choice for transcriptome-wide RNA modification mapping. Here we provide an overview of approaches for RNA modification detection based on their RT-signature, specific chemicals, antibody-dependent (Ab) enrichment, or combinations thereof. We further discuss sources of artifacts in genome-wide modification maps, and experimental and computational considerations to overcome them. The future in this field is tightly linked to the development of new specific chemical reagents, highly specific Ab against RNA modifications and use of single-molecule RNA sequencing techniques.

Sequence-specific and Shape-selective RNA Recognition by the Human RNA 5-Methylcytosine Methyltransferase NSun6

Human NSun6 is an RNA methyltransferase that catalyzes the transfer of the methyl group from S-adenosyl-l-methionine (SAM) to C72 of tRNA^Thr and tRNA^Cys In the current study, we used mass spectrometry to demonstrate that human NSun6 indeed introduces 5-methylcytosine (m⁵C) into tRNA, as expected. To further reveal the tRNA recognition mechanism of human NSun6, we measured the methylation activity of human NSun6 and its kinetic parameters for different tRNA substrates and their mutants. We showed that human NSun6 requires a well folded, full-length tRNA as its substrate. In the acceptor region, the CCA terminus, the target site C72, the discriminator base U73, and the second and third base pairs (2:71 and 3:70) of the acceptor stem are all important RNA recognition elements for human NSun6. In addition, two specific base pairs (11:24 and 12:23) in the D-stem of the tRNA substrate are involved in interacting with human NSun6. Together, our findings suggest that human NSun6 relies on a delicate network for RNA recognition, which involves both the primary sequence and tertiary structure of tRNA substrates.

Diastereoselectivity of 5-Methyluridine Osmylation Is Inverted inside an RNA Chain

In this study, we investigated the reaction of the osmium tetroxide-bipyridine complex with pyrimidines in RNA. This reagent, which reacts with the diastereotopic 5-6 double bond, thus leading to the formation of two diastereomers, was used in the past to label thymidine and 5-methylcytosine in DNA. In light of the growing interest in post-transcriptional RNA modifications, we addressed the question of whether this reagent could be used for labeling of the naturally occurring RNA modifications 5-methylcytosine and 5-methyluridine. On nucleoside level, 5-methylcytosine and 5-methyluridine revealed a 5- and 12-fold preference, respectively, over their nonmethylated equivalents. Performing the reaction on an RNA level, we could show that the steric environment of a pentanucleotide has a major detrimental impact on the reaction rate of osmylation. Interestingly, this drop in reactivity was due to a dramatic change in diastereoselectivity, which in turn resulted from impediment of the preferred attack via the si side. Thus, while on the nucleoside level, the absolute configuration of the major product of osmylation of 5-methyluridine was (5R,6S)-5-methyluridine glycol-dioxoosmium-bipyridine, reaction with an RNA pentanucleotide afforded the corresponding (5S,6R)-diastereomer as the major product. The change in diastereoselectivity lead to an almost complete loss of selectivity toward 5-methylcytosine in a pentanucleotide context, while 5-methyluridine remained about 8 times more reactive than the canonical pyrimidines. On the basis of these findings, we evaluate the usefulness of osmium tetroxide-bipyridine as a potential label for the 5-methyluridine modification in transcriptome-wide studies.

Systematic characterization of novel lncRNAs responding to phosphate starvation in Arabidopsis thaliana

BACKGROUND

Previously, several long non-coding RNAs (lncRNAs) were characterized as regulators in phosphate (Pi) starvation responses. However, systematic studies of novel lncRNAs involved in the Pi starvation signaling pathways have not been reported.

RESULTS

Here, we used a genome-wide sequencing and bioinformatics approach to identify both poly(A) + and poly(A)- lncRNAs that responded to Pi starvation in Arabidopsis thaliana. We sequenced shoot and root transcriptomes of the Arabidopsis seedlings grown under Pi-sufficient and Pi-deficient conditions, and predicted 1212 novel lncRNAs, of which 78 were poly(A)- lncRNAs. By employing strand-specific RNA libraries, we discovered many novel antisense lncRNAs for the first time. We further defined 309 lncRNAs that were differentially expressed between P+ and P- conditions in either shoots or roots. Through Gene Ontology enrichment of the associated protein-coding genes (co-expressed or close on the genome), we found that many lncRNAs were adjacent or co-expressed with the genes involved in several Pi starvation related processes, including cell wall organization and photosynthesis. In total, we identified 104 potential lncRNA targets of PHR1, a key regulator for transcriptional response to Pi starvation. Moreover, we identified 16 candidate lncRNAs as potential targets of miR399, another key regulator of plant Pi homeostasis.

CONCLUSIONS

Altogether, our data provide a rich resource of candidate lncRNAs involved in the Pi starvation regulatory network.

NSUN3 and ABH1 modify the wobble position of mt-tRNAMet to expand codon recognition in mitochondrial translation

Mitochondrial gene expression uses a non‐universal genetic code in mammals. Besides reading the conventional AUG codon, mitochondrial (mt‐)tRNA^Met mediates incorporation of methionine on AUA and AUU codons during translation initiation and on AUA codons during elongation. We show that the RNA methyltransferase NSUN3 localises to mitochondria and interacts with mt‐tRNA^Met to methylate cytosine 34 (C34) at the wobble position. NSUN3 specifically recognises the anticodon stem loop (ASL) of the tRNA, explaining why a mutation that compromises ASL basepairing leads to disease. We further identify ALKBH1/ABH1 as the dioxygenase responsible for oxidising m⁵C34 of mt‐tRNA^Met to generate an f⁵C34 modification. In vitro codon recognition studies with mitochondrial translation factors reveal preferential utilisation of m⁵C34 mt‐tRNA^Met in initiation. Depletion of either NSUN3 or ABH1 strongly affects mitochondrial translation in human cells, implying that modifications generated by both enzymes are necessary for mt‐tRNA^Met function. Together, our data reveal how modifications in mt‐tRNA^Met are generated by the sequential action of NSUN3 and ABH1, allowing the single mitochondrial tRNA^Met to recognise the different codons encoding methionine.

Nucleotide-Level Profiling of m⁵C RNA Methylation

Mapping the position and quantifying the level of 5-methylcytosine (m(5)C) as a modification in different types of cellular RNA is an important objective in the emerging field of epitranscriptomics. Bisulfite conversion has long been the gold standard for detection of m(5)C in DNA but it can also be applied to RNA. Here, we detail methods for bisulfite treatment of RNA, locus-specific PCR amplification and detection of candidate sites by sequencing on the Illumina MiSeq platform.

Translating the epitranscriptome

RNA modifications are indispensable for the translation machinery to provide accurate and efficient protein synthesis. Whereas the importance of transfer RNA (tRNA) and ribosomal RNA (rRNA) modifications has been well described and is unquestioned for decades, the significance of internal messenger RNA (mRNA) modifications has only recently been revealed. Novel experimental methods have enabled the identification of thousands of modified sites within the untranslated and translated regions of mRNAs. Thus far, N⁶‐methyladenosine (m⁶A), pseudouridine (Ψ), 5‐methylcytosine (m⁵C) and N¹‐methyladenosine (m¹A) were identified in eukaryal, and to some extent in prokaryal mRNAs. Several of the functions of these mRNA modifications have previously been reported, but many aspects remain elusive. Modifications can be important factors for the direct regulation of protein synthesis. The potential diversification of genomic information and regulation of RNA expression through editing and modifying mRNAs is versatile and many questions need to be addressed to completely elucidate the role of mRNA modifications. Herein, we summarize and highlight some recent findings on various co‐ and post‐transcriptional modifications, describing the impact of these processes on gene expression, with emphasis on protein synthesis. WIREs RNA 2017, 8:e1375. doi: 10.1002/wrna.1375

For further resources related to this article, please visit the WIREs website.

Messenger RNA modifications: Form, distribution, and function

RNA contains more than 100 distinct modifications that promote the functions of stable noncoding RNAs in translation and splicing. Recent technical advances have revealed widespread and sparse modification of messenger RNAs with N(6)-methyladenosine (m(6)A), 5-methylcytosine (m(5)C), and pseudouridine (Ψ). Here we discuss the rapidly evolving understanding of the location, regulation, and function of these dynamic mRNA marks, collectively termed the epitranscriptome. We highlight differences among modifications and between species that could instruct ongoing efforts to understand how specific mRNA target sites are selected and how their modification is regulated. Diverse molecular consequences of individual m(6)A modifications are beginning to be revealed, but the effects of m(5)C and Ψ remain largely unknown. Future work linking molecular effects to organismal phenotypes will broaden our understanding of mRNA modifications as cell and developmental regulators.

Stem cell function and stress response are controlled by protein synthesis

Whether protein synthesis and cellular stress response pathways interact to control stem cell function is currently unknown. Here we show that mouse skin stem cells synthesize less protein than their immediate progenitors in vivo, even when forced to proliferate. Our analyses reveal that activation of stress response pathways drives both a global reduction of protein synthesis and altered translational programmes that together promote stem cell functions and tumorigenesis. Mechanistically, we show that inhibition of post-transcriptional cytosine-5 methylation locks tumour-initiating cells in this distinct translational inhibition programme. Paradoxically, this inhibition renders stem cells hypersensitive to cytotoxic stress, as tumour regeneration after treatment with 5-fluorouracil is blocked. Thus, stem cells must revoke translation inhibition pathways to regenerate a tissue or tumour.

Mechanism and biological role of Dnmt2 in Nucleic Acid Methylation

A group of homologous nucleic acid modification enzymes called Dnmt2, Trdmt1, Pmt1, DnmA, and Ehmet in different model organisms catalyze the transfer of a methyl group from the cofactor S-adenosyl-methionine (SAM) to the carbon-5 of cytosine residues. Originally considered as DNA MTases, these enzymes were shown to be tRNA methyltransferases about a decade ago. Between the presumed involvement in DNA modification-related epigenetics, and the recent foray into the RNA modification field, significant progress has characterized Dnmt2-related research. Here, we review this progress in its diverse facets including molecular evolution, structural biology, biochemistry, chemical biology, cell biology and epigenetics.

Guitar: An R/Bioconductor Package for Gene Annotation Guided Transcriptomic Analysis of RNA-Related Genomic Features

Biological features, such as genes and transcription factor binding sites, are often denoted with genome-based coordinates as the genomic features. While genome-based representation is usually very effective in correlating various biological features, it can be tedious to examine the relationship between RNA-related genomic features and the landmarks of RNA transcripts with existing tools due to the difficulty in the conversion between genome-based coordinates and RNA-based coordinates. We developed here an open source Guitar R/Bioconductor package for sketching the transcriptomic view of RNA-related biological features represented by genome based coordinates. Internally, Guitar package extracts the standardized RNA coordinates with respect to the landmarks of RNA transcripts, with which hundreds of millions of RNA-related genomic features can then be efficiently analyzed within minutes. We demonstrated the usage of Guitar package in analyzing posttranscriptional RNA modifications (5-methylcytosine and N6-methyladenosine) derived from high-throughput sequencing approaches (MeRIP-Seq and RNA BS-Seq) and show that RNA 5-methylcytosine (m⁵C) is enriched in 5′UTR. The newly developed Guitar R/Bioconductor package achieves stable performance on the data tested and revealed novel biological insights. It will effectively facilitate the analysis of RNA methylation data and other RNA-related biological features in the future.

DRME: Count-based differential RNA methylation analysis at small sample size scenario

Differential methylation, which concerns difference in the degree of epigenetic regulation via methylation between two conditions, has been formulated as a beta or beta-binomial distribution to address the within-group biological variability in sequencing data. However, a beta or beta-binomial model is usually difficult to infer at small sample size scenario with discrete reads count in sequencing data. On the other hand, as an emerging research field, RNA methylation has drawn more and more attention recently, and the differential analysis of RNA methylation is significantly different from that of DNA methylation due to the impact of transcriptional regulation. We developed DRME to better address the differential RNA methylation problem. The proposed model can effectively describe within-group biological variability at small sample size scenario and handles the impact of transcriptional regulation on RNA methylation. We tested the newly developed DRME algorithm on simulated and 4 MeRIP-Seq case-control studies and compared it with Fisher's exact test. It is in principle widely applicable to several other RNA-related data types as well, including RNA Bisulfite sequencing and PAR-CLIP. The code together with an MeRIP-Seq dataset is available online (https://github.com/lzcyzm/DRME) for evaluation and reproduction of the figures shown in this article.

Emerging roles of RNA modifications in bacteria

RNA modifications are known to abound in stable tRNA and rRNA, where they cluster around functionally important regions. However, RNA-seq based techniques profiling entire transcriptomes are now uncovering an abundance of modified ribonucleotides in mRNAs and noncoding RNAs, too. While most of the recent progress in understanding the regulatory influence of these new RNA modifications stems from eukaryotes, there is growing evidence in bacteria for modified nucleotides beyond the stable RNA species, including modifications of small regulatory RNAs. Given their small genome size, good genetic tractability, and ample knowledge of modification enzymes, bacteria offer excellent model systems to decipher cellular functions of RNA modifications in many diverse physiological contexts. This review highlights how new global approaches combining classic analysis with new sequencing techniques may usher in an era of bacterial epitranscriptomics.

Mass spectrometry of modified RNAs: recent developments

A common feature of ribonucleic acids (RNAs) is that they can undergo a variety of chemical modifications. As nearly all of these chemical modifications result in an increase in the mass of the canonical nucleoside, mass spectrometry has long been a powerful approach for identifying and characterizing modified RNAs. Over the past several years, significant advances have been made in method development and software for interpreting tandem mass spectra resulting in approaches that can yield qualitative and quantitative information on RNA modifications, often at the level of sequence specificity. We discuss these advances along with instrumentation developments that have increased our ability to extract such information from relatively complex biological samples. With the increasing interest in how these modifications impact the epitranscriptome, mass spectrometry will continue to play an important role in bioanalytical investigations revolving around RNA.

The Antisense Transcriptome and the Human Brain

The transcriptome of a cell is made up of a varied array of RNA species, including protein-coding RNAs, long non-coding RNAs, short non-coding RNAs, and circular RNAs. The cellular transcriptome is dynamic and can change depending on environmental factors, disease state and cellular context. The human brain has perhaps the most diverse transcriptome profile that is enriched for many species of RNA, including antisense transcripts. Antisense transcripts are produced when both the plus and minus strand of the DNA helix are transcribed at a particular locus. This results in an RNA transcript that has a partial or complete overlap with an intronic or exonic region of the sense transcript. While antisense transcription is known to occur at some level in most organisms, this review focuses specifically on antisense transcription in the brain and how regulation of genes by antisense transcripts can contribute to functional aspects of the healthy and diseased brain. First, we discuss different techniques that can be used in the identification and quantification of antisense transcripts. This is followed by examples of antisense transcription and modes of regulatory function that have been identified in the brain.

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

Mass spectrometry (MS) is a powerful technique for characterizing noncanonical nucleobases and other chemical modifications in small RNAs, yielding rich chemical information that is complementary to high-throughput indirect sequencing. However, mass spectra are often prohibitively complex when fragment ions are analyzed following either solution phase hydrolysis or gas phase fragmentation. For all but the simplest cases, ions arising from multiple fragmentation events, alternative fragmentation pathways, and diverse salt adducts frequently obscure desired single-cut fragment ions. Here we show that it is possible to take advantage of predictable regularities in liquid chromatographic (LC) separation of optimized RNA digests to greatly simplify the interpretation of complex MS data. A two-dimensional analysis of extracted compound chromatograms permits straightforward and robust de novo sequencing, using a novel Monte Carlo algorithm that automatically generates bidirectional paired-end reads, pinpointing the position of modified nucleotides in a sequence. We demonstrate that these advances permit routine LC-MS sequencing of RNAs containing noncanonical nucleotides, and we furthermore examine the applicability of this approach to the study of oligonucleotides containing artificial modifications as well as those commonly observed in post-transcriptionally modified RNAs.

The identification and characterization of non-coding and coding RNAs and their modified nucleosides by mass spectrometry

The analysis of ribonucleic acids (RNA) by mass spectrometry has been a valuable analytical approach for more than 25 years. In fact, mass spectrometry has become a method of choice for the analysis of modified nucleosides from RNA isolated out of biological samples. This review summarizes recent progress that has been made in both nucleoside and oligonucleotide mass spectral analysis. Applications of mass spectrometry in the identification, characterization and quantification of modified nucleosides are discussed. At the oligonucleotide level, advances in modern mass spectrometry approaches combined with the standard RNA modification mapping protocol enable the characterization of RNAs of varying lengths ranging from low molecular weight short interfering RNAs (siRNAs) to the extremely large 23 S rRNAs. New variations and improvements to this protocol are reviewed, including top-down strategies, as these developments now enable qualitative and quantitative measurements of RNA modification patterns in a variety of biological systems.

Dynamic modulation of Dnmt2-dependent tRNA methylation by the micronutrient queuine

Dnmt2 enzymes are cytosine-5 methyltransferases that methylate C38 of several tRNAs. We report here that the activities of two Dnmt2 homologs, Pmt1 from Schizosaccharomyces pombe and DnmA from Dictyostelium discoideum, are strongly stimulated by prior queuosine (Q) modification of the substrate tRNA. In vivo tRNA methylation levels were stimulated by growth of cells in queuine-containing medium; in vitro Pmt1 activity was enhanced on Q-containing RNA; and queuine-stimulated in vivo methylation was abrogated by the absence of the enzyme that inserts queuine into tRNA, eukaryotic tRNA-guanine transglycosylase. Global analysis of tRNA methylation in S. pombe showed a striking selectivity of Pmt1 for tRNA(Asp) methylation, which distinguishes Pmt1 from other Dnmt2 homologs. The present analysis also revealed a novel Pmt1- and Q-independent tRNA methylation site in S. pombe, C34 of tRNA(Pro). Notably, queuine is a micronutrient that is scavenged by higher eukaryotes from the diet and gut microflora. This work therefore reveals an unanticipated route by which the environment can modulate tRNA modification in an organism.

The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent

The combination of Reverse Transcription (RT) and high-throughput sequencing has emerged as a powerful combination to detect modified nucleotides in RNA via analysis of either abortive RT-products or of the incorporation of mismatched dNTPs into cDNA. Here we simultaneously analyze both parameters in detail with respect to the occurrence of N-1-methyladenosine (m¹A) in the template RNA. This naturally occurring modification is associated with structural effects, but it is also known as a mediator of antibiotic resistance in ribosomal RNA. In structural probing experiments with dimethylsulfate, m¹A is routinely detected by RT-arrest. A specifically developed RNA-Seq protocol was tailored to the simultaneous analysis of RT-arrest and misincorporation patterns. By application to a variety of native and synthetic RNA preparations, we found a characteristic signature of m¹A, which, in addition to an arrest rate, features misincorporation as a significant component. Detailed analysis suggests that the signature depends on RNA structure and on the nature of the nucleotide 3′ of m¹A in the template RNA, meaning it is sequence dependent. The RT-signature of m¹A was used for inspection and confirmation of suspected modification sites and resulted in the identification of hitherto unknown m¹A residues in trypanosomal tRNA.

Conservation of tRNA and rRNA 5-methylcytosine in the kingdom Plantae

BACKGROUND

Post-transcriptional methylation of RNA cytosine residues to 5-methylcytosine (m(5)C) is an important modification that regulates RNA metabolism and occurs in both eukaryotes and prokaryotes. Yet, to date, no transcriptome-wide identification of m(5)C sites has been undertaken in plants. Plants provide a unique comparative system for investigating the origin and evolution of m(5)C as they contain three different genomes, the nucleus, mitochondria and chloroplast. Here we use bisulfite conversion of RNA combined with high-throughput IIlumina sequencing (RBS-seq) to identify single-nucleotide resolution of m(5)C sites in non-coding ribosomal RNAs and transfer RNAs of all three sub-cellular transcriptomes across six diverse species that included, the single-celled algae Nannochloropsis oculata, the macro algae Caulerpa taxifolia and multi-cellular higher plants Arabidopsis thaliana, Brassica rapa, Triticum durum and Ginkgo biloba.

RESULTS

Using the plant model Arabidopsis thaliana, we identified a total of 39 highly methylated m(5)C sites in predicted structural positions of nuclear tRNAs and 7 m(5)C sites in rRNAs from nuclear, chloroplast and mitochondrial transcriptomes. Both the nucleotide position and percent methylation of tRNAs and rRNAs m(5)C sites were conserved across all species analysed, from single celled algae N. oculata to multicellular plants. Interestingly the mitochondrial and chloroplast encoded tRNAs were devoid of m(5)C in A. thaliana and this is generally conserved across Plantae. This suggests independent evolution of organelle methylation in animals and plants, as animal mitochondrial tRNAs have m(5)C sites. Here we characterize 5 members of the RNA 5-methylcytosine family in Arabidopsis and extend the functional characterization of TRDMT1 and NOP2A/OLI2. We demonstrate that nuclear tRNA methylation requires two evolutionarily conserved methyltransferases, TRDMT1 and TRM4B. trdmt1 trm4b double mutants are hypersensitive to the antibiotic hygromycin B, demonstrating the function of tRNA methylation in regulating translation. Additionally we demonstrate that nuclear large subunit 25S rRNA methylation requires the conserved RNA methyltransferase NSUN5. Our results also suggest functional redundancy of at least two of the NOP2 paralogs in Arabidopsis.

CONCLUSIONS

Our data demonstrates widespread occurrence and conservation of non-coding RNA methylation in the kingdom Plantae, suggesting important and highly conserved roles of this post-transcriptional modification.

ARM-seq: AlkB-facilitated RNA methylation sequencing reveals a complex landscape of modified tRNA fragments

High throughput RNA sequencing has accelerated discovery of the complex regulatory roles of small RNAs, but RNAs containing modified nucleosides may escape detection when those modifications interfere with reverse transcription during RNA-seq library preparation. Here we describe AlkB-facilitated RNA Methylation sequencing (ARM-Seq) which uses pre-treatment with Escherichia coli AlkB to demethylate 1-methyladenosine, 3-methylcytidine, and 1-methylguanosine, all commonly found in transfer RNAs. Comparative methylation analysis using ARM-Seq provides the first detailed, transcriptome-scale map of these modifications, and reveals an abundance of previously undetected, methylated small RNAs derived from tRNAs. ARM-Seq demonstrates that tRNA-derived small RNAs accurately recapitulate the m¹A modification state for well-characterized yeast tRNAs, and generates new predictions for a large number of human tRNAs, including tRNA precursors and mitochondrial tRNAs. Thus, ARM-Seq provides broad utility for identifying previously overlooked methyl-modified RNAs, can efficiently monitor methylation state, and may reveal new roles for tRNA-derived RNAs as biomarkers or signaling molecules.

Enhanced detection of post-transcriptional modifications using a mass-exclusion list strategy for RNA modification mapping by LC-MS/MS

There has been a renewed appreciation for the dynamic nature of ribonucleic acid (RNA) modifications and for the impact of modified RNAs on organism health resulting in an increased emphasis on developing analytical methods capable of detecting modifications within specific RNA sequence contexts. Here we demonstrate that a DNA-based exclusion list enhances data dependent liquid chromatography tandem mass spectrometry (LC-MS/MS) detection of post-transcriptionally modified nucleosides within specific RNA sequences. This approach is possible because all post-transcriptional modifications of RNA, except pseudouridine, result in a mass increase in the canonical nucleoside undergoing chemical modification. Thus, DNA-based sequences reflect the state of the RNA prior to or in the absence of modification. The utility of this exclusion list strategy is demonstrated through the RNA modification mapping of total tRNAs from the bacteria Escherichia coli, Lactococcus lactis, and Streptomyces griseus. Creation of a DNA-based exclusion list is shown to consistently enhance the number of detected modified ribonuclease (RNase) digestion products by ∼20%. All modified RNase digestion products that were detected during standard data dependent acquisition (DDA) LC-MS/MS were also detected when the DNA-based exclusion list was used. Consequently, the increase in detected modified RNase digestion products is attributed to new experimental information only obtained when using the exclusion list. This exclusion list strategy should be broadly applicable to any class of RNA and improves the utility of mass spectrometry approaches for discovery-based analyses of RNA modifications, such as are required for studies of the epitranscriptome.

Eukaryotic rRNA Modification by Yeast 5-Methylcytosine-Methyltransferases and Human Proliferation-Associated Antigen p120

Modified nucleotide 5-methylcytosine (m5C) is frequently present in various eukaryotic RNAs, including tRNAs, rRNAs and in other non-coding RNAs, as well as in mRNAs. RNA:m5C-methyltranferases (MTases) Nop2 from S. cerevisiae and human proliferation-associated nucleolar antigen p120 are both members of a protein family called Nop2/NSUN/NOL1. Protein p120 is well-known as a tumor marker which is over-expressed in various cancer tissues. Using a combination of RNA bisulfite sequencing and HPLC-MS/MS analysis, we demonstrated here that p120 displays an RNA:m5C- MTase activity, which restores m5C formation at position 2870 in domain V of 25S rRNA in a nop2Δ yeast strain. We also confirm that yeast proteins Nop2p and Rcm1p catalyze the formation of m5C in domains V and IV, respectively. In addition, we do not find any evidence of m5C residues in yeast 18S rRNA. We also performed functional complementation of Nop2-deficient yeasts by human p120 and studied the importance of different sequence and structural domains of Nop2 and p120 for yeast growth and m5C-MTase activity. Chimeric protein formed by Nop2 and p120 fragments revealed the importance of Nop2 N-terminal domain for correct protein localization and its cellular function. We also validated that the presence of Nop2, rather than the m5C modification in rRNA itself, is required for pre-rRNA processing. Our results corroborate that Nop2 belongs to the large family of pre-ribosomal proteins and possesses two related functions in pre-rRNA processing: as an essential factor for cleavages and m5C:RNA:modification. These results support the notion of quality control during ribosome synthesis by such modification enzymes.

Capture, Unfolding, and Detection of Individual tRNA Molecules Using a Nanopore Device

Transfer RNAs (tRNA) are the most common RNA molecules in cells and have critical roles as both translators of the genetic code and regulators of protein synthesis. As such, numerous methods have focused on studying tRNA abundance and regulation, with the most widely used methods being RNA-seq and microarrays. Though revolutionary to transcriptomics, these assays are limited by an inability to encode tRNA modifications in the requisite cDNA. These modifications are abundant in tRNA and critical to their function. Here, we describe proof-of-concept experiments where individual tRNA molecules are examined as linear strands using a biological nanopore. This method utilizes an enzymatically ligated synthetic DNA adapter to concentrate tRNA at the lipid bilayer of the nanopore device and efficiently denature individual tRNA molecules, as they are pulled through the α-hemolysin (α-HL) nanopore. Additionally, the DNA adapter provides a loading site for ϕ29 DNA polymerase (ϕ29 DNAP), which acts as a brake on the translocating tRNA. This increases the dwell time of adapted tRNA in the nanopore, allowing us to identify the region of the nanopore signal that is produced by the translocating tRNA itself. Using adapter-modified Escherichia coli tRNA^fMet and tRNA^Lys, we show that the nanopore signal during controlled translocation is dependent on the identity of the tRNA. This confirms that adapter-modified tRNA can translocate end-to-end through nanopores and provide the foundation for future work in direct sequencing of individual transfer RNA with a nanopore-based device.

Experimental Approaches for Target Profiling of RNA Cytosine Methyltransferases

RNA cytosine methyltransferases (m(5)C-RMTs) constitute an important class of RNA-modifying enzymes, methylating specific cytosines within particular RNA targets in both coding and noncoding RNAs. Almost all organisms express at least one m(5)C-RMT, and vertebrates often express different types or variants of m(5)C-RMTs in different cell types. Deletion or mutation of particular m(5)C-RMTs is connected to severe pathological manifestations ranging from developmental defects to infertility and mental retardation. Some m(5)C-RMTs show spatiotemporal patterns of expression and activity requiring careful experimental design for their analysis in order to capture their context-dependent targets. An essential step for understanding the functions of both the enzymes and the modified cytosines is defining the one-to-one connection between particular m(5)C-RMTs and their target cytosines. Recent technological and methodological advances have provided researchers with new tools to comprehensively explore RNA cytosine methylation and methyltransferases. Here, we describe three complementary approaches applicable for both discovery and validation of candidate target sites of specific m(5)C-RMTs.

Whole transcriptome analysis with sequencing: methods, challenges and potential solutions

Whole transcriptome analysis plays an essential role in deciphering genome structure and function, identifying genetic networks underlying cellular, physiological, biochemical and biological systems and establishing molecular biomarkers that respond to diseases, pathogens and environmental challenges. Here, we review transcriptome analysis methods and technologies that have been used to conduct whole transcriptome shotgun sequencing or whole transcriptome tag/target sequencing analyses. We focus on how adaptors/linkers are added to both 5' and 3' ends of mRNA molecules for cloning or PCR amplification before sequencing. Challenges and potential solutions are also discussed. In brief, next generation sequencing platforms have accelerated releases of the large amounts of gene expression data. It is now time for the genome research community to assemble whole transcriptomes of all species and collect signature targets for each gene/transcript, and thus use known genes/transcripts to determine known transcriptomes directly in the near future.

On the mechanism of RNA phosphodiester backbone cleavage in the absence of solvent

Ribonucleic acid (RNA) modifications play an important role in the regulation of gene expression and the development of RNA-based therapeutics, but their identification, localization and relative quantitation by conventional biochemical methods can be quite challenging. As a promising alternative, mass spectrometry (MS) based approaches that involve RNA dissociation in 'top-down' strategies are currently being developed. For this purpose, it is essential to understand the dissociation mechanisms of unmodified and posttranscriptionally or synthetically modified RNA. Here, we have studied the effect of select nucleobase, ribose and backbone modifications on phosphodiester bond cleavage in collisionally activated dissociation (CAD) of positively and negatively charged RNA. We found that CAD of RNA is a stepwise reaction that is facilitated by, but does not require, the presence of positive charge. Preferred backbone cleavage next to adenosine and guanosine in CAD of (M+nH)(n+) and (M-nH)(n-) ions, respectively, is based on hydrogen bonding between nucleobase and phosphodiester moieties. Moreover, CAD of RNA involves an intermediate that is sufficiently stable to survive extension of the RNA structure and intramolecular proton redistribution according to simple Coulombic repulsion prior to backbone cleavage into C: and Y: ions from phosphodiester bond cleavage.

A mass spectrometry-based method for comprehensive quantitative determination of post-transcriptional RNA modifications: the complete chemical structure of Schizosaccharomyces pombe ribosomal RNAs

We present a liquid chromatography–mass spectrometry (LC-MS)-based method for comprehensive quantitative identification of post-transcriptional modifications (PTMs) of RNA. We incorporated an in vitro-transcribed, heavy isotope-labeled reference RNA into a sample RNA solution, digested the mixture with a number of RNases and detected the post-transcriptionally modified oligonucleotides quantitatively based on shifts in retention time and the MS signal in subsequent LC-MS. This allowed the determination and quantitation of all PTMs in Schizosaccharomyces pombe ribosomal (r)RNAs and generated the first complete PTM maps of eukaryotic rRNAs at single-nucleotide resolution. There were 122 modified sites, most of which appear to locate at the interface of ribosomal subunits where translation takes place. We also identified PTMs at specific locations in rRNAs that were altered in response to growth conditions of yeast cells, suggesting that the cells coordinately regulate the modification levels of RNA.

RNA 5-Methylcytosine Analysis by Bisulfite Sequencing

Cells have developed molecular machineries, which can chemically modify DNA and RNA nucleosides. One particular and chemically simple modification, (cytosine-5) methylation (m(5)C), has been detected both in RNA and DNA suggesting universal use of m(5)C for the function of these nucleotide polymers. m(5)C can be reproducibly mapped to abundant noncoding RNAs (transfer RNA, tRNA and ribosomal RNA, rRNA), and recently, also nonabundant RNAs (including mRNAs) have been reported to carry this modification. Quantification of m(5)C content in total RNA preparations indicates that a limited number of RNAs carry this modification and suggests specific functions for (cytosine-5) RNA methylation. What exactly is the biological function of m(5)C in RNA? Before attempting to address this question, m(5)C needs to be mapped specifically and reproducibly, preferably on a transcriptome-wide scale. To facilitate the detection of m(5)C in its sequence context, RNA bisulfite sequencing (RNA-BisSeq) has been developed. This method relies on the efficient chemical deamination of nonmethylated cytosine, which can be read out as single nucleotide polymorphism (nonmethylated cytosine as thymine vs. methylated cytosine as cytosine), when differentially comparing cDNA libraries to reference sequences after DNA sequencing. Here, the basic protocol of RNA-BisSeq, its current applications and limitations are described.

Site-specific modification of the 6-amino group of adenosine in RNA by an interstrand functionality-transfer reaction with an s-functionalized 4-thiothymidine

Non-natural RNA modifications have been widely used to study the function and structure of RNA. Expanding the study of RNA further requires versatile and efficient tools for site-specific RNA modification. We recently established a new strategy for the site-specific modification of RNA based on a functionality-transfer reaction between an oligodeoxynucleotide (ODN) probe and an RNA substrate. 2'-Deoxy-6-thioguanosine was used to anchor the transfer group, and the 4-amino group of cytosine or the 2-amino group of guanine was specifically modified. In this study, 2'-deoxy-4-thiothymidine was adopted as a new platform to target the 6-amino group of adenosine. The (E)-pyridinyl vinyl keto transfer group was attached to the 4-thioT in the ODN probe, and it was efficiently and specifically transferred to the 6-amino group of the opposing adenosine in RNA in the presence of CuCl2 . This method expands the available RNA target sites for specific modification.

Novel RNA modifications in the nervous system: form and function

Modified RNA molecules have recently been shown to regulate nervous system functions. This mini-review and associated mini-symposium provide an overview of the types and known functions of novel modified RNAs in the nervous system, including covalently modified RNAs, edited RNAs, and circular RNAs. We discuss basic molecular mechanisms involving RNA modifications as well as the impact of modified RNAs and their regulation on neuronal processes and disorders, including neural fate specification, intellectual disability, neurodegeneration, dopamine neuron function, and substance use disorders.

Methylation of ribosomal RNA by NSUN5 is a conserved mechanism modulating organismal lifespan

Several pathways modulating longevity and stress resistance converge on translation by targeting ribosomal proteins or initiation factors, but whether this involves modifications of ribosomal RNA is unclear. Here, we show that reduced levels of the conserved RNA methyltransferase NSUN5 increase the lifespan and stress resistance in yeast, worms and flies. Rcm1, the yeast homologue of NSUN5, methylates C2278 within a conserved region of 25S rRNA. Loss of Rcm1 alters the structural conformation of the ribosome in close proximity to C2278, as well as translational fidelity, and favours recruitment of a distinct subset of oxidative stress-responsive mRNAs into polysomes. Thus, rather than merely being a static molecular machine executing translation, the ribosome exhibits functional diversity by modification of just a single rRNA nucleotide, resulting in an alteration of organismal physiological behaviour, and linking rRNA-mediated translational regulation to modulation of lifespan, and differential stress response.

Sequencing of first-strand cDNA library reveals full-length transcriptomes

Massively parallel strand-specific sequencing of RNA (ssRNA-seq) has emerged as a powerful tool for profiling complex transcriptomes. However, many current methods for ssRNA-seq suffer from the underrepresentation of both the 5' and 3' ends of RNAs, which can be attributed to second-strand cDNA synthesis. The 5' and 3' ends of RNA harbour crucial information for gene regulation; namely, transcription start sites (TSSs) and polyadenylation sites. Here we report a novel ssRNA-seq method that does not involve second-strand cDNA synthesis, as we Directly Ligate sequencing Adaptors to the First-strand cDNA (DLAF). This novel method with fewer enzymatic reactions results in a higher quality of the libraries than the conventional method. Sequencing of DLAF libraries followed by a novel analysis pipeline enables the profiling of both 5' ends and polyadenylation sites at near-base resolution. Therefore, DLAF offers the first genomics tool to obtain the 'full-length' transcriptome with a single library.