Cleavage of methylated CCCGGG sequences containing either N4-methylcytosine or 5-methylcytosine with MspI, HpaII, SmaI, XmaI and Cfr9I restriction endonucleases

Methods for analysis of specific DNA methylation status

Methylation of CpG dinucleotides plays a crucial role in the regulation of gene expression and therefore in the development of different pathologies. Aberrant methylation has been associated to the majority of the diseases, including cancer, neurodegenerative, cardiovascular and autoimmune disorders. Analysis of DNA methylation patterns is crucial to understand the underlying molecular mechanism of these diseases. Moreover, DNA methylation patterns could be used as biomarker for clinical management, such as diagnosis, prognosis and treatment response. Nowadays, a variety of high throughput methods for DNA methylation have been developed to analyze the methylation status of a high number of CpGs at once or even the whole genome. However, identification of specific methylation patterns at specific loci is essential for validation and also as a tool for diagnosis. In this review, we describe the most commonly used approaches to evaluate specific DNA methylation. There are three main groups of techniques that allow the identification of specific regions that are differentially methylated: bisulfite conversion-based methods, restriction enzyme-based approaches, and affinity enrichment-based assays. In the first group, specific restriction enzymes recognize and cleave unmethylated DNA, leaving methylated sequences intact. Bisulfite conversion methods are the most popular approach to distinguish methylated and unmethylated DNA. Unmethylated cytosines are deaminated to uracil by sodium bisulfite treatment, while the methyl cytosines remain unconverted. In the last group, proteins with methylation binding domains or antibodies against methyl cytosines are used to recognize methylated DNA. In this review, we provide the theoretical basis and the framework of each technique as well as the analysis of their strength and the weaknesses.

Transient N4 -Acyl-DNA Protection against Cleavage by Restriction Endonucleases

A set of five N⁴ -acyl-modified 2'-deoxycytidine 5'-triphosphates were incorporated into modified DNA by using phi29 DNA polymerase, and cleavage by selected restriction endonucleases was studied. Modified DNA containing N⁴ -acyl functional groups in either one or both strands of a DNA molecule was resistant to the majority of restriction enzymes tested, whereas modifications outside of the recognition sequences were well tolerated. The N⁴ -acylated cytidine derivatives were subjected to competitive nucleotide incorporation by using phi29 DNA polymerase, showing that a high-fidelity phi29 DNA polymerase efficiently used the modified analogues in the presence of its natural counterpart. These N⁴ modifications were also demonstrated to be easily removed in an aqueous ethanolamine solution, in which all steps, including primer extension, demodification, and cleavage by restriction endonuclease, could be performed in a one-pot procedure that eliminated additional purification stages. It is suggested that N⁴ -modified nucleotides are promising building blocks for a programmable; transient; and, most importantly, straightforward DNA protection against specific endonucleases.

Rapid and affordable genome-wide bisulfite DNA sequencing by XmaI-reduced representation bisulfite sequencing

Aim: To develop a reduced representation bisulfite sequencing (RRBS) approach for rapid and affordable genome-wide DNA methylation analysis. Methods: We have selected restriction endonuclease XmaI to produce RRBS library fragments. After digestion and partial fill-in DNA fragments were ligated to barcoded adapters, bisulfite converted, size-selected, and sequenced on the Ion Torrent Personal Genome Machine. XmaI-RRBS results were compared with the previously published RRBS data. Results: We have developed an XmaI-RRBS method for rapid and affordable genome-wide DNA methylation analysis, with library preparation taking only 4 days and sequencing possible within 4 h. We have also addressed several challenges in order to further improve the RRBS technology. XmaI-RRBS may be performed on degraded DNA samples and is compatible with the bench-top next-generation sequencing machines.

Two Inducible Prophages of an Antarctic Pseudomonas sp. ANT_H14 Use the Same Capsid for Packaging Their Genomes - Characterization of a Novel Phage Helper-Satellite System

Two novel prophages ФAH14a and ФAH14b of a psychrotolerant Antarctic bacterium Pseudomonas sp. ANT_H14 have been characterized. They were simultaneously induced with mitomycin C and packed into capsids of the same size and protein composition. The genome sequences of ФAH14a and ФAH14b have been determined. ФAH14b, the phage with a smaller genome (16,812 bp) seems to parasitize ФAH14a (55,060 bp) and utilizes its capsids, as only the latter encodes a complete set of structural proteins. Both viruses probably constitute a phage helper-satellite system, analogous to the P2-P4 duo. This study describes the architecture and function of the ФAH14a and ФAH14b genomes. Moreover, a functional analysis of a ФAH14a-encoded lytic enzyme and a DNA methyltransferase was performed. In silico analysis revealed the presence of the homologs of ФAH14a and ФAH14b in other Pseudomonas genomes, which may suggest that helper-satellite systems related to the one described in this work are common in pseudomonads.

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.

Epigenetic variability in the genetically uniform forest tree species Pinus pinea L

There is an increasing interest in understanding the role of epigenetic variability in forest species and how it may contribute to their rapid adaptation to changing environments. In this study we have conducted a genome-wide analysis of cytosine methylation pattern in Pinus pinea, a species characterized by very low levels of genetic variation and a remarkable degree of phenotypic plasticity. DNA methylation profiles of different vegetatively propagated trees from representative natural Spanish populations of P. pinea were analyzed with the Methylation Sensitive Amplified Polymorphism (MSAP) technique. A high degree of cytosine methylation was detected (64.36% of all scored DNA fragments). Furthermore, high levels of epigenetic variation were observed among the studied individuals. This high epigenetic variation found in P. pinea contrasted with the lack of genetic variation based on Amplified Fragment Length Polymorphism (AFLP) data. In this manner, variable epigenetic markers clearly discriminate individuals and differentiates two well represented populations while the lack of genetic variation revealed with the AFLP markers fail to differentiate at both, individual or population levels. In addition, the use of different replicated trees allowed identifying common polymorphic methylation sensitive MSAP markers among replicates of a given propagated tree. This set of MSAPs allowed discrimination of the 70% of the analyzed trees.

