High-throughput sequencing for 1-methyladenosine (m(1)A) mapping in RNA.

Detection and mapping of modified nucleotides in RNAs is a difficult and laborious task. Several physico-chemical approaches based on differential properties of modified nucleotides can be used, however, most of these methods do not allow high-throughput analysis. Here we describe in details a method for mapping of rather common 1-methyladenosine (m(1)A) residues using high-throughput next generation sequencing (NGS). Since m(1)A residues block primer extension during reverse transcription (RT), the accumulation of abortive products as well as the nucleotide misincorporation can be detected in the sequencing data. The described library preparation protocol allows to capture both types of cDNA products essential for further bioinformatic analysis. We demonstrate that m(1)A residues produce characteristic arrest and mismatch rates and combination of both can be used for their detection as well as for discrimination of m(1)A from other modified A residues present in RNAs.

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(1)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(1)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(1)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(1)A in the template RNA, meaning it is sequence dependent. The RT-signature of m(1)A was used for inspection and confirmation of suspected modification sites and resulted in the identification of hitherto unknown m(1)A residues in trypanosomal tRNA.

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.

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.