Comprehensive Methylome Characterization of Mycoplasma genitalium and Mycoplasma pneumoniae at Single-Base Resolution

In the bacterial world, methylation is most commonly associated with restriction-modification systems that provide a defense mechanism against invading foreign genomes. In addition, it is known that methylation plays functionally important roles, including timing of DNA replication, chromosome partitioning, DNA repair, and regulation of gene expression. However, full DNA methylome analyses are scarce due to a lack of a simple methodology for rapid and sensitive detection of common epigenetic marks (ie N⁶-methyladenine (6 mA) and N⁴-methylcytosine (4 mC)), in these organisms. Here, we use Single-Molecule Real-Time (SMRT) sequencing to determine the methylomes of two related human pathogen species, Mycoplasma genitalium G-37 and Mycoplasma pneumoniae M129, with single-base resolution. Our analysis identified two new methylation motifs not previously described in bacteria: a widespread 6 mA methylation motif common to both bacteria (5′-CTAT-3′), as well as a more complex Type I m6A sequence motif in M. pneumoniae (5′-GAN₇TAY-3′/3′-CTN₇ ATR-5′). We identify the methyltransferase responsible for the common motif and suggest the one involved in M. pneumoniae only. Analysis of the distribution of methylation sites across the genome of M. pneumoniae suggests a potential role for methylation in regulating the cell cycle, as well as in regulation of gene expression. To our knowledge, this is one of the first direct methylome profiling studies with single-base resolution from a bacterial organism.     

Differential sensitivity of CG and CCG DNA sequences to ethionine-induced hypomethylation of the Nicotiana tabacum genome

Plant DNA is distinguished from the DNA of all other organisms by its high content of 5-methylcytosine (5mC). 5mC levels may amount to 30% of total cytosines, distributed between the sequences CG and CXG. The results presented here show that the methylation status of CXG sequences could be influenced by culturing tobacco tissues on subtoxic concentrations of ethionine. The hypomethylating effect of ethionine, evaluated as the capability of MspI or HpaII to cleave the DNA, proved to be rather specific for CCG and differed from that or 5-azacytidine which did not discriminate between CG and CXG sequences.     

Bisulfighter: accurate detection of methylated cytosines and differentially methylated regions

Analysis of bisulfite sequencing data usually requires two tasks: to call methylated cytosines (mCs) in a sample, and to detect differentially methylated regions (DMRs) between paired samples. Although numerous tools have been proposed for mC calling, methods for DMR detection have been largely limited. Here, we present Bisulfighter, a new software package for detecting mCs and DMRs from bisulfite sequencing data. Bisulfighter combines the LAST alignment tool for mC calling, and a novel framework for DMR detection based on hidden Markov models (HMMs). Unlike previous attempts that depend on empirical parameters, Bisulfighter can use the expectation-maximization algorithm for HMMs to adjust parameters for each data set. We conduct extensive experiments in which accuracy of mC calling and DMR detection is evaluated on simulated data with various mC contexts, read qualities, sequencing depths and DMR lengths, as well as on real data from a wide range of biological processes. We demonstrate that Bisulfighter consistently achieves better accuracy than other published tools, providing greater sensitivity for mCs with fewer false positives, more precise estimates of mC levels, more exact locations of DMRs and better agreement of DMRs with gene expression and DNase I hypersensitivity. The source code is available at http://epigenome.cbrc.jp/bisulfighter.     

Direct detection of methylation in genomic DNA

The identification of methylated sites on bacterial genomic DNA would be a useful tool to study the major roles of DNA methylation in prokaryotes: distinction of self and nonself DNA, direction of post-replicative mismatch repair, control of DNA replication and cell cycle, and regulation of gene expression. Three types of methylated nucleobases are known: N⁶-methyladenine, 5-methylcytosine and N⁴-methylcytosine. The aim of this study was to develop a method to detect all three types of DNA methylation in complete genomic DNA. It was previously shown that N⁶-methyladenine and 5-methylcytosine in plasmid and viral DNA can be detected by intersequence trace comparison of methylated and unmethylated DNA. We extended this method to include N⁴-methylcytosine detection in both in vitro and in vivo methylated DNA. Furthermore, application of intersequence trace comparison was extended to bacterial genomic DNA. Finally, we present evidence that intrasequence comparison suffices to detect methylated sites in genomic DNA. In conclusion, we present a method to detect all three natural types of DNA methylation in bacterial genomic DNA. This provides the possibility to define the complete methylome of any prokaryote.     

The methylomes of six bacteria

