Abstract
Cellular RNAs are subject to a myriad of different chemical modifications that play important roles in controlling RNA expression and function. Dysregulation of certain RNA modifications, the so-called 'epitranscriptome', contributes to human disease. One limitation in studying the functional, physiological, and pathological roles of the epitranscriptome is the availability of methods for the precise mapping of individual RNA modifications throughout the transcriptome. 3-Methylcytidine (m3C) modification of certain tRNAs is well established and was also recently detected in mRNA. However, methods for the specific mapping of m3C throughout the transcriptome are lacking. Here, we developed a m3C-specific technique, Hydrazine-Aniline Cleavage sequencing (HAC-seq), to profile the m3C methylome at single-nucleotide resolution. We applied HAC-seq to analyze ribosomal RNA (rRNA)-depleted total RNAs in human cells. We found that tRNAs are the predominant m3C-modified RNA species, with 17 m3C modification sites on 11 cytoplasmic and 2 mitochondrial tRNA isoacceptors in MCF7 cells. We found no evidence for m3C-modification of mRNA or other non-coding RNAs at comparable levels to tRNAs in these cells. HAC-seq provides a novel method for the unbiased, transcriptome-wide identification of m3C RNA modification at single-nucleotide resolution, and could be widely applied to reveal the m3C methylome in different cells and tissues.