Mammalian DNA N6-methyladenosine: Challenges and new insights.

DNA N6-methyldeoxyadenosine (6mA) modification was first discovered in Bacterium coli in the 1950s. Over the next several decades, 6mA was recognized as a critical DNA modification in the genomes of prokaryotes and protists. While important in prokaryotes, less is known about the presence and functional roles of DNA 6mA in eukaryotes, particularly in mammals. Taking advantage of recent technology advances that made 6mA detection and sequencing possible, studies over the past several years have brought new insights into 6mA biology in mammals. In this perspective, we present recent progress, discuss challenges, and pose four questions for future research regarding mammalian DNA 6mA.

N6-Deoxyadenosine Methylation in Mammalian Mitochondrial DNA.

N6-Methyldeoxyadenosine (6mA) has recently been shown to exist and play regulatory roles in eukaryotic genomic DNA (gDNA). However, the biological functions of 6mA in mammals have yet to be adequately explored, largely due to its low abundance in most mammalian genomes. Here, we report that mammalian mitochondrial DNA (mtDNA) is enriched for 6mA. The level of 6mA in HepG2 mtDNA is at least 1,300-fold higher than that in gDNA under normal growth conditions, corresponding to approximately four 6mA modifications on each mtDNA molecule. METTL4, a putative mammalian methyltransferase, can mediate mtDNA 6mA methylation, which contributes to attenuated mtDNA transcription and a reduced mtDNA copy number. Mechanistically, the presence of 6mA could repress DNA binding and bending by mitochondrial transcription factor (TFAM). Under hypoxia, the 6mA level in mtDNA could be further elevated, suggesting regulatory roles for 6mA in mitochondrial stress response. Our study reveals DNA 6mA as a regulatory mark in mammalian mtDNA.

Randomly incorporated genomic N6-methyldeoxyadenosine delays zygotic transcription initiation in a cnidarian.

N6-methyldeoxyadenosine (6mA) is a chemical alteration of DNA, observed across all realms of life. Although the functions of 6mA are well understood in bacteria and protists, its roles in animal genomes have been controversial. We show that 6mA randomly accumulates in early embryos of the cnidarian Hydractinia symbiolongicarpus, with a peak at the 16-cell stage followed by clearance to background levels two cell cycles later, at the 64-cell stage-the embryonic stage at which zygotic genome activation occurs in this animal. Knocking down Alkbh1, a putative initiator of animal 6mA clearance, resulted in higher levels of 6mA at the 64-cell stage and a delay in the initiation of zygotic transcription. Our data are consistent with 6mA originating from recycled nucleotides of degraded m6A-marked maternal RNA postfertilization. Therefore, while 6mA does not function as an epigenetic mark in Hydractinia, its random incorporation into the early embryonic genome inhibits transcription. In turn, Alkbh1 functions as a genomic 6mA "cleaner," facilitating timely zygotic genome activation. Given the random nature of genomic 6mA accumulation and its ability to interfere with gene expression, defects in 6mA clearance may represent a hitherto unknown cause of various pathologies.

An Adversarial DNA N6-Methyladenine-Sensor Network Preserves Polycomb Silencing.

Adenine N6 methylation in DNA (6mA) is widespread among bacteria and phage and is detected in mammalian genomes, where its function is largely unexplored. Here we show that 6mA deposition and removal are catalyzed by the Mettl4 methyltransferase and Alkbh4 dioxygenase, respectively, and that 6mA accumulation in genic elements corresponds with transcriptional silencing. Inactivation of murine Mettl4 depletes 6mA and causes sublethality and craniofacial dysmorphism in incross progeny. We identify distinct 6mA sensor domains of prokaryotic origin within the MPND deubiquitinase and ASXL1, a component of the Polycomb repressive deubiquitinase (PR-DUB) complex, both of which act to remove monoubiquitin from histone H2A (H2A-K119Ub), a repressive mark. Deposition of 6mA by Mettl4 triggers the proteolytic destruction of both sensor proteins, preserving genome-wide H2A-K119Ub levels. Expression of the bacterial 6mA methyltransferase Dam, in contrast, fails to destroy either sensor. These findings uncover a native, adversarial 6mA network architecture that preserves Polycomb silencing.

N6-Methyladenine DNA Modification in the Human Genome.