Six bacterial genomes, Geobacter metallireducens GS-15, Chromohalobacter salexigens, Vibrio breoganii 1C-10, Bacillus cereus ATCC 10987, Campylobacter jejuni subsp. jejuni 81-176 and C. jejuni NCTC 11168, all of which had previously been sequenced using other platforms were re-sequenced using single-molecule, real-time (SMRT) sequencing specifically to analyze their methylomes. In every case a number of new N⁶-methyladenine (^m6A) and N⁴-methylcytosine (^m4C) methylation patterns were discovered and the DNA methyltransferases (MTases) responsible for those methylation patterns were assigned. In 15 cases, it was possible to match MTase genes with MTase recognition sequences without further sub-cloning. Two Type I restriction systems required sub-cloning to differentiate their recognition sequences, while four MTase genes that were not expressed in the native organism were sub-cloned to test for viability and recognition sequences. Two of these proved active. No attempt was made to detect 5-methylcytosine (^m5C) recognition motifs from the SMRT® sequencing data because this modification produces weaker signals using current methods. However, all predicted ^m6A and ^m4C MTases were detected unambiguously. This study shows that the addition of SMRT sequencing to traditional sequencing approaches gives a wealth of useful functional information about a genome showing not only which MTase genes are active but also revealing their recognition sequences.     

Examination of the DNA substrate selectivity of DNA cytosine methyltransferases using mass tagging

The biological significance of cytosine methylation is as yet incompletely understood, but substantial and growing evidence strongly suggests that perturbation of methylation patterns, resulting from the infidelity of DNA cytosine methyltransferase, is an important component of the development of human cancer. We have developed a novel in vitro assay that allows us to quantitatively determine the DNA substrate preferences of cytosine methylases. This approach, which we call mass tagging, involves the labeling of target cytosine residues in synthetic DNA duplexes with stable isotopes, such as ¹⁵N. Methylation is then measured by the formation of 5-methylcytosine (5mC) by gas chromatography/mass spectrometry. The DNA substrate selectivity is determined from the mass spectrum of the product 5mC. With the non-symmetrical duplex DNA substrate examined in this study we find that the bacterial methyltransferase HpaII (duplex DNA recognition sequence CCGG) methylates the one methylatable cytosine of each strand similarly. Introduction of an A-C mispair at the methylation site shifts methylation exclusively to the mispaired cytosine residue. In direct competition assays with HpaII methylase we observe that the mispaired substrate is methylated more extensively than the fully complementary, normal substrate, although both have one HpaII methylation site. Through the use of this approach we will be able to learn more about the mechanisms by which methylation patterns can become altered.     

Selective Capture of 5-hydroxymethylcytosine from Genomic DNA

5-methylcytosine (5-mC) constitutes ~2-8% of the total cytosines in human genomic DNA and impacts a broad range of biological functions, including gene expression, maintenance of genome integrity, parental imprinting, X-chromosome inactivation, regulation of development, aging, and cancer¹. Recently, the presence of an oxidized 5-mC, 5-hydroxymethylcytosine (5-hmC), was discovered in mammalian cells, in particular in embryonic stem (ES) cells and neuronal cells^2-4. 5-hmC is generated by oxidation of 5-mC catalyzed by TET family iron (II)/α-ketoglutarate-dependent dioxygenases^{2, 3}. 5-hmC is proposed to be involved in the maintenance of embryonic stem (mES) cell, normal hematopoiesis and malignancies, and zygote development^{2, 5-10}. To better understand the function of 5-hmC, a reliable and straightforward sequencing system is essential. Traditional bisulfite sequencing cannot distinguish 5-hmC from 5-mC¹¹. To unravel the biology of 5-hmC, we have developed a highly efficient and selective chemical approach to label and capture 5-hmC, taking advantage of a bacteriophage enzyme that adds a glucose moiety to 5-hmC specifically¹².Here we describe a straightforward two-step procedure for selective chemical labeling of 5-hmC. In the first labeling step, 5-hmC in genomic DNA is labeled with a 6-azide-glucose catalyzed by β-GT, a glucosyltransferase from T4 bacteriophage, in a way that transfers the 6-azide-glucose to 5-hmC from the modified cofactor, UDP-6-N3-Glc (6-N3UDPG). In the second step, biotinylation, a disulfide biotin linker is attached to the azide group by click chemistry. Both steps are highly specific and efficient, leading to complete labeling regardless of the abundance of 5-hmC in genomic regions and giving extremely low background. Following biotinylation of 5-hmC, the 5-hmC-containing DNA fragments are then selectively captured using streptavidin beads in a density-independent manner. The resulting 5-hmC-enriched DNA fragments could be used for downstream analyses, including next-generation sequencing.Our selective labeling and capture protocol confers high sensitivity, applicable to any source of genomic DNA with variable/diverse 5-hmC abundances. Although the main purpose of this protocol is its downstream application (i.e., next-generation sequencing to map out the 5-hmC distribution in genome), it is compatible with single-molecule, real-time SMRT (DNA) sequencing, which is capable of delivering single-base resolution sequencing of 5-hmC.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.     

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.     

