Construction and use of chimeric SPR/phi 3T DNA methyltransferases in the definition of sequence recognizing enzyme regions.

Multispecific DNA methyltransferases (Mtases) of temperate Bacillus subtilis phages SPR and phi 3T methylate the internal cytosine of the sequence GGCC. They differ in their capacity to methylate additional sequences. These are CCGG and CC(A/T)GG in SPR and GCNGC in phi 3T. Introducing unique restriction sites at equivalent locations within the two genes facilitated the construction of chimeric genes. These expressed Mtase activity at a level comparable to that of the parental genes. The methylation specificity of chimeric enzymes was correlated with the location of chimeric fusions. This analysis, which also included the use of mutant genes, showed that domains involved in the recognition of target sequences unique to each enzyme [CCGG, CC(A/T)GG or GCNGC] are represented by the central non-conserved parts of the proteins, whilst recognition of the sequence (GGCC), which is a target for both enzymes, is determined by an adjacent conserved region.     

In vivo DNA protection by relaxed-specificity SinI DNA methyltransferase variants

The SinI DNA methyltransferase, a component of the SinI restriction-modification system, recognizes the sequence GG(A/T)CC and methylates the inner cytosine to produce 5-methylcytosine. Previously isolated relaxed-specificity mutants of the enzyme also methylate, at a lower rate, GG(G/C)CC sites. In this work we tested the capacity of the mutant enzymes to function in vivo as the counterpart of a restriction endonuclease, which can cleave either site. The viability of Escherichia coli cells carrying recombinant plasmids with the mutant methyltransferase genes and expressing the GGNCC-specific Sau96I restriction endonuclease from a compatible plasmid was investigated. The sau96IR gene on the latter plasmid was transcribed from the araBAD promoter, allowing tightly controlled expression of the endonuclease. In the presence of low concentrations of the inducer arabinose, cells synthesizing the N172S or the V173L mutant enzyme displayed increased plating efficiency relative to cells producing the wild-type methyltransferase, indicating enhanced protection of the cell DNA against the Sau96I endonuclease. Nevertheless, this protection was not sufficient to support long-term survival in the presence of the inducer, which is consistent with incomplete methylation of GG(G/C)CC sites in plasmid DNA purified from the N172S and V173L mutants. Elevated DNA ligase activity was shown to further increase viability of cells producing the V173L variant and Sau96I endonuclease.     

Sequential order of target-recognizing domains in multispecific DNA-methyltransferases.

In the multispecific DNA(cytosine-5)-methyltransferases (Mtases) of Bacillus subtilis phages SPR and phi 3T the domains responsible for recognition of DNA methylation targets CCA/TGG, CCGG, GGCC (SPR) and GCNGC, GGCC (phi 3T) represent contiguous sequences of approximately 50 amino acids each. These domains are tandemly arranged and do not overlap. They are part of a 'variable' segment within the enzymes which is flanked by 'conserved' amino acids, which are very similar amongst bacterial monospecific and the multispecific Mtases studied here. These results follow from a mutational analysis of the SPR and phi 3T Mtase genes. They further support our concept of a modular enzyme organization, according to which variability of type II Mtases with respect to target recognition is achieved by a combination of the same enzyme core with a variety of target-recognizing domains.     

Restriction and modification in Bacillus subtilis: nucleotide sequence, functional organization and product of the DNA methyltransferase gene of bacteriophage SPR

Bacillus subtilis phage SPR codes for a DNA methyltransferase (Mtase) which methylates the 5' cytosine in the sequence GGCC and both cytosines in the sequence CCGG. A 2126-bp fragment of SPR DNA containing the Mtase gene has been sequenced. This fragment has only one significant open reading frame of 1347 bp, which corresponds to the Mtase gene. Within the sequence the Mtase promoter has been defined by S1 mapping. The size of the SPR Mtase predicted from the deduced amino acid composition is 49.9 kDal. This is in agreement with both the Mr of the purified enzyme and with that of the SPR Mtase gene product identified here by minicell technique. Base changes leading to mutants affected in Mtase activity were localized within the Mtase gene.     

Specifically unmethylated cytidylic-guanylate sites in Herpesvirus saimiri DNA in tumor cells.

