N6-methyldeoxyadenosine directs nucleosome positioning in Tetrahymena DNA

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.

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.

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.

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.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

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.

Identification and Mechanism of Action of the Global Secondary Metabolism Regulator SaraC in Stereum hirsutum

DNA methylation is an important factor in the regulation of gene expression. In analyzing genomic data of Stereum hirsutum FP-91666, we found a hypothetical bifunctional transcription regulator/O6Meguanine-DNA methyltransferase (named SaraC), which is widely present in both bacteria and fungi, and confirmed that its function in bacteria is mainly for DNA reparation. In this paper, we confirmed that SaraC has the function of DNA binding and demethylation through surface plasma resonance and reaction experiments in vitro. Then, we achieved the overexpression of SaraC (OES) in S. hirsutum, sequenced the methylation and transcription levels of the whole-genome, and further conducted untargeted metabolomics analyses of the OES transformants and the wild type (WT). The results confirmed that the overall-methylation levels of the transformants were significantly downregulated, and various genes related to secondary metabolism were upregulated. Through comparative untargeted metabolomic analyses, it showed that OES SA6 transformant produced a greater number of hybrid polyketides, and we identified 2 novel hybrid polyketides from the fermentation products of SA6. Our results show that overexpression SaraC can effectively stimulate the expression of secondary-metabolism-related genes, which could be a broad-spectrum tool for discovery of metabolites due to its cross-species conservation. IMPORTANCE Fungi are one of the important sources of active compounds. However, in fungi, most of the secondary metabolic biosynthetic gene clusters are weakly expressed or silenced under conventional culture conditions. How to efficiently excavate potential new compounds contained in fungi is becoming a research hot spot in the world. In this study, we found a DNA demethylation protein (SaraC) and confirmed that it is a global secondary metabolism regulator in Stereum hirsutum FP-91666. In the past, SaraC-like proteins were mainly regarded as DNA repair proteins, but our findings proved that it will be a powerful tool for mining secondary metabolites for overexpression of SaraC, which can effectively stimulate the expression of genes related to secondary metabolism.

METTL4-mediated nuclear N6-deoxyadenosine methylation promotes metastasis through activating multiple metastasis-inducing targets

BACKGROUND

DNA N6-methyldeoxyadenosine (6mA) is rarely present in mammalian cells and its nuclear role remains elusive.

RESULTS

Here we show that hypoxia induces nuclear 6mA modification through a DNA methyltransferase, METTL4, in hypoxia-induced epithelial-mesenchymal transition (EMT) and tumor metastasis. Co-expression of METTL4 and 6mA represents a prognosis marker for upper tract urothelial cancer patients. By RNA sequencing and 6mA chromatin immunoprecipitation-exonuclease digestion followed by sequencing, we identify lncRNA RP11-390F4.3 and one novel HIF-1α co-activator, ZMIZ1, that are co-regulated by hypoxia and METTL4. Other genes involved in hypoxia-mediated phenotypes are also regulated by 6mA modification. Quantitative chromatin isolation by RNA purification assay shows the occupancy of lncRNA RP11-390F4.3 on the promoters of multiple EMT regulators, indicating lncRNA-chromatin interaction. Knockdown of lncRNA RP11-390F4.3 abolishes METTL4-mediated tumor metastasis. We demonstrate that ZMIZ1 is an essential co-activator of HIF-1α.

CONCLUSIONS

We show that hypoxia results in enriched 6mA levels in mammalian tumor cells through METTL4. This METTL4-mediated nuclear 6mA deposition induces tumor metastasis through activating multiple metastasis-inducing genes. METTL4 is characterized as a potential therapeutic target in hypoxic tumors.

The SWC4 subunit of the SWR1 chromatin remodeling complex is involved in varying virulence of Metarhizium brunneum isolates offering role of epigenetic regulation of pathogenicity

