HAMR: high-throughput annotation of modified ribonucleotides

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.

Genome-Wide Approaches for RNA Structure Probing

RNA molecules of all types fold into complex secondary and tertiary structures that are important for their function and regulation. Structural and catalytic RNAs such as ribosomal RNA (rRNA) and transfer RNA (tRNA) are central players in protein synthesis, and only function through their proper folding into intricate three-dimensional structures. Studies of messenger RNA (mRNA) regulation have also revealed that structural elements embedded within these RNA species are important for the proper regulation of their total level in the transcriptome. More recently, the discovery of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) has shed light on the importance of RNA structure to genome, transcriptome, and proteome regulation. Due to the relatively small number, high conservation, and importance of structural and catalytic RNAs to all life, much early work in RNA structure analysis mapped out a detailed view of these molecules. Computational and physical methods were used in concert with enzymatic and chemical structure probing to create high-resolution models of these fundamental biological molecules. However, the recent expansion in our knowledge of the importance of RNA structure to coding and regulatory RNAs has left the field in need of faster and scalable methods for high-throughput structural analysis. To address this, nuclease and chemical RNA structure probing methodologies have been adapted for genome-wide analysis. These methods have been deployed to globally characterize thousands of RNA structures in a single experiment. Here, we review these experimental methodologies for high-throughput RNA structure determination and discuss the insights gained from each approach.

RNA Polymerase III Output Is Functionally Linked to tRNA Dimethyl-G26 Modification

Control of the differential abundance or activity of tRNAs can be important determinants of gene regulation. RNA polymerase (RNAP) III synthesizes all tRNAs in eukaryotes and it derepression is associated with cancer. Maf1 is a conserved general repressor of RNAP III under the control of the target of rapamycin (TOR) that acts to integrate transcriptional output and protein synthetic demand toward metabolic economy. Studies in budding yeast have indicated that the global tRNA gene activation that occurs with derepression of RNAP III via maf1-deletion is accompanied by a paradoxical loss of tRNA-mediated nonsense suppressor activity, manifested as an antisuppression phenotype, by an unknown mechanism. We show that maf1-antisuppression also occurs in the fission yeast S. pombe amidst general activation of RNAP III. We used tRNA-HydroSeq to document that little changes occurred in the relative levels of different tRNAs in maf1Δ cells. By contrast, the efficiency of N2,N2-dimethyl G26 (m(2)2G26) modification on certain tRNAs was decreased in response to maf1-deletion and associated with antisuppression, and was validated by other methods. Over-expression of Trm1, which produces m(2)2G26, reversed maf1-antisuppression. A model that emerges is that competition by increased tRNA levels in maf1Δ cells leads to m(2)2G26 hypomodification due to limiting Trm1, reducing the activity of suppressor-tRNASerUCA and accounting for antisuppression. Consistent with this, we show that RNAP III mutations associated with hypomyelinating leukodystrophy decrease tRNA transcription, increase m(2)2G26 efficiency and reverse antisuppression. Extending this more broadly, we show that a decrease in tRNA synthesis by treatment with rapamycin leads to increased m(2)2G26 modification and that this response is conserved among highly divergent yeasts and human cells.

The Epitranscriptome and Innate Immunity

Our knowledge of the variety and abundances of RNA base modifications is rapidly increasing. Modified bases have critical roles in tRNAs, rRNAs, translation, splicing, RNA interference, and other RNA processes, and are now increasingly detected in all types of transcripts. Can new biological principles associated with this diversity of RNA modifications, particularly in mRNAs and long non-coding RNAs, be identified? This review will explore this question by focusing primarily on adenosine to inosine (A-to-I) RNA editing by the adenine deaminase acting on RNA (ADAR) enzymes that have been intensively studied for the past 20 years and have a wide range of effects. Over 100 million adenosine to inosine editing sites have been identified in the human transcriptome, mostly in embedded Alu sequences that form potentially innate immune-stimulating dsRNA hairpins in transcripts. Recent research has demonstrated that inosine in the epitranscriptome and ADAR1 protein establish innate immune tolerance for host dsRNA formed by endogenous sequences. Innate immune sensors that detect viral nucleic acids are among the readers of epitranscriptome RNA modifications, though this does preclude a wide range of other modification effects.