How to interpret methylation sensitive amplified polymorphism (MSAP) profiles?

BACKGROUND

DNA methylation plays a key role in development, contributes to genome stability, and may also respond to external factors supporting adaptation and evolution. To connect different types of stimuli with particular biological processes, identifying genome regions with altered 5-methylcytosine distribution at a genome-wide scale is important. Many researchers are using the simple, reliable, and relatively inexpensive Methylation Sensitive Amplified Polymorphism (MSAP) method that is particularly useful in studies of epigenetic variation. However, electrophoretic patterns produced by the method are rather difficult to interpret, particularly when MspI and HpaII isoschizomers are used because these enzymes are methylation-sensitive, and any C within the CCGG recognition motif can be methylated in plant DNA.

RESULTS

Here, we evaluate MSAP patterns with respect to current knowledge of the enzyme activities and the level and distribution of 5-methylcytosine in plant and vertebrate genomes. We discuss potential caveats related to complex MSAP patterns and provide clues regarding how to interpret them. We further show that addition of combined HpaII + MspI digestion would assist in the interpretation of the most controversial MSAP pattern represented by the signal in the HpaII but not in the MspI profile.

CONCLUSIONS

We recommend modification of the MSAP protocol that definitely discerns between putative hemimethylated mCCGG and internal CmCGG sites. We believe that our view and the simple improvement will assist in correct MSAP data interpretation.

Scoring and analysis of methylation-sensitive amplification polymorphisms for epigenetic population studies

DNA methylation is an important, heritable epigenetic modification in most eukaryotic organisms that is connected with numerous biological processes. To study the impact of natural epigenetic variation in an ecological or evolutionary context, epigenetic studies are increasingly using methylation-sensitive amplification polymorphism (MSAP) for surveys at the population or species level. However, no consensus exists on how to interpret and score the multistate information obtained from the MSAP banding patterns. Here, we review the previously used scoring approaches for population epigenetic studies and develop new alternatives. To assess effects of the different approaches on parameters of epigenetic diversity and differentiation, we applied eight scoring schemes to a case study of three populations of the plant species Viola elatior. For a total number of 168 detected polymorphic MSAP fragments, the number of ultimately scored polymorphic epiloci ranged between 78 and 286 depending on the particular scoring scheme. Both, estimates of epigenetic diversity and differentiation varied strongly between scoring approaches. However, linear regression and PCoA revealed qualitatively similar patterns, suggesting that the scoring approaches are largely consistent. For single-locus analyses of MSAP data, for example the search for loci under selection, we advocate a new scoring approach that separately takes into account different methylation types and thus seems appropriate for drawing more detailed conclusions in ecological or evolutionary contexts. An R script (MSAP_score.r) for scoring and basic data analysis is provided.

High sensitivity 5-hydroxymethylcytosine detection in Balb/C brain tissue

DNA hydroxymethylation is a long known modification of DNA, but has recently become a focus in epigenetic research. Mammalian DNA is enzymatically modified at the 5(th) carbon position of cytosine (C) residues to 5-mC, predominately in the context of CpG dinucleotides. 5-mC is amenable to enzymatic oxidation to 5-hmC by the Tet family of enzymes, which are believed to be involved in development and disease. Currently, the biological role of 5-hmC is not fully understood, but is generating a lot of interest due to its potential as a biomarker. This is due to several groundbreaking studies identifying 5-hydroxymethylcytosine in mouse embryonic stem (ES) and neuronal cells. Research techniques, including bisulfite sequencing methods, are unable to easily distinguish between 5-mC and 5-hmC . A few protocols exist that can measure global amounts of 5-hydroxymethylcytosine in the genome, including liquid chromatography coupled with mass spectrometry analysis or thin layer chromatography of single nucleosides digested from genomic DNA. Antibodies that target 5-hydroxymethylcytosine also exist, which can be used for dot blot analysis, immunofluorescence, or precipitation of hydroxymethylated DNA, but these antibodies do not have single base resolution.In addition, resolution depends on the size of the immunoprecipitated DNA and for microarray experiments, depends on probe design. Since it is unknown exactly where 5-hydroxymethylcytosine exists in the genome or its role in epigenetic regulation, new techniques are required that can identify locus specific hydroxymethylation. The EpiMark 5-hmC and 5-mC Analysis Kit provides a solution for distinguishing between these two modifications at specific loci. The EpiMark 5-hmC and 5-mC Analysis Kit is a simple and robust method for the identification and quantitation of 5-methylcytosine and 5-hydroxymethylcytosine within a specific DNA locus. This enzymatic approach utilizes the differential methylation sensitivity of the isoschizomers MspI and HpaII in a simple 3-step protocol. Genomic DNA of interest is treated with T4-BGT, adding a glucose moeity to 5-hydroxymethylcytosine. This reaction is sequence-independent, therefore all 5-hmC will be glucosylated; unmodified or 5-mC containing DNA will not be affected. This glucosylation is then followed by restriction endonuclease digestion. MspI and HpaII recognize the same sequence (CCGG) but are sensitive to different methylation states. HpaII cleaves only a completely unmodified site: any modification (5-mC, 5-hmC or 5-ghmC) at either cytosine blocks cleavage. MspI recognizes and cleaves 5-mC and 5-hmC, but not 5-ghmC. The third part of the protocol is interrogation of the locus by PCR. As little as 20 ng of input DNA can be used. Amplification of the experimental (glucosylated and digested) and control (mock glucosylated and digested) target DNA with primers flanking a CCGG site of interest (100-200 bp) is performed. If the CpG site contains 5-hydroxymethylcytosine, a band is detected after glucosylation and digestion, but not in the non-glucosylated control reaction. Real time PCR will give an approximation of how much hydroxymethylcytosine is in this particular site. In this experiment, we will analyze the 5-hydroxymethylcytosine amount in a mouse Babl/C brain sample by end point PCR.

