Identification and analysis of adenine N6-methylation sites in the rice genome

Semiconservative transmission of DNA N 6-adenine methylation in a unicellular eukaryote

Although DNA N 6-adenine methylation (6mA) is best known in prokaryotes, its presence in eukaryotes has recently generated great interest. Biochemical and genetic evidence supports that AMT1, an MT-A70 family methyltransferase (MTase), is crucial for 6mA deposition in unicellular eukaryotes. Nonetheless, the 6mA transmission mechanism remains to be elucidated. Taking advantage of single-molecule real-time circular consensus sequencing (SMRT CCS), here we provide definitive evidence for semiconservative transmission of 6mA in Tetrahymena thermophila In wild-type (WT) cells, 6mA occurs at the self-complementary ApT dinucleotide, mostly in full methylation (full-6mApT); after DNA replication, hemi-methylation (hemi-6mApT) is transiently present on the parental strand, opposite to the daughter strand readily labeled by 5-bromo-2'-deoxyuridine (BrdU). In ΔAMT1 cells, 6mA predominantly occurs as hemi-6mApT. Hemi-to-full conversion in WT cells is fast, robust, and processive, whereas de novo methylation in ΔAMT1 cells is slow and sporadic. In Tetrahymena, regularly spaced 6mA clusters coincide with the linker DNA of nucleosomes arrayed in the gene body. Importantly, in vitro methylation of human chromatin by the reconstituted AMT1 complex recapitulates preferential targeting of hemi-6mApT sites in linker DNA, supporting AMT1's intrinsic and autonomous role in maintenance methylation. We conclude that 6mA is transmitted by a semiconservative mechanism: full-6mApT is split by DNA replication into hemi-6mApT, which is restored to full-6mApT by AMT1-dependent maintenance methylation. Our study dissects AMT1-dependent maintenance methylation and AMT1-independent de novo methylation, reveals a 6mA transmission pathway with a striking similarity to 5-methylcytosine (5mC) transmission at the CpG dinucleotide, and establishes 6mA as a bona fide eukaryotic epigenetic mark.

CAPTURE: Comprehensive anti-cancer peptide predictor with a unique amino acid sequence encoder

Anticancer peptides (ACPs) key properties including bioactivity, high efficacy, low toxicity, and lack of drug resistance make them ideal candidates for cancer therapies. To deeply explore the potential of ACPs and accelerate development of cancer therapies, although 53 Artificial Intelligence supported computational predictors have been developed for ACPs and non ACPs classification but only one predictor has been developed for ACPs functional types annotations. Moreover, these predictors extract amino acids distribution patterns to transform peptides sequences into statistical vectors that are further fed to classifiers for discriminating peptides sequences and annotating peptides functional classes. Overall, these predictors remain fail in extracting diverse types of amino acids distribution patterns from peptide sequences. The paper in hand presents a unique CARE encoder that transforms peptides sequences into statistical vectors by extracting 4 different types of distribution patterns including correlation, distribution, composition, and transition. Across public benchmark dataset, proposed encoder potential is explored under two different evaluation settings namely; intrinsic and extrinsic. Extrinsic evaluation indicates that 12 different machine learning classifiers achieve superior performance with the proposed encoder as compared to 55 existing encoders. Furthermore, an intrinsic evaluation reveals that, unlike existing encoders, the proposed encoder generates more discriminative clusters for ACPs and non-ACPs classes. Across 8 public benchmark ACPs and non-ACPs classification datasets, proposed encoder and Adaboost classifier based CAPTURE predictor outperforms existing predictors with an average accuracy, recall and MCC score of 1%, 4%, and 2% respectively. In generalizeability evaluation case study, across 7 benchmark anti-microbial peptides classification datasets, CAPTURE surpasses existing predictors by an average AU-ROC of 2%. CAPTURE predictive pipeline along with label powerset method outperforms state-of-the-art ACPs functional types predictor by 5%, 5%, 5%, 6%, and 3% in terms of average accuracy, subset accuracy, precision, recall, and F1 respectively. CAPTURE web application is available at https://sds_genetic_analysis.opendfki.de/CAPTURE.

Identification of AtALKBH1A and AtALKBH1D as DNA N6-adenine demethylases in Arabidopsis thaliana

DNA N6-methyladenine (6 mA) has recently been discovered as a novel DNA modification in animals and plants. In mammals, AlkB homolog 1 (ALKBH1) has been identified as a DNA 6 mA demethylase. ALKBH1 tightly controls the DNA 6 mA methylation level of mammalian genomes and plays important role in regulating gene expression. DNA 6 mA methylation has also been reported to exist in plant genomes, however, the plant DNA 6 mA demethylases and their function remain largely unknown. Here we identify homologs of ALKBH1 as DNA 6 mA demethylases in Arabidopsis. We discover that there are four homologs of ALKBH1, AtALKBH1A, AtALKBH1B, AtALKBH1C and AtALKBH1D, in Arabidopsis. In vitro enzymatic activity studies reveal that AtALKBH1A and 1D can efficiently erase DNA 6 mA methylation. Loss of function of AtALKBH1A and AtALKBH1D causes elevated DNA 6 mA methylation levels in vivo. atalkbh1a/1d mutant displays delayed seed gemination. Based on our RNA-seq data, we find some regulators of seed gemination are dysregulated in atalkbh1a/1d, and the dysregulation is correlated with changes of DNA 6 mA methylation levels. This study identifies plant DNA 6 mA demethylases and reports the function of DNA 6 mA methylation in regulating seed germination.

Molecular basis of an atypical dsDNA 5mC/6mA bifunctional dioxygenase CcTet from Coprinopsis cinerea in catalyzing dsDNA 5mC demethylation

