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

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.

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.

Dynamic DNA N 6-adenine methylation (6mA) governs the encystment process, showcased in the unicellular eukaryote Pseudocohnilembus persalinus

The formation of resting cysts commonly found in unicellular eukaryotes is a complex and highly regulated survival strategy against environmental stress that involves drastic physiological and biochemical changes. Although most studies have focused on the morphology and structure of cysts, little is known about the molecular mechanisms that control this process. Recent studies indicate that DNA N 6-adenine methylation (6mA) could be dynamically changing in response to external stimuli; however, its potential role in the regulation of cyst formation remains unknown. We used the ciliate Pseudocohnilembus persalinus, which can be easily induced to form cysts to investigate the dynamic pattern of 6mA in trophonts and cysts. Single-molecule real-time (SMRT) sequencing reveals high levels of 6mA in trophonts that decrease in cysts, along with a conversion of symmetric 6mA to asymmetric 6mA. Further analysis shows that 6mA, a mark of active transcription, is involved in altering the expression of encystment-related genes through changes in 6mA levels and 6mA symmetric-to-asymmetric conversion. Most importantly, we show that reducing 6mA levels by knocking down the DNA 6mA methyltransferase PpAMT1 accelerates cyst formation. Taken together, we characterize the genome-wide 6mA landscape in P. persalinus and provide insights into the role of 6mA in gene regulation under environmental stress in eukaryotes. We propose that 6mA acts as a mark of active transcription to regulate the encystment process along with symmetric-to-asymmetric conversion, providing important information for understanding the molecular response to environmental cues from the perspective of 6mA modification.

Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets

N6-methyladenosine (m6A) is a post-transcriptional epigenetic change with transcriptional stability and functionality regulated by specific m6A-modifying enzymes. However, the significance of genes modified by m6A and enzymes specific to m6A regulation in the context of pulmonary arterial hypertension (PAH) remains largely unexplored. MeRIP-seq and RNA-seq were applied to explore variances in m6A and RNA expression within the pulmonary artery tissues of control and monocrotaline-induced PAH rats. Functional enrichments were analyzed using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. To screen candidate m6A-related genes, the STRING and Metascape databases were used to construct a protein-protein interaction network followed by a real-time PCR validation of their expression. The expression level of an m6A regulator was further investigated using immunohistochemical staining, immunofluorescence, and Western blot techniques. Additionally, proliferation assays were conducted on primary rat pulmonary artery smooth muscle cells (PASMCs). We identified forty-two differentially expressed genes that exhibited either hypermethylated or hypomethylated m6A. These genes are predominantly related to the extracellular matrix structure, MAPK, and PI3K/AKT pathways. A candidate gene, centromere protein F (CENPF), was detected with increased expression in the PAH group. Additionally, we first identified an m6A reader, leucine rich pentatricopeptide repeat containing (LRPPRC), which was downregulated in the PAH rat model. The in vitro downregulation of Lrpprc mediated by siRNA resulted in the enhanced proliferation and elevated expression of Cenpf mRNA in primary rat PASMCs. Our study revealed a modified transcriptome-wide m6A landscape and associated regulatory mechanisms in the pulmonary arteries of PAH rats, potentially offering a novel target for therapeutic strategies in the future.

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

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

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

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

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

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

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

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.

Shaping eukaryotic epigenetic systems by horizontal gene transfer

DNA methylation constitutes one of the pillars of epigenetics, relying on covalent bonds for addition and/or removal of chemically distinct marks within the major groove of the double helix. DNA methyltransferases, enzymes which introduce methyl marks, initially evolved in prokaryotes as components of restriction-modification systems protecting host genomes from bacteriophages and other invading foreign DNA. In early eukaryotic evolution, DNA methyltransferases were horizontally transferred from bacteria into eukaryotes several times and independently co-opted into epigenetic regulatory systems, primarily via establishing connections with the chromatin environment. While C5-methylcytosine is the cornerstone of plant and animal epigenetics and has been investigated in much detail, the epigenetic role of other methylated bases is less clear. The recent addition of N4-methylcytosine of bacterial origin as a metazoan DNA modification highlights the prerequisites for foreign gene co-option into the host regulatory networks, and challenges the existing paradigms concerning the origin and evolution of eukaryotic regulatory systems.

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.

