Microarray-based method to evaluate the accuracy of restriction endonucleases HpaII and MspI

A double-strand DNA (ds DNA) microarray was fabricated to analyze the structural perturbations caused by methylation and the different base mismatches in the interaction of the restriction endonucleases HpaII and MspI with DNA. First, a series of synthesized oligonucleotides were arrayed on the aldehyde-coated glass slides. Second, these oligonucleotides were hybridized with target sequences to obtain ds DNA microarray, which includes several types of double strands with the fully methylated, semi-methylated, and unmethylated canonical recognition sequences, semi-methylated and unmethylated base mismatches within the recognition sequences. The cleavage experiments were carried out under normal buffer conditions. The results indicated that MspI could partially cleave methylated and semi-methylated canonical recognition sequences. In contrast, HpaII could not cleave methylated and semi-methylated canonical recognition sequences. HpaII and MspI could both cleave the unmethylated canonical recognition sequence. However, HpaII could partially cleave the sequence containing one GG mismatch and not cleave other base mismatches in the corresponding recognition site. In contrast, MspI could not recognize the base mismatches within the recognition sequence. A good reproducibility was observed in several parallel experiments. The experiment indicates that the microarray technology has great potentials in high-throughput identifying important interactions between protein and DNA.     

Cloning of the MspI modification enzyme. The site of modification and its effects on cleavage by MspI and HpaII

The gene for the MspI modification enzyme from Moraxella was cloned in Escherichia coli using the plasmid vector pBR322. Selection of transformants carrying the gene was based on the resistance of the modified plasmid encoding the enzyme to cleavage by MspI. Both chromosomal and plasmid DNA were modified in the selected clones. None of the clones obtained produced the cognate restriction enzyme which suggests that in this system the genes for the restriction enzyme and methylase are not closely linked. Crude cell extracts prepared from the recombinant strains, but not the host (E. coli HB101), contain an S-adenosylmethionine-dependent methyltransferase specific for the MspI recognition site, CCGG. Production of the enzyme is 3-4-fold greater in the transformants than in the original Moraxella strain. 5-Methylcytosine was identified as the product of the reaction chromatographically. The outer cytosine of the recognition sequence, *CCGG, was shown to be the site of methylation by DNA-sequencing methods. This modification blocks cleavage by both MspI and its isoschizomer HpaII. HpaII, but not MspI, is able to cleave the unmethylated strand of a hemimethylated substrate. The relevance of these results to the use of MspI and HpaII to analyze patterns of methylation in genomic DNA is discussed.     

HhaI and HpaII DNA methyltransferases bind DNA mismatches, methylate uracil and block DNA repair.

The hydrolytic deamination of 5-methylcytosine (5-mC) to thymine (T) is believed to be responsible for the high mutability of the CpG dinucleotide in DNA. We have shown a possible alternate mechanism for mutagenesis at CpG in which HpaII DNA-(cytosine-5) methyltransferase (M.HpaII) can enzymatically deaminate cytosine (C) to uracil (U) in DNA [Shen, J.-C., Rideout, W.M., III and Jones, P.A., Cell, 71, 1073-1080, (1992)]. Both the hydrolytic deamination of 5-mC and enzymatic deamination of C create premutagenic DNA mismatches (G:U and G:T) with the guanine (G) originally paired to the normal C. Surprisingly, we found that DNA-(cytosine-5) methyltransferases have higher affinities for these DNA mismatches than for their normal G:C targets and are capable of transferring a methyl group to the 5-position of U, creating T at low efficiencies. This binding by methyltransferase to mismatches at the recognition site prevented repair of G:U mismatches by uracil DNA glycosylase in vitro.     

How to decide which are the most pertinent overly-represented features during gene set enrichment analysis

Background  