The eukaryotic epigenetic modifications 5-methyldeoxycytosine (5mC) and N6-methyldeoxyadenine (6mA) have indispensable regulatory roles in gene expression and embryonic development. We recently identified an atypical bifunctional dioxygenase CcTet from Coprinopsis cinerea that works on both 5mC and 6mA demethylation. The nonconserved residues Gly331 and Asp337 of CcTet facilitate 6mA accommodation, while D337F unexpectedly abolishes 5mC oxidation activity without interfering 6mA demethylation, indicating a prominent distinct but unclear 5mC oxidation mechanism to the conventional Tet enzymes. Here, we assessed the molecular mechanism of CcTet in catalyzing 5mC oxidation by representing the crystal structure of CcTet-5mC-dsDNA complex. We identified the distinct mechanism by which CcTet recognizes 5mC-dsDNA compared to 6mA-dsDNA substrate. Moreover, Asp337 was found to have a central role in compensating for the loss of a critical 5mC-stablizing H-bond observed in conventional Tet enzymes, and stabilizes 5mC and subsequent intermediates through an H-bond with the N4 atom of the substrates. These findings improve our understanding of Tet enzyme functions in the dsDNA 5mC and 6mA demethylation pathways, and provide useful information for future discovery of small molecular probes targeting Tet enzymes in DNA active demethylation processes.

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.

Identification and quantification of N6-methyladenosine by chemical derivatization coupled with 19F NMR spectroscopy

N 6-Methyladenosine (6mA) is a well-known prokaryotic DNA modification that has been shown to play epigenetic roles in eukaryotic DNA. Accurate detection and quantification of 6mA are prerequisites for molecular understanding of the impact of 6mA modification on DNA. However, the existing methods have several problems, such as high false-positive rate, time-consuming and complex operating procedures. Chemical sensors for the selective detection of 6mA modification are rarely reported in the literature. Fluorinated phenylboronic acid combined with 19F NMR analysis is an effective method for determining DNA or RNA modification. In this study, we presented a simple and fast chemical method for labelling the 6th imino group of 6mA using a boric-acid-derived probe. Besides, the trifluoromethyl group of trifluoromethyl phenylboronic acid (2a) could detect 6mA modification through 19F NMR. Combined with this sensor system, 6mA modification could be detected well and quickly in 6 types of deoxynucleoside mixtures and DNA samples. Taken together, the method developed in the current study has potential for specific detection of 6mA in biological samples.

Genome-wide adenine N6-methylation map reveals epigenomic regulation of lipid accumulation in Nannochloropsis

Epigenetic marks on histones and DNA, such as DNA methylation at N6-adenine (6mA), play crucial roles in gene expression and genome maintenance, but their deposition and function in microalgae remain largely uncharacterized. Here, we report a genome-wide 6mA map for the model industrial oleaginous microalga Nannochloropsis oceanica produced by single-molecule real-time sequencing. Found in 0.1% of adenines, 6mA sites are mostly enriched at the AGGYV motif, more abundant in transposons and 3' untranslated regions, and associated with active transcription. Moreover, 6mA gradually increases in abundance along the direction of gene transcription and shows special positional enrichment near splicing donor and transcription termination sites. Highly expressed genes tend to show greater 6mA abundance in the gene body than do poorly expressed genes, indicating a positive interaction between 6mA and general transcription factors. Furthermore, knockout of the putative 6mA methylase NO08G00280 by genome editing leads to changes in methylation patterns that are correlated with changes in the expression of molybdenum cofactor, sulfate transporter, glycosyl transferase, and lipase genes that underlie reductions in biomass and oil productivity. By contrast, knockout of the candidate demethylase NO06G02500 results in increased 6mA levels and reduced growth. Unraveling the epigenomic players and their roles in biomass productivity and lipid metabolism lays a foundation for epigenetic engineering of industrial microalgae.

SoftVoting6mA: An improved ensemble-based method for predicting DNA N6-methyladenine sites in cross-species genomes

The DNA N6-methyladenine (6mA) is an epigenetic modification, which plays a pivotal role in biological processes encompassing gene expression, DNA replication, repair, and recombination. Therefore, the precise identification of 6mA sites is fundamental for better understanding its function, but challenging. We proposed an improved ensemble-based method for predicting DNA N6-methyladenine sites in cross-species genomes called SoftVoting6mA. The SoftVoting6mA selected four (electron-ion-interaction pseudo potential, One-hot encoding, Kmer, and pseudo dinucleotide composition) codes from 15 types of encoding to represent DNA sequences by comparing their performances. Similarly, the SoftVoting6mA combined four learning algorithms using the soft voting strategy. The 5-fold cross-validation and the independent tests showed that SoftVoting6mA reached the state-of-the-art performance. To enhance accessibility, a user-friendly web server is provided at http://www.biolscience.cn/SoftVoting6mA/.

Significant differences in the degree of genomic DNA N6-methyladenine modifications in Acidithiobacillus ferrooxidans with two different culture substrates