Evaluation of N 6-methyldeoxyadenosine antibody-based genomic profiling in eukaryotes

Low-level DNA N 6-methyldeoxyadenosine (DNA-m6A) has recently been reported across various eukaryotes. Although anti-m6A antibody-based approaches are commonly used to measure DNA-m6A levels, this approach is known to be confounded by DNA secondary structures, RNA contamination, and bacterial contamination. To evaluate for these confounding features, we introduce an approach for systematically validating the selectivity of antibody-based DNA-m6A methods and use a highly selective anti-DNA-m6A antibody to reexamine patterns of DNA-m6A in C. reinhardtii, A. thaliana, and D. melanogaster Our findings raise caution about the use of antibody-based methods for endogenous m6A quantification and mapping in eukaryotes.

Exploration of the Nuclear Proteomes in the Ciliate Oxytricha trifallax

Nuclear dimorphism is a fundamental feature of ciliated protozoa, which have separate somatic and germline genomes in two distinct organelles within a single cell. The transcriptionally active somatic genome, contained within the physically larger macronucleus, is both structurally and functionally different from the silent germline genome housed in the smaller micronucleus. This difference in genome architecture is particularly exaggerated in Oxytricha trifallax, in which the somatic genome comprises tens of thousands of gene-sized nanochromosomes maintained at a high and variable ploidy, while the germline has a diploid set of megabase-scale chromosomes. To examine the compositional differences between the nuclear structures housing the genomes, we performed a proteomic survey of both types of nuclei and of macronuclear histones using quantitative mass spectrometry. We note distinct differences between the somatic and germline nuclei, with many functional proteins being highly enriched in one of the two nuclei. To validate our conclusions and the efficacy of nuclear separation, we used protein localization through a combination of transformations and immunofluorescence. We also note that the macronuclear histones strikingly display only activating marks, consistent with the conclusion that the macronucleus is the hub of transcription. These observations suggest that the compartmentalization of different genome features into separate structures has been accompanied by a similar specialization of nuclear components that maintain and facilitate the functions of the genomes specific to each nucleus.

Structural insights into DNA N6-adenine methylation by the MTA1 complex

N⁶-methyldeoxyadenine (6mA) has recently been reported as a prevalent DNA modification in eukaryotes. The Tetrahymena thermophila MTA1 complex consisting of four subunits, namely MTA1, MTA9, p1, and p2, is the first identified eukaryotic 6mA methyltransferase (MTase) complex. Unlike the prokaryotic 6mA MTases which have been biochemically and structurally characterized, the operation mode of the MTA1 complex remains largely elusive. Here, we report the cryogenic electron microscopy structures of the quaternary MTA1 complex in S-adenosyl methionine (SAM)-bound (2.6 Å) and S-adenosyl homocysteine (SAH)-bound (2.8 Å) states. Using an AI-empowered integrative approach based on AlphaFold prediction and chemical cross-linking mass spectrometry, we further modeled a near-complete structure of the quaternary complex. Coupled with biochemical characterization, we revealed that MTA1 serves as the catalytic core, MTA1, MTA9, and p1 likely accommodate the substrate DNA, and p2 may facilitate the stabilization of MTA1. These results together offer insights into the molecular mechanism underpinning methylation by the MTA1 complex and the potential diversification of MTases for N⁶-adenine methylation.

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.

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

ENet-6mA: Identification of 6mA Modification Sites in Plant Genomes Using ElasticNet and Neural Networks

