Genome-Wide Methylation Patterns in Salmonella enterica Subsp. enterica Serovars

The methylation of DNA bases plays an important role in numerous biological processes including development, gene expression, and DNA replication. Salmonella is an important foodborne pathogen, and methylation in Salmonella is implicated in virulence. Using single molecule real-time (SMRT) DNA-sequencing, we sequenced and assembled the complete genomes of eleven Salmonella enterica isolates from nine different serovars, and analysed the whole-genome methylation patterns of each genome. We describe 16 distinct N⁶-methyladenine (m6A) methylated motifs, one N⁴-methylcytosine (m4C) motif, and one combined m6A-m4C motif. Eight of these motifs are novel, i.e., they have not been previously described. We also identified the methyltransferases (MTases) associated with 13 of the motifs. Some motifs are conserved across all Salmonella serovars tested, while others were found only in a subset of serovars. Eight of the nine serovars contained a unique methylated motif that was not found in any other serovar (most of these motifs were part of Type I restriction modification systems), indicating the high diversity of methylation patterns present in Salmonella.     

Base-resolution detection of N4-methylcytosine in genomic DNA using 4mC-Tet-assisted-bisulfite- sequencing

Restriction-modification (R-M) systems pose a major barrier to DNA transformation and genetic engineering of bacterial species. Systematic identification of DNA methylation in R-M systems, including N⁶-methyladenine (6mA), 5-methylcytosine (5mC) and N⁴-methylcytosine (4mC), will enable strategies to make these species genetically tractable. Although single-molecule, real time (SMRT) sequencing technology is capable of detecting 4mC directly for any bacterial species regardless of whether an assembled genome exists or not, it is not as scalable to profiling hundreds to thousands of samples compared with the commonly used next-generation sequencing technologies. Here, we present 4mC-Tet-assisted bisulfite-sequencing (4mC-TAB-seq), a next-generation sequencing method that rapidly and cost efficiently reveals the genome-wide locations of 4mC for bacterial species with an available assembled reference genome. In 4mC-TAB-seq, both cytosines and 5mCs are read out as thymines, whereas only 4mCs are read out as cytosines, revealing their specific positions throughout the genome. We applied 4mC-TAB-seq to study the methylation of a member of the hyperthermophilc genus, Caldicellulosiruptor, in which 4mC-related restriction is a major barrier to DNA transformation from other species. In combination with MethylC-seq, both 4mC- and 5mC-containing motifs are identified which can assist in rapid and efficient genetic engineering of these bacteria in the future.     

Plasmid-controlled variation in the content of methylated bases in single-stranded DNA bacteriophages M13 and fd

The N⁶-methyladenine and 5-methylcytosine contents in the DNA of bacteriophages M13 and fd have been analyzed. The results are summarized as follows. (1) After growth in bacteria harboring the N-3ft^? drug resistance-factor, fd and M13 are observed to contain approximately 1 to 2 more 5-methylcytosine residues per DNA molecule than after growth in the parental drug-sensitive host; no effect on the N⁶-methyladenine content is produced by the plasmid. (2) After growth in bacteria harboring P1 prophage, fd and M13 are observed to contain approximately 2 to 3 more N⁶-methyladenine residues per DNA molecule than after growth in the parental P1-sensitive host; no apparent effect on the 5-methylcytosine content was produced by the P1 plasmid. (3) In agreement with others, fd carrying B-host specificity (fd·B) is observed to contain 2 more N⁶-methyladenine residues/DNA molecule than fd·K.     

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.     

Determination of trace amounts of 5-methylcytosine in DNA by reverse-phase high-performance liquid chromatography

A high-performance liquid chromatographic method to separate five major bases (cytosine, thymine, guanine, adenine, and uracil) and three minor methylated bases (5-methylcytosine, N6-methyladenine, and 7-methylguanine) has been developed using a volatile mobile phase under isocratic conditions. It is extended to quantitate 5-methylcytosine in trace amounts (1 in 20,000 cytosine residues). The suitability of the method has been verified by estimating 5-methylcytosine in DNAs of phi X174 and pBR322. The method has been applied to quantitate the extent of cytosine methylation in DNA of larval silk glands of Bombyx mori. Our results confirm that the pupal DNA of Drosophila melanogaster does not contain detectable amounts of 5-methylcytosine.     