DASHR: database of small human noncoding RNAs

Small non-coding RNAs (sncRNAs) are highly abundant RNAs, typically <100 nucleotides long, that act as key regulators of diverse cellular processes. Although thousands of sncRNA genes are known to exist in the human genome, no single database provides searchable, unified annotation, and expression information for full sncRNA transcripts and mature RNA products derived from these larger RNAs. Here, we present the Database of small human noncoding RNAs (DASHR). DASHR contains the most comprehensive information to date on human sncRNA genes and mature sncRNA products. DASHR provides a simple user interface for researchers to view sequence and secondary structure, compare expression levels, and evidence of specific processing across all sncRNA genes and mature sncRNA products in various human tissues. DASHR annotation and expression data covers all major classes of sncRNAs including microRNAs (miRNAs), Piwi-interacting (piRNAs), small nuclear, nucleolar, cytoplasmic (sn-, sno-, scRNAs, respectively), transfer (tRNAs), and ribosomal RNAs (rRNAs). Currently, DASHR (v1.0) integrates 187 smRNA high-throughput sequencing (smRNA-seq) datasets with over 2.5 billion reads and annotation data from multiple public sources. DASHR contains annotations for ∼48 000 human sncRNA genes and mature sncRNA products, 82% of which are expressed in one or more of the curated tissues. DASHR is available at http://lisanwanglab.org/DASHR.

tDRmapper: challenges and solutions to mapping, naming, and quantifying tRNA-derived RNAs from human small RNA-sequencing data

Background

Small RNA-sequencing has revealed the diversity and high abundance of small RNAs derived from tRNAs, referred to as tRNA-derived RNAs. However, at present, there is no standardized nomenclature and there are no methods for accurate annotation and quantification of these small RNAs. tRNA-derived RNAs have unique features that limit the utility of conventional alignment tools and quantification methods.

Results

We describe here the challenges of mapping, naming, and quantifying tRNA-derived RNAs and present a novel method that addresses them, called tDRmapper. We then use tDRmapper to perform a comparative analysis of tRNA-derived RNA profiles across different human cell types and diseases. We found that (1) tRNA-derived RNA profiles can differ dramatically across different cell types and disease states, (2) that positions and types of chemical modifications of tRNA-derived RNAs vary by cell type and disease, and (3) that entirely different tRNA-derived RNA species can be produced from the same parental tRNA depending on the cell type.

Conclusion

tDRmappernot only provides a standardized nomenclature and quantification scheme, but also includes graphical visualization that facilitates the discovery of novel tRNA and tRNA-derived RNA biology.

The expanding world of small RNAs in plants

Plant genomes encode various small RNAs that function in distinct, yet overlapping, genetic and epigenetic silencing pathways. However, the abundance and diversity of small-RNA classes varies among plant species, suggesting coevolution between environmental adaptations and gene-silencing mechanisms. Biogenesis of small RNAs in plants is well understood, but we are just beginning to uncover their intricate regulation and activity. Here, we discuss the biogenesis of plant small RNAs, such as microRNAs, secondary siRNAs and heterochromatic siRNAs, and their diverse cellular and developmental functions, including in reproductive transitions, genomic imprinting and paramutation. We also discuss the diversification of small-RNA-directed silencing pathways through the expansion of RNA-dependent RNA polymerases, DICER proteins and ARGONAUTE proteins.

Chemical Modifications Mark Alternatively Spliced and Uncapped Messenger RNAs in Arabidopsis

