6mA-Sniper: Quantifying 6mA sites in eukaryotes at single-nucleotide resolution

While N(6)-methyldeoxyadenine (6mA) modification is a fundamental regulation in prokaryotes, its prevalence and functions in eukaryotes are controversial. Here, we report 6mA-Sniper to quantify 6mA sites in eukaryotes at single-nucleotide resolution, and delineate a 6mA profile in Caenorhabditis elegans with 2034 sites. Twenty-six of 39 events with Mnl I restriction endonuclease sites were verified, demonstrating the feasibility of this method. The levels of 6mA sites pinpointed by 6mA-Sniper are generally increased after Pseudomonas aeruginosa infection, but decreased in strains with the removal of METL-9, the dominant 6mA methyltransferase. The enrichment of these sites on specific motif of [GC]GAG, the selective constrains on them, and their coordinated changes with METL-9 levels thus support an active shaping of the 6mA profile by methyltransferase. Moreover, for regions marked by 6mA sites that emerged after infection, an enrichment of up-regulated genes was detected, possibly mediated through a mutual exclusive cross-talk between 6mA and H3K27me3 modification. We thus highlight 6mA regulation as a previously neglected regulator in eukaryotes.