Detection of internal N7-methylguanosine (m7G) RNA modifications by mutational profiling sequencing

Abstract Methylation of guanosine on position N7 (m7G) on internal RNA positions has been found in all domains of life and have been implicated in human disease. Here, we present m7G Mutational Profiling sequencing (m7G-MaP-seq), which allows high throughput detection of m7G modifications at nucleotide resolution. In our method, m7G modified positions are converted to abasic sites by reduction with sodium borohydride, directly recorded as cDNA mutations through reverse transcription and sequenced. We detect positions with increased mutation rates in the reduced and control samples taking the possibility of sequencing/alignment error into account and use replicates to calculate statistical significance based on log likelihood ratio tests. We show that m7G-MaP-seq efficiently detects known m7G modifications in rRNA with mutational rates up to 25% and we map a previously uncharacterised evolutionarily conserved rRNA modification at position 1581 in Arabidopsis thaliana SSU rRNA. Furthermore, we identify m7G modifications in budding yeast, human and arabidopsis tRNAs and demonstrate that m7G modification occurs before tRNA splicing. We do not find any evidence for internal m7G modifications being present in other small RNA, such as miRNA, snoRNA and sRNA, including human Let-7e. Likewise, high sequencing depth m7G-MaP-seq analysis of mRNA from E. coli or yeast cells did not identify any internal m7G modifications.

[1]  Todd M. Lowe,et al.  ARM-Seq: AlkB-facilitated RNA methylation sequencing reveals a complex landscape of modified tRNA fragments , 2015, Nature Methods.

[2]  Lulu Hu,et al.  Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA. , 2019, Molecular cell.

[3]  Aldema Sas-Chen,et al.  A Chemical Signature for Cytidine Acetylation in RNA. , 2018, Journal of the American Chemical Society.

[4]  V. Heurgué-Hamard,et al.  The human 18S rRNA base methyltransferases DIMT1L and WBSCR22-TRMT112 but not rRNA modification are required for ribosome biogenesis , 2015, Molecular biology of the cell.

[5]  S. Douthwaite,et al.  Identifying modifications in RNA by MALDI mass spectrometry. , 2007, Methods in enzymology.

[6]  J. Bujnicki,et al.  MODOMICS: a database of RNA modification pathways—2013 update , 2012, Nucleic Acids Res..

[7]  Anya V Grozhik,et al.  Distinguishing RNA modifications from noise in epitranscriptome maps. , 2018, Nature chemical biology.

[8]  Anders H. Lund,et al.  Profiling of 2′-O-Me in human rRNA reveals a subset of fractionally modified positions and provides evidence for ribosome heterogeneity , 2016, Nucleic acids research.

[9]  S. Douthwaite,et al.  Structure of the bifunctional methyltransferase YcbY (RlmKL) that adds the m7G2069 and m2G2445 modifications in Escherichia coli 23S rRNA , 2012, Nucleic acids research.

[10]  M. Bohnsack,et al.  WBSCR 22 / Merm 1 is required for late nuclear pre-ribosomal RNA processing and mediates N 7-methylation of G 1639 in human 18 S rRNA , 2014 .

[11]  P. Moore,et al.  The crystal structure of yeast phenylalanine tRNA at 1.93 A resolution: a classic structure revisited. , 2000, RNA.

[12]  Marc Graille,et al.  Structural and functional studies of Bud23–Trm112 reveal 18S rRNA N7-G1575 methylation occurs on late 40S precursor ribosomes , 2014, Proceedings of the National Academy of Sciences.

[13]  Heng Li,et al.  Improving SNP discovery by base alignment quality , 2011, Bioinform..

[14]  Chengqi Yi,et al.  Base-Resolution Mapping Reveals Distinct m1A Methylome in Nuclear- and Mitochondrial-Encoded Transcripts. , 2017, Molecular cell.

[15]  John S. Mattick,et al.  The RNA modification landscape in human disease , 2017, RNA.

[16]  Anna Tramontano,et al.  Accurate energies of hydrogen bonded nucleic acid base pairs and triplets in tRNA tertiary interactions , 2006, Nucleic acids research.

[17]  W. Wintermeyer,et al.  A specific chemical chain scission of tRNA at 7‐methylguanosine , 1970, FEBS letters.