The restriction endonucleases MspI (CCGG), HpaII (CCGG), FnuDII (CGCG), and HaeIII (GGCC) were used to study the methylation of Herpesvirus saimiri DNA in tumor cells taken directly from tumor-bearing animals. No evidence was found for methylation of the 5' terminal C in the sequence CCGG or of the internal C in the sequence GGCC, but extensive methylation of CG was detected. Fifteen HpaII sites and 17 FnuDII sites were detected in the unique DNA region of the H. saimiri strain used. Twenty-eight of the 32 sites were methylated in greater than 90% of the viral DNA molecules in tumor cells, but the remaining 4 sites were unmethylated in greater than 95% of the viral DNA molecules in tumor cells. The locations of the four specifically unmethylated sites were mapped and appeared to be identical in the four different induced leukemias examined (one owl monkey and three white-lipped marmosets). The nonproducer 1670 tumor cell line, in continuous passage for over 7 years, contained four similar specifically unmethylated sites. Possibilities for the physiological significance of the unmethylated sites are discussed.     

Evolutionary transitions to new DNA methyltransferases through target site expansion and shrinkage

DNA-binding and modifying proteins show high specificity but also exhibit a certain level of promiscuity. Such latent promiscuous activities comprise the starting points for new protein functions, but this hypothesis presents a paradox: a new activity can only evolve if it already exists. How then, do novel activities evolve? DNA methyltransferases, for example, are highly divergent in their target sites, but how transitions toward novel sites occur remains unknown. We performed laboratory evolution of the DNA methyltransferase M.HaeIII. We found that new target sites emerged primarily through expansion of the original site, GGCC, and the subsequent shrinkage of evolved expanded sites. Variants evolved for sites that are promiscuously methylated by M.HaeIII [GG(^A/_T)CC and GGCGCC] carried mutations in ‘gate-keeper’ residues. They could thereby methylate novel target sites such as GCGC and GGATCC that were neither selected for nor present in M.HaeIII. These ‘generalist’ intermediates were further evolved to obtain variants with novel target specificities. Our results demonstrate the ease by which new DNA-binding and modifying specificities evolve and the mechanism by which they occur at both the protein and DNA levels.     

Heterodimeric DNA methyltransferases as a platform for creating designer zinc finger methyltransferases for targeted DNA methylation in cells

The ability to target methylation to specific genomic sites would further the study of DNA methylation’s biological role and potentially offer a tool for silencing gene expression and for treating diseases involving abnormal hypomethylation. The end-to-end fusion of DNA methyltransferases to zinc fingers has been shown to bias methylation to desired regions. However, the strategy is inherently limited because the methyltransferase domain remains active regardless of whether the zinc finger domain is bound at its cognate site and can methylate non-target sites. We demonstrate an alternative strategy in which fragments of a DNA methyltransferase, compromised in their ability to methylate DNA, are fused to two zinc fingers designed to bind 9 bp sites flanking a methylation target site. Using the naturally heterodimeric DNA methyltransferase M.EcoHK31I, which methylates the inner cytosine of 5′-YGGCCR-3′, we demonstrate that this strategy can yield a methyltransferase capable of significant levels of methylation at the target site with undetectable levels of methylation at non-target sites in Escherichia coli. However, some non-target methylation could be detected at higher expression levels of the zinc finger methyltransferase indicating that further improvements will be necessary to attain the desired exclusive target specificity.     

Presence of 5-methylcytosine in CC(A/T)GG sequences (Dcm methylation) in DNAs from different bacteria.

The presence of CC(A/T)GG sequences with methylated internal cytosine (Dcm methylation) was determined in DNA from different genera of eubacteria. This methylation was studied by using restriction enzymes EcoRII and BstNI, which cleave unmethylated or methylated CC(A/T)GG sequences. Dcm methylation was only detected in genera of the family Enterobacteriaceae closely related to Escherichia: Shigella, Citrobacter, Salmonella, and Klebsiella.     

Methylation and expression of a metallothionein promoter ovine growth hormone fusion gene (MToGH1) in transgenic mice

