Characterization of eukaryotic DNA N(6)-methyladenine by a highly sensitive restriction enzyme-assisted sequencing

Ense-i6mA: Identification of DNA N6-Methyladenine Sites Using XGB-RFE Feature Selection and Ensemble Machine Learning

DNA N6-methyladenine (6mA) is an important epigenetic modification that plays a vital role in various cellular processes. Accurate identification of the 6mA sites is fundamental to elucidate the biological functions and mechanisms of modification. However, experimental methods for detecting 6mA sites are high-priced and time-consuming. In this study, we propose a novel computational method, called Ense-i6mA, to predict 6mA sites. Firstly, five encoding schemes, i.e., one-hot encoding, gcContent, Z-Curve, K-mer nucleotide frequency, and K-mer nucleotide frequency with gap, are employed to extract DNA sequence features. Secondly, eXtreme gradient boosting coupled with recursive feature elimination is applied to remove noisy features for avoiding over-fitting, reducing computing time and complexity. Then, the best subset of features is fed into base-classifiers composed of Extra Trees, eXtreme Gradient Boosting, Light Gradient Boosting Machine, and Support Vector Machine. Finally, to minimize generalization errors, the prediction probabilities of the base-classifiers are aggregated by averaging for inferring the final 6mA sites results. We conduct experiments on two species, i.e., Arabidopsis thaliana and Drosophila melanogaster, to compare the performance of Ense-i6mA against the recent 6mA sites prediction methods. The experimental results demonstrate that the proposed Ense-i6mA achieves area under the receiver operating characteristic curve values of 0.967 and 0.968, accuracies of 91.4% and 92.0%, and Mathew's correlation coefficient values of 0.829 and 0.842 on two benchmark datasets, respectively, and outperforms several existing state-of-the-art methods.

G-Quadruplex-Mediated Transcriptional Regulation of SYT7: Implications for Tumor Progression and Therapeutic Strategies

Synaptotagmin 7 (SYT7), a member of the synaptotagmin family, exhibits high expression in various tumors and is closely associated with patient prognosis. The tight regulation of SYT7 expression assumes paramount significance in the progression of tumorigenesis. In this study, we detected a high GC content in the first 1000 bp of the promoter region of SYT7, suggesting a potential role of the G-quadruplex in its transcriptional regulation. Circular dichroism spectroscopy results showed that -187 to -172 bp sequence can form a typical parallel G-quadruplex structure, and site mutation revealed the critical role of the ninth guanine in its formation. Then, treatment of two ligands of G-quadruplex (TMPyP4 and Pyridostatin) reduced both the expression of SYT7 and subsequent tumor proliferation, demonstrating the potential of the G-quadruplex as a targeted therapy for tumors. By shedding light on the pivotal role of the G-quadruplex in regulating SYT7 transcription, our study not only advances our comprehension of this intricate regulatory mechanism but also emphasizes the significance of SYT7 in tumor proliferation. These findings collectively contribute to a more comprehensive understanding of the interplay between G-quadruplex regulation and SYT7 function in tumor development.

Deep learning based method for predicting DNA N6-methyladenosine sites

DNA N6 methyladenine (6mA) plays an important role in many biological processes, and accurately identifying its sites helps one to understand its biological effects more comprehensively. Previous traditional experimental methods are very labor-intensive and traditional machine learning methods also seem to be somewhat insufficient as the database of 6mA methylation groups becomes progressively larger, so we propose a deep learning-based method called multi-scale convolutional model based on global response normalization (CG6mA) to solve the prediction problem of 6mA site. This method is tested with other methods on three different kinds of benchmark datasets, and the results show that our model can get more excellent prediction results.

Development of SPQC sensor based on the specific recognition of TAL-effectors for locus-specific detection of 6-methyladenine in DNA

N6-methyladenosine (6mA) plays a pivotal role in diverse biological processes, including cancer, bacterial toxin secretion, and bacterial drug resistance. However, to date there has not been a selective, sensitive, and simple method for quantitative detection of 6mA at single base resolution. Herein, we present a series piezoelectric quartz crystal (SPQC) sensor based on the specific recognition of transcription-activator-like effectors (TALEs) for locus-specific detection of 6mA. Detection sensitivity is enhanced through the use of a hybridization chain reaction (HCR) in conjunction with silver staining. The limit of detection (LOD) of the sensor was 0.63 pM and can distinguish single base mismatches. We demonstrate the applicability of the sensor platform by quantitating 6mA DNA at a specific site in biological matrix. The SPQC sensor presented herein offers a promising platform for in-depth study of cancer, bacterial toxin secretion, and bacterial drug resistance.

Ready for renascence in mosquito: The regulation of gene expression in Plasmodium sexual development

Malaria remains one of the most perilous vector-borne infectious diseases for humans globally. Sexual gametocyte represents the exclusive stage at which malaria parasites are transmitted from the vertebrate to the Anopheles host. The feasible and effective approach to prevent malaria transmission is by addressing the sexual developmental processes, that is, gametocytogenesis and gametogenesis. Thus, this review will comprehensively cover advances in the regulation of gene expression surrounding the transmissible stages, including epigenetic, transcriptional, and post-transcriptional control.

Silver-Coordinated Watson-Crick Pairing-Driven Three-Dimensional DNA Walker for Locus-Specific Detection of Genomic N6-Methyladenine and N4-Methylcytosine at the Single-Molecule Level