Transformation of cytosine to uracil in single-stranded DNA via their oxime sulfonates

The novel transformation of cytosine oxime sulfonate (C*) to uracil sulfonate (U*) proceeded by a Co(II)-assisted benzoyl peroxide reaction, eventually leading to the conversion of cytosine to uracil.

Chemical approach toward efficient DNA methylation analysis

The development of a reaction for the detection of the presence/absence of one methyl group in a very long DNA strand is a chemically and biologically challenging research subject. Several newly designed chemical assays for the typing of DNA methylation are reported and discussed in this paper. A new concept of sequence-specific short-term methylation analysis, supported by a chemical basis, is the starting point for a novel methylation-typing assay, which will supersede conventional methods.

DNA methylase activity as a marker for the presence of a family of phage-like elements conferring efflux-mediated macrolide resistance in streptococci

Recently, two related chimeric genetic elements (Tn1207.3 and Phi10394.4) were shown to carry the macrolide efflux gene mef in Streptococcus pyogenes (group A streptococci [GAS]). The dissemination of elements belonging to the Tn1207.3/Phi10394.4 family in recent isolates of GAS, Streptococcus dysgalactiae subsp. equisimilis, Streptococcus pneumoniae, and Streptococcus agalactiae recovered in Portugal was surveyed. In total, 149 GAS, 18 S. pneumoniae, 4 S. dysgalactiae subsp. equisimilis, and 5 S. agalactiae isolates from infections, presenting the M phenotype of macrolide resistance and containing the mef gene, were screened for the presence of Tn1207.3/Phi10394.4 by PCR targeting open reading frames (ORFs) specific for these related elements. All the GAS isolates tested and one of the S. dysgalactiae subsp. equisimilis isolates carried Tn1207.3. However, neither of these elements was found in the isolates of the other streptococcal species. It was also noted that the DNAs of the isolates carrying Tn1207.3 were resistant to cleavage by the endonuclease SmaI. Cloning and expression of ORF12 of Tn1207.3 in Escherichia coli showed that it encoded a methyltransferase that rendered DNA refractory to cleavage by SmaI (M.Spy10394I). Using this characteristic as a marker for the presence of the Tn1207.3/Phi10394.4 family, we reviewed the literature and concluded that these genetic elements are widely distributed among tetracycline-susceptible GAS isolates presenting the M phenotype from diverse geographic origins and may have played an important role in the dissemination of macrolide resistance in this species.

TspMI, a thermostable isoschizomer of XmaI (5′C/CCGGG3′): characterization and single molecule imaging with DNA

TspMI, a thermostable isoschizomer of XmaI from a Thermus sp., has been characterized. The enzyme was purified to homogeneity using Cibacron-Blue 3GA agarose, Heparin agarose, SP sephadex C50, and Mono-Q fast protein liquid chromatography and was found to be a homodimer of 40 kDa. Restriction mapping and run-off sequencing of TspMI-cleaved DNA ends depicted that it cleaved at 5'C/CCGGG3' to generate a four-base, 5'-CCGG overhang. The enzyme was sensitive to methylation of second and third cytosines in its recognition sequence. TspMI worked optimally at 60 degrees C with 6 mM Mg(2+), no Na(+)/K(+), and showed no star activity in the presence of 25% glycerol. The enzyme could efficiently digest the DNA labeled with a higher concentration of YOYO-I (one dye molecule to one nucleotide), making it a useful candidate for real-time imaging experiments. Single molecule interaction between TspMI and lambda DNA was studied using total internal reflection fluorescence microscopy. The enzyme survived 30 polymerase chain reaction (PCR) cycles in the presence of 10% glycerol and 0.5 M trehalose without any activity loss and, hence, is suitable for incorporation in restriction-endonuclease-mediated selective-PCR for various applications.

Sulfolobus tengchongensis spindle-shaped virus STSV1: virus-host interactions and genomic features

A virus infecting the hyperthermophilic archaeon Sulfolobus tengchongensis has been isolated from a field sample from Tengchong, China, and characterized. The virus, denoted STSV1 (Sulfolobus tengchongensis spindle-shaped virus 1), has the morphology of a spindle (230 by 107 nm) with a tail of variable length (68 nm on average) at one end and is the largest of the known spindle-shaped viruses. After infecting its host, the virus multiplied rapidly to high titers (>10(10) PFU/ml). Replication of the virus retarded host growth but did not cause lysis of the host cells. STSV1 did not integrate into the host chromosome and existed in a carrier state. The STSV1 DNA was modified in an unusual fashion, presumably by virally encoded modification systems. STSV1 harbors a double-stranded DNA genome of 75,294 bp, which shares no significant sequence similarity to those of fuselloviruses. The viral genome contains a total of 74 open reading frames (ORFs), among which 14 have a putative function. Five ORFs encode viral structural proteins, including a putative coat protein of high abundance. The products of the other nine ORFs are probably involved in polysaccharide biosynthesis, nucleotide metabolism, and DNA modification. The viral genome divides into two nearly equal halves of opposite gene orientation. This observation as well as a GC-skew analysis point to the presence of a putative viral origin of replication in the 1.4-kb intergenic region between ORF1 and ORF74. Both morphological and genomic features identify STSV1 as a novel virus infecting the genus Sulfolobus.