DNA N6-methyladenine (6mA) modification is the most prevalent DNA modification in prokaryotes, but whether it exists in human cells and whether it plays a role in human diseases remain enigmatic. Here, we showed that 6mA is extensively present in the human genome, and we cataloged 881,240 6mA sites accounting for ∼0.051% of the total adenines. [G/C]AGG[C/T] was the most significantly associated motif with 6mA modification. 6mA sites were enriched in the coding regions and mark actively transcribed genes in human cells. DNA 6mA and N6-demethyladenine modification in the human genome were mediated by methyltransferase N6AMT1 and demethylase ALKBH1, respectively. The abundance of 6mA was significantly lower in cancers, accompanied by decreased N6AMT1 and increased ALKBH1 levels, and downregulation of 6mA modification levels promoted tumorigenesis. Collectively, our results demonstrate that DNA 6mA modification is extensively present in human cells and the decrease of genomic DNA 6mA promotes human tumorigenesis.

Epigenetic marks or not? The discovery of novel DNA modifications in eukaryotes.

DNA modifications add another layer of complexity to the eukaryotic genome to regulate gene expression, playing critical roles as epigenetic marks. In eukaryotes, the study of DNA epigenetic modifications has been confined to 5mC and its derivatives for decades. However, rapid developing approaches have witnessed the expansion of DNA modification reservoirs during the past several years, including the identification of 6mA, 5gmC, 4mC, and 4acC in diverse organisms. However, whether these DNA modifications function as epigenetic marks requires careful consideration. In this review, we try to present a panorama of all the DNA epigenetic modifications in eukaryotes, emphasizing recent breakthroughs in the identification of novel DNA modifications. The characterization of their roles in transcriptional regulation as potential epigenetic marks is summarized. More importantly, the pathways for generating or eliminating these DNA modifications, as well as the proteins involved are comprehensively dissected. Furthermore, we briefly discuss the potential challenges and perspectives, which should be taken into account while investigating novel DNA modifications.

DNA N6-methyldeoxyadenosine in mammals and human disease.

N6-methyladenine (6mA) is the most prevalent DNA modification in prokaryotes. However, its presence and significance in eukaryotes remain elusive. Recently, with methodology advances in detection and sequencing of 6mA in eukaryotes, 6mA is back in the spotlight. Although multiple studies have reported that 6mA is an important epigenetic mark in eukaryotes and plays a regulatory role in DNA transcription, transposon activation, stress response, and other bioprocesses, there are some discrepancies in the current literature. We review the recent advances in 6mA research in eukaryotes, especially in mammals. In particular, we describe the abundance/distribution of 6mA, its potential role in regulating gene expression, identified regulators, and pathological roles in human diseases, especially in cancer. The limitations faced by the field and future perspectives in 6mA research are also discussed.

DNA-MP: a generalized DNA modifications predictor for multiple species based on powerful sequence encoding method.

Accurate prediction of deoxyribonucleic acid (DNA) modifications is essential to explore and discern the process of cell differentiation, gene expression and epigenetic regulation. Several computational approaches have been proposed for particular type-specific DNA modification prediction. Two recent generalized computational predictors are capable of detecting three different types of DNA modifications; however, type-specific and generalized modifications predictors produce limited performance across multiple species mainly due to the use of ineffective sequence encoding methods. The paper in hand presents a generalized computational approach "DNA-MP" that is competent to more precisely predict three different DNA modifications across multiple species. Proposed DNA-MP approach makes use of a powerful encoding method "position specific nucleotides occurrence based 117 on modification and non-modification class densities normalized difference" (POCD-ND) to generate the statistical representations of DNA sequences and a deep forest classifier for modifications prediction. POCD-ND encoder generates statistical representations by extracting position specific distributional information of nucleotides in the DNA sequences. We perform a comprehensive intrinsic and extrinsic evaluation of the proposed encoder and compare its performance with 32 most widely used encoding methods on $17$ benchmark DNA modifications prediction datasets of $12$ different species using $10$ different machine learning classifiers. Overall, with all classifiers, the proposed POCD-ND encoder outperforms existing $32$ different encoders. Furthermore, combinedly over 5-fold cross validation benchmark datasets and independent test sets, proposed DNA-MP predictor outperforms state-of-the-art type-specific and generalized modifications predictors by an average accuracy of 7% across 4mc datasets, 1.35% across 5hmc datasets and 10% for 6ma datasets. To facilitate the scientific community, the DNA-MP web application is available at https://sds_genetic_analysis.opendfki.de/DNA_Modifications/.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut.

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

BERT6mA: prediction of DNA N6-methyladenine site using deep learning-based approaches.

