6mA-Sniper: Quantifying 6mA Sites in Eukaryotes at Single-Nucleotide Resolution

While N6-methyldeoxyadenine (6mA) modification has been linked to fundamental regulatory processes in prokaryotes, its prevalence and functional implications in eukaryotes are controversial. Here, we report 6mA-Sniper to quantify 6mA sites in eukaryotes at single-nucleotide resolution. With 6mA-Sniper, we delineated an accurate 6mA profile in C. elegans with 2,034 sites, significantly enriched on sequences of [GC]GAG motif. Twenty-six of 39 6mA events with MnlI restriction endonuclease sites were experimentally verified, demonstrating the feasibility of this method. Notably, the enrichment of these 6mA sites on a specific sequence motif, their within-population conservation and the combinatorial patterns, and the selective constrains on them jointly support an active model for the shaping of the profile by some undiscovered methyltransferases. In a joint study (Cell Research, in revision), Ma et al. reported METL-9 as a new methyltransferase in shaping the basal and stress-dependent 6mA profile in C. elegans. Notably, with the 6mA profile identified by 6mA-Sniper at single-nucleotide resolution, we found that the levels of 6mAs are significantly decreased in strains with the removal of METL-9 (METL-9 KO-OP50), while generally increased after P. aeruginosa infection, further verified the efficiency of 6mA-Sniper in accurately pinpointing 6mA sites. Moreover, for the regions marked by 998 6mA sites emerged specifically after the infection, we identified an enrichment of the upregulated genes after the infection. The gene upregulations are likely mediated through a mutual exclusive crosstalk between 6mA and H3K27me3 modification, as supported by their co-occurrence, and the signal of increased H3K27me3 at regions marked by 6mAs depleted in METL-9 KO-OP50 strains. Notably, in different C. elegans strains, the cross-strain genetic variants removing 6mA sites are associated with the decreased expression of their host genes, and the removal of two randomly-selected 6mA events with genome editing directly decreased the expression of their host genes. We thus highlight 6mA regulation as a previously-neglected regulator of transcriptome in eukaryotes.