qDNAmod: a statistical model-based tool to reveal intercellular heterogeneity of DNA modification from SMRT sequencing data

In an isogenic cell population, phenotypic heterogeneity among individual cells is common and critical for survival of the population under different environment conditions. DNA modification is an important epigenetic factor that can regulate phenotypic heterogeneity. The single molecule real-time (SMRT) sequencing technology provides a unique platform for detecting a wide range of DNA modifications, including N6-methyladenine (6-mA), N4-methylcytosine (4-mC) and 5-methylcytosine (5-mC). Here we present qDNAmod, a novel bioinformatic tool for genome-wide quantitative profiling of intercellular heterogeneity of DNA modification from SMRT sequencing data. It is capable of estimating proportion of isogenic haploid cells, in which the same loci of the genome are differentially modified. We tested the reliability of qDNAmod with the SMRT sequencing data of Streptococcus pneumoniae strain ST556. qDNAmod detected extensive intercellular heterogeneity of DNA methylation (6-mA) in a clonal population of ST556. Subsequent biochemical analyses revealed that the recognition sequences of two type I restriction–modification (R-M) systems are responsible for the intercellular heterogeneity of DNA methylation initially identified by qDNAmod. qDNAmod thus represents a valuable tool for studying intercellular phenotypic heterogeneity from genome-wide DNA modification.     

Differential sensitivity of CG and CCG DNA sequences to ethionine-induced hypomethylation of the Nicotiana tabacum genome.

Plant DNA is distinguished from the DNA of all other organisms by its high content of 5-methylcytosine (5mC). 5mC levels may amount to 30% of total cytosines, distributed between the sequences CG and CXG. The results presented here show that the methylation status of CXG sequences could be influenced by culturing tobacco tissues on subtoxic concentrations of ethionine. The hypomethylating effect of ethionine, evaluated as the capability of MspI or HpaII to cleave the DNA, proved to be rather specific for CCG and differed from that or 5-azacytidine which did not discriminate between CG and CXG sequences.     

Engineering selectivity of Cutibacterium acnes phages by epigenetic imprinting

Cutibacterium acnes (C. acnes) is a gram-positive bacterium and a member of the human skin microbiome. Despite being the most abundant skin commensal, certain members have been associated with common inflammatory disorders such as acne vulgaris. The availability of the complete genome sequences from various C. acnes clades have enabled the identification of putative methyltransferases, some of them potentially belonging to restriction-modification (R-M) systems which protect the host of invading DNA. However, little is known on whether these systems are functional in the different C. acnes strains. To investigate the activity of these putative R-M and their relevance in host protective mechanisms, we analyzed the methylome of six representative C. acnes strains by Oxford Nanopore Technologies (ONT) sequencing. We detected the presence of a 6-methyladenine modification at a defined DNA consensus sequence in strain KPA171202 and recombinant expression of this R-M system confirmed its methylation activity. Additionally, a R-M knockout mutant verified the loss of methylation properties of the strain. We studied the potential of one C. acnes bacteriophage (PAD20) in killing various C. acnes strains and linked an increase in its specificity to phage DNA methylation acquired upon infection of a methylation competent strain. We demonstrate a therapeutic application of this mechanism where phages propagated in R-M deficient strains selectively kill R-M deficient acne-prone clades while probiotic ones remain resistant to phage infection.     

Methylation of parental and progeny DNA strands in Physarum polycephalum

Although 5-methylcytosine comprises 4 to 8% of the cytosine residues in the major nuclear DNA of Physarum polycephalum (Evans &; Evans, 1970), only 1 % of the cytosine residues of progeny DNA become methylated during replication. Further methylation occurs during the same and subsequent mitotic cycles, so that 6 to 7 cycles after its synthesis, 5-methylcytosine comprises 5 to 7% of the DNA-cytosine residues of a single generation of DNA. The extent of methylation occurring during the S period has been measured by the determination of the specific activity of the precursor (S-adenosylmethionine) and the product (DNA-5-methylcytosine) and by comparison of the radioactivity in DNA-cytosine and DNA-5-methylcytosine after incorporation of [¹⁴C]deoxycytidine. Continuing methylation of parental DNA has been shown, by density shift experiments and by the conversion of prelabeled DNA-cytosine to DNA-5-methylcytosine. The DNA-5-methylcytosine once formed was found to be stable.     

Analysis of TET Expression/Activity and 5mC Oxidation during Normal and Malignant Germ Cell Development