We have examined transgene methylation in the DNA from the livers of a pedigree of mice carrying three copies of an integrated MToGH1 transgene. Utilizing the methylation-sensitive isoschizomersMsp I andHpa II, Southern blot analysis revealed that all second generation animals derived from a transgenic female had hypermethylated DNA, whereas first generation animals sired by a transgenic male displayed a range of methylation phenotypes ranging from no methylation to hypermethylation of the transgene sequences. Of the mice that exhibited hypermethylation of the transgene in CpG dinucleotides (CmCGG), a minority of these animals also exhibited apparent CpC methylation (i.e. inhibition ofMsp I cutting, presumably blocked by methylation of the outer C of CCGG). Methylation was also examined in the inner C of CC(A/T)GG sequences in the MToGH1 transgene using the isoschizomer pairBstN I andEcoR II. A minority of MToGH1 animals in the F₁ generation showed clear evidence of methylation in these sites as well as in the inner and outer Cs of CCGG sites. An examination of MToGH1 expression in terms of oGH levels in serum revealed that there was a high degree of variation in the levels of circulating oGH between animals of this pedigree. There was a weak inverse relationship between the serum level of oGH and the extent of methylation of the transgene. In particular, mice exhibiting CpC together with CpG methylation were found to have very low levels of circulating oGH. Our results highlight the nature and complexity of epigenetic factors associated with transgene sequences which may ultimately influence expression of introduced genes in the mammalian genome.     

Statistical evidence for a biochemical pathway of natural, sequence-targeted G/C to C/G transversion mutagenesis in Haemophilus influenzae Rd.

Markov chain analysis of the Haemophilus influenzae Rd genome reveals striking under-representation of three palindromic tetranucleotide strings (CCGG, GGCC and CATG), accompanied by over-representation of six tetranucleotide strings that are derived from the former by exchanging strand location of the two residues making up a G/C nucleotide pair at the terminal palindrome position. Constraints are outlined for a molecular model able to explain the phenomenon as the result of sequence-targeted, enzyme-driven G/C to C/G transversion mutagenesis. Possible participation in the process by components of known DNA mismatch repair or restriction/modification systems (in particular, cytosine methylation) is discussed. The effect widens the spectrum of enzyme-driven, specific mutagenesis beyond the formerly described C/G to T/A transition (VSP repair of Escherichia coli). Potential evolutionary benefits of enzymatic pathways of specific mutagenesis can be envisioned.     

Exact size and organization of DNA target-recognizing domains of multispecific DNA-(cytosine-C5)-methyltransferases.

A large portion of the sequences of type II DNA-(cytosine-C5)-methyltransferases (C5-MTases) represent highly conserved blocks of amino acids. General steps in the methylation reaction performed by C5-MTases have been found to be mediated by some of these domains. C5-MTases carry, in addition at the same relative location, a region variable in size and amino acid composition, part of which is associated with the capacity of each C5-MTase to recognize its characteristic target. Individual target-recognizing domains (TRDs) for the targets CCGG (M), CC(A/T)GG (E), GGCC (H), GCNGC (F) and G(G/A/T)GC(C/A/T)C (B) could be identified in the C-terminal part of the variable region of multispecific C5-MTases. With experiments reported here, we have established the organization of the variable regions of the multispecific MTases M.SPRI, M.phi3TI, M.H2I and M.rho 11SI at the resolution of individual amino acids. These regions comprise 204, 175, 268 and 268 amino acids, respectively. All variable regions are bipartite. They contain at their N-terminal side a very similar sequence of 71 amino acids. The integrity of this sequence must be assured to provide enzyme activity. Bracketed by 6-10 'linker' amino acids, they have, depending on the enzyme studied, towards their C-terminal end ensembles of individual TRDs of 38 (M), 39 (E), 40 (H), 44 (F) and 54 (B) amino acids. TRDs of different enzymes with equal specificity have the same size. TRDs do not overlap but are either separated by linker amino acids or abut each other.     

Organization of multispecific DNA methyltransferases encoded by temperate Bacillus subtilis phages.