The search for enriched features has become widely used to characterize a set of genes or proteins. A key aspect of this technique is its ability to identify correlations amongst heterogeneous data such as Gene Ontology annotations, gene expression data and genome location of genes. Despite the rapid growth of available data, very little has been proposed in terms of formalization and optimization. Additionally, current methods mainly ignore the structure of the data which causes results redundancy. For example, when searching for enrichment in GO terms, genes can be annotated with multiple GO terms and should be propagated to the more general terms in the Gene Ontology. Consequently, the gene sets often overlap partially or totally, and this causes the reported enriched GO terms to be both numerous and redundant, hence, overwhelming the researcher with non-pertinent information. This situation is not unique, it arises whenever some hierarchical clustering is performed (e.g. based on the gene expression profiles), the extreme case being when genes that are neighbors on the chromosomes are considered.     

In situ nick translation of human metaphase chromosomes with the restriction enzymes MspI and HpaII reveals an R-band pattern.

Several restriction enzymes (HindIII, HaeIII, MspI, HpaII, EcoRI, KpnI, and NotI) were evaluated for their ability to induce bands in human metaphase chromosomes during in situ nick translation. MspI and HpaII were able to induce a completely developed R-band pattern. Preferential cleavage of R-band chromatin is due to the presence of unmethylated CpG-residues present in CpG-rich islands, which are apparently unevenly distributed and mainly concentrated in R-bands.     

Changes in methylation of tissue cultured soybean cells detected by digestion with the restriction enzymes HpaII and MspI

Each of the tandemly arranged 5S RNA genes of soybean contain two CCGG sites which, if unmethylated, can be digested by both MspI and HpaII. Methylation of the internal cytosine (CmeCGG) prevents digestion by HpaII but allows digestions by MspI.Suspension cultures were prepared from soybean plants and the DNA from these cultures was examined for the susceptibility of 5S RNA genes to digestion by MspI and HpaII. 5S genes from DNA extracted from intact plants can be partially digested with MspI but not at all by HpaII. In contrast, shortly after cells were cultured the 5S RNA could be hydrolyzed by both HpaII and MspI. After prolonged cell culture, the 5S genes from some cell lines were found to have become partially or even completely resistant to HpaII digestion. The results suggest that lack of methylation can occur when cells are cultured and that such methylation may play a role in the heritable changes observed in cell culture.Research supported by Grant 01498 from the National Institutes of Environmental Health Sciences     

How to decide whether small samples comply with an equidistribution

The decision whether a measured distribution complies with an equidistribution is a central element of many biostatistical methods. High throughput differential expression measurements, for instance, necessitate to judge possible over-representation of genes. The reliability of this judgement, however, is strongly affected when rarely expressed genes are pooled. We propose a method that can be applied to frequency ranked distributions and that yields a simple but efficient criterion to assess the hypothesis of equiprobable expression levels. By applying our technique to surrogate data we exemplify how the decision criterion can differentiate between a true equidistribution and a triangular distribution. The distinction succeeds even for small sample sizes where standard tests of significance (e.g. chi(2)) fail. Our method will have a major impact on several problems of computational biology where rare events baffle a reliable assessment of frequency distributions. The program package is available upon request from the authors.     

MspI, an isoschizomer of hpaII which cleaves both unmethylated and methylated hpaII sites.

The cleavage of DNA by restriction endonucleases HpaII and HapII is prevented by the presence of a 5-methyl group at the internal C residue of its recognition sequence CCGG. MspI, an isoschizomer of HpaII available from New England Biolabs, cleaves DNA irrespective of the presence of a methyl group at this position. This enzyme cleaves DNA from Haemophilus parainfluenzae and Haemophilus aphrophilus readily while HpaII and HapII cannot degrade these DNAs. Practically all HpaII sites in mammalian sperm DNA are also protected by methylation at the internal C position since HpaII and HapII barely cleave this DNA (average molecular weight 40 kb). MspI, however, cleaves the DNA to an average size of about 5 kb.     

The 5'-GGATCC-3' cleavage specificity of BamHI is increased to 5'-CCGGATCCGG-3' by sequential double methylation with M.HpaII and M.BamHI