N6-Methyladenine (6mdA) and N4-methylcytosine (4mdC) are the two most dominant DNA modifications in both prokaryotes and eukaryotes, but standard hybridization-based techniques cannot be applied for the 6mdA/4mdC assay. Herein, we demonstrate the silver-coordinated Watson-Crick pairing-driven three-dimensional (3D) DNA walker for locus-specific detection of genomic 6mdA/4mdC at the single-molecule level. 6mdA-DNA and 4mdC-DNA can selectively hybridize with the binding probes (BP1 and BP2) to form 6mdA-DNA-BP1 and 4mdC-DNA-BP2 duplexes. The 6mdA-C/4mdC-A mismatches cannot be stabilized by AgI, and thus, 18-nt BP1/BP2 cannot be extended by the catalysis of KF exonuclease. Through toehold-mediated strand displacement (TMSD), the signal probe (SP1/SP2) functionalized on the gold nanoparticles (AuNPs) can competitively bind to BP1/BP2 in 6mdA-DNA-BP1/4mdC-DNA-BP2 duplex to obtain SP1-18-nt BP1 and SP2-18-nt BP2 duplexes. The resulting DNA duplexes can act as the substrates of lambda exonuclease, leading to the cleavage of SP1/SP2 and the release of Cy3/Cy5 and 18-nt BP1/BP2. The released 18-nt BP1/BP2 can subsequently serve as the walker DNA, moving along the surface of the AuNP to activate dynamic 3D DNA walking and releasing abundant Cy3/Cy5. The released Cy3/Cy5 can be quantified by single-molecule imaging. This nanosensor exhibits high sensitivity with a limit of detection (LOD) of 9.80 × 10-15 M for 6mdA-DNA and 9.97 × 10-15 M for 4mdC-DNA. It can discriminate 6mdA-/4mdC-DNA from unmodified genomic DNAs, distinguish 0.01% 6mdA-/4mdC-DNA from excess unmethylated DNAs, and quantify 6mdA-/4mdC-DNA at specific sites in genomic DNAs of liver cancer cells and Escherichia coli plasmid cloning vector, providing a new platform for locus-specific analysis of 6mdA/4mdC in genomic DNAs.

Sequencing of N6-methyl-deoxyadenosine at single-base resolution across the mammalian genome

Although DNA N6-methyl-deoxyadenosine (6mA) is abundant in bacteria and protists, its presence and function in mammalian genomes have been less clear. We present Direct-Read 6mA sequencing (DR-6mA-seq), an antibody-independent method, to measure 6mA at base resolution. DR-6mA-seq employs a unique mutation-based strategy to reveal 6mA sites as misincorporation signatures without any chemical or enzymatic modulation of 6mA. We validated DR-6mA-seq through the successful mapping of the well-characterized G(6mA)TC motif in the E. coli DNA. As expected, when applying DR-6mA-seq to mammalian systems, we found that genomic DNA (gDNA) 6mA abundance is generally low in most mammalian tissues and cells; however, we did observe distinct gDNA 6mA sites in mouse testis and glioblastoma cells. DR-6mA-seq provides an enabling tool to detect 6mA at single-base resolution for a comprehensive understanding of DNA 6mA in eukaryotes.

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

While N6-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.

Single-Nucleotide Resolution Mapping of N6-Methyladenine in Genomic DNA

N6-Methyladenine (6mA) is a naturally occurring DNA modification in both prokaryotes and eukaryotes. Herein, we developed a deaminase-mediated sequencing (DM-seq) method for genome-wide mapping of 6mA at single-nucleotide resolution. The method capitalizes on the selective deamination of adenine, but not 6mA, in DNA mediated by an evolved adenine deaminase, ABE8e. By employing this method, we achieved genome-wide mapping of 6mA in Escherichia coli and in mammalian mitochondrial DNA (mtDNA) at single-nucleotide resolution. We found that the 6mA sites are mainly located in the GATC motif in the E. coli genome. We also identified 17 6mA sites in mtDNA of HepG2 cells, where all of the 6mA sites are distributed in the heavy strand of mtDNA. We envision that DM-seq will be a valuable tool for uncovering new functions of 6mA in DNA and for exploring its potential roles in mitochondria-related human diseases.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

Force-Induced Visualization of Nucleic Acid Functions with Single-Nucleotide Resolution

Nucleic acids are major targets for molecular sensing because of their wide involvement in biological functions. Determining their presence, movement, and binding specificity is thus well pursued. However, many current techniques are usually sophisticated, expensive, and often lack single-nucleotide resolution. In this paper, we report the force-induced visualization method that relies on the novel concept of mechanical force to determine the functional positions of nucleic acids with single-nucleotide resolution. The use of an adjustable mechanical force overcomes the variation of analyte concentration and differences in buffer conditions that are common in biological settings. Two examples are described to validate the method: one is probing the mRNA movement during ribosomal translocation, and the other is revealing the interacting sites and strengths of DNA-binding drugs based on the force amplitude. The flexibility of the method, simplicity of the associated device, and capability of multiplexed detection will potentially enable a broad range of biomedical applications.

The Global Changes of N6-methyldeoxyadenosine in Response to Low Temperature in Arabidopsis thaliana and Rice

N6-methyldeoxyadenosine (6mA) is a recently discovered DNA modification involved in regulating plant adaptation to abiotic stresses. However, the mechanisms and changes of 6mA under cold stress in plants are not yet fully understood. Here, we conducted a genome-wide analysis of 6mA and observed that 6mA peaks were predominantly present within the gene body regions under both normal and cold conditions. In addition, the global level of 6mA increased both in Arabidopsis and rice after the cold treatment. The genes that exhibited an up-methylation showed enrichment in various biological processes, whereas there was no significant enrichment observed among the down-methylated genes. The association analysis revealed a positive correlation between the 6mA level and the gene expression level. Joint analysis of the 6mA methylome and transcriptome of Arabidopsis and rice unraveled that fluctuations in 6mA levels caused by cold exposure were not correlated to changes in transcript levels. Furthermore, we discovered that orthologous genes modified by 6mA showed high expression levels; however, only a minor amount of differentially 6mA-methylated orthologous genes were shared between Arabidopsis and rice under low-temperature conditions. In conclusion, our study provides information on the role of 6mA in response to cold stress and reveals its potential for regulating the expression of stress-related genes.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

A Unique m6A-Dependent Restriction Endonuclease from an Archaeal Virus