N6-methyladenine (6mA) is associated with important roles in DNA replication, DNA repair, transcription, regulation of gene expression. Several experimental methods were used to identify DNA modifications. However, these experimental methods are costly and time-consuming. To detect the 6mA and complement these shortcomings of experimental methods, we proposed a novel, deep leaning approach called BERT6mA. To compare the BERT6mA with other deep learning approaches, we used the benchmark datasets including 11 species. The BERT6mA presented the highest AUCs in eight species in independent tests. Furthermore, BERT6mA showed higher and comparable performance with the state-of-the-art models while the BERT6mA showed poor performances in a few species with a small sample size. To overcome this issue, pretraining and fine-tuning between two species were applied to the BERT6mA. The pretrained and fine-tuned models on specific species presented higher performances than other models even for the species with a small sample size. In addition to the prediction, we analyzed the attention weights generated by BERT6mA to reveal how the BERT6mA model extracts critical features responsible for the 6mA prediction. To facilitate biological sciences, the BERT6mA online web server and its source codes are freely accessible at https://github.com/kuratahiroyuki/BERT6mA.git, respectively.

Structural Basis of Nucleic Acid Recognition and 6mA Demethylation by Caenorhabditis elegans NMAD-1A.

N6-methyladenine (6mA) of DNA is an emerging epigenetic mark in the genomes of Chlamydomonas, Caenorhabditis elegans, and mammals recently. Levels of 6mA undergo drastic fluctuation and thus affect fertility during meiosis and early embryogenesis. Here, we showed three complex structures of 6mA demethylase C. elegans NMAD-1A, a canonical isoform of NMAD-1 (F09F7.7). Biochemical results revealed that NMAD-1A prefers 6mA Bubble or Bulge DNAs. Structural studies of NMAD-1A revealed an unexpected "stretch-out" conformation of its Flip2 region, a conserved element that is usually bent over the catalytic center to facilitate substrate base flipping in other DNA demethylases. Moreover, the wide channel between the Flip1 and Flip2 of the NMAD-1A explained the observed preference of NMAD-1A for unpairing substrates, of which the flipped 6mA was primed for catalysis. Structural analysis and mutagenesis studies confirmed that key elements such as carboxy-terminal domain (CTD) and hypothetical zinc finger domain (ZFD) critically contributed to structural integrity, catalytic activity, and nucleosome binding. Collectively, our biochemical and structural studies suggest that NMAD-1A prefers to regulate 6mA in the unpairing regions and is thus possibly associated with dynamic chromosome regulation and meiosis regulation.

DNA methylation on C5-Cytosine and N6-Adenine in the Bursaphelenchus xylophilus genome.

BACKGROUND: The pinewood nematode is the causal agent of the pine wilt disease, which causes severe ecological and economic losses in coniferous forests. The invasion of pine wood nematode has undergone various rapid adaptations to a wide range of temperatures and to new hosts and vector insects. DNA methylation may play crucial roles in the rapid adaptation of PWN during invasion. However, whether the PWN genome contins functional DNA modifications remains elusive. RESULTS: Here, we detected the extensive presence of 5-methylcytosine (5mC) and N6-methyladenine (6mA) in the B. xylophilus genome, with low methylation levels at most positions. Cytosines were methylated in the CpG, CHG. and CHH sequence contexts, with the lowest methylation levels at CpG sites. The methylation levels of CpG and 6mA in gene regions showed opposite trends. The changes in the abundance of 5mC and 6mA showed the same trends in response to temperature change, but opposite trends during development. Sequence and phylogenetic analyses showed that the proteins BxDAMT and BxNMAD have typical characteristics of a methylase and demethylase, respectively, and are conserved among species. CONCLUSIONS: These findings shed light on the epigenetic modifications present in the genome of PWN, and will improve our understanding of its invasiveness and evolution.

Leveraging the attention mechanism to improve the identification of DNA N6-methyladenine sites.

DNA N6-methyladenine is an important type of DNA modification that plays important roles in multiple biological processes. Despite the recent progress in developing DNA 6mA site prediction methods, several challenges remain to be addressed. For example, although the hand-crafted features are interpretable, they contain redundant information that may bias the model training and have a negative impact on the trained model. Furthermore, although deep learning (DL)-based models can perform feature extraction and classification automatically, they lack the interpretability of the crucial features learned by those models. As such, considerable research efforts have been focused on achieving the trade-off between the interpretability and straightforwardness of DL neural networks. In this study, we develop two new DL-based models for improving the prediction of N6-methyladenine sites, termed LA6mA and AL6mA, which use bidirectional long short-term memory to respectively capture the long-range information and self-attention mechanism to extract the key position information from DNA sequences. The performance of the two proposed methods is benchmarked and evaluated on the two model organisms Arabidopsis thaliana and Drosophila melanogaster. On the two benchmark datasets, LA6mA achieves an area under the receiver operating characteristic curve (AUROC) value of 0.962 and 0.966, whereas AL6mA achieves an AUROC value of 0.945 and 0.941, respectively. Moreover, an in-depth analysis of the attention matrix is conducted to interpret the important information, which is hidden in the sequence and relevant for 6mA site prediction. The two novel pipelines developed for DNA 6mA site prediction in this work will facilitate a better understanding of the underlying principle of DL-based DNA methylation site prediction and its future applications.