Site-specific DNA methylation is known to block cleavage by a number of restriction endonucleases. We show that methylation at 'non-canonical' DNA modification sites can also block methylation by five of 13 DNA methyltransferases (MTases) tested. Furthermore, MTases and endonucleases that recognize the same nucleotide sequence can differ in their sensitivity to non-canonical methylation. In particular, BamHI endonuclease can cut 5'-GGATCm5C efficiently, whereas M.BamHI cannot methylate this modified sequence. Methyltransferase/endonuclease pairs which differ in their sensitivity to non-canonical methylation can be exploited to generate rare DNA cleavage sites. For example, we show that M.HpaII, M.BamHI, and BamHI can be used sequentially in a three-step procedure to specifically cleave DNA at the 10-bp sequence 5'-CCGGATCCGG. Several highly selective DNA cutting strategies are made possible by these sequential double methylation-blocking reactions.     

Cloning and characterization of the HpaII methylase gene.

The HpaII restriction-modification system from Haemophilus parainfluenzae recognizes the DNA sequence CCGG. The gene for the HpaII methylase has been cloned into E. coli and its nucleotide sequence has been determined. The DNA of the clones is fully protected against cleavage by the HpaII restriction enzyme in vitro, indicating that the methylase gene is active in E. coli. The clones were isolated in an McrA-strain of E. coli; attempts to isolate them in an McrA+ strain were unsuccessful. The clones do not express detectable HpaII restriction endonuclease activity, suggesting that either the endonuclease gene is not expressed well in E. coli, or that it is not present in its entirety in any of the clones that we have isolated. The derived amino acid sequence of the HpaII methylase shows overall similarity to other cytosine methylases. It bears a particularly close resemblance to the sequences of the HhaI, BsuFI and MspI methylases. When compared with three other methylases that recognize CCGG, the variable region of the HpaII methylase, which is believed to be responsible for sequence specific recognition, shows some similarity to the corresponding regions of the BsuFI and MspI methylases, but is rather dissimilar to that of the SPR methylase.     

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

The cleavage specificity of R.Cfr9I was determined to be C decreases CCGGG whereas the methylation specificity of M.Cfr9I was C4mCCGGG. The action of MspI, HpaII, SmaI, XmaI and Cfr9I restriction endonucleases on an unmethylated parent d(GGACCCGGGTCC) dodecanucleotide duplex and a set of oligonucleotide duplexes, containing all possible substitutions of either 4mC or 5mC for C in the CCCGGG sequence, was investigated. It was found that 4mC methylation, in contrast to 5mC, renders the CCCGGG site resistant to practically all the investigated endonucleases. The cleavage of methylated substrates with restriction endonucleases is discussed.     

The sequence GGCmCGG is resistant to MspI cleavage.

MspI essentially fails to cut the sequence GGCmCGG at enzyme concentrations which give total digestion of CCGG, CmCGG and GGCCGG sites. This result explains why certain sites in mammalian DNA are resistant to both MspI and HpaII and shows that this results from an idiosynchracy of MspI rather than a novel form of DNA methylation at this site in mammalian cells.     

M.HhaI binds tightly to substrates containing mismatches at the target base.

The (cytosine-5) DNA methyltransferase M.HhaI causes its target cytosine base to be flipped completely out of the DNA helix upon binding. We have investigated the effects of replacing the target cytosine by other, mismatched bases, including adenine, guanine, thymine and uracil. We find that M.HhaI binds more tightly to such mismatched substrates and can even transfer a methyl group to uracil if a G:U mismatch is present. Other mismatched substrates in which the orphan guanine is changed exhibit similar behavior. Overall, the affinity of DNA binding correlates inversely with the stability of the target base pair, while the nature of the target base appears irrelevant for complex formation. The presence of a cofactor analog. S-adenosyl-L-homocysteine, greatly enhances the selectivity of the methyltransferase for cytosine at the target site. We propose that the DNA methyltransferases have evolved from mismatch binding proteins and that base flipping was, and still is, a key element in many DNA-enzyme interactions.