Identification of differential RNA modifications from nanopore direct RNA sequencing with xPore.

RNA modifications, such as N6-methyladenosine (m6A), modulate functions of cellular RNA species. However, quantifying differences in RNA modifications has been challenging. Here we develop a computational method, xPore, to identify differential RNA modifications from nanopore direct RNA sequencing (RNA-seq) data. We evaluate our method on transcriptome-wide m6A profiling data, demonstrating that xPore identifies positions of m6A sites at single-base resolution, estimates the fraction of modified RNA species in the cell and quantifies the differential modification rate across conditions. We apply xPore to direct RNA-seq data from six cell lines and multiple myeloma patient samples without a matched control sample and find that many m6A sites are preserved across cell types, whereas a subset exhibit significant differences in their modification rates. Our results show that RNA modifications can be identified from direct RNA-seq data with high accuracy, enabling analysis of differential modifications and expression from a single high-throughput experiment.

METTL4 catalyzes m6Am methylation in U2 snRNA to regulate pre-mRNA splicing.

N 6-methylation of 2'-O-methyladenosine (Am) in RNA occurs in eukaryotic cells to generate N6,2'-O-dimethyladenosine (m6Am). Identification of the methyltransferase responsible for m6Am catalysis has accelerated studies on the function of m6Am in RNA processing. While m6Am is generally found in the first transcribed nucleotide of mRNAs, the modification is also found internally within U2 snRNA. However, the writer required for catalyzing internal m6Am formation had remained elusive. By sequencing transcriptome-wide RNA methylation at single-base-resolution, we identified human METTL4 as the writer that directly methylates Am at U2 snRNA position 30 into m6Am. We found that METTL4 localizes to the nucleus and its conserved methyltransferase catalytic site is required for U2 snRNA methylation. By sequencing human cells with overexpressed Mettl4, we determined METTL4's in vivo target RNA motif specificity. In the absence of Mettl4 in human cells, U2 snRNA lacks m6Am thereby affecting a subset of splicing events that exhibit specific features such as 3' splice-site weakness and an increase in exon inclusion. These findings suggest that METTL4 methylation of U2 snRNA regulates splicing of specific pre-mRNA transcripts.

Atlas of quantitative single-base-resolution N6-methyl-adenine methylomes.

Various methyltransferases and demethylases catalyse methylation and demethylation of N6-methyladenosine (m6A) and N6,2'-O-dimethyladenosine (m6Am) but precise methylomes uniquely mediated by each methyltransferase/demethylase are still lacking. Here, we develop m6A-Crosslinking-Exonuclease-sequencing (m6ACE-seq) to map transcriptome-wide m6A and m6Am at quantitative single-base-resolution. This allows for the generation of a comprehensive atlas of distinct methylomes uniquely mediated by every individual known methyltransferase or demethylase. Our atlas reveals METTL16 to indirectly impact manifold methylation targets beyond its consensus target motif and highlights the importance of precision in mapping PCIF1-dependent m6Am. Rather than reverse RNA methylation, we find that both ALKBH5 and FTO instead maintain their regulated sites in an unmethylated steady-state. In FTO's absence, anomalous m6Am disrupts snRNA interaction with nuclear export machinery, potentially causing aberrant pre-mRNA splicing events.

Single-nucleotide-resolution sequencing of human N6-methyldeoxyadenosine reveals strand-asymmetric clusters associated with SSBP1 on the mitochondrial genome.

N6-methyldeoxyadenosine (6mA) is a well-characterized DNA modification in prokaryotes but reports on its presence and function in mammals have been controversial. To address this issue, we established the capacity of 6mA-Crosslinking-Exonuclease-sequencing (6mACE-seq) to detect genome-wide 6mA at single-nucleotide-resolution, demonstrating this by accurately mapping 6mA in synthesized DNA and bacterial genomes. Using 6mACE-seq, we generated a human-genome-wide 6mA map that accurately reproduced known 6mA enrichment at active retrotransposons and revealed mitochondrial chromosome-wide 6mA clusters asymmetrically enriched on the heavy-strand. We identified a novel putative 6mA-binding protein in single-stranded DNA-binding protein 1 (SSBP1), a mitochondrial DNA (mtDNA) replication factor known to coat the heavy-strand, linking 6mA with the regulation of mtDNA replication. Finally, we characterized AlkB homologue 1 (ALKBH1) as a mitochondrial protein with 6mA demethylase activity and showed that its loss decreases mitochondrial oxidative phosphorylation. Our results show that 6mA clusters play a previously unappreciated role in regulating human mitochondrial function, despite 6mA being an uncommon DNA modification in the human genome.