Typing of nonencapsulated haemophilus strains by repetitive-element sequence-based PCR using intergenic dyad sequences

Intergenic dyad sequences (IDS) are short repeated elements that have been described for several Haemophilus genomes and for only two other bacterial genera. We developed a repetitive-element sequence-based PCR using an IDS-specific primer as a typing method (IDS-PCR) for nonencapsulated Haemophilus strains and compared this technique with pulsed-field gel electrophoresis (PFGE) of DNA restricted with SmaI. IDS-PCR was rapid, easy to perform, and reproducible, with a high discriminatory capacity for nontypeable Haemophilus influenzae (NTHI) strains. The 69 NTHI strains tested generated 65 different banding patterns. Epidemiologically related strains gave similar or identical fingerprints, and all of the unrelated strains except two showed different patterns. These results were in agreement with those obtained by PFGE. For 20 genital strains usually identified as being biotype IV NTHI and belonging to a cryptic genospecies of Haemophilus with remarkable genetic homogeneity, four bands were significantly present and six bands were significantly absent from the fingerprints. The 20 strains were gathered in 11 closely related profiles, whereas PFGE provided no band when DNA was treated with SmaI. IDS-PCR improved the differentiation previously obtained within this species by ribotyping and multilocus enzyme electrophoresis. Our findings suggest that IDS-PCR is a rapid, reliable, and discriminatory method for typing NTHI strains and is currently the most efficient method for distinguishing strains within the cryptic genospecies of HAEMOPHILUS:

Circular permutation of DNA cytosine-N4 methyltransferases: in vivo coexistence in the BcnI system and in vitro probing by hybrid formation

Sequence analysis of the BcnI restriction-modification system from Bacillus centrosporus revealed four open reading frames (bcnIC, bcnIR, bcnIB and bcnIA) that are arranged as two converging collinear pairs. One pair encodes a putative small regulatory protein, C.BcnI, and the restriction endonuclease R.BcnI. The other two gene products are the DNA cytosine-N4 methyltransferases M.BcnIA and M.BcnIB, which differ by circular permutation of conserved sequence motifs. The BcnI methyltransferases are isospecific on double-stranded DNA [methylation specificity CC(C/G)GG], but M.BcnIA can also methylate the target sites in single-stranded DNA. Functional analysis shows that bcnIA is dispensable (bcnIB is capable of protecting the DNA against the in vivo activity of bcnIR); in contrast, no stable clones were obtained if bcnIB alone was deleted from the system. By analogy with the DpnII system, the second methylase M.BcnIA may play a role in the transformation proficiency of its gram-positive host. The interchangeability of homologous elements in the beta class of cytosine-N4 methylases was probed by hybrid formation between M.BcnIB and its closest homolog M.Cfr9I (CCCGGG) employing a novel semi-random strategy combined with selection for catalytic activity. The fusion points in the active hybrids mapped in a narrow region located between sequence motifs X and I. Our data illustrate that recombination of two related sequences by circular permutation may serve as an evolutionary mechanism for creating new specificities of amino MTases.

Crystal structure of the Bse634I restriction endonuclease: comparison of two enzymes recognizing the same DNA sequence

Crystal structures of Type II restriction endonucleases demonstrate a conserved common core and active site residues but diverse structural elements involved in DNA sequence discrimination. Comparative structural analysis of restriction enzymes recognizing the same nucleotide sequence might therefore contribute to our understanding of the structural diversity of specificity determinants within restriction enzymes. We have solved the crystal structure of the Bacillus stearothermophilus restriction endonuclease Bse634I by the multiple isomorphous replacement technique to 2.17 A resolution. Bse634I is an isoschisomer of the Cfr10I restriction enzyme whose crystal structure has been reported previously. Comparative structural analysis of the first pair of isoschisomeric enzymes revealed conserved structural determinants of sequence recognition and catalysis. However, conformations of the N-terminal subdomains differed between Bse634I/Cfr10I, suggesting a rigid body movement that might couple DNA recognition and catalysis. Structural similarities extend to the quaternary structure level: crystal contacts suggest that Bse634I similarly to Cfr10I is arranged as a tetramer. Kinetic analysis reveals that Bse634I is able to interact simultaneously with two recognition sites supporting the tetrameric architecture of the protein. Thus, restriction enzymes Bse634I, Cfr10I and NgoMIV, recognizing overlapping nucleotide sequences, exhibit a conserved tetrameric architecture that is of functional importance.

Recognition of native DNA methylation by the PvuII restriction endonuclease

Recognizing the methylation status of specific DNA sequences is central to the function of many systems in eukaryotes and prokaryotes. Restriction-modification systems have to distinguish between 'self' and 'non-self' DNA and depend on the inability of restriction endonucleases to cleave their DNA substrates when the DNA is appropriately methylated. These endonucleases thus provide a model system for studying the recognition of DNA methylation by proteins. We have characterized the interaction of R.PVU:II with DNA containing the physiologically relevant N4-methylcytosine modification. R.PVU:II binds (N4m)C-modified DNA and cleaves it very slowly. Methylated strands in hemimethylated duplexes were cleaved at a higher rate than in fully methylated duplexes, in parallel with a higher binding affinity for hemimethylated DNA. The co-crystal structures of R.PVU:II-DNA, together with a mutagenesis study, have implicated specific amino acids in recognition of the methylatable base; one of these is His84. We report that replacing His84 with Ala reduced the rate of cleavage of unmodified DNA but, in contrast, slightly increased the cleavage of (N4m)C-modified DNA.

Molecular evolution of DNA-(cytosine-N4) methyltransferases: evidence for their polyphyletic origin