DNA N6-methyladenine (6mA) modification is widespread in organisms and plays an important functional role in the regulation of cellular processes. As a model organism in biohydrometallurgy, Acidithiobacillus ferrooxidans can obtain energy from the oxidation of ferrous iron (Fe2+) and various reduced inorganic sulfides (RISCs) under acidic conditions. To determine the linkage between genomic DNA methylation and the switching between the two oxidative metabolic pathways in A. ferrooxidans, the 6mA landscape in the genome of A. ferrooxidans cultured under different conditions was evaluated by using 6mA-IP-seq. A total of 214 and 47 high-confidence peaks of 6mA were identified under the Fe2+ and RISCs oxidizing conditions, respectively (P<10-5), suggesting that genomic methylation was greater under Fe2+ oxidizing conditions. 6mA experienced a decline at the transcription start site (TSS) and occurs frequently in gene bodies under both oxidizing conditions. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that 7 KEGG pathways were mapped into and most of the differentially methylated genes were enriched in oxidative phosphorylation and metabolic pathways. Fourteen genes were selected for studying the effect of differences in methylation on mRNA expression. Thirteen genes, excluding petA-1, demonstrated a decrease in mRNA expression as methylation levels increased. Overall, the 6mA methylation enrichment patterns are similar under two conditions but show differences in the enriched pathways. The phenomenon of upregulated gene methylation levels coupled with downregulated expression suggests a potential association between the regulation mechanisms of 6mA and the Fe2+ and RISCs oxidation pathways.

Exploring N6-methyladenosine (m6A) modification in tree species: opportunities and challenges

N 6-methyladenosine (m6A) in eukaryotes is the most common and widespread internal modification in mRNA. The modification regulates mRNA stability, translation efficiency, and splicing, thereby fine-tuning gene regulation. In plants, m6A is dynamic and critical for various growth stages, embryonic development, morphogenesis, flowering, stress response, crop yield, and biomass. Although recent high-throughput sequencing approaches have enabled the rapid identification of m6A modification sites, the site-specific mechanism of this modification remains unclear in trees. In this review, we discuss the functional significance of m6A in trees under different stress conditions and discuss recent advancements in the quantification of m6A. Quantitative and functional insights into the dynamic aspect of m6A modification could assist researchers in engineering tree crops for better productivity and resistance to various stress conditions.

The Epigenomic Features and Potential Functions of PEG- and PDS-Favorable DNA G-Quadruplexes in Rice

A G-quadruplex (G4) is a typical non-B DNA structure and involved in various DNA-templated events in eukaryotic genomes. PEG and PDS chemicals have been widely applied for promoting the folding of in vivo or in vitro G4s. However, how PEG and PDS preferentially affect a subset of G4 formation genome-wide is still largely unknown. We here conducted a BG4-based IP-seq in vitro under K++PEG or K++PDS conditions in the rice genome. We found that PEG-favored IP-G4s+ have distinct sequence features, distinct genomic distributions and distinct associations with TEGs, non-TEGs and subtypes of TEs compared to PDS-favored ones. Strikingly, PEG-specific IP-G4s+ are associated with euchromatin with less enrichment levels of DNA methylation but with more enriched active histone marks, while PDS-specific IP-G4s+ are associated with heterochromatin with higher enrichment levels of DNA methylation and repressive marks. Moreover, we found that genes with PEG-specific IP-G4s+ are more expressed than those with PDS-specific IP-G4s+, suggesting that PEG/PDS-specific IP-G4s+ alone or coordinating with epigenetic marks are involved in the regulation of the differential expression of related genes, therefore functioning in distinct biological processes. Thus, our study provides new insights into differential impacts of PEG and PDS on G4 formation, thereby advancing our understanding of G4 biology.

METTL4-mediated N6-methyladenine DNA modification regulates thermotolerance in Arabidopsis thaliana

DNA N6-methyladenine (6 mA) is an evolutionarily conserved DNA modification in procaryotes and eukaryotes. The DNA 6 mA methylation is tightly controlled by 6 mA regulatory proteins. DNA N6-adenine methyltransferase 1 (DAMT-1) has been identified as a DNA 6 mA methyltransferase in animals. In plants, DNA 6 mA methylation has been found, however, the DNA 6 mA methyltransferases and their function in plants are largely unknown. In our study, we find METTL4 is a DNA 6 mA methyltransferase in Arabidopsis thaliana. Both in vitro and in vivo evidences support the DNA 6 mA methyltransferase activity of METTL4. mettl4 mutant is hypersensitive to heat stress, suggesting DNA 6 mA methylation plays important role in heat stress adaption. RNA-seq and 6 mA IP-qPCR analysis show that METTL4 participates in heat stress tolerance by regulating expression of heat responsive genes. Our study find METTL4 is a plant DNA 6 mA methyltransferase and illustrates its function in regulating heat stress response.

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.

MethSemble-6mA: an ensemble-based 6mA prediction server and its application on promoter region of LBD gene family in Poaceae