N6-methyladenine (6mA) has been recognized as a key epigenetic alteration that affects a variety of biological activities. Precise prediction of 6mA modification sites is essential for understanding the logical consistency of biological activity. There are various experimental methods for identifying 6mA modification sites, but in silico prediction has emerged as a potential option due to the very high cost and labor-intensive nature of experimental procedures. Taking this into consideration, developing an efficient and accurate model for identifying N6-methyladenine is one of the top objectives in the field of bioinformatics. Therefore, we have created an in silico model for the classification of 6mA modifications in plant genomes. ENet-6mA uses three encoding methods, including one-hot, nucleotide chemical properties (NCP), and electron–ion interaction potential (EIIP), which are concatenated and fed as input to ElasticNet for feature reduction, and then the optimized features are given directly to the neural network to get classified. We used a benchmark dataset of rice for five-fold cross-validation testing and three other datasets from plant genomes for cross-species testing purposes. The results show that the model can predict the N6-methyladenine sites very well, even cross-species. Additionally, we separated the datasets into different ratios and calculated the performance using the area under the precision–recall curve (AUPRC), achieving 0.81, 0.79, and 0.50 with 1:10 (positive:negative) samples for F. vesca, R. chinensis, and A. thaliana, respectively.

Evolutionary bioenergetics of ciliates

Understanding why various organisms evolve alternative ways of living requires information on both the fitness advantages of phenotypic modifications and the costs of constructing and operating cellular features. Although the former has been the subject of a myriad of ecological studies, almost no attention has been given to how organisms allocate resources to alternative structures and functions. We address these matters by capitalizing on an array of observations on diverse ciliate species and from the emerging field of evolutionary bioenergetics. A relatively robust and general estimator for the total cost of a cell per cell cycle (in units of ATP equivalents) is provided, and this is then used to understand how the magnitudes of various investments scale with cell size. Among other things, we examine the costs associated with the large macronuclear genomes of ciliates, as well as ribosomes, various internal membranes, osmoregulation, cilia, and swimming activities. Although a number of uncertainties remain, the general approach taken may serve as blueprint for expanding this line of work to additional traits and phylogenetic lineages.

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.

Means, mechanisms and consequences of adenine methylation in DNA

N⁶-methyl-2'-deoxyadenosine (6mA or m⁶dA) has been reported in the DNA of prokaryotes and eukaryotes ranging from unicellular protozoa and algae to multicellular plants and mammals. It has been proposed to modulate DNA structure and transcription, transmit information across generations and have a role in disease, among other functions. However, its existence in more recently evolved eukaryotes remains a topic of debate. Recent technological advancements have facilitated the identification and quantification of 6mA even when the modification is exceptionally rare, but each approach has limitations. Critical assessment of existing data, rigorous design of future studies and further development of methods will be required to confirm the presence and biological functions of 6mA in multicellular eukaryotes.

Detecting DNA Methylations in the Hyperthermoacidophilic Crenarchaeon Sulfolobus acidocaldarius Using SMRT Sequencing

DNA methylations are one of the most well-known epigenetic modifications along with histone modifications and noncoding RNAs. They are found at specific sites along the DNA in all domains of life, with 5-mC and 6-mA/4-mC being well-characterized in eukaryotes and bacteria respectively, and they have not only been described as contributing to the structure of the double helix itself but also as regulators of DNA-based processes such as replication, transcription, and recombination. Different methods have been developed to accurately identify and/or map methylated motifs to decipher the involvement of DNA methylations in regulatory networks that affect the cellular state.Although DNA methylations have been detected along archaeal genomes, their involvement as regulators of DNA-based processes remains the least known. To highlight the importance of DNA methylations in the control of key cellular mechanisms and their dynamics in archaea cells, we have used single-molecule real-time (SMRT) sequencing. This sequencing technology allows the identification and direct mapping of the methylated motifs along the genome of an organism. In this chapter, we present a step-by-step protocol for detecting DNA methylations in the hyperthermophilic crenarchaeon Sulfolobus acidocaldarius using SMRT sequencing. This protocol can easily be adapted to other prokaryotes.