Prokaryotes possess numerous diverse defense systems to resist viral infections, while some viruses have also evolved antiviral defense systems to exclude other viruses in cases of multiple infections. Here, we report the first virus-derived modification-dependent restriction endonuclease (HHPV4I) from the archaeal virus HHPV4 (Haloarcula hispanica pleomorphic virus 4). HHPV4I contains an SRA domain, a winged helix (wH) domain, and an HNH domain; recognizes the Gm6ATC site; and specifically binds to Gm6ATC site-containing DNA. Both the wH domain and the HNH domain are responsible for DNA binding. Unlike the well-known m6A-specific restriction enzyme DpnI, HHPV4I only efficiently cleaves DNA with a fully methylated Gm6ATC site and cleaves DNA both upstream and downstream of the Gm6ATC sites on both DNA strands. Furthermore, HHPV4I preferentially cleaves DNA between VR bases (V = A/G/C, R = A/G) 4 to 20 nt away from the Gm6ATC site. Thus, the cleavage pattern of HHPV4I is distinct from those of all of the presently characterized restriction endonucleases. Mutations in the wH domain of HHPV4I do not alter m6A-dependent endonuclease activity, but they decrease recognition sequence specificity, thus expanding the cleaving capacity to more m6A-containing DNA sequences. The wH domain provides a target for searching, developing, and engineering novel m6A-dependent endonucleases. IMPORTANCE Many modification-dependent restriction endonucleases (MDREs) were identified in prokaryotes and recognized modified cytosine bases, such as 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), and glucosyl-5-hydroxymethylcytosine (g5hmC). The first virus-derived MDRE (HHPV4I) from the archaeal virus HHPV4 was identified in this study. The viral MDRE suggested a new strategy employed by the virus to exclude other viruses in the case of multiple replications. HHPV4I is a novel N6-methyladenine (m6A)-dependent restriction endonuclease, while the cleavage pattern of HHPV4I is distinct from the well-known m6A-dependent restriction endonuclease DpnI. HHPV4I recognizes Gm6ATC sites and cleaves DNA both upstream and downstream of the Gm6ATC sites on both DNA strands. It preferentially cleaves DNA between VR bases (V = A/G/C, R = A/G) 4 to 20 nt away from the Gm6ATC sites. Furthermore, mutations in the HHPV4I wH domain can alter the sequence specificity without impeding the m6A-dependent DNA cleavage activity, providing a target for engineering more m6A-dependent endonucleases with different sequence specificities.

Legionellapneumophila induces methylomic changes in ten-eleven translocation to ensure bacterial reproduction in human lung epithelial cells

Introduction. Legionella pneumophila is a Gram-negative flagellated bacteria that can infect human lungs and cause a severe form of pneumonia named Legionnaires' disease.Hypothesis. We hypothesize that L. pneumophila infection induces methylomic changes in methylcytosine dioxygenases, ten-eleven translocation (TET) genes, and controls DNA methylation following infection.Aim. In the current research, we sought to further investigate DNA methylation changes in human lung epithelial cells upon L. pneumophila infection and determine how methylation inhibitor agents disturb L. pneumophila reproduction.Methodology. A549 cell line was used in L. pneumophila infection and inhibitors' treatment, including 5-azacytidine (5-AZA) and (-)-epigallocatechin-3-O-gallate (EGCG).Results. Interestingly, DNA methylation analysis of infected A549 using sodium bisulfite PCR and the methylation-sensitive HpaII enzyme showed potential methylation activity within the promoter regions of ten-eleven translocation (TET) genes located on CpG/397-8 and CpG/385-6 of TET1 and TET3, respectively. Such methylation changes in TET effectors decreased their expression profile following infection, indicated by quantitative real-time PCR (RT-qPCR), immunoblotting and flow cytometry. Furthermore, pre-treatment of A549 cells with 5-AZA or EGCG significantly decreased the bacterial reproduction characterized by the expression of L. pneumophila 16S ribosomal RNA and the c.f.u. ml^-1 of bacterial particles. Moreover, both methylation inhibitors showed potent inhibition of methionine synthase (MS) expression, which was further confirmed by the docking analysis of inhibitor ligands and crystal structure of MS protein.Conclusion. These data provide evidence for the methylomic changes in the promoter region of TET1 and TET3 by L. pneumophila infection in the A549 cell line and suggest the anti-bacterial properties of 5-AZA and EGCG, as methylation inhibitors, are due to targeting the epigenetic effector methionine synthase.

I-DNAN6mA: Accurate Identification of DNA N6-Methyladenine Sites Using the Base-Pairing Map and Deep Learning

The recent discovery of numerous DNA N⁶-methyladenine (6mA) sites has transformed our perception about the roles of 6mA in living organisms. However, our ability to understand them is hampered by our inability to identify 6mA sites rapidly and cost-efficiently by existing experimental methods. Developing a novel method to quickly and accurately identify 6mA sites is critical for speeding up the progress of its function detection and understanding. In this study, we propose a novel computational method, called I-DNAN6mA, to identify 6mA sites and complement experimental methods well, by leveraging the base-pairing rules and a well-designed three-stage deep learning model with pairwise inputs. The performance of our proposed method is benchmarked and evaluated on four species, i.e., Arabidopsis thaliana, Drosophila melanogaster, Rice, and Rosaceae. The experimental results demonstrate that I-DNAN6mA achieves area under the receiver operating characteristic curve values of 0.967, 0.963, 0.947, 0.976, and 0.990, accuracies of 91.5, 92.7, 88.2, 0.938, and 96.2%, and Mathew's correlation coefficient values of 0.855, 0.831, 0.763, 0.877, and 0.924 on five benchmark data sets, respectively, and outperforms several existing state-of-the-art methods. To our knowledge, I-DNAN6mA is the first approach to identify 6mA sites using a novel image-like representation of DNA sequences and a deep learning model with pairwise inputs. I-DNAN6mA is expected to be useful for locating functional regions of DNA.

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.

A review of methods for predicting DNA N6-methyladenine sites

Deoxyribonucleic acid(DNA) N6-methyladenine plays a vital role in various biological processes, and the accurate identification of its site can provide a more comprehensive understanding of its biological effects. There are several methods for 6mA site prediction. With the continuous development of technology, traditional techniques with the high costs and low efficiencies are gradually being replaced by computer methods. Computer methods that are widely used can be divided into two categories: traditional machine learning and deep learning methods. We first list some existing experimental methods for predicting the 6mA site, then analyze the general process from sequence input to results in computer methods and review existing model architectures. Finally, the results were summarized and compared to facilitate subsequent researchers in choosing the most suitable method for their work.

High-precision mapping reveals rare N6-deoxyadenosine methylation in the mammalian genome