The host – pathogen interaction is a multifactorial process subject to a co-evolutionary arms race consisting of rapid changes in both host and pathogen, controlled at the genetic and epigenetic levels. Previously, we showed intra-species variation in disease progression and pathogenicity in aphids for Metarhizium brunneum isolates MbK and Mb7. Herein, we compared genomic, epigenetic, and metabolomic variations between these isolates and their effects on pathogenicity. Genomic variation could not completely explain the observed differences between the isolates. However, differential N6-adenine methylation (6 mA) and its correlation to reduced expression of the essential SWC4 subunit of SWR1 chromatin-remodelling complex (SWR1-C) led us to hypothesize a role for swc4 in the varying pathogenicity. Mutagenesis of the essential swc4 gene in MbKisolate resulted in reduction of secondary-metabolite (SM) secretion and impaired virulence in Galleria mellonella. Our results suggest the role of SWC4 in the regulation of SMs and the role of both SWC4 and SWR1-C in virulence of M. brunneum isolates. A better understanding of epigenetic regulation of SM production and secretion in entomopathogenic fungi may enable theirmanipulation for better biocontrol performance, and expand possibilities for environmentally friendly pest control.

ALKBH4 Stabilization Is Required for Arsenic-Induced 6mA DNA Methylation Inhibition, Keratinocyte Malignant Transformation, and Tumorigenicity

Inorganic arsenic is one of the well-known human skin carcinogens. However, the molecular mechanism by which arsenic promotes carcinogenesis remains unclear. Previous studies have established that epigenetic changes, including changes in DNA methylation, are among the critical mechanisms that drive carcinogenesis. N⁶-methyladenine (6mA) methylation on DNA is a widespread epigenetic modification that was initially found on bacterial and phage DNA. Only recently has 6mA been identified in mammalian genomes. However, the function of 6mA in gene expression and cancer development is not well understood. Here, we show that chronic low doses of arsenic induce malignant transformation and tumorigenesis in keratinocytes and lead to the upregulation of ALKBH4 and downregulation of 6mA on DNA. We found that reduced 6mA levels in response to low levels of arsenic were mediated by the upregulation of the 6mA DNA demethylase ALKBH4. Moreover, we found that arsenic increased ALKBH4 protein levels and that ALKBH4 deletion impaired arsenic-induced tumorigenicity in vitro and in mice. Mechanistically, we found that arsenic promoted ALKBH4 protein stability through reduced autophagy. Together, our findings reveal that the DNA 6mA demethylaseALKBH4 promotes arsenic tumorigenicity and establishes ALKBH4 as a promising target for arsenic-induced tumorigenesis.

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.

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.

Enhanced nucleosome assembly at CpG sites containing an extended 5-methylcytosine analogue

Methylation of cytosine to 5-methylcytosine (mC) at CpG sites is a prevalent reversible epigenetic mark in vertebrates established by DNA methyltransferases (MTases); the attached methyl groups can alter local structure of DNA and chromatin as well as binding of dedicated proteins. Nucleosome assembly on methylated DNA has been studied extensively, however little is known how the chromatin structure is affected by larger chemical variations in the major groove of DNA. Here, we studied the nucleosome formation in vitro on DNA containing an extended 5mC analog, 5-(6-azidohex-2-ynyl)cytosine (ahyC) installed at biological relevant CpG sites. We found that multiple ahyC residues on 80-Widom and Hsp70 promoter DNA fragments proved compatible with nucleosome assembly. Moreover, unlike mC, ahyC increases the affinity of histones to the DNA, partially altering nucleosome positioning, stability, and the action of chromatin remodelers. Based on molecular dynamics calculations, we suggest that these new features are due to increased DNA flexibility at ahyC-modified sites. Our findings provide new insights into the biophysical behavior of modified DNA and open new ways for directed design of synthetic nucleosomes.

Broad domains of histone marks in the highly compact Paramecium macronuclear genome

The unicellular ciliate Paramecium contains a large vegetative macronucleus with several unusual characteristics, including an extremely high coding density and high polyploidy. As macronculear chromatin is devoid of heterochromatin, our study characterizes the functional epigenomic organization necessary for gene regulation and proper Pol II activity. Histone marks (H3K4me3, H3K9ac, H3K27me3) reveal no narrow peaks but broad domains along gene bodies, whereas intergenic regions are devoid of nucleosomes. Our data implicate H3K4me3 levels inside ORFs to be the main factor associated with gene expression, and H3K27me3 appears in association with H3K4me3 in plastic genes. Silent and lowly expressed genes show low nucleosome occupancy, suggesting that gene inactivation does not involve increased nucleosome occupancy and chromatin condensation. Because of a high occupancy of Pol II along highly expressed ORFs, transcriptional elongation appears to be quite different from that of other species. This is supported by missing heptameric repeats in the C-terminal domain of Pol II and a divergent elongation system. Our data imply that unoccupied DNA is the default state, whereas gene activation requires nucleosome recruitment together with broad domains of H3K4me3. In summary, gene activation and silencing in Paramecium run counter to the current understanding of chromatin biology.

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.

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles' local geometry consists of straight segments connected by sharp bends compressing the DNA's inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles' inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