Deep6mAPred: A CNN and Bi-LSTM-based deep learning method for predicting DNA N6-methyladenosine sites across plant species.

DNA N6-methyladenine (6mA) is a key DNA modification, which plays versatile roles in the cellular processes, including regulation of gene expression, DNA repair, and DNA replication. DNA 6mA is closely associated with many diseases in the mammals and with growth as well as development of plants. Precisely detecting DNA 6mA sites is of great importance to exploration of 6mA functions. Although many computational methods have been presented for DNA 6mA prediction, there is still a wide gap in the practical application. We presented a convolution neural network (CNN) and bi-directional long-short term memory (Bi-LSTM)-based deep learning method (Deep6mAPred) for predicting DNA 6mA sites across plant species. The Deep6mAPred stacked the CNNs and the Bi-LSTMs in a paralleling manner instead of a series-connection manner. The Deep6mAPred also employed the attention mechanism for improving the representations of sequences. The Deep6mAPred reached an accuracy of 0.9556 over the independent rice dataset, far outperforming the state-of-the-art methods. The tests across plant species showed that the Deep6mAPred is of a remarkable advantage over the state of the art methods. We developed a user-friendly web application for DNA 6mA prediction, which is freely available at http://106.13.196.152:7001/ for all the scientific researchers. The Deep6mAPred would enrich tools to predict DNA 6mA sites and speed up the exploration of DNA modification.

DNA adenine methylation in eukaryotes: Enzymatic mark or a form of DNA damage?

6-methyladenine (6mA) is fairly abundant in nuclear DNA of basal fungi, ciliates and green algae. In these organisms, 6mA is maintained near transcription start sites in ApT context by a parental-strand instruction dependent maintenance methyltransferase and is positively associated with transcription. In animals and plants, 6mA levels are high only in organellar DNA. The 6mA levels in nuclear DNA are very low. They are attributable to nucleotide salvage and the activity of otherwise mitochondrial METTL4, and may be considered as a price that cells pay for adenine methylation in RNA and/or organellar DNA. Cells minimize this price by sanitizing dNTP pools to limit 6mA incorporation, and by converting 6mA that has been incorporated into DNA back to adenine. Hence, 6mA in nuclear DNA should be described as an epigenetic mark only in basal fungi, ciliates and green algae, but not in animals and plants.

N6-methyladenine: A Rare and Dynamic DNA Mark.

Chromatin, consisting of deoxyribonucleic acid (DNA) wrapped around histone proteins, facilitates DNA compaction and allows identical DNA code to confer many different cellular phenotypes. This biological versatility is accomplished in large part by post-translational modifications to histones and chemical modifications to DNA. These modifications direct the cellular machinery to expand or compact specific chromatin regions and mark certain regions of the DNA as important for cellular functions. While each of the four bases that make up DNA can be modified (Iyer et al., Prog Mol Biol Transl Sci. 101:25-104, 2011), this chapter will focus on methylation of the 6th position on adenines (6mA). 6mA is a prevalent modification in unicellular organisms and until recently was thought to be restricted to them. A flurry of conflicting studies have proposed that 6mA either does not exist, is present at low levels, or is present at relatively high levels and regulates complex processes in different multicellular eukaryotes. Here, we will briefly describe the history of 6mA, examine its evolutionary conservation, and evaluate the current methods for detecting 6mA. We will discuss the proteins that have been reported to bind and regulate 6mA and examine the known and potential functions of this modification in eukaryotes. Finally, we will close with a discussion of the ongoing debate about whether 6mA exists as a directed DNA modification in multicellular eukaryotes.

A Mutation-Based Method for Pinpointing a DNA N6 -Methyladenine Methyltransferase Modification Site at Single Base Resolution.