During mammalian development the fertilized zygote and primordial germ cells lose their DNA methylation within one cell cycle leading to the concept of active DNA demethylation. Recent studies identified the TET hydroxylases as key enzymes responsible for active DNA demethylation, catalyzing the oxidation of 5-methylcytosine to 5-hydroxymethylcytosine. Further oxidation and activation of the base excision repair mechanism leads to replacement of a modified cytosine by an unmodified one. In this study, we analyzed the expression/activity of TET1-3 and screened for the presence of 5mC oxidation products in adult human testis and in germ cell cancers. By analyzing human testis sections, we show that levels of 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine are decreasing as spermatogenesis proceeds, while 5-methylcytosine levels remain constant. These data indicate that during spermatogenesis active DNA demethylation becomes downregulated leading to a conservation of the methylation marks in mature sperm. We demonstrate that all carcinoma in situ and the majority of seminomas are hypomethylated and hypohydroxymethylated compared to non-seminomas. Interestingly, 5-formylcytosine and 5-carboxylcytosine were detectable in all germ cell cancer entities analyzed, but levels did not correlate to the 5-methylcytosine or 5-hydroxymethylcytosine status. A meta-analysis of gene expression data of germ cell cancer tissues and corresponding cell lines demonstrates high expression of TET1 and the DNA glycosylase TDG, suggesting that germ cell cancers utilize the oxidation pathway for active DNA demethylation. During xenograft experiments, where seminoma-like TCam-2 cells transit to an embryonal carcinoma-like state DNMT3B and DNMT3L where strongly upregulated, which correlated to increasing 5-methylcytosine levels. Additionally, 5-hydroxymethylcytosine levels were elevated, demonstrating that de novo methylation and active demethylation accompanies this transition process. Finally, mutations of IDH1 (IDH1 ^R132) and IDH2 (IDH2 ^R172) leading to production of the TET inhibiting oncometabolite 2-hydroxyglutarate in germ cell cancer cell lines were not detected.     

Single-base resolution analysis of DNA epigenome via high-throughput sequencing

Epigenetic changes caused by DNA methylation and histone modifications play important roles in the regulation of various cellular processes and development. Recent discoveries of 5-methylcytosine (5mC) oxidation derivatives including 5-hydroxymethylcytosine (5hmC), 5-formylcytsine (5fC) and 5-carboxycytosine (5caC) in mammalian genome further expand our understanding of the epigenetic regulation. Analysis of DNA modification patterns relies increasingly on sequencing-based profiling methods. A number of different approaches have been established to map the DNA epigenomes with single-base resolution, as represented by the bisulfite-based methods, such as classical bisulfite sequencing (BS-seq), TAB-seq (TET-assisted bisulfite sequencing), oxBS-seq (oxidative bisulfite sequencing) and etc. These methods have been used to generate base-resolution maps of 5mC and its oxidation derivatives in genomic samples. The focus of this review will be to discuss the chemical methodologies that have been developed to detect the cytosine derivatives in the genomic DNA.     

GSE Is a Maternal Factor Involved in Active DNA Demethylation in Zygotes

After fertilization, the sperm and oocyte genomes undergo extensive epigenetic reprogramming to form a totipotent zygote. The dynamic epigenetic changes during early embryo development primarily involve DNA methylation and demethylation. We have previously identified Gse (gonad-specific expression gene) to be expressed specifically in germ cells and early embryos. Its encoded protein GSE is predominantly localized in the nuclei of cells from the zygote to blastocyst stages, suggesting possible roles in the epigenetic changes occurring during early embryo development. Here, we report the involvement of GSE in epigenetic reprogramming of the paternal genome during mouse zygote development. Preferential binding of GSE to the paternal chromatin was observed from pronuclear stage 2 (PN2) onward. A knockdown of GSE by antisense RNA in oocytes produced no apparent effect on the first and second cell cycles in preimplantation embryos, but caused a significant reduction in the loss of 5-methylcytosine (5mC) and the accumulation of 5-hydroxymethylcytosine (5hmC) in the paternal pronucleus. Furthermore, DNA methylation levels in CpG sites of LINE1 transposable elements, Lemd1, Nanog and the upstream regulatory region of the Oct4 (also known as Pou5f1) gene were clearly increased in GSE-knockdown zygotes at mid-pronuclear stages (PN3-4), but the imprinted H19-differential methylated region was not affected. Importantly, DNA immunoprecipitation of 5mC and 5hmC also indicates that knockdown of GSE in zygotes resulted in a significant reduction of the conversion of 5mC to 5hmC on LINE1. Therefore, our results suggest an important role of maternal GSE for mediating active DNA demethylation in the zygote.