N6-deoxyadenosine methylation (6mA) is the most widespread type of DNA modification in prokaryotes and is also abundantly distributed in some unicellular eukaryotes. However, 6mA levels are remarkably low in mammals. The lack of a precise and comprehensive mapping method has hindered more advanced investigations of 6mA. Here, we report a new method MM-seq (modification-induced mismatch sequencing) for genome-wide 6mA mapping based on a novel detection principle. We found that modified DNA bases are prone to form a local open region that allows capture by antibody, for example, via a DNA breathing or base-flipping mechanism. Specified endonuclease or exonuclease can recognize the antibody-stabilized mismatch-like structure and mark the exact modified sites for sequencing readout. Using this method, we examined the genomic positions of 6mA in bacteria (E. coli), green algae (C. reinhardtii), and mammalian cells (HEK239T, Huh7, and HeLa cells). In contrast to bacteria and green algae, human cells possess a very limited number of 6mA sites which are sporadically distributed across the genome of different cell types. After knocking out the RNA m6A methyltransferase METTL3 in mouse ES cells, 6mA becomes mostly diminished. Our results imply that rare 6mA in the mammalian genome is introduced by RNA m6A machinery via a non-targeted mechanism.

Light Dependent Changes in Adenylate Methylation of the Promoter of the Mitochondrial Citrate Synthase Gene in Maize (Zea mays L.) Leaves

Limited methyl-specific restriction of genomic DNA by endonuclease MAL1 revealed the changes in its methyl status caused by adenine modification in maize (Zea mays L.) leaves under different light conditions (dark, light, irradiation by red and far-red light). Incubation in the light and irradiation by red light exhibited an activating effect on DNA adenine methylase activity, which was reflected in an increase in the number of methylated adenines in GATC sites. Far-red light and darkness exhibited an opposite effect. The use of nitrite conversion of DNA followed by methyladenine-dependent restriction by MboI nuclease revealed a phytochrome B-dependent mechanism of regulation of the methyl status of adenine in the GATC sites in the promoter of the gene encoding the mitochondrial isoform of citrate synthase. Irradiation of plants with red light caused changes in the adenine methyl status of the analyzed amplicon, as evidenced by the presence of restriction products of 290, 254, and 121 nucleotides. Adenine methylation occurred at all three GATC sites in the analyzed DNA sequence. It is concluded that adenylate methylation is controlled by phytochrome B via the transcription factor PIF4 and represents an important mechanism for the tricarboxylic acid cycle regulation by light.

Association of vitamin B12/ferritin deficiency in cancer patients with methylomic changes at promotors of TET methylcytosine dioxygenases

Aim: To investigate potential DNA methylation in methylcytosine dioxygenases and correlation of TET genes with vitamin B12/ferritin levels in cancer patients. Materials & methods: 200 blood samples were obtained from both cancer patients and healthy individuals. Results: The expression of DNMT1, DNMT3a and DNMT3b was increased in patients with low vitamin B12 and ferritin levels, while the expression of MTR, TET1 and TET3 significantly decreased. DNA methylation analysis in patients with deficient vitamin B12/ferritin levels showed methylomic changes within the location 318/CG and 385/CG in the promoter region of TET1 and TET3, respectively. Conclusion: Vitamin B12/ferritin deficiency contributes to DNA methylation progress in cancer patients.

Mapping single-nucleotide m6A by m6A-REF-seq

The past few years have witnessed rapid progress in the field of RNA modifications. As the most prevailing modification on eukaryotic mRNA, m6A is characterized to play a vital role in various cellular activities. However, limitations of the detection method impede functional studies of m6A. Here we introduce m6A-REF-seq, a powerful and straightforward method to identify m6A at single-nucleotide resolution. m6A-REF-seq relies on the recognition of RNA endonuclease MazF towards m6A at the ACA motif, providing an orthogonal method independent of the m6A antibody being adopted by most of current methods. We describe a detailed protocol to perform m6A-REF-seq, including NGS library construction and sequencing data analysis. In particular, we describe an optimized assay to validate individual m6A sites identified by m6A-REF-seq, which can also be applied to detect any candidate m6A sites.

Same modification, different location: the mythical role of N6-adenine methylation in plant genomes

The present review summarizes recent advances in the understanding of 6mA in DNA as an emergent epigenetic mark with distinctive characteristics, discusses its importance in plant genomes, and highlights its chemical nature and functions. Adenine methylation is an epigenetic modification present in DNA (6mA) and RNA (m6A) that has a regulatory function in many cellular processes. This modification occurs through a reversible reaction that covalently binds a methyl group, usually at the N⁶ position of the purine ring. This modification carries biophysical properties that affect the stability of nucleic acids as well as their binding affinity with other molecules. DNA 6mA has been related to genome stability, gene expression, DNA replication, and repair mechanisms. Recent advances have shown that 6mA in plant genomes is related to development and stress response. In this review, we present recent advances in the understanding of 6mA in DNA as an emergent epigenetic mark with distinctive characteristics. We discuss the key elements of this modification, focusing mainly on its importance in plant genomes. Furthermore, we highlight its chemical nature and the regulatory effects that it exerts on gene expression and plant development. Finally, we emphasize the functions of 6mA in photosynthesis, stress, and flowering.

Genome-wide DNA N6-adenine methylation in sea buckthorn (Hippophae rhamnoides L.) fruit development

As a new epigenetic mark, DNA N6-adenine (6mA) methylation plays an important role in various biological processes and has been reported in many prokaryotic organisms in recent years. However, the distribution patterns and functions of DNA 6mA modification have been poorly studied in non-model crops. In this study, we observed that the methylation ratio of 6mA was about 0.016% in the sea buckthorn (Hippophae rhamnoides L.) genome using mass spectrometry. We first constructed a comprehensive 6mA landscape in sea buckthorn genome using nanopore sequencing at single-base resolution. Distribution analysis suggested that 6mA methylated sites were widely distributed in the sea buckthorn chromosomes, which were similar to those in Arabidopsis and rice. Furthermore, reduced 6mA DNA methylation is associated with different expression of genes related to the fruit-ripening process in sea buckthorn. Our results revealed that 6mA DNA modification could be considered an important epigenomic mark and contributes to the fruit ripening process in plants.

NT-seq: a chemical-based sequencing method for genomic methylome profiling