DNA N6 -methyladenine (6mA) has recently received notable attention due to an increased finding of its functional roles in higher eukaryotes. Here we report an enzyme-assisted chemical labeling method to pinpoint the DNA 6mA methyltransferase (MTase) substrate modification site at single base resolution. A designed allyl-substituted MTase cofactor was applied in the catalytic transfer reaction, and the allyl group was installed to the N6 -position of adenine within a specific DNA sequence to form N6 -allyladenine (6aA). The iodination of 6aA allyl group induced the formation of 1, N6 -cyclized adenine which caused mutations during DNA replication by a polymerase. Thus the modification site could be precisely detected by a mutation signal. We synthesized 6aA deoxynucleoside and deoxynucleotide model compounds and a 6aA-containing DNA probe, and screened nine DNA polymerases to define an optimal system capable of detecting the substrate modification site of a DNA 6mA MTase at single-base resolution.

Interplay of DNA and RNA N 6 -methyladenine with R-loops in regulating gene transcription in Arabidopsis.

R-loops and covalent modifications of N 6 -methyladenine on DNA (D-6 mA) or RNA (R-m6A) have been documented to function in various cellular processes in eukaryotes. However, the relationships between R-loops and both covalent modifications are still elusive in plants. Here, we integrated existing ssDRIP-seq with D-6 mA and R-m6A data from Arabidopsis thaliana. We found that the presence of either of both modifications facilitates R-loop formation and transcription of overlapping genes. Interestingly, our study suggests that the presence of R-m6A is key to affect R-loop intensity and positively regulate gene transcription. Moreover, the presence of D-6 mA plays an additive role to facilitate the effect of R-m6A on R-loop intensity, however, D-6 mA may negatively regulate gene transcription when coexisted with R-m6A. Our analyses indicate that D-6 mA, R-m6A, or histone marks may act individually and cooperatively with R-loops in regulating gene transcription. Our study is the first to link R-loops with D-6 mA and R-m6A in plants, thereby providing new insights into interactions between R-loops with D-6 mA, R-m6A, and histone marks for regulating gene transcription. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-01010-5.

i6mA-Fuse: improved and robust prediction of DNA 6 mA sites in the Rosaceae genome by fusing multiple feature representation.

DNA N6-methyladenine (6 mA) is one of the most vital epigenetic modifications and involved in controlling the various gene expression levels. With the avalanche of DNA sequences generated in numerous databases, the accurate identification of 6 mA plays an essential role for understanding molecular mechanisms. Because the experimental approaches are time-consuming and costly, it is desirable to develop a computation model for rapidly and accurately identifying 6 mA. To the best of our knowledge, we first proposed a computational model named i6mA-Fuse to predict 6 mA sites from the Rosaceae genomes, especially in Rosa chinensis and Fragaria vesca. We implemented the five encoding schemes, i.e., mononucleotide binary, dinucleotide binary, k-space spectral nucleotide, k-mer, and electron-ion interaction pseudo potential compositions, to build the five, single-encoding random forest (RF) models. The i6mA-Fuse uses a linear regression model to combine the predicted probability scores of the five, single encoding-based RF models. The resultant species-specific i6mA-Fuse achieved remarkably high performances with AUCs of 0.982 and 0.978 and with MCCs of 0.869 and 0.858 on the independent datasets of Rosa chinensis and Fragaria vesca, respectively. In the F. vesca-specific i6mA-Fuse, the MBE and EIIP contributed to 75% and 25% of the total prediction; in the R. chinensis-specific i6mA-Fuse, Kmer, MBE, and EIIP contribute to 15%, 65%, and 20% of the total prediction. To assist high-throughput prediction for DNA 6 mA identification, the i6mA-Fuse is publicly accessible at https://kurata14.bio.kyutech.ac.jp/i6mA-Fuse/.

The complex phylogenetic relationships of a 4mC/6mA DNA methyltransferase in prokaryotes.

DNA methyltransferases are proteins that modify DNA via attachment of methyl groups to nucleobases and are ubiquitous across the bacterial, archaeal, and eukaryotic domains of life. Here, we investigated the complex evolutionary history of the large and consequential 4mC/6mA DNA methyltransferase protein family using phylogenetic reconstruction of amino acid sequences. We present a well-supported phylogeny of this family based on systematic sampling of taxa across superphyla of bacteria and archaea. We compared the phylogeny to a current representation of the species tree of life and found that the 4mC/6mA methyltransferase family has a strikingly complex evolutionary history that likely began sometime after the last universal common ancestor of life diverged into the bacterial and archaeal lineages and probably involved many horizontal gene transfers within and between domains. Despite the complexity of its evolutionary history, we inferred that only one significant shift in molecular evolutionary rate characterizes the diversification of this protein family.