Identifying disease-causing variants among a large number of single nucleotide variants (SNVs) is still a major challenge. Recently, N6-methyladenosine (m6A) has become a research hotspot because of its critical roles in many fundamental biological processes and a variety of diseases. Therefore, it is important to evaluate the effect of variants on m6A modification, in order to gain a better understanding of them. Here, we report m6AVar (http://m6avar.renlab.org), a comprehensive database of m6A-associated variants that potentially influence m6A modification, which will help to interpret variants by m6A function. The m6A-associated variants were derived from three different m6A sources including miCLIP/PA-m6A-seq experiments (high confidence), MeRIP-Seq experiments (medium confidence) and transcriptome-wide predictions (low confidence). Currently, m6AVar contains 16 132 high, 71 321 medium and 326 915 low confidence level m6A-associated variants. We also integrated the RBP-binding regions, miRNA-targets and splicing sites associated with variants to help users investigate the effect of m6A-associated variants on post-transcriptional regulation. Because it integrates the data from genome-wide association studies (GWAS) and ClinVar, m6AVar is also a useful resource for investigating the relationship between the m6A-associated variants and disease. Overall, m6AVar will serve as a useful resource for annotating variants and identifying disease-causing variants.

The N 6 -methyladenosine (m 6 A) modification influences various mRNA metabolic events and tumorigenesis, however, its functions in nonsense-mediated mRNA decay (NMD) and whether NMD detects induced carcinogenesis pathways remain undefined. Here, we showed that the m 6 A methyltransferase METTL3 sustained its oncogenic role by modulating NMD of splicing factors and alternative splicing isoform switches in glioblastoma (GBM). Methylated RNA immunoprecipitation-seq (MeRIP-seq) analyses showed that m 6 A modification peaks were enriched at metabolic pathwayrelated transcripts in glioma stem cells (GSC) compared with neural progenitor cells. In addition, the clinical aggressiveness of malignant gliomas was associated with elevated expression of METTL3. Furthermore, silencing METTL3 or overexpressing dominant-negative mutant METTL3 suppressed the growth and self-renewal of GSCs. Integrated transcriptome and MeRIP-seq analyses revealed that downregulating the expression of METTL3 decreased m 6 A modification levels of serineand arginine-rich splicing factors (SRSF), which led to YTHDC1-dependent NMD of SRSF transcripts and decreased SRSF protein expression. Reduced expression of SRSFs led to larger changes in alternative splicing isoform switches. Importantly, the phenotypes mediated by METTL3 deficiency could be rescued by downregulating BCL-X or NCOR2 isoforms. Overall, these results establish a novel function of m 6 A in modulating NMD and uncover the mechanism by which METTL3 promotes GBM tumor growth and progression.Significance: These findings establish the oncogenic role of m 6 A writer METTL3 in glioblastoma stem cells.

No abstract