Bioinformatics and functional analysis define four distinct groups of AlkB DNA-dioxygenases in bacteria

The iron(II)- and 2-oxoglutarate (2OG)-dependent dioxygenase AlkB from Escherichia coli (EcAlkB) repairs alkylation damage in DNA by direct reversal. EcAlkB substrates include methylated bases, such as 1-methyladenine (m¹A) and 3-methylcytosine (m³C), as well as certain bulkier lesions, for example the exocyclic adduct 1,N⁶-ethenoadenine (εA). EcAlkB is the only bacterial AlkB protein characterized to date, and we here present an extensive bioinformatics and functional analysis of bacterial AlkB proteins. Based on sequence phylogeny, we show that these proteins can be subdivided into four groups: denoted 1A, 1B, 2A and 2B; each characterized by the presence of specific conserved amino acid residues in the putative nucleotide-recognizing domain. A scattered distribution of AlkB proteins from the four different groups across the bacterial kingdom indicates a substantial degree of horizontal transfer of AlkB genes. DNA repair activity was associated with all tested recombinant AlkB proteins. Notably, both a group 2B protein from Xanthomonas campestris and a group 2A protein from Rhizobium etli repaired etheno adducts, but had negligible activity on methylated bases. Our data indicate that the majority, if not all, of the bacterial AlkB proteins are DNA repair enzymes, and that some of these proteins do not primarily target methylated bases.     

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.     

Effect of enzymatic methylation on the chemical,physical, and biological properties of DNA

DNA methylase was partially purified from Escherichia coli W and used to methylate DNA from Bacillus subtilis and bacteriophage φ105. The former DNA was methylated 1.17% and the latter 0.87%. The products were 6-methyladenine (85%) and 5-methylcytosine (15%) in both cases. The methylated DNA was stable toward depurination and viscosity loss at elevated temperatures. Methylation led to a 50% decrease in transforming activity in two strains of B. subtilis and no change in a third strain. The ability of phage φ105 DNA to rescue a defective phage strain was decreased 50% by methylation. No changes were observed in the ability of methylated DNA to serve as a template for DNA polymerase or RNA polymerase. The pattern of cleavage of DNA by a variety of restriction endonucleases was not affected by methylation. There were no changes in the physicochemical properties of DNA on methylation as measured by hyperchromicity on heating, formaldehyde denaturation, viscosity, and sedimentation.     

DNA methylation in mycobacteria: Absence of methylation at GATC (Dam) and CCA/TGG (Dcm) sequences

Abstract The presence of 6-methyladenine and 5-methylcytosine at Dam (GATC) and Dcm (CCA/TGG) sites in DNA of mycobacterial species was investigated using isoschizomer restriction enzymes. In all species examined, Dam and Dcm recognition sequences were not methylated indicating the absence of these methyltransferases. On the other hand, high performance liquid chromatographic analysis of genomic DNA from Mycobacterium smegmatis and Mycobacterium tuberculosis showed significant levels of 6-methyladenine and 5-methylcytosine suggesting the presence of DNA methyltransferases other than Dam and Dcm. Occurrence of methylation was also established by a sensitive genetic assay.     

Characterization of DNA methyltransferase specificities using single-molecule, real-time DNA sequencing

DNA methylation is the most common form of DNA modification in prokaryotic and eukaryotic genomes. We have applied the method of single-molecule, real-time (SMRT®) DNA sequencing that is capable of direct detection of modified bases at single-nucleotide resolution to characterize the specificity of several bacterial DNA methyltransferases (MTases). In addition to previously described SMRT sequencing of N6-methyladenine and 5-methylcytosine, we show that N4-methylcytosine also has a specific kinetic signature and is therefore identifiable using this approach. We demonstrate for all three prokaryotic methylation types that SMRT sequencing confirms the identity and position of the methylated base in cases where the MTase specificity was previously established by other methods. We then applied the method to determine the sequence context and methylated base identity for three MTases with unknown specificities. In addition, we also find evidence of unanticipated MTase promiscuity with some enzymes apparently also modifying sequences that are related, but not identical, to the cognate site.     

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.     

DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates

The Gram-negative bacterium Neisseria meningitidis features extensive genetic variability. To present, proposed virulence genotypes are also detected in isolates from asymptomatic carriers, indicating more complex mechanisms underlying variable colonization modes of N. meningitidis.We applied the Single Molecule, Real-Time (SMRT) sequencing method from Pacific Biosciences to assess the genome-wide DNA modification profiles of two genetically related N. meningitidis strains, both of serogroup A. The resulting DNA methylomes revealed clear divergences, represented by the detection of shared and of strain-specific DNA methylation target motifs. The positional distribution of these methylated target sites within the genomic sequences displayed clear biases, which suggest a functional role of DNA methylation related to the regulation of genes.DNA methylation in N. meningitidis has a likely underestimated potential for variability, as evidenced by a careful analysis of the ORF status of a panel of confirmed and predicted DNA methyltransferase genes in an extended collection of N. meningitidis strains of serogroup A. Based on high coverage short sequence reads, we find phase variability as a major contributor to the variability in DNA methylation. Taking into account the phase variable loci, the inferred functional status of DNA methyltransferase genes matched the observed methylation profiles.Towards an elucidation of presently incompletely characterized functional consequences of DNA methylation in N. meningitidis, we reveal a prominent colocalization of methylated bases with Single Nucleotide Polymorphisms (SNPs) detected within our genomic sequence collection. As a novel observation we report increased mutability also at 6mA methylated nucleotides, complementing mutational hotspots previously described at 5mC methylated nucleotides.These findings suggest a more diverse role of DNA methylation and Restriction-Modification (RM) systems in the evolution of prokaryotic genomes.     

Methylation of adenine in the nuclear DNA of Tetrahymena is internucleosomal and independent of histone H1

There are about 50 copies of each chromosome in the somatic macronucleus of the ciliated protozoan Tetrahymena. Approximately 0.8% of the adenine residues in the macronuclear DNA of Tetrahymena are methylated to N⁶-methyladenine. The degree of methylation varies between sites from a very low percentage to >90%. In this study a correlation was found between nucleosome positioning and DNA methylation. Eight GATC sites with different levels of methylation were examined. There was a direct correlation between the degree of methylation and proximity to linker DNA at these sites. Although methylation occurs preferentially in linker DNA, the patterns and extent of methylation in a histone H1 knockout strain were virtually indistinguishable from those in wild-type cells.     

The methylation of tRNA in KB cells after treatment with actinomycin D

Extensive shifts in the distribution of labeled methylated constituents of tRNA were observed in KB cells treated with actinomycin D for 30 min prior to a 90-min pulse with ³H-CH₃-methionine. Although this treatment completely blocked the synthesis of tRNA, methylation continued to the extent of 12–15% of controls (pulsed without antibiotic). Under this condition the relative proportion of radioactivity incorporated into 3-methylcytosine, N²-methylguanine and 2′-O-methylribose was markedly increased (170–235% of control values), it was moderately reduced in 1-methyladenine, 5-methyl-cytosine and 5-methyluracil (35–70% of controls) and markedly reduced in 1-, 7- and N²N²-methylguanines (15–30% of controls). These data suggest that specific types of methylations occur at particular times during the processing of pre-tRNAs.     

Differential methylation of mitochondrial and nuclear DNA in cultured mouse, hamster and virus-transformed hamster cells. In vivo and in vitro methylation