B. subtilis phage rho 11s codes for a multispecific DNA methyltransferase (Mtase) which methylates cytosine within the sequences GGCC and GAGCTC. The Mtase gene of rho 11s was isolated and sequenced. It has 1509 bp, corresponding to 503 amino acids (aa). The enzyme's Mr of 57.2 kd predicted from the nucleotide sequence was verified by direct Mr determinations of the Mtase. A comparison of the aa sequence of the rho 11s Mtase with those of related phages SPR and phi 3%, which differ in their methylation potential, revealed generalities in the building plan of such enzymes. At least 70% of the aa of each enzyme are contained in two regions of 243 and 109 aa at the N and C termini respectively, which are highly conserved among the three enzymes. In each enzyme, variable sequences separate the conserved regions. Variability is generated through the single or multiple use of related and unrelated sequence motifs. We propose that the recognition of those DNA target sequences, which are unique for each of the three enzymes, is determined by these variable regions. Evolutionary relationships between the three enzymes are discussed.     

Multispecific DNA methyltransferases from Bacillus subtilis phages. Properties of wild-type and various mutant enzymes with altered DNA affinity

Temperate Bacillus subtilis phages SPR, phi 3T, rho 11 and SP beta code for DNA methyltransferases, each having multiple sequence specificities. The SPR wild-type and various mutant methyltransferases were overproduced 1000-fold in Escherichia coli and were purified by three consecutive chromatographic steps. The stable form of these multispecific enzymes in solution are monomers with a relative molecular mass (Mr) of about 50,000. The methyl-transfer kinetics of the SPR wild-type and mutant enzymes were determined with DNA substrates carrying either none or one of the three recognition sequences (GGCC, CCGG, CCATGG). Evaluation of the catalytic properties for DNA and S-adenosylmethionine binding suggested that the NH2-terminal part of the protein is important for both non-sequence-specific DNA binding and S-adenosylmethionine binding as well as transfer of methyl groups. On the other hand, mutations in the COOH-terminal part lead to weaker site-specific interactions of the enzyme. Antibodies raised against the purified SPR enzyme specifically immunoprecipitated the phi 3T, rho 11 and SP beta methyltransferases, bu failed to precipitate the chromosomally coded enzymes from B. subtilis (BsuRI) and B. sphaericus (BspRI). Immunoaffinity chromatography is an efficient purification step for the related phage methyltransferases.     

Promiscuous methylation of non-canonical DNA sites by HaeIII methyltransferase

The cytosine C5 methyltransferase M.HaeIII recognises and methylates the central cytosine of its canonical site GGCC. Here we report that M.HaeIII can also, with lower efficiency, methylate cytosines located in a wide range of non-canonical sequences. Using bisulphite sequencing we mapped the methyl- cytosine residues in DNA methylated in vitro and in vivo by M.HaeIII. Methyl-cytosine residues were observed in multiple sequence contexts, most commonly, but not exclusively, at star sites (sites differing by a single base from the canonical sequence). The most frequently used star sites had changes at positions 1 and 4, but there is little or no methylation at star sites changed at position 2. The rate of methylation of non-canonical sites can be quite significant: a DNA substrate lacking a canonical site was methylated by M.HaeIII in vitro at a rate only an order of magnitude slower than an otherwise identical substrate containing the canonical site. In vivo methylation of non-canonical sites may therefore be significant and may have provided the starting point for the evolution of restriction–modification systems with novel sequence specificities.     

Different methylation pattern of melon satellite DNA sequences in hypocotyl and callus tissues

Satellite DNA sequences from Cucumis melo have been examined with respect to modification at CCGG sequences in hypocotyls and in callus tissues. For this purpose, restriction fragments given by HpaII and MspI were compared (both enzymes recognize CCGG sequences but have different sensitivity to methylation at this site). Whereas the methylation level of satellite DNA sequences is on average higher in hypocotyls than in callus tissues, the comparison of partially methylated repeat units of satellite DNA reveals that in callus tissues, all methylated restriction sites are doubly methylated.     

Methylation and repeats in silent and nonsense mutations of p53