OUP accepted manuscript

Bacteriophage exclusion (‘BREX’) phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation, and restriction. We identified a previously uncharacterized gene in the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR (‘BREX Regulator’), demonstrate that it forms a homodimer and specifically binds a DNA target site upstream of its transcription start site. Deletion of the BrxR gene causes cell toxicity, reduces restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene, or a point mutation blocking its DNA binding ability, has little effect on restriction, implying that the BrxR coding sequence and BrxR protein play independent functional roles. We speculate that elements within the BrxR coding sequence are involved in cis regulation of anti-phage activity, while the BrxR protein itself plays an additional regulatory role, perhaps during horizontal transfer.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

In Vitro Evaluation of DNA Damage Effect Markers toward Five Nitrogen Mustards Based on Liquid Chromatography-Tandem Mass Spectrometry

Endogenous DNA lesions frequently occur due to internal effects such as oxidative stress, inflammation, endogenous alkylation, and epigenetic modifications. However, exposure to chemical toxicants from the environment, diet, or drugs can also induce significant endogenous DNA damage. The quantification of endogenous DNA damage effect markers might reflect the actual DNA damage level of chemical toxicants. Herein, we report a liquid chromatography-triple quadrupole tandem mass spectrometry (LC-QqQ MS/MS) method for simultaneous determination of eight representative endogenous DNA damage biomarkers, including five endogenous DNA damage effect markers (oxidative damage, 8-oxo-dG; lipid peroxidation, εdA and N²-Et-dG; inflammation, 5-Cl-dC; and endogenous alkylation, O⁶-Me-dG), and three epigenetic modifications (5-m-dC, 5-hm-dC, and N⁶-Me-dA). The method validation was performed, and the linear range was 0.05 pg to 2 ng (on-column), the limit of detection was 0.02 pg (on-column), and the precision, accuracy, matrix effect, and recovery were all between 85 and 115%. We then applied this method to evaluate endogenous DNA damage to human embryonic lung fibroblast cells exposed to five nitrogen mustards [NMs, i.e., HN1, HN2, HN3, chlorambucil (CB), and cyclophosphamide (CTX)], where curcumin exposure was used as a control due to its inability to induce the formation of endogenous DNA adducts. The amounts of eight DNA adducts in the low-, middle-, and high-concentration exposure groups of five NMs were almost all significantly different from those in the blank group (P < 0.05). We obtained a positive correlation between the contents of eight DNA damage biomarkers and the inhibition dose of five NMs, except for N²-Et-dG and 5-Cl-dC. Via further principal component analysis and partial least squares discriminant analysis, we clustered all NMs into three units with different cytotoxicity levels, that is, HN2 and HN1 (highly toxic), HN3 and CB (moderately toxic), and CTX (less toxic). Moreover, for the same concentration of HN1/2/3 exposure groups, as the cytotoxicity increased according to the order of HN3 < HN1 < HN2, the contents of 8-oxo-dG, 5-m-dC, 5-hm-dC, and N⁶-Me-dA increased, whereas the content of O⁶-Me-dG decreased. Therefore, the contents of these DNA damage effect markers were somewhat related to the cytotoxicity and concentration of NMs. We hope that this method will provide an alternative evaluation approach for the toxicological effects of NMs and the safety of the medication.

HILIC-MS/MS for the Determination of Methylated Adenine Nucleosides in Human Urine