DNA Methylation on N6-Adenine Regulates the Hyphal Development during Dimorphism in the Early-Diverging Fungus Mucor lusitanicus

The epigenetic modifications control the pathogenicity of human pathogenic fungi, which have been poorly studied in Mucorales, causative agents of mucormycosis. This order belongs to a group referred to as early-diverging fungi that are characterized by high levels of N6-methyldeoxy adenine (6mA) in their genome with dense 6mA clusters associated with actively expressed genes. AlkB enzymes can act as demethylases of 6mA in DNA, with the most remarkable eukaryotic examples being mammalian ALKBH1 and Caenorhabditis elegans NMAD-1. The Mucor lusitanicus (formerly M. circinelloides f. lusitanicus) genome contains one gene, dmt1, and two genes, dmt2 and dmt3, encoding proteins similar to C. elegans NMAD-1 and ALKBH1, respectively. The function of these three genes was analyzed by the generation of single and double deletion mutants for each gene. Multiple processes were studied in the mutants, but defects were only found in single and double deletion mutants for dmt1. In contrast to the wild-type strain, dmt1 mutants showed an increase in 6mA levels during the dimorphic transition, suggesting that 6mA is associated with dimorphism in M. lusitanicus. Furthermore, the spores of dmt1 mutants challenged with macrophages underwent a reduction in polar growth, suggesting that 6mA also has a role during the spore–macrophage interaction that could be important in the infection process.

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.

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.

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

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

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

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

A distinct class of eukaryotic MT-A70 methyltransferases maintain symmetric DNA N6-adenine methylation at the ApT dinucleotides as an epigenetic mark associated with transcription

Rediscovered as a potential eukaryotic epigenetic mark, DNA N⁶-adenine methylation (6mA) varies across species in abundance and its relationships with transcription. Here we characterize AMT1—representing a distinct MT-A70 family methyltransferase—in the ciliate Tetrahymena thermophila. AMT1 loss-of-function leads to severe defects in growth and development. Single Molecule, Real-Time (SMRT) sequencing reveals that AMT1 is required for the bulk of 6mA and all symmetric methylation at the ApT dinucleotides. The detection of hemi-methylated ApT sites suggests a semi-conservative mechanism for maintaining symmetric methylation. AMT1 affects expression of many genes; in particular, RAB46, encoding a Rab family GTPase involved in contractile vacuole function, is likely a direct target. The distribution of 6mA resembles H3K4 methylation and H2A.Z, two conserved epigenetic marks associated with RNA polymerase II transcription. Furthermore, strong 6mA and nucleosome positioning in wild-type cells is attenuated in ΔAMT1 cells. Our results support that AMT1-catalyzed 6mA is an integral part of the transcription-associated epigenetic landscape. AMT1 homologues are generally found in protists and basal fungi featuring ApT hyper-methylation associated with transcription, which are missing in animals, plants, and true fungi. This dichotomy of 6mA functions and the underlying molecular mechanisms may have implications in eukaryotic diversification.

Ten-eleven translocation proteins and their role beyond DNA demethylation - what we can learn from the fly

Ten-eleven Translocation (TET) proteins have emerged as a family of epigenetic regulators that are important during development and have been implicated in various types of cancers. TET is a highly conserved protein that has orthologues in almost all multicellular organisms. Here, we review recent literature on the novel substrate specificity of this family of DNA 5-methylcytosine demethylases on DNA 6-methyladenine and RNA 5-methylcytosine that were first identified in the invertebrate model Drosophila. We focus on the biological role of these novel epigenetic marks in the fruit fly and mammals and highlight TET proteins' critical function during development specifically in brain development.

N6-Deoxyadenosine Methylation in Mammalian Mitochondrial DNA

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

Long-range single-molecule mapping of chromatin accessibility in eukaryotes