DNA N4-cytosine methyltransferases (N4mC MTases) are a family of S-adenosyl-L-methionine (AdoMet)-dependent MTases. Members of this family were previously found to share nine conserved sequence motifs, but the evolutionary basis of these similarities has never been studied in detail. We performed phylogenetic analysis of 37 known and potential new family members from the multiple sequence alignment using distance matrix, parsimony and maximum likelihood approaches to infer the evolutionary relationship among the N4mC MTases and classify them into groups of orthologs. All the treeing algorithms employed as well as results of exhaustive sequence database searching support a scenario, in which the majority of N4mC MTases, except for M. Bal I and M. Bam HI, arose by divergence from a common ancestor. Interestingly, MTases M. Bal I and M. Bam HI apparently originated from N6-adenine MTases and represent the most recent addendum to the N4mC MTase family. In addition to the previously reported nine sequence motifs, two more conserved sequence patches were detected. Phylogenetic analysis also provided the evidence for massive horizontal transfer of MTase genes, presumably with the whole restriction-modification systems, between Bacteria and Archaea.

Prophage carriage as a molecular epidemiological marker in Streptococcus pneumoniae

The great majority of clinical isolates of Streptococcus pneumoniae carry prophages that may be identified through their hybridization with a DNA probe specific for the pneumococcal lytA gene (M. Ramirez, E. Severina, and A. Tomasz, J. Bacteriol. 181:3618-3625, 1999). We now show that the lytA hybridization pattern of chromosomal SmaI digests is stable for a given strain during extensive serial culturing in the laboratory; the pattern is specific for the strain's clonal type, as defined by pulsed-field gel electrophoretis (PFGE) pattern, and variations in PFGE subtypes may be explained by changes in the number and chromosomal localization of this prophage(s). These observations indicate that the lytA hybridization pattern may be used as a molecular epidemiological marker that offers additional resolution of the genetic background of S. pneumoniae strains.

Perinatal lethality in H19 enhancers-Igf2 transgenic mice

The insulin-like growth factor II (IGFII) is a mitogen for a number of cell types in vitro and is required for normal embryonic growth. It has been hypothesized that overexpression of IGF2 is responsible for the increased growth and tumor predisposition in patients with Beckwith-Wiedemann syndrome. Association of increased levels of IGFII with increased growth is also incorporated in a current model for the evolution of Igf2 imprinting. Different experimental approaches to increasing IGFII levels in the mouse have yielded different results with respect to its effects on growth, viability, and tumor development. To investigate the consequences of IGf2 overexpression in the embryonic period, without alterations in the activity of other genes, we produced transgenic mice that express the Igf2 gene under the control of the H19 enhancers. Transgene expression in the embryonic period had no significant effect on the overall size of the embryos, but was associated with perinatal lethality in homozygous, and some heterozygous, mice. A large fraction of homozygous mice also developed a cleft palate. These findings indicate that overexpression of Igf2 can have an adverse effect on viability in the absence of a pronounced effect on overall body growth. The results are consistent with the view that growth and perinatal viability are affected differently by Igf2 overexpression in endodermal and mesodermal tissues.

Solid phase synthesis and restriction endonuclease cleavage of oligodeoxynucleotides containing 5-(hydroxymethyl)-cytosine

Emerging data suggest an important role for cytosine methylation in tumorigenesis. Simultaneously, recent studies indicate a significant contribution of endogenous oxidative DNA damage to the development of human disease. Oxidation of the 5-methyl group of 5-methylcytosine (5mC) residues in DNA results in the formation of 5-(hydroxymethyl)cytosine (hmC). The biological consequences ofhmC residues in vertebrate DNA are as yet unknown; however, conversion of the hydrophobic methyl group to the hydrophilic hydroxymethyl group may substantially alter the interaction of sequence-specific binding proteins with DNA. Central to both biophysical and biochemical studies on the potential consequences of specific DNA damage products such as hmC are efficient methods for the synthesis of oligodeoxynucleotides containing such modified bases at selected positions. In this paper, we describe a method for the placement of hmC residues in oligodeoxynucleotides using established phosphoramidite chemistry. In addition, we have examined the influence of specific hmC residues on enzymatic cleavage of oligodeoxynucleotides by the methylation-sensitive restriction endonucleases MspI and HpaII.

Reduced DNA methylation in Arabidopsis thaliana results in abnormal plant development

Arabidopsis plants transformed with an antisense construct of an Arabidopsis methyltransferase cDNA (METI) have reduced cytosine methylation in CG dinucleotides. Methylation levels in progeny of five independent transformants ranged from 10% to 100% of the wild type. Removal of the antisense construct by segregation in sexual crosses did not fully restore methylation patterns in the progeny, indicating that methylation patterns are subject to meiotic inheritance in Arabidopsis. Plants with decreased methylation displayed a number of phenotypic and developmental abnormalities, including reduced apical dominance, smaller plant size, altered leaf size and shape, decreased fertility, and altered flowering time. Floral organs showed homeotic transformations that were associated with ectopic expression of the floral homeotic genes AGAMOUS and APETALA3 in leaf tissue. These observations suggest that DNA methylation plays an important role in regulating many developmental pathways in plants and that the developmental abnormalities seen in the methyltransferase antisense plants may be due to dysregulation of gene expression.

Overproduction of DNA cytosine methyltransferases causes methylation and C --> T mutations at non-canonical sites