The Lateral Organ Boundaries Domain (LBD) containing genes are a set of plant-specific transcription factors and are crucial for controlling both organ development and defense mechanisms as well as anthocyanin synthesis and nitrogen metabolism. It is imperative to understand how methylation regulates gene expression, through predicting methylation sites of their promoters particularly in major crop species. In this study, we developed a user-friendly prediction server for accurate prediction of 6mA sites by incorporating a robust feature set, viz., Binary Encoding of Mono-nucleotide DNA. Our model,MethSemble-6mA, outperformed other state-of-the-art tools in terms of accuracy (93.12%). Furthermore, we investigated the pattern of probable 6mA sites at the upstream promoter regions of the LBD-containing genes in Triticum aestivum and its allied species using the developed tool. On average, each selected species had four 6mA sites, and it was found that with speciation and due course of evolution in wheat, the frequency of methylation have reduced, and a few sites remain conserved. This obviously cues gene birth and gene expression alteration through methylation over time in a species and reflects functional conservation throughout evolution. Since DNA methylation is a vital event in almost all plant developmental processes (e.g., genomic imprinting and gametogenesis) along with other life processes, our findings on epigenetic regulation of LBD-containing genes have dynamic implications in basic and applied research. Additionally, MethSemble-6mA (http://cabgrid.res.in:5799/) will serve as a useful resource for a plant breeders who are interested to pursue epigenetic-based crop improvement research.

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.

A DNA adenine demethylase impairs PRC2-mediated repression of genes marked by a specific chromatin signature

BACKGROUND

The Fe (II)- and α-ketoglutarate-dependent AlkB family dioxygenases are implicated in nucleotide demethylation. AlkB homolog1 (ALKBH1) is shown to demethylate DNA adenine methylation (6mA) preferentially from single-stranded or unpaired DNA, while its demethylase activity and function in the chromatin context are unclear.

RESULTS

Here, we find that loss-of-function of the rice ALKBH1 gene leads to increased 6mA in the R-loop regions of the genome but has a limited effect on the overall 6mA level. However, in the context of mixed tissues, rather than on individual loci, the ALKBH1 mutation or overexpression mainly affects the expression of genes with a specific combination of chromatin modifications in the body region marked with H3K4me3 and H3K27me3 but depleted of DNA CG methylation. In the similar context of mixed tissues, further analysis reveals that the ALKBH1 protein preferentially binds to genes marked by the chromatin signature and has a function to maintain a high H3K4me3/H3K27me3 ratio by impairing the binding of Polycomb repressive complex 2 (PRC2) to the targets, which is required for both the basal and stress-induced expression of the genes.

CONCLUSION

Our findings unravel a function of ALKBH1 to control the balance between the antagonistic histone methylations for gene activity and provide insight into the regulatory mechanism of PRC2-mediated H3K27me3 deposition within the gene body region.

Aberrant DNA N6 -methyladenine incorporation via adenylate kinase 1 is suppressed by ADAL deaminase-dependent 2'-deoxynucleotide pool sanitation

Intracellular decay of N6 -methyladenine (m6A)-containing RNA potentially induces aberrant N6 -methyl-2'-adenine (6mdA) misincorporation into DNA. Biophysically, misincorporated 6mdA may destabilize the DNA duplex in a manner similar to bona fide methylated 6mdA DNA, thereby affecting DNA replication and transcription. Utilizing heavy stable isotope labeling and ultrasensitive UHPLC-MS/MS assay, we demonstrate that intracellular m6A-RNA decay does not generate free 6mdA species, nor lead to any misincorporated DNA 6mdA in most mammalian cell lines tested, unveiling the existence of a sanitation mechanism that prevents 6mdA misincorporation. Depletion of deaminase ADAL increases the levels of free 6mdA species, concomitant with the presence of DNA-misincorporated 6mdA resulting from intracellular RNA m6A decay, suggesting that ADAL catabolizes 6mdAMP in vivo. Furthermore, we show that the overexpression of adenylate kinase 1 (AK1) promotes 6mdA misincorporation, while AK1 knockdown diminishes 6mdA incorporation, in ADAL-deficient cells. We conclude that ADAL together with other factors (such as MTH1) contributes to 2'-deoxynucleotide pool sanitation in most cells but compromised sanitation (e.g., in NIH3T3 cells) and increased AK1 expression may facilitate aberrant 6mdA incorporation. This sanitation mechanism may provide a framework for the maintenance of the epigenetic 6mdA landscape.

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.

N6-methyladenosine RNA modification regulates cotton drought response in a Ca2+ and ABA-dependent manner

N⁶ -methyladenosine (m⁶ A) is the most prevalent internal modification present in mRNAs, and is considered to participate in a range of developmental and biological processes. Drought response is highly regulated at the genomic, transcriptional and post-transcriptional levels. However, the biological function and regulatory mechanism of m⁶ A modification in the drought stress response is still poorly understood. We generated a transcriptome-wide m⁶ A map using drought-resistant and drought-sensitive varieties of cotton under different water deficient conditions to uncover patterns of m⁶ A methylation in cotton response to drought stress. The results reveal that m⁶ A represents a common modification and exhibit dramatic changes in distribution during drought stress. More 5'UTR m⁶ A was deposited in the drought-resistant variety and was associated with a positive effect on drought resistance by regulating mRNA abundance. Interestingly, we observed that increased m⁶ A abundance was associated with increased mRNA abundance under drought, contributing to drought resistance, and vice versa. The demethylase GhALKBH10B was found to decrease m⁶ A levels, facilitating the mRNA decay of ABA signal-related genes (GhZEP, GhNCED4 and GhPP2CA) and Ca²⁺ signal-related genes (GhECA1, GhCNGC4, GhANN1 and GhCML13), and mutation of GhALKBH10B enhanced drought resistance at seedling stage in cotton. Virus-induced gene silencing (VIGS) of two Ca²⁺ -related genes, GhECA1 and GhCNGC4, reduced drought resistance with the decreased m⁶ A enrichment on silenced genes in cotton. Collectively, we reveal a novel mechanism of post-transcriptional modification involved in affecting drought response in cotton, by mediating m⁶ A methylation on targeted transcripts in the ABA and Ca²⁺ signalling transduction pathways.

Identification and quantification of DNA N6-methyladenine modification in mammals: A challenge to modern analytical technologies

DNA N⁶-methyladenine modification (6mA) is a predominant epigenetic mark in prokaryotes but rarely present in multicellular metazoa. The analytical technologies have been developed for sensitive detection of 6mA, including ultra-high performance liquid chromatography coupled with mass spectrometry (UHPLC-MS/MS) and single molecule real-time sequencing (SMRTseq). However, it remains challenging to detect 6mA at global level and/or in the context of sequence in multicellular metazoa (including mammals). This mini-review brings insights into current dilemma and potential solutions for the identification and quantifications of 6mA in mammals.

ProFAB-open protein functional annotation benchmark

As the number of protein sequences increases in biological databases, computational methods are required to provide accurate functional annotation with high coverage. Although several machine learning methods have been proposed for this purpose, there are still two main issues: (i) construction of reliable positive and negative training and validation datasets, and (ii) fair evaluation of their performances based on predefined experimental settings. To address these issues, we have developed ProFAB: Open Protein Functional Annotation Benchmark, which is a platform providing an infrastructure for a fair comparison of protein function prediction methods. ProFAB provides filtered and preprocessed protein annotation datasets and enables the training and evaluation of function prediction methods via several options. We believe that ProFAB will be useful for both computational and experimental researchers by enabling the utilization of ready-to-use datasets and machine learning algorithms for protein function prediction based on Gene Ontology terms and Enzyme Commission numbers. ProFAB is available at https://github.com/kansil/ProFAB and https://profab.kansil.org.

Beyond transcription factors: more regulatory layers affecting soybean gene expression under abiotic stress

Abiotic stresses such as nutritional imbalance, salt, light intensity, and high and low temperatures negatively affect plant growth and development. Through the course of evolution, plants developed multiple mechanisms to cope with environmental variations, such as physiological, morphological, and molecular adaptations. Epigenetic regulation, transcription factor activity, and post-transcriptional regulation operated by RNA molecules are mechanisms associated with gene expression regulation under stress. Epigenetic regulation, including histone and DNA covalent modifications, triggers chromatin remodeling and changes the accessibility of transcription machinery leading to alterations in gene activity and plant homeostasis responses. Soybean is a legume widely produced and whose productivity is deeply affected by abiotic stresses. Many studies explored how soybean faces stress to identify key elements and improve productivity through breeding and genetic engineering. This review summarizes recent progress in soybean gene expression regulation through epigenetic modifications and circRNAs pathways, and points out the knowledge gaps that are important to study by the scientific community. It focuses on epigenetic factors participating in soybean abiotic stress responses, and chromatin modifications in response to stressful environments and draws attention to the regulatory potential of circular RNA in post-transcriptional processing.

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.

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.

Plant synthetic epigenomic engineering for crop improvement

Efforts have been directed to redesign crops with increased yield, stress adaptability, and nutritional value through synthetic biology-the application of engineering principles to biology. A recent expansion in our understanding of how epigenetic mechanisms regulate plant development and stress responses has unveiled a new set of resources that can be harnessed to develop improved crops, thus heralding the promise of "synthetic epigenetics." In this review, we summarize the latest advances in epigenetic regulation and highlight how innovative sequencing techniques, epigenetic editing, and deep learning-driven predictive tools can rapidly extend these insights. We also proposed the future directions of synthetic epigenetics for the development of engineered smart crops that can actively monitor and respond to internal and external cues throughout their life cycles.

i6mA-Caps: a CapsuleNet-based framework for identifying DNA N6-methyladenine sites

MOTIVATION

DNA N6-methyladenine (6mA) has been demonstrated to have an essential function in epigenetic modification in eukaryotic species in recent research. 6mA has been linked to various biological processes. It's critical to create a new algorithm that can rapidly and reliably detect 6mA sites in genomes to investigate their biological roles. The identification of 6mA marks in the genome is the first and most important step in understanding the underlying molecular processes, as well as their regulatory functions.

RESULTS

In this article, we proposed a novel computational tool called i6mA-Caps which CapsuleNet based a framework for identifying the DNA N6-methyladenine sites. The proposed framework uses a single encoding scheme for numerical representation of the DNA sequence. The numerical data is then used by the set of convolution layers to extract low-level features. These features are then used by the capsule network to extract intermediate-level and later high-level features to classify the 6mA sites. The proposed network is evaluated on three datasets belonging to three genomes which are Rosaceae, Rice and Arabidopsis thaliana. Proposed method has attained an accuracy of 96.71%, 94% and 86.83% for independent Rosaceae dataset, Rice dataset and A.thaliana dataset respectively. The proposed framework has exhibited improved results when compared with the existing top-of-the-line methods.

AVAILABILITY AND IMPLEMENTATION

A user-friendly web-server is made available for the biological experts which can be accessed at: http://nsclbio.jbnu.ac.kr/tools/i6mA-Caps/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Epigenomic Features and Potential Functions of K+ and Na+ Favorable DNA G-Quadruplexes in Rice

DNA G-quadruplexes (G4s) are non-canonical four-stranded DNA structures involved in various biological processes in eukaryotes. Molecularly crowded solutions and monovalent cations have been reported to stabilize in vitro and in vivo G4 formation. However, how K⁺ and Na⁺ affect G4 formation genome-wide is still unclear in plants. Here, we conducted BG4-DNA-IP-seq, DNA immunoprecipitation with anti-BG4 antibody coupled with sequencing, under K⁺ and Na⁺ + PEG conditions in vitro. We found that K⁺-specific IP-G4s had a longer peak size, more GC and PQS content, and distinct AT and GC skews compared to Na⁺-specific IP-G4s. Moreover, K⁺- and Na⁺-specific IP-G4s exhibited differential subgenomic enrichment and distinct putative functional motifs for the binding of certain trans-factors. More importantly, we found that K⁺-specific IP-G4s were more associated with active marks, such as active histone marks, and low DNA methylation levels, as compared to Na⁺-specific IP-G4s; thus, K⁺-specific IP-G4s in combination with active chromatin features facilitate the expression of overlapping genes. In addition, K⁺- and Na⁺-specific IP-G4 overlapping genes exhibited differential GO (gene ontology) terms, suggesting they may have distinct biological relevance in rice. Thus, our study, for the first time, explores the effects of K⁺ and Na⁺ on global G4 formation in vitro, thereby providing valuable resources for functional G4 studies in rice. It will provide certain G4 loci for the biotechnological engineering of rice in the future.

Activating stress memory: eustressors as potential tools for plant breeding

Plants are continuously exposed to stress conditions, such that they have developed sophisticated and elegant survival strategies, which are reflected in their phenotypic plasticity, priming capacity, and memory acquisition. Epigenetic mechanisms play a critical role in modulating gene expression and stress responses, allowing malleability, reversibility, stability, and heritability of favourable phenotypes to enhance plant performance. Considering the urgency to improve our agricultural system because of going impacting climate change, potential and sustainable strategies rely on the controlled use of eustressors, enhancing desired characteristics and yield and shaping stress tolerance in crops. However, for plant breeding purposes is necessary to focus on the use of eustressors capable of establishing stable epigenetic marks to generate a transgenerational memory to stimulate a priming state in plants to face the changing environment.

BcMettl4-Mediated DNA Adenine N6-Methylation Is Critical for Virulence of Botrytis cinerea

DNA adenine N⁶-methylation (6mA) plays a critical role in various biological functions, but its occurrence and functions in filamentous plant pathogens are largely unexplored. Botrytis cinerea is an important pathogenic fungus worldwide. A systematic analysis of 6mA in B. cinerea was performed in this study, revealing that 6mA is widely distributed in the genome of this fungus. The 2 kb regions flanking many genes, particularly the upstream promoter regions, were susceptible to methylation. The role of BcMettl4, a 6mA methyltransferase, in the virulence of B. cinerea was investigated. BcMETTL4 disruption and point mutations of its catalytic motif “DPPW” both resulted in significant 6mA reduction in the genomic DNA and in reduced virulence of B. cinerea. RNA-Seq analysis revealed a total of 13 downregulated genes in the disruption mutant ΔBcMettl4 in which methylation occurred at the promoter sites. These were involved in oxidoreduction, secretory pathways, autophagy and carbohydrate metabolism. Two of these genes, BcFDH and BcMFS2, were independently disrupted. Knockout of BcFDH led to reduced sclerotium formation, while disruption of BcMFS2 resulted in dramatically decreased conidium formation and pathogenicity. These observations indicated that 6mA provides potential epigenetic markers in B. cinerea and that BcMettl4 regulates virulence in this important plant pathogen.

DNA N6-Methyladenine Modification in Eukaryotic Genome

DNA methylation is treated as an important epigenetic mark in various biological activities. In the past, a large number of articles focused on 5 mC while lacking attention to N6-methyladenine (6 mA). The presence of 6 mA modification was previously discovered only in prokaryotes. Recently, with the development of detection technologies, 6 mA has been found in several eukaryotes, including protozoans, metazoans, plants, and fungi. The importance of 6 mA in prokaryotes and single-celled eukaryotes has been widely accepted. However, due to the incredibly low density of 6 mA and restrictions on detection technologies, the prevalence of 6 mA and its role in biological processes in eukaryotic organisms are highly debated. In this review, we first summarize the advantages and disadvantages of 6 mA detection methods. Then, we conclude existing reports on the prevalence of 6 mA in eukaryotic organisms. Next, we highlight possible methyltransferases, demethylases, and the recognition proteins of 6 mA. In addition, we summarize the functions of 6 mA in eukaryotes. Last but not least, we summarize our point of view and put forward the problems that need further research.

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.

Structural basis for MTA1c-mediated DNA N6-adenine methylation

DNA N6-adenine methylation (6 mA) has recently been found to play a crucial role in epigenetic regulation in eukaryotes. MTA1c, a newly discovered 6 mA methyltransferase complex in ciliates, is composed of MTA1, MTA9, p1 and p2 subunits and specifically methylates ApT dinucleotides, yet its mechanism of action remains unknown. Here, we report the structures of Tetrahymena thermophila MTA1 (TthMTA1), Paramecium tetraurelia MTA9 (PteMTA9)-TthMTA1 binary complex, as well as the structures of TthMTA1-p1-p2 and TthMTA1-p2 complexes in apo, S-adenosyl methionine-bound and S-adenosyl homocysteine-bound states. We show that MTA1 is the catalytically active subunit, p1 and p2 are involved in the formation of substrate DNA-binding channel, and MTA9 plays a structural role in the stabilization of substrate binding. We identify that MTA1 is a cofactor-dependent catalytically active subunit, which exhibits stable SAM-binding activity only after assembly with p2. Our structures and corresponding functional studies provide a more detailed mechanistic understanding of 6 mA methylation.

A fungal dioxygenase CcTet serves as a eukaryotic 6mA demethylase on duplex DNA

N⁶-methyladenosine (6mA) is a DNA modification that has recently been found to play regulatory roles during mammalian early embryo development and mitochondrial transcription. We found that a dioxygenase CcTet from the fungus Coprinopsis cinerea is also a dsDNA 6mA demethylase. It oxidizes 6mA to the intermediate N⁶-hydroxymethyladenosine (6hmA) with robust activity of 6mA-containing duplex DNA (dsDNA) as well as isolated genomics DNA. Structural characterization revealed that CcTet utilizes three flexible loop regions and two key residues-D337 and G331-in the active pocket to preferentially recognize substrates on dsDNA. A CcTet D337F mutant protein retained the catalytic activity on 6mA but lost activity on 5-methylcytosine. Our findings uncovered a 6mA demethylase that works on dsDNA, suggesting potential 6mA demethylation in fungi and elucidating 6mA recognition and the catalytic mechanism of CcTet. The CcTet D337F mutant protein also provides a chemical biology tool for future functional manipulation of DNA 6mA in vivo.

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.

PredMHC: An Effective Predictor of Major Histocompatibility Complex Using Mixed Features

The major histocompatibility complex (MHC) is a large locus on vertebrate DNA that contains a tightly linked set of polymorphic genes encoding cell surface proteins essential for the adaptive immune system. The groups of proteins encoded in the MHC play an important role in the adaptive immune system. Therefore, the accurate identification of the MHC is necessary to understand its role in the adaptive immune system. An effective predictor called PredMHC is established in this study to identify the MHC from protein sequences. Firstly, PredMHC encoded a protein sequence with mixed features including 188D, APAAC, KSCTriad, CKSAAGP, and PAAC. Secondly, three classifiers including SGD, SMO, and random forest were trained on the mixed features of the protein sequence. Finally, the prediction result was obtained by the voting of the three classifiers. The experimental results of the 10-fold cross-validation test in the training dataset showed that PredMHC can obtain 91.69% accuracy. Experimental results on comparison with other features, classifiers, and existing methods showed the effectiveness of PredMHC in predicting the MHC.

DNA N6-methyldeoxyadenosine in mammals and human disease

N⁶-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.

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.

Context-dependent DNA polymerization effects can masquerade as DNA modification signals

Background

Single molecule measurements of DNA polymerization kinetics provide a sensitive means to detect both secondary structures in DNA and deviations from primary chemical structure as a result of modified bases. In one approach to such analysis, deviations can be inferred by monitoring the behavior of DNA polymerase using single-molecule, real-time sequencing with zero-mode waveguide. This approach uses a Single Molecule Real Time (SMRT)-sequencing measurement of time between fluorescence pulse signals from consecutive nucleosides incorporated during DNA replication, called the interpulse duration (IPD).

Results

In this paper we present an analysis of loci with high IPDs in two genomes, a bacterial genome (E. coli) and a eukaryotic genome (C. elegans). To distinguish the potential effects of DNA modification on DNA polymerization speed, we paired an analysis of native genomic DNA with whole-genome amplified (WGA) material in which DNA modifications were effectively removed. Adenine modification sites for E. coli are known and we observed the expected IPD shifts at these sites in the native but not WGA samples. For C. elegans, such differences were not observed. Instead, we found a number of novel sequence contexts where IPDs were raised relative to the average IPDs for each of the four nucleotides, but for which the raised IPD was present in both native and WGA samples.

Conclusion

The latter results argue strongly against DNA modification as the underlying driver for high IPD segments for C. elegans, and provide a framework for separating effects of DNA modification from context-dependent DNA polymerase kinetic patterns inherent in underlying DNA sequence for a complex eukaryotic genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08471-2.

MGF6mARice: prediction of DNA N6-methyladenine sites in rice by exploiting molecular graph feature and residual block

DNA N6-methyladenine (6mA) is produced by the N6 position of the adenine being methylated, which occurs at the molecular level, and is involved in numerous vital biological processes in the rice genome. Given the shortcomings of biological experiments, researchers have developed many computational methods to predict 6mA sites and achieved good performance. However, the existing methods do not consider the occurrence mechanism of 6mA to extract features from the molecular structure. In this paper, a novel deep learning method is proposed by devising DNA molecular graph feature and residual block structure for 6mA sites prediction in rice, named MGF6mARice. Firstly, the DNA sequence is changed into a simplified molecular input line entry system (SMILES) format, which reflects chemical molecular structure. Secondly, for the molecular structure data, we construct the DNA molecular graph feature based on the principle of graph convolutional network. Then, the residual block is designed to extract higher level, distinguishable features from molecular graph features. Finally, the prediction module is used to obtain the result of whether it is a 6mA site. By means of 10-fold cross-validation, MGF6mARice outperforms the state-of-the-art approaches. Multiple experiments have shown that the molecular graph feature and residual block can promote the performance of MGF6mARice in 6mA prediction. To the best of our knowledge, it is the first time to derive a feature of DNA sequence by considering the chemical molecular structure. We hope that MGF6mARice will be helpful for researchers to analyze 6mA sites in rice.

Comprehensive Analysis of DNA 5-Methylcytosine and N6-Adenine Methylation by Nanopore Sequencing in Hepatocellular Carcinoma

DNA methylation is a widespread epigenetic signal in human genome. With Nanopore technology, differential methylation modifications including 5-methylcytosine (5mC) and 6-methyladenine (6mA) can be identified. 5mC is the most important modification in mammals, although 6mA may also function in growth and development as well as in pathogenesis. While the role of 5mC at CpG islands in promoter regions associated with transcriptional regulation has been well studied, but the relationship between 6mA and transcription is still unclear. Thus, we collected two pairs of tumor tissues and adjacent normal tissues from hepatocellular carcinoma (HCC) surgical samples for Nanopore sequencing and transcriptome sequencing. It was found that 2,373 genes had both 5mC and 6mA, along with up- and down-regulated methylation sites. These genes were regarded as unstable methylation genes. Compared with 6mA, 5mC had more inclined distribution of unstable methylation sites. Chi-square test showed that the levels of 5mC were consistent with both up- and down-regulated genes, but 6mA was not significant. Moreover, the top three unstable methylation genes, TBC1D3H, CSMD1, and ROBO2, were all related to cancer. Transcriptome and survival analyses revealed four potential tumor suppressor genes including KCNIP4, CACNA1C, PACRG, and ST6GALNAC3. In this study, we firstly proposed to combine 5mC and 6mA methylation sites to explore functional genes, and further research found top of these unstable methylation genes might be functional and some of them could serve as potential tumor suppressor genes. Our study provided a new solution for epigenetic regulation research and therapy of HCC.

Epigenomic features of DNA G-quadruplexes and their roles in regulating rice gene transcription

A DNA G-quadruplex (G4) is a non-canonical four-stranded nucleic acid structure involved in many biological processes in mammals. The current knowledge on plant DNA G4s, however, is limited; whether and how DNA G4s impact gene expression in plants is still largely unknown. Here, we applied a protocol referred to as BG4-DNA-IP-seq followed by a comprehensive characterization of DNA G4s in rice (Oryza sativa L.); we next integrated dG4s (experimentally detectable G4s) with existing omics data and found that dG4s exhibited differential DNA methylation between transposable element (TE) and non-TE genes. dG4 regions displayed genic-dependent enrichment of epigenomic signatures; finally, we showed that these sites displayed a positive association with expression of DNA G4-containing genes when located at promoters, and a negative association when located in the gene body, suggesting localization-dependent promotional/repressive roles of DNA G4s in regulating gene transcription. This study reveals interrelations between DNA G4s and epigenomic signatures, as well as implicates DNA G4s in modulating gene transcription in rice. Our study provides valuable resources for the functional characterization or bioengineering of some of key DNA G4s in rice.

N4-acetyldeoxycytosine DNA modification marks euchromatin regions in Arabidopsis thaliana

BACKGROUND

Direct analogs of chemically modified bases that carry important epigenetic information, such as 5-methylcytosine (m5C)/5-methyldeoxycytosine (5mC), 5-hydroxymethylcytosine (hm5C)/5-hydroxymethyldeoxycytosine (5hmC), and N6-methyladenosine (m6A)/N6-methyldeoxyadenosine (6mA), are detected in both RNA and DNA, respectively. The modified base N4-acetylcytosine (ac4C) is well studied in RNAs, but its presence and epigenetic roles in cellular DNA have not been explored.

RESULTS

Here, we demonstrate the existence of N4-acetyldeoxycytosine (4acC) in genomic DNA of Arabidopsis with multiple detection methods. Genome-wide profiling of 4acC modification reveals that 4acC peaks are mostly distributed in euchromatin regions and present in nearly half of the expressed protein-coding genes in Arabidopsis. 4acC is mainly located around transcription start sites and positively correlates with gene expression levels. Imbalance of 5mC does not directly affect 4acC modification. We also characterize the associations of 4acC with 5mC and histone modifications that cooperatively regulate gene expression. Moreover, 4acC is also detected in genomic DNA of rice, maize, mouse, and human by mass spectrometry.

CONCLUSIONS

Our findings reveal 4acC as a hitherto unknown DNA modification in higher eukaryotes. We identify potential interactions of this mark with other epigenetic marks in gene expression regulation.

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.

RNA demethylation increases the yield and biomass of rice and potato plants in field trials

RNA N⁶-methyladenosine (m⁶A) modifications are essential in plants. Here, we show that transgenic expression of the human RNA demethylase FTO in rice caused a more than threefold increase in grain yield under greenhouse conditions. In field trials, transgenic expression of FTO in rice and potato caused ~50% increases in yield and biomass. We demonstrate that the presence of FTO stimulates root meristem cell proliferation and tiller bud formation and promotes photosynthetic efficiency and drought tolerance but has no effect on mature cell size, shoot meristem cell proliferation, root diameter, plant height or ploidy. FTO mediates substantial m⁶A demethylation (around 7% of demethylation in poly(A) RNA and around 35% decrease of m⁶A in non-ribosomal nuclear RNA) in plant RNA, inducing chromatin openness and transcriptional activation. Therefore, modulation of plant RNA m⁶A methylation is a promising strategy to dramatically improve plant growth and crop yield.

Using k-mer embeddings learned from a Skip-gram based neural network for building a cross-species DNA N6-methyladenine site prediction model

This study used k-mer embeddings as effective feature to identify DNA N6-Methyladenine sites in plant genomes and obtained improved performance without substantial effort in feature extraction, combination and selection. Identification of DNA N6-methyladenine sites has been a very active topic of computational biology due to the unavailability of suitable methods to identify them accurately, especially in plants. Substantial results were obtained with a great effort put in extracting, heuristic searching, or fusing a diverse types of features, not to mention a feature selection step. In this study, we regarded DNA sequences as textual information and employed natural language processing techniques to decipher hidden biological meanings from those sequences. In other words, we considered DNA, the human life book, as a book corpus for training DNA language models. K-mer embeddings then were generated from these language models to be used in machine learning prediction models. Skip-gram neural networks were the base of the language models and ensemble tree-based algorithms were the machine learning algorithms for prediction models. We trained the prediction model on Rosaceae genome dataset and performed a comprehensive test on 3 plant genome datasets. Our proposed method shows promising performance with AUC performance approaching an ideal value on Rosaceae dataset (0.99), a high score on Rice dataset (0.95) and improved performance on Rice dataset while enjoying an elegant, yet efficient feature extraction process.

Epigenome and Epitranscriptome: Potential Resources for Crop Improvement

Crop breeding faces the challenge of increasing food demand, especially under climatic changes. Conventional breeding has relied on genetic diversity by combining alleles to obtain desired traits. In recent years, research on epigenetics and epitranscriptomics has shown that epigenetic and epitranscriptomic diversity provides additional sources for crop breeding and harnessing epigenetic and epitranscriptomic regulation through biotechnologies has great potential for crop improvement. Here, we review epigenome and epitranscriptome variations during plant development and in response to environmental stress as well as the available sources for epiallele formation. We also discuss the possible strategies for applying epialleles and epitranscriptome engineering in crop breeding.