All exonic CG sequences in p53 are methylated; this epigenetic modification is correlated with frequent G:C-->A:T transitions in p53. Recent reports reveal the presence in p53 of non-CG methylation in CC and CCC sequences, complementary to sites of selective guanosine adduct formation (GG and GGG), and the association of genetic instability with methylation at repetitive sequences. We presently investigated the distribution of methylation sites and repetitive elements in silent and nonsense p53 mutations (2051) among the IARC's TP53 somatic mutation database for exons 5-8. Silent mutations are nonrandom, but mostly involve G:C-->A:T transitions (62%); in particular C-->T mutations (39% of all silent mutations) are mostly correlated with CC and CCC sequences, while G-->A mutations with GG sequences. Sequence analysis of all non-G:C-->A:T silent mutations reveals the frequent formation of new methylation sites (CG), new CCC and GGG sequences in the resulting sequence, refinement of symmetry elements at interrupted microsatellite-like sequences and formation of small repeats (55.3%). The G:C-->A:T silent mutations characterize cancers associated with cigarette smoking (e.g. bladder or lung and bronchus cancer versus colorectal cancer); on the contrary, non-G:C-->A:T silent mutations have similar frequencies in most cancers. Nonsense mutations in exons 5-8, all resulting in mutants lacking amino acids 307-393, which are crucial for p53 activity, were also analyzed. The frequency of nonsense mutations is higher at methylated sites or repeats 1-2 nucleotides removed from methylation sites. Frameshift mutations are also more frequent at repeated sequences. The frequent G:C-->A:T silent mutations could indicate that CC and CCC sequences of exons 5-8 are occasionally targets of non-CpG methylation of cytosine. This process of de novo methylation in the presence of microsatellite-like sequences and small repeats might influence the genetic stability of a variety of genes.     

DNA methylation in mouse A-repeats in DNA methyltransferase-knockout ES cells and in normal cells determined by bisulfite genomic sequencing

Mouse ES cells with a null mutation of the known DNA methyltransferase retain some residual DNA methylation and can methylate foreign sequences de novo. We have used bisulfite genomic sequencing to examine the sequence specificity and distributions of methylation of a hypermethylated CG island sequence, mouse A-repeats. There were 13 CG dinucleotides in the region examined, 12 of which were methylated to variable extents in all DNAs. We found that: (1) there is considerable residual DNA methylation in ES cells lacking the known DNA methyltransferase (29% of normal methylation in the complete knockout ES DNA); (2) this other activity methylates at exactly the same CG sites as the major methyltransferase; and (3) differences in the distribution of methylated sites between A-repeats in these DNAs are consistent with this other activity methylating in a random de novo fashion. Also, the lack of any methylation in non-CG sites argues that, in other studies where non-CG methylation sites have been found by bisulfite sequencing, detection of such sites of non-CG methylation is not an inherent artifact in this methodology.     

Accuracy of DNA methylation pattern preservation by the Dnmt1 methyltransferase

DNA methyltransferase 1 (Dnmt1) has a central role in copying the pattern of DNA methylation after replication which is one manifestation of epigenetic inheritance. With oligonculeotide substrates we show that mouse Dnmt1 has a 30- to 40-fold preference for hemimethylated DNA that is almost lost after addition of fully methylated oligonucleotides. Using long hemimethylated DNA substrates that carry defined methylation patterns and bisulfite analysis of the methylation reaction products, we show a 15-fold preference for hemimethylated CG sites. Dnmt1 moves along the DNA in a random walk methylating hemimethylated substrates with high processivity (>50 sites are visited on average which corresponds to linear diffusion over 6000 bp). The frequency of skipping sites is very low (<0.3%) and there is no detectable flanking sequence preference. CGCTC sites tend to terminate the processive methylation of DNA by Dnmt1. Unmethylated DNA is modified non-processively with a preference for methylation at CCGG sites. We simulate the propagation of methylation patterns using a stochastic model with the specificity of Dnmt1 observed here and conclude that either methylation of several sites is required to propagate the methylation information over several cellular generations or additional epigenetic information must be used.     

Construction of physical and genetic maps of Chlamydia trachomatis serovar L2 by pulsed-field gel electrophoresis.

We constructed the physical map of Chlamydia trachomatis serovar L2 by using three restriction endonucleases, NotI (GC[GGCCGC), SgrAI (C(A/G)[CCGG(T/G)G), and Sse8387I (CCTGCA[GG), and we analyzed the fragments by pulsed-field gel electrophoresis. A total of 25 restriction endonuclease sites and 13 genes and/or operons were located on the map. The genome size was determined to be 1,045 kb. Neither highly transcribed chlamydia genes nor developmental cycle-specific genes were clustered on the genome.