Multicopy clones of Escherichia coli cytosine methyltransferases Dcm and EcoRII methylase (M. EcoRII) cause an approximately 50-fold increase in C --> T mutations at their canonical site of methylation, 5'-CmeCAGG (meC is 5-methylcytosine). These plasmids also cause transition mutations at the second cytosine in the sequences CCGGG at approximately 10-fold lower frequency. Similarly, M. HpaII was found to cause a significant increase in C --> T mutations at a CCAG site, in addition to causing mutations at its canonical site of methylation, CCGG. Using a plasmid that substantially overproduces M. EcoRII, in vivo methylation at CCSGG (S is C or G) and other non-canonical sites could be detected using a gel electrophoretic assay. There is a direct correlation between the level of M. EcoRII activity in cells, the extent of methylation at non-canonical sites and frequency of mutations at these same sites. Overproduction of M. EcoRII in cells also causes degradation of DNA and induction of the SOS response. In vitro, M. EcoRII methylates an oligonucleotide duplex containing a CCGGG site at a slow rate, suggesting that overproduction of the enzyme is essential for significant amounts of such methylation to occur. Together these results show that cytosine methyltransferases occasionally methylate cellular DNA at non-canonical sites and suggest that in E. coli, methylation-specific restriction systems and sequence specificity of the DNA mismatch correction systems may have evolved to accommodate this fact. These results also suggest that mutational effects of cytosine methyltransferases may be much broader than previously imagined.

Binding of histone H1 to DNA is indifferent to methylation at CpG sequences

The possibility that histone H1 binds preferentially to DNA containing 5-methylcytosine in the dinucleotide CpG is appealing, as it could help to explain the repressive effects of methylation on gene activity. In this study, the affinity of purified H1 for methylated and non-methylated DNA sequences has been tested using both naked DNA and chromatin. Based on a variety of assays (bandshifts, filter-binding assays, Southwestern blots, and nuclease sensitivity assays), we conclude that H1 has no significant preference for binding to naked methylated DNA. Similarly, H1 showed the same affinities for methylated and non-methylated DNA when assembled into chromatin in a Xenopus oocyte extract. Thus potential cooperative interaction of H1 with polynucleosomal complexes is not enhanced by the presence of DNA methylation.

Determination of methylation specificity of DsaV methyltransferase by a simple biochemical method

We have developed a simple new method that can identify the base methylated by a sequence-specific DNA methyltransferase and have used it to identify the cytosine that is methylated by DsaV methyltransferase (M. DsaV) within its recognition sequence 5'-CCNGG. The method utilizes the fact that exonuclease III of E. coli does not degrade DNA ends with 3' overhangs and cannot hydrolyze a phosphorothioate linkage. DNA duplexes containing phosphorothioate linkages at specific positions were methylated with M. DsaV in the presence of [methyl-3H] S-adenosylmethionine and were subjected to exonuclease III digestion. The pattern of [methyl-3H] dCMP release from the duplexes was consistent with the methylation of the internal cytosine in CCNGG, but not of the outer cytosine. To establish the accuracy of this method, we confirmed the known specificity of EcoRII methyltransferase by the method. We also confirmed the specificity of M. DsaV using an established biochemical method that involves the use of a type IIS restriction enzyme. Methylation of CCWGG (W = A or T) sequences at the internal cytosines is native to E. coli and is not restricted by the modified cytosine restriction (Mcr) systems. Surprisingly, the gene for M. DsaV was significantly restricted by the McrBC system. We interpret this to mean that M. DsaV may occasionally methylate at sequences other than CCNGG or may occasionally methylate the outer cytosine in its recognition sequence.

The bulk chromatin structure of a murine transgene does not vary with its transcriptional or DNA methylation status

The DNA methylation status of HRD, a murine transgene, can be controlled by the genetic background upon which it is carried. We found the transgene to be transcribed in competent tissues only when undermethylated. Chromatin structure over the transgene was assayed by nuclear accessibility with DNase I, MspI, and PstI. While the transgene was up to fivefold more resistant to MspI when methylated than when not methylated, we observed no such difference with DNase I or PstI. We suggest that methyl-CpG-binding proteins are responsible for the difference observed with MspI, but that the chromatin structures are otherwise similarly compacted. Methylation could, therefore, play a regulatory role in gene expression beyond that which can be accomplished by bulk chromatin structure alone.

Organization and sequence of the HpaII restriction-modification system and adjacent genes