DNA methylation plays vital roles in both prokaryotes and eukaryotes. There are three forms of DNA methylation in prokaryotes: N⁶-methyladenine (6mA), N⁴-methylcytosine (4mC), and 5-methylcytosine (5mC). Although many sequencing methods have been developed to sequence specific types of methylation, few technologies can be used for efficiently mapping multiple types of methylation. Here, we present NT-seq for mapping all three types of methylation simultaneously. NT-seq reliably detects all known methylation motifs in two bacterial genomes and can be used for identifying de novo methylation motifs. NT-seq provides a simple and efficient solution for detecting multiple types of DNA methylation.

Critical assessment of DNA adenine methylation in eukaryotes using quantitative deconvolution

The discovery of N6-methyldeoxyadenine (6mA) across eukaryotes led to a search for additional epigenetic mechanisms. However, some studies have highlighted confounding factors that challenge the prevalence of 6mA in eukaryotes. We developed a metagenomic method to quantitatively deconvolve 6mA events from a genomic DNA sample into species of interest, genomic regions, and sources of contamination. Applying this method, we observed high-resolution 6mA deposition in two protozoa. We found that commensal or soil bacteria explained the vast majority of 6mA in insect and plant samples. We found no evidence of high abundance of 6mA in Drosophila, Arabidopsis, or humans. Plasmids used for genetic manipulation, even those from Dam methyltransferase mutant Escherichia coli, could carry abundant 6mA, confounding the evaluation of candidate 6mA methyltransferases and demethylases. On the basis of this work, we advocate for a reassessment of 6mA in eukaryotes.

Deep transformers and convolutional neural network in identifying DNA N6-methyladenine sites in cross-species genomes

As one of the most common post-transcriptional epigenetic modifications, N6-methyladenine (6 mA), plays an essential role in various cellular processes and disease pathogenesis. Therefore, accurately identifying 6 mA modifications is necessary for a deep understanding of cellular processes and other possible functional mechanisms. Although a few computational methods have been proposed, their respective models were developed with small training datasets. Hence, their practical application is quite limited in genome-wide detection. To overcome the existing limitations, we present a novel model based on transformer architecture and deep learning to identify DNA 6 mA sites from the cross-species genome. The model is constructed on a benchmark dataset and explored a feature derived from pre-trained transformer word embedding approaches. Subsequently, a convolutional neural network was employed to learn the generated features and generate the prediction outcomes. As a result, our predictor achieved excellent performance during independent test with the accuracy and Matthews correlation coefficient (MCC) of 79.3% and 0.58, respectively. Overall, its performance achieved better accuracy than the baseline models and significantly outperformed the existing predictors, demonstrating the effectiveness of our proposed hybrid framework. Furthermore, our model is expected to assist biologists in accurately identifying 6mAs and formulate the novel testable biological hypothesis. We also release source codes and datasets freely at https://github.com/khanhlee/bert-dna for front-end users.

Femtomolar and locus-specific detection of N6-methyladenine in DNA by integrating double-hindered replication and nucleic acid-functionalized MB@Zr-MOF

In this study, a novel electrochemical biosensor was constructed for ultrasensitive and locus-specific detection of N⁶-Methyladenine (m⁶A) in DNA using double-hindered replication and nucleic acid-coated methylene blue (MB)@Zr-MOF. Based on the combination of m⁶A-impeded replication and Ag^I-mediated mismatch replication, this mode could effectively stop the extension of the strand once DNA polymerase encountered m⁶A site, which specifically distinguish the m⁶A site from natural A site in DNA. Also, Zr-MOF with high porosity and negative surface potential features was carefully chose to load cationic MB, resulting a stable and robust MB@Zr-MOF electrochemical tag. As a result, the developed biosensor exhibited a wide linear range from 1 fM to 1 nM with detection limit down to 0.89 fM. Profiting from the high sensitivity and selectivity, the biosensing strategy revealed good applicability, which had been demonstrated by quantitating m⁶A DNA at specific site in biological matrix. Thus, the biosensor provides a promising platform for locus-specific m⁶A DNA analysis.

The exploration of N6-deoxyadenosine methylation in mammalian genomes

N⁶-methyladenine (N⁶-mA, m⁶dA, or 6mA), a prevalent DNA modification in prokaryotes, has recently been identified in higher eukaryotes, including mammals. Although 6mA has been well-studied in prokaryotes, the function and regulatory mechanism of 6mA in eukaryotes are still poorly understood. Recent studies indicate that 6mA can serve as an epigenetic mark and play critical roles in various biological processes, from transposable-element suppression to environmental stress response. Here, we review the significant advances in methodology for 6mA detection and major progress in understanding the regulation and function of this non-canonical DNA methylation in eukaryotes, predominantly mammals.

Advances in mapping the epigenetic modifications of 5-methylcytosine (5mC), N6-methyladenine (6mA), and N4-methylcytosine (4mC)

DNA modification plays a pivotal role in regulating gene expression in cell development. As prevalent markers on DNA, 5-methylcytosine (5mC), N6-methyladenine (6mA), and N4-methylcytosine (4mC) can be recognized by specific methyltransferases, facilitating cellular defense and the versatile regulation of gene expression in eukaryotes and prokaryotes. Recent advances in DNA sequencing technology have permitted the positions of different modifications to be resolved at the genome-wide scale, which has led to the discovery of several novel insights into the complexity and functions of multiple methylations. In this review, we summarize differences in the various mapping approaches and discuss their pros and cons with respect to their relative read depths, speeds, and costs. We also discuss the development of future sequencing technologies and strategies for improving the detection resolution of current sequencing technologies. Lastly, we speculate on the potentially instrumental role that these sequencing technologies might play in future research.

Quantification and mapping of DNA modifications

Apart from the four canonical nucleobases, DNA molecules carry a number of natural modifications. Substantial evidence shows that DNA modifications can regulate diverse biological processes. Dynamic and reversible modifications of DNA are critical for cell differentiation and development. Dysregulation of DNA modifications is closely related to many human diseases. The research of DNA modifications is a rapidly expanding area and has been significantly stimulated by the innovations of analytical methods. With the recent advances in methods and techniques, a series of new DNA modifications have been discovered in the genomes of prokaryotes and eukaryotes. Deciphering the biological roles of DNA modifications depends on the sensitive detection, accurate quantification, and genome-wide mapping of modifications in genomic DNA. This review provides an overview of the recent advances in analytical methods and techniques for both the quantification and genome-wide mapping of natural DNA modifications. We discuss the principles, advantages, and limitations of these developed methods. It is anticipated that new methods and techniques will resolve the current challenges in this burgeoning research field and expedite the elucidation of the functions of DNA modifications.