Information has been lacking as to whether mitochondrial DNA of animal cells is methylated. The methylation patterns of mitochondrial and nuclear DNAs of several mammalian cell lines have therefore been compared by four methods: (1) in vivo transfer of the methyl group from [methyl-³H]methionine; (2) in vivo incorporation of [³²P]orthophosphate and a combination of (1) and (2); (3) in vivo incorporation of [³H]deoxycytidine; (4) in vitro methylation of DNAs with ³H-labeled S-adenosylmethionine as methyl donor and DNA methylase preparations from L cell nuclei. The cell lines were mouse L cells, $\text{BHK21}\text{C13}$, $\text{C13}\text{B4}$ (baby hamster kidney cells transformed by the Bryan strain of Rouse sarcoma virus), and PyY (BHK cells transformed by polyoma virus). DNA bases were separated chromatographically, using 5-methylcytosine, 6-methylaminopurine and, in some cases, 7-methylguanine as markers.Mitochondrial DNA was found to be significantly less methylated than nuclear DNA with respect to 5-methylcytosine in all cell types studied and by all methods used. The relative advantages and disadvantages of each method have been discussed. The level of 5-methylcytosine in mitochondrial DNA as compared with that in nuclear DNA was estimated as one-fourth to one-fourteenth in various cell lines. The estimated 5-methylcytosine content per circular mitochondrial DNA molecule (mol. wt 10 × 10⁶) was about 12 methylcytosine residues for L cells and 24, 30 and 36 methylcytosine residues for BHK, B4 and PyY cells, respectively. Relative to cytosine residues, the estimate was one 5-methylcytosine per 500 cytosine residues of mitochondrial DNA and one 5-methylcytosine per 36 cytosine residues of nuclear DNA from L-cells. The values for methylcytosine of mitochondrial DNA are presumed to be maximal. PyY cells as compared with other cells had the highest methylcytosine content of both mitochondrial and nuclear DNA as estimated by method (3). No methylation of nuclear DNA was observed in confluent L cells.Evidence for the presence of DNA methylase activity associated with mitochondrial fractions was obtained. This activity could be distinguished from other cellular DNA methylase activity by differential response to mercaptoethanol. Radioactivity from ³H-labeled S-adenosylmethionine was found only in 5-methyl-cytosine of DNA.     

Depletion of CG-Specific Methylation in Mycoplasma hyorhinis Genomic DNA after Host Cell Invasion

Adaptation to the environment requires pathogenic bacteria to alter their gene expression in order to increase long-term survival in the host. Here, we present the first experimental evidence that bacterial DNA methylation affects the intracellular survival of pathogenic Mycoplasma hyorhinis. Using bisulfite sequencing, we identified that the M. hyorhinis DNA methylation landscape was distinct in free-living M. hyorhinis relative to the internalized bacteria surviving in the infected human cells. We determined that genomic GATC sites were consistently highly methylated in the bacterial chromosome suggesting that the bacterial GATC-specific 5-methylcytosine DNA methyltransferase was fully functional both pre- and post-infection. In contrast, only the low CG methylation pattern was observed in the mycoplasma genome in the infective bacteria that invaded and then survived in the host cells. In turn, two distinct populations, with either high or low CG methylation, were detected in the M. hyorhinis cultures continually grown in the rich medium independently of host cells. We also identified that M. hyorhinis efficiently evaded endosomal degradation and uses exocytosis to exit infected human cells enabling re-infection of additional cells. The well-orchestrated changes in the chromosome methylation landscape play a major regulatory role in the mycoplasma life cycle.     

Host specificity of DNA produced by Escherichia coli : XV. The role of nucleotide methylation in in vitro B-specific modification

It is shown that in vitro Escherichia coli strain B-specific modification of the replicative form of bacteriophage fd DNA is accompanied by the methylation of certain adenine moieties to form N-6-methyladenine. The reaction follows first order kinetics and saturation is reached when about four adenines are methylated per replicative form. No methyl groups are transferred to B-modified DNA.     

Sequence specificity of DNA methylases from Bacillus amyloliquefaciens and Bacillus brevis