Posttranscriptional chemical modification of RNA bases is a widespread and physiologically relevant regulator of RNA maturation, stability, and function. While modifications are best characterized in short, noncoding RNAs such as tRNAs, growing evidence indicates that mRNAs and long noncoding RNAs (lncRNAs) are likewise modified. Here, we apply our high-throughput annotation of modified ribonucleotides (HAMR) pipeline to identify and classify modifications that affect Watson-Crick base pairing at three different levels of the Arabidopsis thaliana transcriptome (polyadenylated, small, and degrading RNAs). We find this type of modifications primarily within uncapped, degrading mRNAs and lncRNAs, suggesting they are the cause or consequence of RNA turnover. Additionally, modifications within stable mRNAs tend to occur in alternatively spliced introns, suggesting they regulate splicing. Furthermore, these modifications target mRNAs with coherent functions, including stress responses. Thus, our comprehensive analysis across multiple RNA classes yields insights into the functions of covalent RNA modifications in plant transcriptomes.

The Arabidopsis epitranscriptome

The most prevalent internal modification of plant messenger RNAs, N(6)-methyladenosine (m(6)A), was first discovered in the 1970s, then largely forgotten. However, the impact of modifications to eukaryote mRNA, collectively known as the epitranscriptome, has recently attracted renewed attention. mRNA methylation is required for normal Arabidopsis development and the first methylation maps reveal that thousands of Arabidopsis mRNAs are methylated. Arabidopsis is likely to be a model of wide utility in understanding the biological impacts of the epitranscriptome. We review recent progress and look ahead with questions awaiting answers to reveal an entire layer of gene regulation that has until recently been overlooked.

Transcriptome-wide measurement of plant RNA secondary structure

RNAs fold into intricate and precise secondary structures. These structural patterns regulate multiple steps of the RNA lifecycle, while also conferring catalytic and scaffolding functions to certain transcripts. Therefore, a full understanding of RNA posttranscriptional regulation requires a comprehensive picture of secondary structure. Here, we review several high throughput sequencing-based methods to globally survey plant RNA secondary structure. These methods are more accurate than computational prediction, and more scalable than physical techniques such as crystallography. We note hurdles to reliably measuring secondary structure, including RNA-binding proteins, RNA base modifications, and intramolecular duplexes. Finally, we survey the functional knowledge that has been gleaned from each of these methods, and identify some unanswered questions that remain.

Transfer RNA detection by small RNA deep sequencing and disease association with myelodysplastic syndromes

Background

Although advances in sequencing technologies have popularized the use of microRNA (miRNA) sequencing (miRNA-seq) for the quantification of miRNA expression, questions remain concerning the optimal methodologies for analysis and utilization of the data. The construction of a miRNA sequencing library selects RNA by length rather than type. However, as we have previously described, miRNAs represent only a subset of the species obtained by size selection. Consequently, the libraries obtained for miRNA sequencing also contain a variety of additional species of small RNAs. This study looks at the prevalence of these other species obtained from bone marrow aspirate specimens and explores the predictive value of these small RNAs in the determination of response to therapy in myelodysplastic syndromes (MDS).

Methods

Paired pre and post treatment bone marrow aspirate specimens were obtained from patients with MDS who were treated with either azacytidine or decitabine (24 pre-treatment specimens, 23 post-treatment specimens) with 22 additional non-MDS control specimens. Total RNA was extracted from these specimens and submitted for next generation sequencing after an additional size exclusion step to enrich for small RNAs. The species of small RNAs were enumerated, single nucleotide variants (SNVs) identified, and finally the differential expression of tRNA-derived species (tDRs) in the specimens correlated with diseasestatus and response to therapy.

Results

Using miRNA sequencing data generated from bone marrow aspirate samples of patients with known MDS (N = 47) and controls (N = 23), we demonstrated that transfer RNA (tRNA) fragments (specifically tRNA halves, tRHs) are one of the most common species of small RNA isolated from size selection. Using tRNA expression values extracted from miRNA sequencing data, we identified six tRNA fragments that are differentially expressed between MDS and normal samples. Using the elastic net method, we identified four tRNAs-derived small RNAs (tDRs) that together can explain 67 % of the variation in treatment response for MDS patients. Similar analysis of specifically mitochondrial tDRs (mt-tDRs) identified 13 mt-tDRs which distinguished disease status in the samples and a single mt-tDR which predited response. Finally, 14 SNVs within the tDRs were found in at least 20 % of the MDS samples and were not observed in any of the control specimens.