Highly sensitive detection of 6mA at single-base resolution based on A-C mismatch

We first demonstrated that 6mA can be selectively recognized based on the selective ligation reaction of DNA ligase toward A-C mismatch and 6mA-C mismatch. This method, when further combined with amplification using RCA, achieved highly sensitive identification of 6mA in dsDNA at single-base resolution.

Case Study of the Response of N6-Methyladenine DNA Modification to Environmental Stressors in the Unicellular Eukaryote Tetrahymena thermophila

Rediscovered as a potential epigenetic mark, N⁶-methyladenine DNA modification (6mA) was recently reported to be sensitive to environmental stressors in several multicellular eukaryotes. As 6mA distribution and function differ significantly in multicellular and unicellular organisms, whether and how 6mA in unicellular eukaryotes responds to environmental stress remains elusive. Here, we characterized the dynamic changes of 6mA under starvation in the unicellular model organism Tetrahymena thermophila. Single-molecule, real-time (SMRT) sequencing reveals that DNA 6mA levels in starved cells are significantly reduced, especially symmetric 6mA, compared to those in vegetatively growing cells. Despite a global 6mA reduction, the fraction of asymmetric 6mA with a high methylation level was increased, which might be the driving force for stronger nucleosome positioning in starved cells. Starvation affects expression of many metabolism-related genes, the expression level change of which is associated with the amount of 6mA change, thereby linking 6mA with global transcription and starvation adaptation. The reduction of symmetric 6mA and the increase of asymmetric 6mA coincide with the downregulation of AMT1 and upregulation of AMT2 and AMT5, which are supposedly the MT-A70 methyltransferases required for symmetric and asymmetric 6mA, respectively. These results demonstrated that a regulated 6mA response to environmental cues is evolutionarily conserved in eukaryotes.

IMPORTANCE Increasing evidence indicated that 6mA could respond to environmental stressors in multicellular eukaryotes. As 6mA distribution and function differ significantly in multicellular and unicellular organisms, whether and how 6mA in unicellular eukaryotes responds to environmental stress remains elusive. In the present work, we characterized the dynamic changes of 6mA under starvation in the unicellular model organism Tetrahymena thermophila. Our results provide insights into how Tetrahymena fine-tunes its 6mA level and composition upon starvation, suggesting that a regulated 6mA response to environmental cues is evolutionarily conserved in eukaryotes.

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

DNA N⁶ -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 N⁶ -position of adenine within a specific DNA sequence to form N⁶ -allyladenine (6aA). The iodination of 6aA allyl group induced the formation of 1, N⁶ -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.

Toward incorporating epigenetics into regulation of gene expression in the parasite Trichomonas vaginalis

Trichomonas vaginalis is an extracellular parasite that colonizes the human urogenital tract, causing a highly prevalent sexually transmitted infection. The parasite must change its transcriptional profile in order to establish and maintain infection. However, few core regulatory elements and transcription factors have been identified to date and little is known about other mechanisms that may control these rapid changes in gene expression during parasite infection. In the last years, epigenetic mechanisms involved in the regulation of gene expression have been gaining major attention in this parasite. In this review, we summarize and discuss the major advances of the last few years with regard to epigenetics (DNA methylation, post-translational histone modifications, and histone variants) in the parasite T. vaginalis. These studies can shed light into our current understanding of this parasite's biology with far-reaching implications for the prognosis and treatment of trichomoniasis.

Detection of N6-Methyladenine in Eukaryotes

DNA N⁶-methyladenine (6mA) is a chemical modification at the N⁶-positon of adenine. In the last decades, 6mA had been found in genome from numerous prokaryotic species, but only existed in a few lower eukaryotes. In prokaryotes, 6mA plays an important role in restriction-modification, DNA replication, and DNA mismatch repair. Because of the too low abundance of 6mA, it was long-stalled whether 6mA existed in multicellular eukaryotes and playing any functions, particularly in mammals. In recent years, partially benefitting from the advances in analytical methods, 6mA was found in the genomes from Drosophila melanogaster, Chlamydomonas algae, Caenorhabditis elegans, zebrafish, Xenopus laevis and mouse embryonic stem cells and even in the human genome. The 6mA was dynamic changed in early embryonic development of fly and zebrafish and much more enriched in gene body of transposons in fly, repetitive regions in zebrafish, around the transcription start sites in Chlamydomonas, and widespread distribution in C. elegans, indicating 6mA probably playing different functions in different species. Meanwhile, 6mA methylases and demethylases were found in fly, worm, and Chlamydomonas. In this chapter, we will briefly review the distribution, regulation, and function of 6mA in eukaryotes and focus on the advances of 6mA analysis methods, especially LC-MS/MS, immunoprecipitation, next-generation sequencing, and single-molecule real-time sequencing technology.

Fluorescence assay based on the thioflavin T-induced conformation switch of G-quadruplexes for TET1 detection

A simple and label-free fluorescence method is developed for the highly sensitive detection of TET1 based on ThT/G-quadruplexes in combination with the specific design of oligonucleotides.

Advances in the profiling of N6-methyladenosine (m6A) modifications

Over 160 RNA modifications have been identified, including N⁷-methylguanine (m⁷G), N⁶-methyladenosine (m⁶A), and 5-methylcytosine (m⁵C). These modifications play key roles in regulating the fate of RNA. In eukaryotes, m⁶A is the most abundant mRNA modification, accounting for over 80% of all RNA methylation modifications. Highly dynamic m⁶A modification may exert important effects on organismal reproduction and development. Significant advances in understanding the mechanism of m⁶A modification have been made using immunoprecipitation, chemical labeling, and site-directed mutagenesis, combined with next-generation sequencing. Single-molecule real-time and nanopore direct RNA sequencing (DRS) approaches provide additional ways to study RNA modifications at the cellular level. In this review, we explore the technical history of identifying m⁶A RNA modifications, emphasizing technological advances in detecting m⁶A modification. In particular, we discuss the challenge of generating accurate dynamic single-base resolution m⁶A maps and also strategies for improving detection specificity. Finally, we outline a roadmap for future research in this area, focusing on the application of RNA epigenetic modification, represented by m⁶A modification.