Mapping open chromatin regions has emerged as a widely used tool for identifying active regulatory elements in eukaryotes. However, existing approaches, limited by reliance on DNA fragmentation and short-read sequencing, cannot provide information about large-scale chromatin states or reveal coordination between the states of distal regulatory elements. We have developed a method for profiling the accessibility of individual chromatin fibers, a single-molecule long-read accessible chromatin mapping sequencing assay (SMAC-seq), enabling the simultaneous, high-resolution, single-molecule assessment of chromatin states at multikilobase length scales. Our strategy is based on combining the preferential methylation of open chromatin regions by DNA methyltransferases with low sequence specificity, in this case EcoGII, an N6-methyladenosine (m6A) methyltransferase, and the ability of nanopore sequencing to directly read DNA modifications. We demonstrate that aggregate SMAC-seq signals match bulk-level accessibility measurements, observe single-molecule nucleosome and transcription factor protection footprints, and quantify the correlation between chromatin states of distal genomic elements.

Epigenetic Modifications of mRNA and DNA in Plants

Advances in the detection and mapping of messenger RNA (mRNA) N⁶-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C), and DNA N⁶-methyldeoxyadenosine (6mA) redefined our understanding of these modifications as additional tiers of epigenetic regulation. In plants, the most prevalent internal mRNA modifications, m⁶A and m⁵C, play crucial and dynamic roles in many processes, including embryo development, stem cell fate determination, trichome branching, leaf morphogenesis, floral transition, stress responses, fruit ripening, and root development. The newly identified and widespread epigenetic marker 6mA DNA methylation is associated with gene expression, plant development, and stress responses. Here, we review the latest research progress on mRNA and DNA epigenetic modifications, including the detection, dynamics, distribution, functions, regulatory proteins, and evolution, with a focus on m⁶A, m⁵C, and 6mA. We also provide some perspectives on future research of the newly identified and unknown epigenetic modifications of mRNA and DNA in plants.

DNA Modification Readers and Writers and Their Interplay

Genomic DNA is modified in a postreplicative manner and several modifications, the enzymes responsible for their deposition as well as proteins that read these modifications, have been described. Here, we focus on the impact of DNA modifications on the DNA helix and review the writers and readers of cytosine modifications and how they interplay to shape genome composition, stability, and function.

N(6)-Methyladenine in eukaryotes

DNA modifications are a major form of epigenetic regulation that eukaryotic cells utilize in concert with histone modifications. While much work has been done elucidating the role of 5-methylcytosine over the past several decades, only recently has it been recognized that N(6)-methyladenine (N6-mA) is present in quantifiable and biologically active levels in the DNA of eukaryotic cells. Unlike prokaryotes which utilize N6-mA to recognize "self" from "foreign" DNA, eukaryotes have been found to use N6-mA in varying ways, from regulating transposable elements to gene regulation in response to hypoxia and stress. In this review, we examine the current state of the N6-mA in research field, and the current understanding of the biochemical mechanisms which deposit and remove N6-mA from the eukaryotic genome.

Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures

Web 3DNA (w3DNA) 2.0 is a significantly enhanced version of the widely used w3DNA server for the analysis, visualization, and modeling of 3D nucleic-acid-containing structures. Since its initial release in 2009, the w3DNA server has continuously served the community by making commonly-used features of the 3DNA suite of command-line programs readily accessible. However, due to the lack of updates, w3DNA has clearly shown its age in terms of modern web technologies and it has long lagged behind further developments of 3DNA per se. The w3DNA 2.0 server presented here overcomes all known shortcomings of w3DNA while maintaining its battle-tested characteristics. Technically, w3DNA 2.0 implements a simple and intuitive interface (with sensible defaults) for increased usability, and it complies with HTML5 web standards for broad accessibility. Featurewise, w3DNA 2.0 employs the most recent version of 3DNA, enhanced with many new functionalities, including: the automatic handling of modified nucleotides; a set of 'simple' base-pair and step parameters for qualitative characterization of non-Watson-Crick double-helical structures; new structural parameters that integrate the rigid base plane and the backbone phosphate group, the two nucleic acid components most reliably determined with X-ray crystallography; in silico base mutations that preserve the backbone geometry; and a notably improved module for building models of single-stranded RNA, double-helical DNA, Pauling triplex, G-quadruplex, or DNA structures 'decorated' with proteins. The w3DNA 2.0 server is freely available, without registration, at http://web.x3dna.org.

Epigenetic Methylations on N6-Adenine and N6-Adenosine with the same Input but Different Output