Discussion

This study highlights the prevalence of tDRs in RNA-seq studies focused on small RNAs. The potential etiologies of these species, both technical and biologic, are discussed as well as important challenges in the interpretation of tDR data.

Conclusions

Our analysis results suggest that tRNA fragments can be accurately detected through miRNA sequencing data and that the expression of these species may be useful in the diagnosis of MDS and the prediction of response to therapy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1929-y) contains supplementary material, which is available to authorized users.

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.

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.

Single-nucleotide-resolution mapping of m6A and m6Am throughout the transcriptome

N⁶-methyladenosine (m6A) is the most abundant modified base in eukaryotic mRNA and has been linked to diverse effects on mRNA fate. Current m6A mapping approaches localize m6A residues to 100–200 nt-long regions of transcripts. The precise position of m6A in mRNAs cannot be identified on a transcriptome-wide level because there are no chemical methods to distinguish between m6A and adenosine. Here we show that anti-m6A antibodies can induce specific mutational signatures at m6A residues after ultraviolet light-induced antibody-RNA crosslinking and reverse transcription. We find these antibodies similarly induce mutational signatures at N⁶,2′-O-dimethyladenosine (m6Am), a nucleotide found at the first encoded position of certain mRNAs. Using these mutational signatures, we map m6A and m6Am at single-nucleotide resolution in human and mouse mRNA and identify snoRNAs as a novel class of m6A-containing ncRNAs.

Inosine modifications in human tRNAs are incorporated at the precursor tRNA level

Transfer RNAs (tRNAs) are key adaptor molecules of the genetic code that are heavily modified post-transcriptionally. Inosine at the first residue of the anticodon (position 34; I34) is an essential widespread tRNA modification that has been poorly studied thus far. The modification in eukaryotes results from a deamination reaction of adenine that is catalyzed by the heterodimeric enzyme adenosine deaminase acting on tRNA (hetADAT), composed of two subunits: ADAT2 and ADAT3. Using high-throughput small RNA sequencing (RNAseq), we show that this modification is incorporated to human tRNAs at the precursor tRNA level and during maturation. We also functionally validated the human genes encoding for hetADAT and show that the subunits of this enzyme co-localize in nucleus in an ADAT2-dependent manner. Finally, by knocking down HsADAT2, we demonstrate that variations in the cellular levels of hetADAT will result in changes in the levels of I34 modification in all its potential substrates. Altogether, we present RNAseq as a powerful tool to study post-transcriptional tRNA modifications at the precursor tRNA level and give the first insights on the biology of I34 tRNA modification in metazoans.

Diverse cell stresses induce unique patterns of tRNA up- and down-regulation: tRNA-seq for quantifying changes in tRNA copy number

Emerging evidence points to roles for tRNA modifications and tRNA abundance in cellular stress responses. While isolated instances of stress-induced tRNA degradation have been reported, we sought to assess the effects of stress on tRNA levels at a systems level. To this end, we developed a next-generation sequencing method that exploits the paucity of ribonucleoside modifications at the 3′-end of tRNAs to quantify changes in all cellular tRNA molecules. Application of this tRNA-seq method to Saccharomyces cerevisiae identified all 76 expressed unique tRNA species out of 295 coded in the yeast genome, including all isoacceptor variants, with highly precise relative (fold-change) quantification of tRNAs. In studies of stress-induced changes in tRNA levels, we found that oxidation (H₂O₂) and alkylation (methylmethane sulfonate, MMS) stresses induced nearly identical patterns of up- and down-regulation for 58 tRNAs. However, 18 tRNAs showed opposing changes for the stresses, which parallels our observation of signature reprogramming of tRNA modifications caused by H₂O₂ and MMS. Further, stress-induced degradation was limited to only a small proportion of a few tRNA species. With tRNA-seq applicable to any organism, these results suggest that translational control of stress response involves a contribution from tRNA abundance.