N6-methyl-2'-deoxyadenosine (m6dA) is a newly discovered DNA epigenetic mark in mammals. N6-methyladenosine (m6A), 2'-O-methyladenosine (Am), N6,2'-O-dimethyladenosine (m6Am), and N6,N6-dimethyladenosine (m62A) are common RNA modifications. Previous studies illustrated the associations between the aberrations of these methylated adenosines in nucleic acids and cancer. Herein, we developed Fe3O4/graphene-based magnetic dispersive solid-phase extraction for the enrichment and hydrophilic interaction liquid chromatography-mass spectrometry (HILIC-MS/MS) for the measurements of m6dA, m6A, Am, m6Am, and m62A in human urine samples. We found that malic acid could improve the HILIC-based separation of these modified nucleosides and markedly enhance the sensitivity of their MS detection. With this method, we accurately quantified the contents of these modified adenine nucleosides in urine samples collected from gastric and colorectal cancer patients as well as healthy controls. We found that, relative to healthy controls, urinary m6dA and Am levels are significantly lower for gastric and colorectal cancer patients; while gastric cancer patients also exhibited lower levels of urinary m6A, the trend was opposite for colorectal cancer patients. Together, we developed a robust analytical method for simultaneous measurements of five methylated adenine nucleosides in human urine, and our results revealed an association between the levels of urinary methylated adenine nucleosides and the occurrence of gastric as well as colorectal cancers, suggesting the potential applications of these modified nucleosides as biomarkers for the early detection of these cancers.

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.

Epigenetic Mechanisms Mediate Nicotine-Induced Reward and Behaviour in Zebrafish

Nicotine induces long-term changes in the neural activity of the mesocorticolimbic reward pathway structures. The mechanisms involved in this process have not been fully characterized. The hypothesis discussed here proposed that epigenetic regulation participates in the installation of persistent adaptations and long-lasting synaptic plasticity generated by nicotine action on the mesolimbic dopamine neurons of zebrafish. The epigenetic mechanisms induced by nicotine entail histone and DNA chemical modifications, which have been described to lead to changes in gene expression. Among the enzymes that catalyze epigenetic chemical modifications, histone deacetylases (HDACs) remove acetyl groups from histones, thereby facilitating DNA relaxation and making DNA more accessible to gene transcription. DNA methylation, which is dependent on DNA methyltransferase (DNMTs) activity, inhibits gene expression by recruiting several methyl binding proteins that prevent RNA polymerase binding to DNA. In zebrafish, phenylbutyrate (PhB), an HDAC inhibitor, abolishes nicotine rewarding properties together with a series of typical reward-associated behaviors. Furthermore, PhB and nicotine alter long- and short-term object recognition memory in zebrafish, respectively. Regarding DNA methylation effects, a methyl group donor L-methionine (L-met) was found to dramatically reduce nicotine-induced conditioned place preference (CPP) in zebrafish. Simultaneous treatment with DNMT inhibitor 5-aza-2’-deoxycytidine (AZA) was found to reverse the L-met effect on nicotine-induced CPP as well as nicotine reward-specific effects on genetic expression in zebrafish. Therefore, pharmacological interventions that modulate epigenetic regulation of gene expression should be considered as a potential therapeutic method to treat nicotine addiction.

Molecular Epigenetics: Chemical Biology Tools Come of Age

The field of epigenetics has exploded over the last two decades, revealing an astonishing level of complexity in the way genetic information is stored and accessed in eukaryotes. This expansion of knowledge, which is very much ongoing, has been made possible by the availability of evermore sensitive and precise molecular tools. This review focuses on the increasingly important role that chemistry plays in this burgeoning field. In an effort to make these contributions more accessible to the nonspecialist, we group available chemical approaches into those that allow the covalent structure of the protein and DNA components of chromatin to be manipulated, those that allow the activity of myriad factors that act on chromatin to be controlled, and those that allow the covalent structure and folding of chromatin to be characterized. The application of these tools is illustrated through a series of case studies that highlight how the molecular precision afforded by chemistry is being used to establish causal biochemical relationships at the heart of epigenetic regulation.

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.

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

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

Transcribed germline-limited coding sequences in Oxytricha trifallax