Meta-i6mA: an interspecies predictor for identifying DNA N6-methyladenine sites of plant genomes by exploiting informative features in an integrative machine-learning framework

DNA N6-methyladenine (6mA) represents important epigenetic modifications, which are responsible for various cellular processes. The accurate identification of 6mA sites is one of the challenging tasks in genome analysis, which leads to an understanding of their biological functions. To date, several species-specific machine learning (ML)-based models have been proposed, but majority of them did not test their model to other species. Hence, their practical application to other plant species is quite limited. In this study, we explored 10 different feature encoding schemes, with the goal of capturing key characteristics around 6mA sites. We selected five feature encoding schemes based on physicochemical and position-specific information that possesses high discriminative capability. The resultant feature sets were inputted to six commonly used ML methods (random forest, support vector machine, extremely randomized tree, logistic regression, naïve Bayes and AdaBoost). The Rosaceae genome was employed to train the above classifiers, which generated 30 baseline models. To integrate their individual strength, Meta-i6mA was proposed that combined the baseline models using the meta-predictor approach. In extensive independent test, Meta-i6mA showed high Matthews correlation coefficient values of 0.918, 0.827 and 0.635 on Rosaceae, rice and Arabidopsis thaliana, respectively and outperformed the existing predictors. We anticipate that the Meta-i6mA can be applied across different plant species. Furthermore, we developed an online user-friendly web server, which is available at http://kurata14.bio.kyutech.ac.jp/Meta-i6mA/.

Mechanisms of mutagenesis induced by DNA lesions: multiple factors affect mutations in translesion DNA synthesis

Environmental mutagens lead to mutagenesis. However, the mechanisms are very complicated and not fully understood. Environmental mutagens produce various DNA lesions, including base-damaged or sugar-modified DNA lesions, as well as epigenetically modified DNA. DNA polymerases produce mutation spectra in translesion DNA synthesis (TLS) through misincorporation of incorrect nucleotides, frameshift deletions, blockage of DNA replication, imbalance of leading- and lagging-strand DNA synthesis, and genome instability. Motif or subunit in DNA polymerases further affects the mutations in TLS. Moreover, protein interactions and accessory proteins in DNA replisome also alter mutations in TLS, demonstrated by several representative DNA replisomes. Finally, in cells, multiple DNA polymerases or cellular proteins collaborate in TLS and reduce in vivo mutagenesis. Summaries and perspectives were listed. This review shows mechanisms of mutagenesis induced by DNA lesions and the effects of multiple factors on mutations in TLS in vitro and in vivo.

SMRT Sequencing of Paramecium Bursaria Chlorella Virus-1 Reveals Diverse Methylation Stability in Adenines Targeted by Restriction Modification Systems

Chloroviruses (family Phycodnaviridae) infect eukaryotic, freshwater, unicellular green algae. A unique feature of these viruses is an abundance of DNA methyltransferases, with isolates dedicating up to 4.5% of their protein coding potential to these genes. This diversity highlights just one of the long-standing values of the chlorovirus model system; where group-wide epigenomic characterization might begin to elucidate the function(s) of DNA methylation in large dsDNA viruses. We characterized DNA modifications in the prototype chlorovirus, PBCV-1, using single-molecule real time (SMRT) sequencing (aka PacBio). Results were compared to total available sites predicted in silico based on DNA sequence alone. SMRT-software detected N6-methyl-adenine (m6A) at GATC and CATG recognition sites, motifs previously shown to be targeted by PBCV-1 DNA methyltransferases M.CviAI and M. CviAII, respectively. At the same time, PacBio analyses indicated that 10.9% of the PBCV-1 genome had large interpulse duration ratio (ipdRatio) values, the primary metric for DNA modification identification. These events represent 20.6x more sites than can be accounted for by all available adenines in GATC and CATG motifs, suggesting base or backbone modifications other than methylation might be present. To define methylation stability, we cross-compared methylation status of each GATC and CATG sequence in three biological replicates and found ∼81% of sites were stably methylated, while ∼2% consistently lack methylation. The remaining 17% of sites were stochastically methylated. When methylation status was analyzed for both strands of each target, we show that palindromes existed in completely non-methylated states, fully-methylated states, or hemi-methylated states, though GATC sites more often lack methylation than CATG sequences. Given that both sequences are targeted by not just methyltransferases, but by restriction endonucleases that are together encoded by PBCV-1 as virus-originating restriction modification (RM) systems, there is strong selective pressure to modify all target sites. The finding that most instances of non-methylation are associated with hemi-methylation is congruent with observations that hemi-methylated palindromes are resistant to cleavage by restriction endonucleases. However, sites where hemi-methylation is conserved might represent a unique regulatory function for PBCV-1. This study serves as a baseline for future investigation into the epigenomics of chloroviruses and their giant virus relatives.

iDNA-MS: An Integrated Computational Tool for Detecting DNA Modification Sites in Multiple Genomes

5hmC, 6mA, and 4mC are three common DNA modifications and are involved in various of biological processes. Accurate genome-wide identification of these sites is invaluable for better understanding their biological functions. Owing to the labor-intensive and expensive nature of experimental methods, it is urgent to develop computational methods for the genome-wide detection of these sites. Keeping this in mind, the current study was devoted to construct a computational method to identify 5hmC, 6mA, and 4mC. We initially used K-tuple nucleotide component, nucleotide chemical property and nucleotide frequency, and mono-nucleotide binary encoding scheme to formulate samples. Subsequently, random forest was utilized to identify 5hmC, 6mA, and 4mC sites. Cross-validated results showed that the proposed method could produce the excellent generalization ability in the identification of the three modification sites. Based on the proposed model, a web-server called iDNA-MS was established and is freely accessible at http://lin-group.cn/server/iDNA-MS.

MASQC: Next Generation Sequencing Assists Third Generation Sequencing for Quality Control in N6-Methyladenine DNA Identification

