- Open Access
N6-methyldeoxyadenosine directs nucleosome positioning in Tetrahymena DNA
- Guan-Zheng Luo†1, 2, 3Email authorView ORCID ID profile,
- Ziyang Hao†2, 3,
- Liangzhi Luo2, 3, 4,
- Mingren Shen5,
- Daniela Sparvoli6,
- Yuqing Zheng5,
- Zijie Zhang2, 3,
- Xiaocheng Weng2, 3,
- Kai Chen2, 3,
- Qiang Cui5,
- Aaron P. Turkewitz6 and
- Chuan He2, 3Email author
© The Author(s). 2018
- Received: 8 February 2018
- Accepted: 22 October 2018
- Published: 19 November 2018
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.
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.
These findings uncover a mechanism by which DNA 6mA assists to shape the nucleosome positioning and potentially affects gene expression.
- N 6-methyldeoxyadenosine
- DNA methylation
DNA modifications play a pivotal role in epigenetic regulation. 5-Methylcytosine (5mC) has been shown to participate in various biological processes as the most characterized DNA epigenetic mark in animals and plants . Another form of DNA modification, N6-methyldeoxyadenosine (6mA or m6dA), was discovered in the genomes of both prokaryotes and eukaryotes more than 40 years ago [2, 3]. In prokaryotes, 6mA is involved in numerous processes such as virus defense, DNA replication, DNA repair, transcription regulation, and transposition of DNA . However, the biological significance of 6mA in eukaryotes remained largely unknown. With newly developed high-throughput sequencing technology and sensitive mass spectrometry , 6mA has been recently found to be a potential epigenetic mark in the genomes of both unicellular and multicellular organisms [6–10] and was suggested to exist in vertebrates [11–13].
Distinct distribution patterns of 6mA have been revealed in various organisms, with versatile functional roles proposed . In our previous work, we found a periodic distribution of 6mA in the green alga Chlamydomonas, suggesting a potential role of 6mA affecting nucleosome positioning and gene expression . The periodic distribution pattern has also been discovered in another eukaryotic organism . Specifically, the genomic locations of the nucleosome and 6mA are anti-correlated, with 6mA marking linker regions between nucleosomes. Although the presence of 6mA has been reported in diverse eukaryotes, genomic distributions and functional implications vary in different species [14, 16]. So far, the periodic pattern around transcription start site (TSS) and the association with nucleosome positioning are exclusively reported in two unicellular organisms: Chlamydomonas and Tetrahymena [2, 15, 17], raising the question of how pervasive this interesting feature is in eukaryotes. In this study, we used Tetrahymena as our model system and explored the periodic distribution pattern of 6mA in vivo and in vitro, suggesting a more conserved functionality of 6mA on nucleosome positioning in eukaryotes.
Organized nucleosome positioning is crucial for gene expression . Constitution of a nucleosome array is directed by both the underlying DNA and chromatin remodelers [19, 20]. Nucleosome formation intrinsically disfavors certain DNA sequences in vitro, especially poly(dA-dT) sequences, suggesting the positioning of nucleosomes could be significantly affected by DNA sequences [21, 22]. The nucleosome locations can even be predicted based on genomic DNA sequences [22, 23]. In vitro nucleosome assembly successfully rebuilt the nucleosome-free region (NFR) around TSS in multiple systems, including yeast and human; however, the phased positioning of nucleosome with periodic pattern can barely be recapitulated . Based on the anti-correlated 6mA nucleosome pattern in both green alga and Tetrahymena, we hypothesize that the DNA 6mA methylation could also direct nucleosome positioning. Indeed, in vitro nucleosome assembly assay faithfully rebuilds the nucleosome arrays on native Tetrahymena genomic DNA but not on unmethylated DNA. We further used model DNA bearing 6mA modification to perform in vitro nucleosome assembly and measured the 6mA abundance of nucleosome-protected regions versus unprotected regions. We found the nucleosome-protected regions contain much less 6mA than unprotected regions, reinforcing our hypothesis that nucleosome wrapping avoids 6mA-containing DNA in a species-independent manner. Computational simulation indicates DNA flexibility can be modulated by a single 6mA site, which can further affect nucleosome positioning. By homology searching, we identified a MTA70 family protein, TAMT-1 (Tetrahymena deoxyadenosine methyltransferase-1), which methylates adenine in the ApT sequence context as demonstrated through an in vitro protein reactivity assay. Knockout of TAMT-1 reduced the total 6mA level and notably altered the transcriptome pattern. Together, we propose that 6mA can direct nucleosome positioning as a DNA marker, which subsequently affects gene expression in an epigenetic way.
Determining the 6mA methylome in Tetrahymena
6mA and nucleosomes are anti-correlated in vivo
6mA directs nucleosome positioning in vitro
To rule out the potential DNA sequence contributions to nucleosome positioning, we computationally predicted the nucleosome occupancy based on the genomic DNA sequence. As expected, the prediction recapitulated the NFR upstream of TSS and enrichment of nucleosome downstream of TSS but did not show any phased pattern (Additional file 1: Figure S3b). Next, we used Phi29 DNA polymerase to do unbiased whole-genome amplification (WGA) with a small amount of native genomic DNA (~ 5 ng) (Additional file 1: Figure S3c), which was followed by nucleosome assembly using the amplified DNA without methylation due to the dramatic dilution of the native genomic DNA (Fig. 3c, Additional file 1: Figure S3d). Although we could obtain clear ~ 150 bp bands after MNase digestion (Additional file 1: Figure S3e), many fewer nucleosome peaks could be identified by sequencing, suggesting that the nucleosome positioning is randomized on the unmethylated DNA substrate. By aligning the predicted nucleosome locations to gene promoters, we only found NFRs upstream of TSS, without detecting uniform arrays downstream of TSS (Fig. 3d). To compare the nucleosome positions relative to single 6mA sites, we depicted the accumulated nucleosome distribution around specific loci which had been identified to be methylated inside cells. Similar to the in vivo distribution, nucleosomes assembled on extracted genomic DNA exhibit the same periodic pattern, with a significant reduction at 6mA sites. Furthermore, nucleosomes assembled on unmethylated DNA distribute equally along the 6mA flanking regions, even with a slight enrichment at the adjacent area (Fig. 3e), presumably reflecting the sequence composition which in turn partially affects nucleosome positioning [22, 23].
To generalize our observation from Tetrahymena genomic DNA, we performed nucleosome assembly on DNA encoding a nucleosome positioning element. This 208 bp synthesized sequence is derived from the Lytechinus variegatus (sea urchin) 5S rDNA, which contains one possible binding site for a nucleosome octamer. A single GATC motif is present in the sequence (see the “Methods” section). After the treatment of Dam methylase, we obtained methylation at the GATC site as validated by UHPLC-QQQ-MS/MS (Additional file 1: Figure S3f). Then, we mixed unmethylated and methylated DNA at 1:1 ratio (Fig. 3f). After nucleosome assembly and MNase digestion, we used UHPLC-QQQ-MS/MS to measure the 6mA level of nucleosome-protected DNA versus the flow-through. The protected portion shows much lower 6mA level than the input or the flow-through portion (Fig. 3g), indicating a single 6mA site is sufficient to alter the preference of DNA wrapping on the nucleosome.
6mA modulates DNA flexibility
TAMT-1 is a 6mA methyltransferase in Tetrahymena
We next investigated potential methyltransferase(s) which could be responsible for 6mA methylation in Tetrahymena. Previous studies have predicted the existence of multiple DNA 6mA methyltransferases in several eukaryotic organisms, most of which can be attributed to three major families of Dam-like, DNMT, and MTA70 . Among them, the MTA70 family is evolved from MunI-like bacterial DNA 6mA methyltransferase and is widely conserved as RNA m6A methyltransferases, including IME4 in yeast, METTL3 in human, and the reported DNA 6mA methyltransferase DAMT-1 in C. elegans. The Tetrahymena genome encodes two genes belonging to this family, TTHERM_00136470 and THERM_00388490. Both of these genes consist of a MTA70 domain and a ZZ-type zinc finger domain for potential DNA binding. Sequence alignment and phylogenetic distribution analysis of the MTA70 family suggested that these two proteins might be primary candidates for DNA 6mA methylation in Tetrahymena (Additional file 1: Figure S6a, b).
We then sought to determine whether TAMT-1 could catalyze 6mA methylation in vitro. The low expression level of recombinant TAMT-1 and degradation of purified protein hindered us to get enough amounts of proteins to perform the entire enzymatic kinetic assay. Instead, we cloned and purified the catalytic domain of TAMT-1 consist of the MTA70 and the ZZ-type zinc finger domain (Additional file 1: Figure S7c). DNA probes with consensus sequences of CATG (Probe-1), GATC (Probe-2), and random AT (Probe-3) in double- and single-stranded forms were incubated with purified TMAT-1 separately. S-(5′-adenosyl)-l-methionine-d3 (d3-SAM) was used as the cofactor for accurate UPLC-QQQ-MS/MS quantification. We calculated the methylation yields by the molar ratio of newly formed d3-6mA to each digested DNA probe (Additional file 1: Figure S7d). The result showed that the catalytic domain of TAMT-1 exhibited no activity with Probe-1 and Probe-3, whereas it showed considerable methyltransferase activity towards Probe-2 with a GATC motif in double- and single-stranded DNA. We further measured the time course of the reaction between the catalytic domain of TAMT-1 and Probe-2 to further confirm the enzyme activity (Additional file 1: Figure S7e). In particular, the enzyme showed a strong preference for the single-stranded DNA probe (tenfold) compared to the corresponding duplex DNA (Fig. 5b). To exclude the possibility that the methylation activity was caused by the contamination of Escherichia coli methyltransferase(s) during purification, and further confirm that the methylation activity was directly mediated by TAMT-1, we mutated the methylation signature motif DPPW and highly conserved amino acid residues E111 and K183 which play important roles in catalysis . The result showed that mutations of DPPW to APPA, together with E111A and K183A, significantly decreased the in vitro methylation activity, suggesting that the methylation activity was mediated by TMAT-1 (Fig. 5c). Taken together, these results confirm that TAMT-1 is a 6mA methyltransferase in Tetrahymena.
Deficiency of 6mA perturbs gene expression
As a ciliated protozoan, Tetrahymena diverged from vertebrates or green alga more than two billion years ago [35–37]. The analogous distribution patterns of 6mA in green alga and Tetrahymena as well as the most recent discoveries in fungi suggest conserved biogenesis pathways and functions of this DNA modification inherited from an ancient common ancestor . The methylation machinery could also be acquired independently by convergent evolution.
We performed in vitro nucleosome array assembly using recombinant histones and model DNA mixed with Dam-methylated DNA from an unrelated organism. Mass spectrometry showed that nucleosomes intrinsically disfavor methylated DNA. Though in vitro experiments showed that 6mA directs nucleosome assembly without any other auxiliary proteins, 6mA may recruit partner proteins in vivo which could reinforce the chromatin architecture. We further showed that the mutual repulsion effect between nucleosome and 6mA is independent of species or DNA sequence. Because 6mA appears to be a prevalent DNA modification in many species including mammals, the effect of 6mA on directing nucleosome positioning in Tetrahymena may also applies to other organisms as this appears to be an intrinsic biophysical property associated with 6mA in DNA.
How could 6mA repel nucleosome wrapping? It has been shown that the presence of methylation could rigidify the double-stranded DNA structure . The methylation can reinforce the base stacking along the duplex. Another source of energy comes from reduced solvation penalty when 6mA is packed along the lineal duplex. A potential bending of the 6mA-containing duplex DNA such as wrapping around the nucleosome would expose the hydrophobic methyl group to the solvent water, inducing negative solvation penalty. The fully methylated 6mA on the opposite strands and clustering of 6mA methylation further amplify these effects, which collectively rigidify duplex DNA.
We identified a potential 6mA methyltransferase, TAMT-1, which demonstrates in vitro activity on ApT dinucleotides. TAMT-1 knockout suppresses 6mA levels in vivo. Only one 6mA methyltransferase has been found in C. elegans in the past . Intriguingly, these two candidates belong to the well-conserved MTA70 protein family, which also includes mRNA methyltransferases . The widespread distribution of this protein family may suggest a more extensive presence of 6mA in uncharacterized organisms [5, 16]. By disrupting TAMT-1, we observed a dramatic decline of 6mA levels in living cells. Additional methyltransferases may exist which could explain the remnant 6mA in the KO strain. Nevertheless, the discoveries of 6mA methyltransferases provided us with a useful target to manipulate 6mA and perform more thorough functional studies in the future.
It is worth noting that in certain organisms, Dnmt5 generates dense 5mC clusters which directly disfavor nucleosomes . A similar mechanistic argument with better packing and reduced solvation penalty could explain the effect of 5mC. It is thus very likely that 5mC may also affect nucleosome positioning in more widespread organisms. Interestingly, in our previous study, we discovered that 5mC and 6mA generally tend to avoid each other in the Chlamydomonas genome, suggesting complementary roles of different DNA modifications in one organism. Tetrahymena provides us with a unique model system to study the interaction between 6mA and nucleosome as it possesses only 6mA but not 5mC, reinforcing our notion that 6mA plays a key role in nucleosome positioning.
Until recently, 5mC and its derivatives were the well-established epigenetic marks in eukaryotic genomes. The recently discovered prevalence and functional implications of 6mA open a new avenue in epigenetic research. In this study, we found that the anti-correlation to nucleosome is a well-conserved feature of 6mA in evolutionary distinct species. Species-independent, 6mA-modified DNA is intrinsically disfavored by a nucleosome. In higher eukaryotes in which 5mC is recognized by reader proteins to suppress gene expression, the presence of 6mA could provide additional tuning of transcription by affecting nucleosome positioning or by resisting bending incurred by various DNA-binding proteins such as transcriptional factors. Overall, our results show the intrinsic properties of the DNA 6mA methylation can significantly impact gene expression and biological outcome in eukaryotes.
In this study, we reveal the intrinsic repulsion between 6mA and nucleosomes, which contributes to shape nucleosome positioning and affect gene expression. By disrupting the newly identified potential methyltransferase, we observed significant transcriptome changes. We propose that disordered nucleosome positioning caused by 6mA depletion is one of the main reasons underlying the observed phenotypes.
Cell strains and DNA/RNA collection
Cell strain SB210 was obtained from Tetrahymena Stock Center (https://tetrahymena.vet.cornell.edu/strains.php) and cultured in proteose peptone (PP) medium (2% proteose peptone, 10 μM FeCl3, or 90 μM sequestrene (Fe-EDTA)) at 31 °C. Culture densities were measured using Z1 Coulter Counter (Beckman Coulter). The genomic DNA was extracted by Quick-DNA™ Miniprep Kit (Zymo Research, Cat. No. D3024). Total RNA was extracted by Direct-zol™ RNA MiniPrep (Zymo Research, Cat. No. R2050) and further purified to mRNA by Dynabeads® mRNA Purification Kit (Thermo Fisher Scientific, Cat. No. 61006).
Generation of TAMT1 knockout strains
Eight hundred sixty-two base pairs of sequence from upstream of tamt-1 (TTHERM_00388490) coding region and a portion of the open reading frame (ORF) (1,139 bp) were amplified and subcloned into the pNEO4 vector flanking the neo4 cassette with SacI/PstI and XhoI/KpnI restriction cutting sites, respectively. The sequences of primers are listed in Additional file 2: Table S2. The construct resulted in the deletion of tamt-1 genomic region from 0 to 907 bp. The knockout vector was linearized by digestion with SacI and KpnI then transformed into CU428.1 cells by biolistic transformation following the reported method . To assess the gene disruption, RT-qPCR was performed.
Total RNA was isolated by using RNeasy Plus Mini kit (Qiagen). Five hundred nanograms total RNA were reverse-transcribed into cDNA with PrimeScript™ RT reagent Kit (Takara) and then subjected to qPCR analysis with FastStart SYBR Green Master Mix (Roche) in a Roche LightCycler 96. ATU1 were used as internal control. Relative changes in expression were calculated using the ΔΔCt method. All RT-qPCR primers are listed in Additional file 2: Table S2.
One microgram purified genomic DNA was segmented to ~ 250 bp by sonicator (Bioruptor, Diagenode). Then, DNA segments were end-repaired, 3′-adenylated, and ligated to Illumina adaptors by using NEBNext® DNA Library Prep Kit (NEB, Cat. No. E6040S). 6mA-containing DNA was enriched by antibody immunoprecipitation (SYSY, Cat. No. 202 003) then subjected to NGS (Illumina, HiSeq 2500). 6mA peaks were called by MACS2 by comparing the pull-down reads versus input DNA. Cutoff FDR < 0.01 was set to increase the reliability of 6mA candidates. Detailed protocol can be found in Reference .
Genomic DNA containing 6mA was first immunoprecipitated following a procedure similar to previously described with 6mA-IP-seq. Target DNA was then covalently cross-linked to the antibody by UV 254 nm irradiation, followed by lambda exonuclease digestion (NEB, Cat. No. M0262S). After adaptor ligation and PCR amplification (18 cycles), DNA library was constructed for NGS (Illumina, HiSeq 2500). Detailed protocol can be found in Reference .
Purified genomic DNA was digested by three different restriction enzymes separately (DpnI, DpnII, and CviAII). Digested DNA segments were further sheared to ~ 250 bp by sonicator (Bioruptor, Diagenode). DNA segments were end-repaired, 3′-adenylated, and ligated to sequencing adaptor by using NEBNext® DNA Library Prep Kit (NEB, Cat. No. E6040S). After PCR amplification for 15 cycles, purified DNA library was constructed for NGS (Illumina, HiSeq 2500). Detailed protocol can be found in References [6, 9].
Nucleosome positioning analysis
The nucleosome positioning is mainly determined by MNase-seq. Cells were lysed, and nuclei were isolated by Nuclei Isolation Kit (Sigma-Aldrich, Cat. No. NUC101). After MNase treatment for 10 min, DNA fragments were purified by DNA Clean & Concentrator Kit (Zymo Research, Cat. No. D4013) and loaded into agarose gel. The ~ 150-bp band was isolated and extracted by Zymoclean™ Gel DNA Recovery Kit (Zymo Research, Cat. No. D4007). DNA library was constructed by NEBNext® DNA Library Prep Kit (NEB, Cat. No. E6040S) according to the standard Illumina DNA library preparation procedures.
The 208-bp model DNA was purchased from NEB (Cat. No. N1202S), which was originated from Lytechinus variegatus (sea urchin) 5S rDNA and supposed to contain one nucleosome binding site. The entire sequence has only one GATC site (see below) where the A can be methylated to 6mA by Dam methyltransferase (NEB, Cat. No. M0222S), ACTTCCAGGGATTTATAAGCCGATGACGTCATAACATCCCTGACCCTTTAAATAGCTTAACTTTCATCAAGCAAGAGCCTACGACCATACCATGCTGAATATACCGGTTCTCGTCCG(A)TCACCGAAGTCAAGCAGCATAGGGCTCGGTTAGTACTTGGATGGGAGACCGCCTGGGAATACCGAATTCCCCGAGGAATTCCAACGAATA
In in vitro nucleosome assembly experiments, purified recombinant human histone H2A and H2B were mixed to generate dimer in advance, as well as histone H3.1/H4 tetramer. Then, 100 pm dimers and 50 pm tetramers were mixed with 50 pm DNA at 2 M NaCl. The salt concentration was sequentially diluted to allow the nucleosome generation on the DNA. Detailed procedures can be found in the EpiMark® Nucleosome Assembly Kit (NEB, Cat. No. E5350S). Assembled DNA/nucleosome component was then digested by MNase, and the protected region was purified for further analysis.
Nucleosome fuzziness is defined as the deviation of nucleosome positions within each unit in a cell population. The loci and fuzziness of nucleosome were calculated by software DANPOS .
Molecular dynamics simulations
The initial atomistic model of an unmodified 33-bp B-DNA fragment was built using w3DNA web interface  to the 3DNA software package . The sequence of the fragment was GCTCACCCGCGCCCATGGTGGGAGCCGGAGACG, where the positions of potential N6-methylation are in italics. The structure of the modified fragment was obtained by modifying the initial atomistic model using PyMol (http://pymol.org/). Molecular dynamics simulations were carried out using AMBER16 package with graphics processing units . Parmbsc1 force field was used in the simulations . Parameters for N6-methyladenosine was adapted from those of RNA N6-methyladenosine . The structures were solvated in TIP3P water  such that the closest distance between DNA atoms and the truncated octahedral water box edge was greater than 10 Å. The system also contained 150 mM NaCl after neutralizing with Na+ ions. Particle mesh Ewald was used to calculate the electrostatic interactions with a grid spacing of about 1.0 Å. The non-bonded cutoff was 12 Å with the missing long-range van der Waals interactions approximated with a long-range continuum correction . Energy minimization was done by 3000 steps of steepest descent followed by 97,000 steps of conjugate gradient. The DNA was fixed initially and then relaxed to minimize the energy of the entire system. One nanosecond of NPT simulation was used to equilibrate the density of the system. Production simulations were carried out in NVT ensemble at 300 K using Langevin dynamics with a collision frequency of 1 ps−1. SHAKE algorithm was used to constrain all bonds involving hydrogen atoms, which allowed an integration step of 2 fs. The trajectories were about 500 ns long for each construct. Conformation snapshots were saved every 20 ps. Structural parameters of DNA were computed using the 3DNA program . The error bars in Fig. 1 and Additional file 1: Figures S1 and S2 are standard errors of the mean calculated using block averaging and conventional error propagation rules.
Validation of 6mA sites
In this work, the fraction of input saved is 10 μl and the fraction for each 6mA-IP is 70 μl. The IP fraction is seven times to the input fraction. The equation above is as follows: ΔCt [normalized 6mA-IP] = (Ct [6mA-IP] − (Ct [Input] − Log2 (7)). The percentage (Input %) value for each sample is calculated as follows: Input % = 100/2 ΔCt [normalized 6mA-IP]. The “Input %” value represents the enrichment of certain 6mA modification on a specific region. Two replicates are performed and averaged to represent the relative abundance.
Sequencing was performed at the University of Chicago Genomics Facility. Imaging was performed at the University of Chicago Integrated Light Microscopy Facility.
This work was supported by the National Institutes of Health HG006827 (to CH), GM105783 (to APT), GM106443 (to QC), and NSF-DMS-1160360 (to QC), and by the National Natural Science Foundation of China (Nos. 91753129 and 31870808 to G-ZL). CH is an investigator of the Howard Hughes Medical Institute (HHMI). LL is supported by the China Scholarship Council (CSC, 201706670015). The Mass Spectrometry Facility of the University of Chicago is funded by the National Science Foundation (CHE-1048528).
Availability of data and materials
The high-throughput data used in this study are deposited in the NCBI GEO database with accession number GSE104699 .
G-ZL and CH designed the experiments. G-ZL and ZH performed most of the experiments with help from LL, ZZ, XW, and KC. ZH and DS constructed the KO cell line. MS and YZ did the molecular simulation. G-ZL, ZH, and CH wrote the manuscript with inputs from QC and APT. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
CH is a scientific founder of Accent Therapeutics and a member of the scientific advisory board.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Smith ZD, Meissner A. DNA methylation: roles in mammalian development. Nat Rev Genet. 2013;14(3):204–20.View ArticleGoogle Scholar
- Gorovsky MA, Hattman S, Pleger GL. (6N)methyl adenine in the nuclear DNA of a eucaryote, Tetrahymena pyriformis. J Cell Biol. 1973;56(3):697–701.View ArticleGoogle Scholar
- Hattman S, et al. Comparative study of DNA methylation in three unicellular eucaryotes. J Bacteriol. 1978;135(3):1156–7.PubMedPubMed CentralGoogle Scholar
- Wion D, Casadesus J. N6-methyl-adenine: an epigenetic signal for DNA-protein interactions. Nat Rev Microbiol. 2006;4(3):183–92.View ArticleGoogle Scholar
- Luo GZ, et al. DNA N(6)-methyladenine: a new epigenetic mark in eukaryotes? Nat Rev Mol Cell Biol. 2015;16(12):705–10.View ArticleGoogle Scholar
- Fu Y, et al. N6-methyldeoxyadenosine marks active transcription start sites in chlamydomonas. Cell. 2015;161(4):879–92.View ArticleGoogle Scholar
- Greer EL, et al. DNA methylation on N6-adenine in C. elegans. Cell. 2015;161(4):868–78.View ArticleGoogle Scholar
- Zhang G, et al. N6-methyladenine DNA modification in Drosophila. Cell. 2015;161(4):893–906.View ArticleGoogle Scholar
- Luo GZ, et al. Characterization of eukaryotic DNA N(6)-methyladenine by a highly sensitive restriction enzyme-assisted sequencing. Nat Commun. 2016;7:11301.View ArticleGoogle Scholar
- Mondo SJ, et al. Widespread adenine N6-methylation of active genes in fungi. Nat Genet. 2017;49(6):964–8.View ArticleGoogle Scholar
- Koziol MJ, et al. Identification of methylated deoxyadenosines in vertebrates reveals diversity in DNA modifications. Nat Struct Mol Biol. 2016;23(1):24–30.View ArticleGoogle Scholar
- Wu TP, et al. DNA methylation on N(6)-adenine in mammalian embryonic stem cells. Nature. 2016;532(7599):329–33.View ArticleGoogle Scholar
- Liu J, et al. Abundant DNA 6mA methylation during early embryogenesis of zebrafish and pig. Nat Commun. 2016;7:13052.View ArticleGoogle Scholar
- Pfeifer GP. Epigenetics: an elusive DNA base in mammals. Nature. 2016;532(7599):319–20.View ArticleGoogle Scholar
- Wang Y, et al. N6-adenine DNA methylation is associated with the linker DNA of H2A.Z-containing well-positioned nucleosomes in Pol II-transcribed genes in Tetrahymena. Nucleic Acids Res. 2017;45(20):11594–606.View ArticleGoogle Scholar
- Iyer LM, Zhang D, Aravind L. Adenine methylation in eukaryotes: apprehending the complex evolutionary history and functional potential of an epigenetic modification. Bioessays. 2016;38(1):27–40.View ArticleGoogle Scholar
- Wang Y, et al. N6-methyladenine DNA modification in the unicellular eukaryotic organism Tetrahymena thermophila. Eur J Protistol. 2016;58:94–102.View ArticleGoogle Scholar
- Kornberg RD, Lorch Y. Twenty-five years of the nucleosome, fundamental particle of the eukaryote chromosome. Cell. 1999;98(3):285–94.View ArticleGoogle Scholar
- Vignali M, et al. ATP-dependent chromatin-remodeling complexes. Mol Cell Biol. 2000;20(6):1899–910.View ArticleGoogle Scholar
- Yen K, et al. Genome-wide nucleosome specificity and directionality of chromatin remodelers. Cell. 2012;149(7):1461–73.View ArticleGoogle Scholar
- Sekinger EA, Moqtaderi Z, Struhl K. Intrinsic histone-DNA interactions and low nucleosome density are important for preferential accessibility of promoter regions in yeast. Mol Cell. 2005;18(6):735–48.View ArticleGoogle Scholar
- Kaplan N, et al. The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009;458(7236):362–U129.View ArticleGoogle Scholar
- Segal E, et al. A genomic code for nucleosome positioning. Nature. 2006;442(7104):772–8.View ArticleGoogle Scholar
- Zhang Y, et al. Intrinsic histone-DNA interactions are not the major determinant of nucleosome positions in vivo. Nat Struct Mol Biol. 2009;16(8):847–52.View ArticleGoogle Scholar
- Chen K, Luo GZ, He C. High-resolution mapping of N(6)-methyladenosine in transcriptome and genome using a photo-crosslinking-assisted strategy. Methods Enzymol. 2015;560:161–85.View ArticleGoogle Scholar
- Bromberg S, Pratt K, Hattman S. Sequence specificity of DNA adenine methylase in the protozoan Tetrahymena-Thermophila. J Bacteriol. 1982;150(2):993–6.PubMedPubMed CentralGoogle Scholar
- Jiang C, Pugh BF. Nucleosome positioning and gene regulation: advances through genomics. Nat Rev Genet. 2009;10(3):161–72.View ArticleGoogle Scholar
- Zhang Z, Pugh BF. High-resolution genome-wide mapping of the primary structure of chromatin. Cell. 2011;144(2):175–86.View ArticleGoogle Scholar
- Czapla L, Swigon D, Olson WK. Sequence-dependent effects in the cyclization of short DNA. J Chem Theory Comput. 2006;2(3):685–95.View ArticleGoogle Scholar
- Ngo TT, et al. Effects of cytosine modifications on DNA flexibility and nucleosome mechanical stability. Nat Commun. 2016;7:10813.View ArticleGoogle Scholar
- Bujnicki JM, et al. Structure prediction and phylogenetic analysis of a functionally diverse family of proteins homologous to the MT-A70 subunit of the human mRNA:m(6)A methyltransferase. J Mol Evol. 2002;55(4):431–44.View ArticleGoogle Scholar
- Chen K, et al. DANPOS: dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 2013;23(2):341–51.View ArticleGoogle Scholar
- Workman JL. Nucleosome displacement in transcription. Genes Dev. 2006;20(15):2009–17.View ArticleGoogle Scholar
- Bai L, Morozov AV. Gene regulation by nucleosome positioning. Trends Genet. 2010;26(11):476–83.View ArticleGoogle Scholar
- Parfrey LW, et al. Estimating the timing of early eukaryotic diversification with multigene molecular clocks. Proc Natl Acad Sci U S A. 2011;108(33):13624–9.View ArticleGoogle Scholar
- Minge MA, et al. Evolutionary position of breviate amoebae and the primary eukaryote divergence. Proc Biol Sci. 2009;276(1657):597–604.View ArticleGoogle Scholar
- Eisen JA, et al. Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote. PLoS Biol. 2006;4(9):e286.View ArticleGoogle Scholar
- Liu J, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10(2):93–5.View ArticleGoogle Scholar
- Huff JT, Zilberman D. Dnmt1-independent CG methylation contributes to nucleosome positioning in diverse eukaryotes. Cell. 2014;156(6):1286–97.View ArticleGoogle Scholar
- Cassidy-Hanley D, et al. Germline and somatic transformation of mating Tetrahymena thermophila by particle bombardment. Genetics. 1997;146(1):135–47.PubMedPubMed CentralGoogle Scholar
- Zheng G, Lu XJ, Olson WK. Web 3DNA--a web server for the analysis, reconstruction, and visualization of three-dimensional nucleic-acid structures. Nucleic Acids Res. 2009;37(Web Server issue):W240–6.View ArticleGoogle Scholar
- Lu XJ, Olson WK. 3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures. Nat Protoc. 2008;3(7):1213–27.View ArticleGoogle Scholar
- Salomon-Ferrer R, et al. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh Ewald. J Chem Theory Comput. 2013;9(9):3878–88.View ArticleGoogle Scholar
- Ivani I, et al. Parmbsc1: a refined force field for DNA simulations. Nat Methods. 2016;13(1):55–8.View ArticleGoogle Scholar
- Aduri R, et al. AMBER force field parameters for the naturally occurring modified nucleosides in RNA. J Chem Theory Comput. 2007;3(4):1464–75.View ArticleGoogle Scholar
- Price DJ, Brooks CL 3rd. A modified TIP3P water potential for simulation with Ewald summation. J Chem Phys. 2004;121(20):10096–103.View ArticleGoogle Scholar
- Case DA, et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26(16):1668–88.View ArticleGoogle Scholar
- Luo G, et al. N6-methyldeoxyadenosine directs nucleosome positioning in Tetrahymena DNA. Datasets. Gene expression omnibus. 2018. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104699. Accessed 17 Oct 2018.