We report the organization of the HpaII restriction and modification (R-M) system from Haemophilus parainfluenzae (recognition sequence: 5'...CCGG...3'), the sequence of the gene coding for the HpaII restriction endonuclease, and the sequence of the upstream flanking DNA. The HpaII system comprises two genes, hpaIIM, coding for the methyltransferase (MTase; 358 amino acids (aa), 40.4 kDa: product, Cm5CGG), and hpaIIR, coding for the restriction endonuclease (ENase; 358 aa, 40.9 kDa: product, C'CGG). The genes are adjacent, they have the same orientation, and they occur in the order hpaIIM then hpaIIR. The ENase bears little as sequence similarity to the isoschizomeric R.BsuFI and R.MspI ENases. Upstream of, and partly overlapping hpaIIM is the coding sequence for a 141-aa protein that resembles the very-short-patch-repair endonuclease (Vsr) of Escherichia coli. Upstream of that is the coding sequence for a protein that resembles valyl-tRNA synthetase (ValS).

Cloning and analysis of translational control for genes encoding the Cfr9I restriction-modification system

The complete type-II Cfr9I restriction-modification (R-M) system of Citrobacter freundii strain RFL9, recognizing the DNA sequence CCCGGG, has been cloned and expressed, and functionally active enzymes have been produced in Escherichia coli. Both the methyltransferase (MTase; M.Cfr9I) and restriction endonuclease (ENase; R.Cfr9I) were found to be encoded on a 2.3-kb cloned fragment in the same transcriptional orientation, but differing in translational phases. The last codon (underlined) (ATGA) of the MTase-encoding gene (Cfr9IM) overlaps with the start codon for the ENase-encoding gene (overlined) (cfr9IR). A nucleotide sequence complementary to a predicted Shine-Dalgarno sequence preceding cfr9IR is within this gene. Predicted free energy (delta G) for formation of the mRNA secondary structure involving these complementary sequences was found to be -16.1 kcal/mol. Amino-acid sequence homology of 80% was found between R.Cfr9I and R.XcyI.

The endonuclease isoschizomers, SmaI and XmaI, bend DNA in opposite orientations

The SmaI and XmaI endonucleases are imperfect isoschizomers that recognize the sequence CCCGGG. SmaI cleaves between the internal CpG to produce blunt end scissions whereas XmaI cleaves between the external cytosines to produce a four base, five prime overhang. Each of the endonucleases forms stable, specific complexes with DNA in the absence of magnesium. Circular permutation analyses of the protein-DNA complexes revealed that each of the endonucleases induces bending of the DNA. Phase sensitive detection analyses verified the existence of the SmaI and XmaI induced bends. Furthermore, bending of the helix axis by the endonucleases appeared to be directed in opposite orientations. The orientation of the SmaI induced bend appeared to be towards the major groove and is reminiscent of the direction of the bend induced by EcoRV which similarly induces blunt end scissions. Conversely, XmaI appeared to bend the DNA towards the minor groove.

Restriction endonuclease cleavage of 5-methyl-deoxycytosine hemimethylated DNA at high enzyme-to-substrate ratios

We have investigated the ability of a large number of restriction enzymes to digest non-canonically hemimethylated DNA at high enzyme-to-substrate ratios. A single-stranded unmethylated phagemid was used as a template to complete synthesis of the second strand using 5-methyl-dCTP to substitute for all the deoxycytosine residues. A fragment of this double-stranded hemimethylated DNA which contains the multiple cloning site region was used as a substrate. For all the enzymes tested, at least some degree of protection from digestion is observed. Sites completely protected from digestion by their cognate enzymes are SalI, BstXI, SacI, SacII, SmaI, SstI, XhoI, PstI, HinfI, BamHI and AccI. Sites partially protected from digestion by their cognate enzymes are XbaI, HindIII, KpnI, SpeI, ClaI, EcoRI and PvuII. Knowledge of the sensitivity of commonly used restriction enzymes to hemimethylated substrates is useful for several applications, which will be discussed.

Structure and evolution of the XcyI restriction-modification system

The XcyI restriction-modification system from Xanthomonas cyanopsidis recognizes the sequence, CCCGGG. The XcyI endonuclease and methylase genes have been cloned and sequenced and were found to be aligned in a head to tail orientation with the methylase preceding and overlapping the endonuclease by one base pair. The nucleotide sequence codes for an N4 cytosine methyltransferase with a predicted molecular weight of 33,500 and an endonuclease comprised of 333 codons and a molecular weight of 36,600. Sequence comparisons revealed significant similarity between the XcyI, CfrI and SmaI methylisomers. In contrast, no similarity was detected between the primary structures of the XcyI and SmaI endonucleases. The XcyI restriction-modification system is highly homologous to the XmaI genes, although the DNA sequences flanking the genes rapidly diverge. The sequence of the XcyI endonuclease contains two motifs which have recently been identified as essential to the activity of the EcoRV endonuclease.

Alw26I, Eco31I and Esp3I--type IIs methyltransferases modifying cytosine and adenine in complementary strands of the target DNA

The specificity of three DNA methyltransferases M.Alw26I, M.Eco31I and M.Esp3I, isolated from Acinetobacter Iwoffi RFL26, Escherichia coli RFL31 and Hafnia alvei RFL3+, respectively, was determined. All the enzymes methylate both strands of asymmetric recognition sites yielding m5C in the top-strand and m6A in the bottom-strand, as below: 5'-GTm5CTC 5'-GGTm5CTC 5'-CGTm5CTC 3'-Cm6AGAG 3'-CCm6AGAG 3'-GCm6AGAG (M.Alw26I) (M.Eco31I) (M.Esp3I) They are the first members of type IIs methyltransferases that modify different types of nucleotides in the recognition sequence.

Sequence-specific DNA binding by the MspI DNA methyltransferase

The MspI methyltransferase (M.MspI) recognizes the sequence CCGG and catalyzes the formation of 5-methylcytosine at the fist C-residue. We have investigated the sequence-specific DNA-binding properties of M.MspI under equilibrium conditions, using gel-mobility shift assays and DNasel footprinting. M.MspI binds to DNA in a sequence-specific manner either alone or in the presence of the normal methyl donor S-adenosyl-L-methionine as well as the analogues, sinefungin and S-adenosyl-L-homocysteine. In the presence of S-adenosyl-L-homocysteine, M.MspI shows the highest binding affinity to DNA containing a hemimethylated recognition sequence (Kd = 3.6 x 10(-7) M), but binds less well to unmethylated DNA (Kd = 8.3 x 10(-7) M). Surprisingly it shows specific, although poor, binding to fully methylated DNA (Kd = 4.2 x 10(-6) M). M.MspI binds approximately 5-fold more tightly to DNA containing its recognition sequence, CCGG, than to nonspecific sequences in the absence of cofactors. In the presence of S-adenosyl-L-methionine, S-adenosyl-L-homocysteine or sinefungin the discrimination between specific and non-specific sequences increases up to 100-fold. DNasel footprinting studies indicate that 16 base pairs of DNA are covered by M.MspI, with the recognition sequence CCGG located asymmetrically within the footprint.

The sequence specificity domain of cytosine-C5 methylases

Prokaryotic DNA[cytosine-C5]methyltransferases (m5C-methylases) share a common architectural arrangement of ten conserved sequence motifs. A series of eleven hybrids have been constructed between the HpaII (recognition sequence: Cm5CGG) and HhaI (recognition sequence: Gm5CGC) DNA-methylases. The hybrids were over-expressed in E.coli and their in vivo methylation phenotypes investigated. Six were inactive by our assay while five of them retained partial methylation activity and full specificity. In all five cases the specificity matched that of the parent methylase which contributed the so-called variable region, located between conserved motifs VIII and IX. This was the only sequence held in common between the active hybrids and for the first time provides unequivocal evidence that the specificity determinants of the mono-specific m5C-methylases are located within the variable region. Correlation of the hybrid methylase structure with the efficiency of methylation suggests that conserved motif IX may interact with the variable region whereas motif X most probably interacts with the N-terminal half of the molecule.

M.Smal is an N4-methylcytosine specific DNA-methylase

An enzymatic activity rendering DNA immune to the action of the Smal restriction endonuclease in the presence of S-adenosyl-L-methionine has been detected in Serratia marcescens Sb. This methylase, M.Smal, modifies the second cytosine residue of the substrate sequence CCCGGG yielding N4-methylcytosine.

Specificity of restriction endonucleases and DNA modification methyltransferases a review (Edition 3)

The properties and sources of all known class-I, class-II and class-III restriction endonucleases (ENases) and DNA modification methyltransferases (MTases) are listed and newly subclassified according to their sequence specificity. In addition, the enzymes are distinguished in a novel manner according to sequence specificity, cleavage position and methylation sensitivity. Furthermore, new nomenclature rules are proposed for unambiguously defined enzyme names. In the various Tables, the enzymes are cross-indexed alphabetically according to their names (Table I), classified according to their recognition sequence homologies (Table II), and characterized within Table II by the cleavage and methylation positions, the number of recognition sites on the DNA of the bacteriophages lambda, phi X174, and M13mp7, the viruses Ad2 and SV40, the plasmids pBR322 and pBR328, and the microorganisms from which they originate. Other tabulated properties of the ENases include relaxed specificities (integrated within Table II), the structure of the generated fragment ends (Table III), interconversion of restriction sites (Table IV) and the sensitivity to different kinds of DNA methylation (Table V). Table VI shows the influence of class-II MTases on the activity of class-II ENases with at least partially overlapping recognition sequences. Table VII lists all class-II restriction endonucleases and MTases which are commercially available. The information given in Table V focuses on the influence of methylation of the recognition sequences on the activity of ENases. This information might be useful for the design of cloning experiments especially in Escherichia coli containing M.EcodamI and M.EcodcmI [H16, M21, U3] or for studying the level and distribution of site-specific methylation in cellular DNA, e.g., 5'- (M)CpG-3' in mammals, 5'-(M)CpNpG-3' in plants or 5'-GpA(M)pTpC-3' in enterobacteria [B29, E4, M30, V4, V13, W24]. In Table IV a cross index for the interconversion of two- and four-nt 5'-protruding ends into new recognition sequences is complied. This was obtained by the fill-in reaction with the Klenow (large) fragment of the E. coli DNA polymerase I (PolIk), or additional nuclease S1 treatment followed by ligation of the modified fragment termini [P3]. Interconversion of restriction sites generates novel cloning sites without the need of linkers. This should improve the flexibility of genetic engineering experiments [K56, P3].(ABSTRACT TRUNCATED AT 400 WORDS)

Cloning and nucleotide sequence of the gene encoding the Ecal DNA methyltransferase

The gene coding for the GGTNACC specific Ecal DNA methyltransferase (M.Ecal) has been cloned in E. coli from Enterobacter cloacae and its nucleotide sequence has been determined. The ecalM gene codes for a protein of 452 amino acids (Mr: 51,111). It was determined that M.Ecal is an adenine methyltransferase. M.Ecal shows limited amino acid sequence similarity to other adenine methyltransferases. A clone that expresses Ecal methyltransferase at high level was constructed.

Sequence motifs characteristic of DNA[cytosine-N4]methyltransferases: similarity to adenine and cytosine-C5 DNA-methylases

The sequences coding for DNA[cytosine-N4]methyltransferases MvaI (from Micrococcus varians RFL19) and Cfr9I (from Citrobacter freundii RFL9) have been determined. The predicted methylases are proteins of 454 and 300 amino acids, respectively. Primary structure comparison of M.Cfr9I and another m4C-forming methylase, M.Pvu II, revealed extended regions of homology. The sequence comparison of the three DNA[cytosine-N4]-methylases using originally developed software revealed two conserved patterns, DPF-GSGT and TSPPY, which were found similar also to those of adenine and DNA[cytosine-C5]-methylases. These data provided a basis for global alignment and classification of DNA-methylase sequences. Structural considerations led us to suggest that the first region could be the binding site of AdoMet, while the second is thought to be directly involved in the modification of the exocyclic amino group.

Cloning, characterization and heterologous expression of the SmaI restriction-modification system

The genes coding for the class-II Serratia marcescens restriction-modification system have been cloned and expressed in E. coli. Recombinant clones, restricted incoming phage only poorly; the recombinant plasmids, however, became fully modified in vivo, i.e. completely resistant against digestion with R.SmaI. The determined nucleotide sequence of the cloned system revealed three open reading frames with lengths of 252 bp, 741 bp, and 876 bp. Through various deletion experiments and an insertion-mutation experiment the 876 bp open reading frame could be assigned to the SmaI DNA modification enzyme and the 741 bp open reading frame to the SmaI restriction endonuclease. Mapping of the transcription start sites of the genes revealed that the SmaI endonuclease is transcribed as polycistronic mRNA together with a 252 bp long preceding open reading frame of unknown function. No homology was found when comparing the amino acid sequence of M.SmaI with the published sequences of m5C-specific DNA modification methyltransferases. On the other hand, a stretch of 14 amino acids in the C-proximal region of M.SmaI shows a significant homology to the C-proximal amino acid sequences of the N6A-methyltransferases M.HinfI and M.DpnIIA and the N4C-methyltransferase M.PvuII.