The germline-soma divide is a fundamental distinction in developmental biology, and different genes are expressed in germline and somatic cells throughout metazoan life cycles. Ciliates, a group of microbial eukaryotes, exhibit germline-somatic nuclear dimorphism within a single cell with two different genomes. The ciliate Oxytricha trifallax undergoes massive RNA-guided DNA elimination and genome rearrangement to produce a new somatic macronucleus (MAC) from a copy of the germline micronucleus (MIC). This process eliminates noncoding DNA sequences that interrupt genes and also deletes hundreds of germline-limited open reading frames (ORFs) that are transcribed during genome rearrangement. Here, we update the set of transcribed germline-limited ORFs (TGLOs) in O. trifallax. We show that TGLOs tend to be expressed during nuclear development and then are absent from the somatic MAC. We also demonstrate that exposure to synthetic RNA can reprogram TGLO retention in the somatic MAC and that TGLO retention leads to transcription outside the normal developmental program. These data suggest that TGLOs represent a group of developmentally regulated protein-coding sequences whose gene expression is terminated by DNA elimination.

Nucleosome Positioning and Spacing: From Mechanism to Function

Eukaryotes associate their genomes with histone proteins, forming nucleosome particles. Nucleosomes regulate and protect the genetic information. They often assemble into evenly spaced arrays of nucleosomes. These regular nucleosome arrays cover significant portions of the genome, in particular over genes. The presence of these evenly spaced nucleosome arrays is highly conserved throughout the entire eukaryotic domain. Here, we review the mechanisms behind the establishment of this primary structure of chromatin with special emphasis on the biogenesis of evenly spaced nucleosome arrays. We highlight the roles that transcription, nucleosome remodelers, DNA sequence, and histone density play towards the formation of evenly spaced nucleosome arrays and summarize our current understanding of their cellular functions. We end with key unanswered questions that remain to be explored to obtain an in-depth understanding of the biogenesis and function of the nucleosome landscape.

The Compact Macronuclear Genome of the Ciliate Halteria grandinella: A Transcriptome-Like Genome with 23,000 Nanochromosomes

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions.

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions. This challenges our usual understanding of chromosomal structure and suggests the possibility of novel mechvanisms in transcriptional regulation. Comprehensive analysis of multiple data sets reveals that Halteria transcription dynamics are influenced by: (i) nonuniform nanochromosome copy numbers correlated with gene-expression level; (ii) dynamic alterations at both the DNA and RNA levels, including alternative internal eliminated sequence (IES) deletions during macronucleus formation and large-scale alternative splicing in transcript maturation; and (iii) extremely short 5′ and 3′ untranslated regions (UTRs) and universal TATA box-like motifs in the compact 5′ subtelomeric regions of most chromosomes. This study broadens the view of ciliate biology and the evolution of unicellular eukaryotes, and identifies Halteria as one of the most compact known eukaryotic genomes, indicating that complex cell structure does not require complex gene architecture.

DNAModAnnot: a R toolbox for DNA modification filtering and annotation

MOTIVATION

Long-read sequencing technologies can be employed to detect and map DNA modifications at the nucleotide resolution on a genome-wide scale. However, published software packages neglect the integration of genomic annotation and comprehensive filtering when analyzing patterns of modified bases detected using Pacific Biosciences (PacBio) or Oxford Nanopore Technologies (ONT) data. Here, we present DNAModAnnot, a R package designed for the global analysis of DNA modification patterns using adapted filtering and visualization tools.

RESULTS

We tested our package using PacBio sequencing data to analyze patterns of the 6-methyladenine (6 mA) in the ciliate Paramecium tetraurelia, in which high 6 mA amounts were previously reported. We found Paramecium tetraurelia 6 mA genome-wide distribution to be similar to other ciliates. We also performed 5-methylcytosine (5mC) analysis in human lymphoblastoid cells using ONT data and confirmed previously known patterns of 5mC. DNAModAnnot provides a toolbox for the genome-wide analysis of different DNA modifications using PacBio and ONT long-read sequencing data.

AVAILABILITY

DNAModAnnot is distributed as a R package available via GitHub (https://github.com/AlexisHardy/DNAModAnnot).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.