Epigenetic modifications on individual bases in DNA and RNA can encode inheritable genetic information beyond the canonical bases. Among the nucleic acid modifications, DNA N6-methadenine (6mA) and RNA N6-methyladenosine (m6A) have recently been well-studied due to the technological development of detection strategies and the functional identification of modification enzymes. The current findings demonstrate a wide spectrum of 6mA and m6A distributions from prokaryotes to eukaryotes and critical roles in multiple cellular processes. It is interesting that the processes of modification in which the methyl group is added to adenine and adenosine are the same, but the outcomes of these modifications in terms of their physiological impacts in organisms are quite different. In this review, we summarize the latest progress in the study of enzymes involved in the 6mA and m6A methylation machinery, including methyltransferases and demethylases, and their functions in various biological pathways. In particular, we focus on the mechanisms by which 6mA and m6A regulate the expression of target genes, and we highlight the future challenges in epigenetic regulation.

Adenine Methylation in Drosophila Is Associated with the Tissue-Specific Expression of Developmental and Regulatory Genes

N6-methyladenine (6mA or m6dA) is a DNA modification that has long been known to play an important role in a variety of biological functions in prokaryotes. This modification has only recently been described in eukaryotes, where it seems to have evolved species-specific functions ranging from nucleosome positioning to transposon repression. In Drosophila, 6mA has been shown to be important for enforcing the tissue specificity of neuronal genes in the brain and suppressing transposable element expression in the ovaries. In this study, we have analyzed the raw signal data from nanopore sequencing to identify 6mA positions in the D. melanogaster genome at single-base resolution. We find that this modification is enriched upstream from transcription start sites, within the introns and 3′ UTRs of genes, as well as in simple repeats. These 6mA positions are enriched for sequence motifs that are recognized by known transcriptional activators involved in development, such as Bicoid and Caudal, and the genes that carry this modification are enriched for functions involved in development, regulation of transcription, and neuronal activity. These genes show high expression specificity in a variety of tissues besides the brain, suggesting that this modification may play a more general role in enforcing the specificity of gene expression across many tissues, throughout development, and between the sexes.

Identification of a DNA N6-Adenine Methyltransferase Complex and Its Impact on Chromatin Organization

DNA N6-adenine methylation (6mA) has recently been described in diverse eukaryotes, spanning unicellular organisms to metazoa. Here, we report a DNA 6mA methyltransferase complex in ciliates, termed MTA1c. It consists of two MT-A70 proteins and two homeobox-like DNA-binding proteins and specifically methylates dsDNA. Disruption of the catalytic subunit, MTA1, in the ciliate Oxytricha leads to genome-wide loss of 6mA and abolishment of the consensus ApT dimethylated motif. Mutants fail to complete the sexual cycle, which normally coincides with peak MTA1 expression. We investigate the impact of 6mA on nucleosome occupancy in vitro by reconstructing complete, full-length Oxytricha chromosomes harboring 6mA in native or ectopic positions. We show that 6mA directly disfavors nucleosomes in vitro in a local, quantitative manner, independent of DNA sequence. Furthermore, the chromatin remodeler ACF can overcome this effect. Our study identifies a diverged DNA N6-adenine methyltransferase and defines the role of 6mA in chromatin organization.

N6-methyldeoxyadenosine directs nucleosome positioning in Tetrahymena DNA

BACKGROUND

N6-methyldeoxyadenosine (6mA or m6dA) was shown more than 40 years ago in simple eukaryotes. Recent studies revealed the presence of 6mA in more prevalent eukaryotes, even in vertebrates. However, functional characterizations have been limited.

RESULTS

We use Tetrahymena thermophila as a model organism to examine the effects of 6mA on nucleosome positioning. Independent methods reveal the enrichment of 6mA near and after transcription start sites with a periodic pattern and anti-correlation relationship with the positions of nucleosomes. The distribution pattern can be recapitulated by in vitro nucleosome assembly on native Tetrahymena genomic DNA but not on DNA without 6mA. Model DNA containing artificially installed 6mA resists nucleosome assembling compared to unmodified DNA in vitro. Computational simulation indicates that 6mA increases dsDNA rigidity, which disfavors nucleosome wrapping. Knockout of a potential 6mA methyltransferase leads to a transcriptome-wide change of gene expression.

CONCLUSIONS

These findings uncover a mechanism by which DNA 6mA assists to shape the nucleosome positioning and potentially affects gene expression.