DNA N6-methyladenine (6mA) modification has been discovered as the most prevalent DNA modification in prokaryotes and eukaryotes, involving gene expression, DNA replication and repair, and host-pathogen interactions. Single-molecule real-time sequencing (SMRT-seq) can detect 6mA events in prokaryotic and eukaryotic genomes at the single-nucleotide level. However, there are no strict and economical quality control methods for high false-positive 6mA events in eukaryotic genomes. Therefore, by analyzing the distribution of 6mA in eukaryotic and prokaryotes, we proposed a method named MASQC (MeDIP-seq assists SMRT-seq for quality control in 6mA identification), which can identify 6mA events without doing the whole genome amplification (WGA) sequencing. The proposed MASQC method was assessed on two eukaryotic genomes and six bacterial genomes, our results demonstrate that MASQC performs well in quality control of false positive 6mA identification for both eukaryotic and prokaryotic genomes.

Mapping and editing of nucleic acid modifications

Modification on nucleic acid plays a pivotal role in controlling gene expression. Various kinds of modifications greatly increase the information-encoding capacity of DNA and RNA by introducing extra chemical group to existing bases instead of altering the genetic sequences. As a marker on DNA or RNA, nucleic acid modification can be recognized by specific proteins, leading to versatile regulation of gene expression. However, modified and regular bases are often indistinguishable by most conventional molecular methods, impeding detailed functional studies that require the information of genomic location. Recently, new technologies are emerging to resolve the positions of varied modifications on both DNA and RNA. Intriguingly, by integrating regional targeting tools and effector proteins, researchers begin to actively control the modification status of desired gene in vivo. In this review, we summarize the characteristics of DNA and RNA modifications, the available mapping and editing tools, and the potential application as well as deficiency of these technologies in basic and translational researches.

Discovery of a new predominant cytosine DNA modification that is linked to gene expression in malaria parasites

DNA cytosine modifications are key epigenetic regulators of cellular processes in mammalian cells, with their misregulation leading to varied disease states. In the human malaria parasite Plasmodium falciparum, a unicellular eukaryotic pathogen, little is known about the predominant cytosine modifications, cytosine methylation (5mC) and hydroxymethylation (5hmC). Here, we report the first identification of a hydroxymethylcytosine-like (5hmC-like) modification in P. falciparum asexual blood stages using a suite of biochemical methods. In contrast to mammalian cells, we report 5hmC-like levels in the P. falciparum genome of 0.2–0.4%, which are significantly higher than the methylated cytosine (mC) levels of 0.01–0.05%. Immunoprecipitation of hydroxymethylated DNA followed by next generation sequencing (hmeDIP-seq) revealed that 5hmC-like modifications are enriched in gene bodies with minimal dynamic changes during asexual development. Moreover, levels of the 5hmC-like base in gene bodies positively correlated to transcript levels, with more than 2000 genes stably marked with this modification throughout asexual development. Our work highlights the existence of a new predominant cytosine DNA modification pathway in P. falciparum and opens up exciting avenues for gene regulation research and the development of antimalarials.

Methylomic Changes of Autophagy-Related Genes by Legionella Effector Lpg2936 in Infected Macrophages

Legionella pneumophila (L. pneumophila) is a Gram-negative bacterium that infects the human respiratory tract causing Legionnaires' disease, a severe form of pneumonia. Recently, rising evidence indicated the ability of Legionella to regulate host defense via its type 4 secretion system including hundreds of effectors that promote intracellular bacterial replication. The host defense against such invaders includes autophagic machinery that is responsible for degradation events of invading pathogens and recycling of cell components. The interplay between host autophagy and Legionella infection has been reported, indicating the role of bacterial effectors in the regulation of autophagy during intracellular replication. Here, we investigated the potential impact of Legionella effector Lpg2936 in the regulation of host autophagy and its role in bacterial replication using mice-derived macrophages and human lung epithelial cells (A549 cells). First, monitoring of autophagic flux following infection revealed a marked reduction of Atg7 and LC3B expression profile and low accumulation levels of autophagy-related LC3-I, LC3-II, and the Atg12-Atg5 protein complex. A novel methyladenine alteration was observed due to irreversible changes of GATC motif to G(6 mA) TC in the promoter region of Atg7 and LC3B indicated by cleaved genomic-DNA using the N6 methyladenine-sensitive restriction enzyme DpnI. Interestingly, RNA interference (RNAi) of Lpg2936 in infected macrophages showed dramatic inhibition of bacterial replication by restoring the expression of autophagy-related proteins. This is accompanied by low production levels of bacterial-associated pro-inflammatory cytokines. Furthermore, a constructed Lpg2936 segment in the GFP expression vector was translocated in the host nucleus and successfully induced methyladenine changes in Atg7 and LC3B promoter region and subsequently regulated autophagy in A549 cells independent of infection. Finally, treatment with methylation inhibitors 5-AZA and (2)-Epigallocatechin-3-gallate (EGCG) was able to restore autophagy-related gene expression and to disrupt bacterial replication in infected macrophages. This cumulative evidence indicates the methylation effect of Legionella effector Lpg2936 on the host autophagy-related molecules Atg7 and LC3B and subsequent reduction in the expression levels of autophagy effectors during intracellular replication of L. pneumophila.

SDM6A: A Web-Based Integrative Machine-Learning Framework for Predicting 6mA Sites in the Rice Genome

DNA N6-adenine methylation (6mA) is an epigenetic modification in prokaryotes and eukaryotes. Identifying 6mA sites in rice genome is important in rice epigenetics and breeding, but non-random distribution and biological functions of these sites remain unclear. Several machine-learning tools can identify 6mA sites but show limited prediction accuracy, which limits their usability in epigenetic research. Here, we developed a novel computational predictor, called the Sequence-based DNA N6-methyladenine predictor (SDM6A), which is a two-layer ensemble approach for identifying 6mA sites in the rice genome. Unlike existing methods, which are based on single models with basic features, SDM6A explores various features, and five encoding methods were identified as appropriate for this problem. Subsequently, an optimal feature set was identified from encodings, and corresponding models were developed individually using support vector machine and extremely randomized tree. First, all five single models were integrated via ensemble approach to define the class for each classifier. Second, two classifiers were integrated to generate a final prediction. SDM6A achieved robust performance on cross-validation and independent evaluation, with average accuracy and Matthews correlation coefficient (MCC) of 88.2% and 0.764, respectively. Corresponding metrics were 4.7%-11.0% and 2.3%-5.5% higher than those of existing methods, respectively. A user-friendly, publicly accessible web server (http://thegleelab.org/SDM6A) was implemented to predict novel putative 6mA sites in rice genome.