[18]  Gary L. Glish,et al.  Tandem Mass Spectrometry of Small, Multiply Charged Oligonucleotides , 1992, Journal of the American Society for Mass Spectrometry.

[19]  Nan Yu,et al.  The Comparative RNA Web (CRW) Site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs , 2002, BMC Bioinformatics.

[20]  Isabelle Hazemann,et al.  Visualization of chemical modifications in the human 80S ribosome structure. , 2017 .

[21]  F. Kirpekar,et al.  Detection of pseudouridine and other modifications in tRNA by cyanoethylation and MALDI mass spectrometry. , 2002, Nucleic acids research.

[22]  M. Taoka,et al.  Landscape of the complete RNA chemical modifications in the human 80S ribosome , 2018, Nucleic acids research.

[23]  Janusz M. Bujnicki,et al.  Bud23 Methylates G1575 of 18S rRNA and Is Required for Efficient Nuclear Export of Pre-40S Subunits , 2008, Molecular and Cellular Biology.

[24]  B. Klaholz,et al.  Visualization of chemical modifications in the human 80S ribosome structure , 2017, Nature.

[25]  Andrew J. Bannister,et al.  METTL1 Promotes let-7 MicroRNA Processing via m7G Methylation , 2019, Molecular cell.

[26]  W. Wintermeyer,et al.  Tertiary structure interactions of 7‐methylguanosine in yeast tRNAPhe as studied by borohydride reduction , 1975, FEBS letters.

[27]  Victor S Lelyveld,et al.  Mettl1/Wdr4-Mediated m7G tRNA Methylome Is Required for Normal mRNA Translation and Embryonic Stem Cell Self-Renewal and Differentiation. , 2018, Molecular cell.

[28]  K. Ochi,et al.  Loss of a conserved 7‐methylguanosine modification in 16S rRNA confers low‐level streptomycin resistance in bacteria , 2007, Molecular microbiology.

[29]  Chengqi Yi,et al.  Transcriptome-wide mapping reveals reversible and dynamic N(1)-methyladenosine methylome. , 2016, Nature chemical biology.

[30]  B. Porse,et al.  A novel partial modification at C2501 in Escherichia coli 23S ribosomal RNA. , 2004, RNA.

[31]  Nikolaos Sidiropoulos,et al.  Reproducible Analysis of Sequencing-Based RNA Structure Probing Data with User-Friendly Tools. , 2015, Methods in enzymology.

[32]  Schraga Schwartz,et al.  The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution , 2017, Nature.

[33]  Y. Takano,et al.  RNA recognition mechanism of eukaryote tRNA (m7G46) methyltransferase (Trm8–Trm82 complex) , 2007, FEBS letters.

[34]  Yuri Motorin,et al.  AlkAniline-Seq: Profiling of m7 G and m3 C RNA Modifications at Single Nucleotide Resolution. , 2018, Angewandte Chemie.

[35]  Andreas Hildebrandt,et al.  The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent , 2015, Nucleic acids research.

[36]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[37]  Mineko Konishi,et al.  Genetic identification of Arabidopsis RID2 as an essential factor involved in pre-rRNA processing. , 2011, The Plant journal : for cell and molecular biology.

[38]  Anders Albrechtsen,et al.  ANGSD: Analysis of Next Generation Sequencing Data , 2014, BMC Bioinformatics.

[39]  Steven Busan,et al.  RNA motif discovery by SHAPE and mutational profiling (SHAPE-MaP) , 2014, Nature Methods.

[40]  Chengqi Yi,et al.  Epitranscriptome sequencing technologies: decoding RNA modifications , 2016, Nature Methods.

[41]  Isabelle Behm-Ansmant,et al.  Use of Specific Chemical Reagents for Detection of Modified Nucleotides in RNA , 2011, Journal of nucleic acids.

[42]  Gideon Rechavi,et al.  The dynamic N1-methyladenosine methylome in eukaryotic messenger RNA , 2016, Nature.

[43]  M. Bohnsack,et al.  WBSCR22/Merm1 is required for late nuclear pre-ribosomal RNA processing and mediates N7-methylation of G1639 in human 18S rRNA , 2015, RNA.

[44]  R. Wolfenden,et al.  1-Methyladenosine. Dimroth rearrangement and reversible reduction. , 1968, Biochemistry.