DNA methylation in Bacillus amyloliquefaciens strain H (Bam)² and Bacillus brevis (Bbv) has been examined by a variety of techniques. In vivo labelling studies revealed that Bam DNA contains no N⁶-methyladenine (MeAde), but contains 5-methylcytosine (MeCyt); approximately 0·7% of the cytosine residues are methylated.DNA methylase activity was partially purified from both Bam and Bbv; the Bam enzyme preparation transferred methyl groups from S-adenosyl-l-[methyl-³H]methionine ([³H]AdoMet) to specific DNA cytosine residues only; in agreement with Vanyushin & Dobritsa (1975), the Bbv enzyme preparation methylated both DNA adenine and cytosine residues. The (partial) sequence specificity of the methylases was determined by analyzing [³H]methyl-labelled dinucleotides obtained from enzymatic digests of DNA methylated in vitro. Bam and Bbv each contain a DNA-cytosine methylase with overlapping sequence specificity; e.g. both enzymes produce G-C1, C1-A and C1-T. This is consistent with a single, twofold symmetrical methylation sequence of 5′ … G-C1-(A or T)-G-C … 3′; this was observed by Vanyushin & Dobritsa (1975) for a different Bbv strain. Bam contains a second DNA-cytosine methylase (not present in Bbv), which produces T-C1 and C1-T. We propose that this methylase is the BamI modification enzyme, and that the modified sequence is 5′ … G-G-A-T-C1-C … 3′.Bbv appears to contain two DNA-adenine methylases which produce the (partial) methylated sequences, 5′ … G-A1-T … 3′ and 5′ … A-A1-G … 3′, respectively; in the former case, all the G-A-T-C sites on Bbv DNA appear to be methylated.     

Rational development of transformation in Clostridium thermocellum ATCC 27405 via complete methylome analysis and evasion of native restriction–modification systems

A major barrier to both metabolic engineering and fundamental biological studies is the lack of genetic tools in most microorganisms. One example is Clostridium thermocellum ATCC 27405^T, where genetic tools are not available to help validate decades of hypotheses. A significant barrier to DNA transformation is restriction–modification systems, which defend against foreign DNA methylated differently than the host. To determine the active restriction–modification systems in this strain, we performed complete methylome analysis via single-molecule, real-time sequencing to detect 6-methyladenine and 4-methylcytosine and the rarely used whole-genome bisulfite sequencing to detect 5-methylcytosine. Multiple active systems were identified, and corresponding DNA methyltransferases were expressed from the Escherichia coli chromosome to mimic the C. thermocellum methylome. Plasmid methylation was experimentally validated and successfully electroporated into C. thermocellum ATCC 27405. This combined approach enabled genetic modification of the C. thermocellum-type strain and acts as a blueprint for transformation of other non-model microorganisms.

     

Hydroquinone Increases 5-Hydroxymethylcytosine Formation through Ten Eleven Translocation 1 (TET1) 5-Methylcytosine Dioxygenase

DNA methylation regulates gene expression throughout development and in a wide range of pathologies such as cancer and neurological disorders. Pathways controlling the dynamic levels and targets of methylation are known to be disrupted by chemicals and are therefore of great interest in both prevention and clinical contexts. Benzene and its metabolite hydroquinone have been shown to lead to decreased levels of DNA methylation, although the mechanism is not known. This study employs a cell culture model to investigate the mechanism of hydroquinone-mediated changes in DNA methylation. Exposures that do not affect HEK293 cell viability led to genomic and methylated reporter DNA demethylation. Hydroquinone caused reactivation of a methylated reporter plasmid that was prevented by the addition of N-acetylcysteine. Hydroquinone also caused an increase in Ten Eleven Translocation 1 activity and global levels of 5-hydroxymethylcytosine. 5-Hydroxymethylcytosine was found enriched at LINE-1 prior to a decrease in both 5-hydroxymethylcytosine and 5-methylcytosine. Ten Eleven Translocation-1 knockdown decreased 5-hydroxymethylcytosine formation following hydroquinone exposure as well as the induction of glutamate-cysteine ligase catalytic subunit and 14-3-3σ. Finally, Ten Eleven Translocation 1 knockdown decreased the percentage of cells accumulating in G₂+M following hydroquinone exposure, indicating that it may have a role in cell cycle changes in response to toxicants. This work demonstrates that hydroquinone exposure leads to active and functional DNA demethylation in HEK293 cells in a mechanism involving reactive oxygen species and Ten Eleven Translocation 1 5-methylcytosine dioxygenase.