Epigenetic Segregation of Microbial Genomes from Complex Samples Using Restriction Endonucleases HpaII and McrB

We describe continuing work to develop restriction endonucleases as tools to enrich targeted genomes of interest from diverse populations. Two approaches were developed in parallel to segregate genomic DNA based on cytosine methylation. First, the methyl-sensitive endonuclease HpaII was used to bind non-CG methylated DNA. Second, a truncated fragment of McrB was used to bind CpG methylated DNA. Enrichment levels of microbial genomes can exceed 100-fold with HpaII allowing improved genomic detection and coverage of otherwise trace microbial genomes from sputum. Additionally, we observe interesting enrichment results that correlate with the methylation states not only of bacteria, but of fungi, viruses, a protist and plants. The methods presented here offer promise for testing biological samples for pathogens and global analysis of population methylomes.     

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.     

Methylation pattern of mouse mitochondrial DNA.

The pattern of methylation of mouse mitochondrial DNA (mtDNA) was studied using several techniques. By employing a sensitive analytical procedure it was possible to show that this DNA contains the modified base 5-methylcytosine (m5Cyt). This residue occurred exclusively at the dinucleotide sequence CpG at a frequency of 3 to 5%. The pattern of methylation was further investigated by determining the state of methylation of several MspI (HpaII) sites. Different sites were found to be methylated to a different extent, implying that methylation of mtDNA is nonrandom. Based on the known base composition and nucleotide sequence of mouse mtDNA, the dinucleotide sequence CpG was found to be underrepresented in this DNA. The features of mtDNA methylation (CpG methylation, partial methylation of specific sites and CpG underrepresentation) are also characteristic of vertebrate nuclear DNA. This resemblance may reflect functional relationship between the mitochondrial and nuclear genomes.     

Identification of differentially methylated sites within unmethylated DNA domains in normal and cancer cells

Altered DNA methylation has been linked to neoplastic cell transformation and is a hallmark of cancer progression. Therefore, the screening for differentially methylated sequences as tumor biomarkers has a significant implication in the clinical setting. To determine the cancer-linked alterations in DNA methylation pattern, we have applied an endonuclease, McrBC, to the existing methylation-sensitive arbitrarily primed polymerase chain reaction (msAP-PCR) method and developed McrBC-msAP-PCR. This modified approach allows detection of differentially methylated sites within unmethylated DNA domains enriched by regulatory sequences and CpG islands. In this method, we used digestion of DNA with the McrBC methylation-sensitive endonuclease to selectively exclude the methylated fraction of DNA, which comprises interspersed and tandem-repeated sequences and exons other than first exons, from analysis. The subsequent digestion of unmethylated DNA fragments with SmaI and HpaII methylation-sensitive restriction endonucleases followed by AP-PCR amplification resulted in the detection of unknown unique sequences associated with cancer-linked methylation changes in genomic DNA. Hypermethylation and hypomethylation are visualized by the increase or decrease in the band intensity of DNA fingerprints. By using this technique, we were able to differentiate clearly, identify, and characterize a number of novel unique DNA sequences with differentially methylated sites in normal and breast cancer cell lines and in normal and rat tumor liver tissues.     

Differential methylation of the c-H-ras gene in normal mouse cells and during skin tumour progression.

We have previously shown that the mouse c-H-ras gene acquires transforming activity in chemically induced skin tumours. We have now investigated the pattern of DNA methylation at HpaII and XhoI sites around the c-H-ras locus in various tissues and stages of epidermal tumour progression. The results of this study suggest a correlation between the methylation state of the c-H-ras gene and its susceptibility to oncogenic conversion by a point-mutation. The locus is substantially undermethylated in normal epidermis in comparison with NIH/3T3 fibroblasts. Intermediate levels of methylation were observed in the other tissues investigated. The undermethylation at HpaII sites in epidermal DNA persists through the morphologically distinct phases of hyperplasia, benign papilloma and malignant carcinoma. Methylation at a specific XhoI site close to the c-H-ras gene is significantly reduced with respect to normal epidermis in some, but not all epidermal tumours. The methylation state of the c-H-ras locus in specific tumours is stably maintained following transfection of these DNAs into NIH/3T3 cells and selection of transformed foci. Demethylation of the locus is not essential in vitro for the transforming activity of DNA from epidermal tumours. The significance of changes in the methylation pattern of the c-H-ras gene in different tissues and during tumour progression is discussed.     

Assessing the DNA methylation status of single cells with the comet assay

The comet assay (single cell gel electrophoresis) is a cost-effective, sensitive, and simple technique that is traditionally used for analyzing and quantifying DNA damage in individual cells. The aim of this study was to determine whether the comet assay could be modified to detect changes in the levels of DNA methylation in single cells. We used the difference in methylation sensitivity of the isoschizomeric restriction endonucleases HpaII and MspI to demonstrate the feasibility of the comet assay to measure the global DNA methylation level of individual cells. The results were verified with the well-established cytosine extension assay. We were able to show variations in DNA methylation after treatment of cultured cells with 5-azacytidine and succinylacetone, an accumulating metabolite in human tyrosinemia type I.     

Caenorhabditis elegans DNA does not contain 5-methylcytosine at any time during development or aging.

DNA, isolated from age-synchronous senescent populations of Caenorhabditis elegans has been quantitatively and qualitatively analyzed for the presence of 5-methylcytosine. High performance liquid chromatography on two wild-type and several mutant strains of C. elegans failed to detect any 5-methylcytosine. The restriction endonuclease isoschizomers, HpaII and MspI, were used to digest genomic DNA after CsCl purification and failed to detect any 5' cytosine methylation at any age. We conclude that C. elegans does not contain detectable (0.01 mole percent) levels of 5-methylcytosine.     

DNA methylation in mouse cells in culture as measured by restriction enzymes

Methylation of DNA in normal mouse cultured 3T3 cells and in their virally or chemically transformed derivatives was studied. DNA methylation was studied by restriction with HpaII, MspI, or HpaII plus MspI. DNA from the chemically transformed cells was cleaved about twice as often with HpaII than was the DNA of normal and virally transformed cells. Digests with MspI and HpaII plus MspI were identical in all cell lines studied. Densitometry of the restriction patterns allowed an estimate of total DNA methylation from the weight average lengths. The chemically transformed cell line showed 25% reduction in methylation compared to the other cell lines. Southern blot hybridization using satellite DNA showed that these sequences followed a pattern of modification similar to that of total DNA.     

Analysis of methylation and distribution of CpG sequences on human active and inactive X chromosomes by in situ nick translation

In situ nick translation of fixed mitotic chromosomes after HpaII or MspI digestion allows us to detect different DNA methylation levels along chromosomes. We used this technique to analyse the methylation levels of CCGG sites in the active and inactive X chromosomes of female human cells. In addition, we analysed the distribution of these sites with respect to the banding pattern. Our data show that the inactive X, as a whole, is more methylated than the active one and that CCGG sequences are preferentially located on R-positive bands.     

The adenovirus type 12 - mouse cell system: permissivity and analysis of integration patterns of viral DNA in tumor cells.

The integration patterns of persisting adenovirus type 12 (Ad12) DNA were analyzed in two Ad12-induced tumors of Balb/c and CBA/J mice and in one tumor cell line derived from an Ad12-induced retinoblastoma of C3H origin. In all three tumors the Ad12 genome was integrated colinearly and various copy numbers of viral DNA were found. Analysis of the Ad12 integration patterns revealed relatively simple offsize band patterns regardless of Ad12 copy numbers. The degree of methylation at the 5''-CCGG-3'' sites in the inserted Ad12 genome was determined using the isoschizomeric restriction endonuclease pair HpaII and MspI. Methylation was rather incomplete in the primary tumor tissues but almost complete in the retinoblastoma line carried in culture for many passages. The levels of expression of the viral genome in the Balb/c tumor and in the retinoblastoma line were determined by in vitro translation of RNA isolated from these cells and selected with appropriate restriction endonuclease fragments of Ad12 DNA. In both instances the 59 K, 19 K, and 17 K proteins of the E1b region were expressed. Proteins of the E1a region appeared very faint in the size class between 22 K and 42 K. The permissivity of Ad12 and the replication of Ad12 DNA in mouse cells were investigated by blotting restricted DNA from cells soon after, and a long time after, infection and by hybridization with 32P-labeled Ad12 DNA. Neither primary mouse kidney cells nor the established L929 mouse cell line supported viral DNA replication. These results raise the question to what extent host cell factors determine Ad12 DNA replication in mammalian cells.     

Characterization of DNA methylation in the rat

In the rat, differentiation and cell proliferation both affect DNA methylation. We studied 5-methylcytosine at the inner cytosine of the sequence C-C-G-G, a common methylation site, using endonuclease MspI (which cleaves C-C-G-G- and C-mC-G-G), and its isoschizomer HpaII (which cleaves only C-C-G-G). DNA from all tissues and cell lines studied was methylated at C-C-G-G, at levels ranging from 45 to 80%, but the methylation sites were not distributed uniformly. Our analysis suggests a model in which cells contain variable amounts of three DNA methylation states, averaging 30-40, 70-80 and 95-100% methylation, respectively. One biological parameter that alters methylation is the proliferative state of the cell. We observed that NRK, a non-transformed cell line, increased its DNA methylation from 45 to 67% when monolayer cultures became confluent and non-dividing. We also observed that a class of repetitive DNA was completely methylated in DNA from all sources except a transformed cell line.     

Analysis of DNA methylation in various swine tissues

DNA methylation is known to play an important role in regulating gene expression during biological development and tissue differentiation in eukaryotes. In this study, we used the fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP) method to assess the extent and pattern of cytosine methylation in muscle, heart, liver, spleen, lung, kidney and stomach from the swine strain Laiwu, and we also examined specific methylation patterns in the seven tissues. In total, 96,371 fragments, each representing a recognition site cleaved by either or both EcoRI + HpaII and EcoRI + MspI, the HpaII and MspI are isoschizomeric enzymes, were amplified using 16 pairs of selective primers. A total of 50,094 sites were found to be methylated at cytosines in seven tissues. The incidence of DNA methylation was approximately 53.99% in muscle, 51.24% in the heart, 50.18% in the liver, 53.31% in the spleen, 51.97% in the lung, 51.15% in the kidney and 53.39% in the stomach, as revealed by the incidence of differential digestion. Additionally, differences in DNA methylation levels imply that such variations may be related to specific gene expression during tissue differentiation, growth and development. Three types of bands were generated in the F-MSAP profile, the total numbers of these three types of bands in the seven tissues were 46,277, 24,801 and 25,293, respectively.In addition, different methylation patterns were observed in seven tissues from pig, and almost all of the methylation patterns detected by F-MSAP could be confirmed by Southern analysis using the isolated amplified fragments as probes. The results clearly demonstrated that the F-MSAP technique can be adapted for use in large-scale DNA methylation detection in the pig genome.     

Genome-wide DNA methylation analysis reveals phytoestrogen modification of promoter methylation patterns during embryonic stem cell differentiation

Background

Environmental challenges during development affect the fetal epigenome, but the period(s) vulnerable to epigenetic dysregulation is(are) not clear. By employing a soy phytoestrogen, genistein, that is known to alter the epigenetic states of the A^vy allele during embryogenesis, we have explored the sensitive period for epigenetic regulation. The post-implantation period, when de novo DNA methylation actively proceeds, is amenable to in vitro analysis using a mouse embryonic stem (ES) cell differentiation system.

Methods and Findings

Mouse ES cells were differentiated in the presence or absence of genistein, and DNA methylation patterns on day 10 were compared by microarray-based promoter methylation analysis coupled with a methylation-sensitive endonuclease (HpaII/McrBC)-dependent enrichment procedure. Moderate changes in methylation levels were observed in a subset of promoters following genistein treatment. Detailed investigation of the Ucp1 and Sytl1 promoters further revealed that genistein does not affect de novo methylation occurring between day 0 and day 4, but interferes with subsequent regulatory processes and leads to a decrease in methylation level for both promoters.

Conclusion

Genistein perturbed the methylation pattern of differentiated ES cells after de novo methylation. Our observations suggest that, for a subset of genes, regulation after de novo DNA methylation in the early embryo may be sensitive to genistein.     

DNA methylation: organ specific variations in the methylation pattern within and around ovalbumin and other chicken genes.

The restriction enzymes HhaI and HpaII, whose activity is inhibited by cytosine methylation within their recognition sites, have been utilised as probes to study methylation in the vicinity of the ovalbumin gene in DNA from various chicken tissues. This was complemented by a preliminary study of methylation in the regions of chicken ovotransferrin (conalbumin), ovomucoid and beta-globin genes. From our data we conclude that HaI or HpaII sites can be divided in 3 classes according to their pattern of methylation in different tissues. In the first class of sites (mV class) the extent of methylation varies in different tissues. The patterns obtained show that methylation at the sites located within and around the 3 genes which code for egg white proteins is in general lowest in oviduct of laying hen, where these genes are expressed. However some sites are not methylated (m- class) and others are 95 to 100% resistant (m+ class) to digestion by HhaI or HpaII in the DNAs of all the tissues which were tested. Our study has also revealed a remarkable number of allelic variants for the presence of HhaI or HpaII sites in the region of the ovalbumin gene.     

In situ methylation of insect chromosomes with methylase HpaII

A method of in situ DNA methylation with the prokaryotic methylase HpaII has been developed on fixed mitotic and meiotic chromosomes of the insect species Baetica ustulata. Incorporation of methyl groups into the chromosomal DNA is revealed by autoradiography using a labelled substrate and by its ability to prevent endonuclease digestion. The method allows direct visualization of clusters of methylatable CCGG sites. The distribution of these clusters in the chromosome complement of Baetica shows two separate domains of heterochromatic DNA which differ in their methylation patterns. Each is distributed at equivalent locations in both homologous and nonhomologous chromosomes. The existence of two compartments, one methylated and the other unmethylated, in the heterochromatic DNA could be interpreted as a remnant of the ancestral echinoderm-like pattern of methylation.     

Functional consequences of haem orientational disorder in sperm-whale and yellow-fin-tuna myoglobins.

With the use of the isoschizomeric restriction endonucleases HpaII and MspI, we found that mouse tumour ornithine decarboxylase (ODC; EC 4.1.1.17) genes are extensively methylated. ODC genes in L1210 mouse leukaemia cells were apparently more methylated than in Ehrlich ascites carcinoma, as revealed by the use of HpaII endonuclease, yet the digestion of genomic DNA isolated from these two murine tumour cell lines with MspI, which cleaves at a CCGG sequence, also with internally methylated cytosine, resulted in an apparently identical restriction pattern. It is possible that the amplification of ODC genes in Ehrlich ascites-carcinoma cells in response to 2-difluoromethylornithine (DFMO) was associated with hypomethylation, or that less-methylated genes were amplified. A human myeloma (Sultan) cell line only revealed three separate hybridization signals when cleaved with HpaII. One of these signals was amplified under the pressure of DFMO. When cleaved with MspI, these three HpaII fragments disappeared and were replaced by a double signal of 2.3-2.4 kilobase-pairs (kbp) in size. The amplified ODC sequences in the Sultan myeloma cell line apparently originated from chromosome 2, as indicated by a unique hybridization signal in a 5.8 kbp HindIII fragment specific for the human ODC locus on chromosome 2. A comparison of different human cells, the Sultan myeloma, a lymphocytic B-cell leukaemia (Ball), normal mononuclear leucocytes and leucocytes obtained from leukaemia patients, revealed interesting differences in the methylation of ODC genes. The use of two restriction endonucleases (HpaII and CfoI), the cleavage site for both of which contains a CG sequence and which only cleave when cytosine is unmethylated, indicated that ODC genes in the lymphocytic leukaemia cells were much less methylated than those in the normal leucocytes or in the Sultan cells.     

In vitro methylation of specific regions of the cloned Moloney sarcoma virus genome inhibits its transforming activity.

The transforming activity of cloned Moloney sarcoma virus (MSV) proviral DNA was inhibited by in vitro methylation of the DNA at cytosine residues, using HpaII and HhaI methylases before transfection into NIH 3T3 cells. The inhibition of transforming activity due to HpaII methylation was reversed by treatment of the transfected cells with 5-azacytidine, a specific inhibitor of methylation. Analysis of the genomic DNA from the transformed cells which resulted from the transfection of methylated MSV DNA revealed that the integrated MSV proviral DNA was sensitive to HpaII digestion in all cell lines examined, suggesting that loss of methyl groups was necessary for transformation. When cells were infected with Moloney murine leukemia virus at various times after transfection with methylated MSV DNA, the amount of transforming virus produced indicated that the loss of methyl groups occurred within 24 h. Methylation of MSV DNA at HhaI sites was as inhibitory to transforming activity as methylation at HpaII sites. In addition, methylation at both HpaII and HhaI sites did not further reduce the transforming activity of the DNA. These results suggested that; whereas methylation of specific sites on the provirus may not be essential for inhibiting the transforming activity of MSV DNA, methylation of specific regions may be necessary. Thus, by cotransfection of plasmids containing only specific regions of the MSV provirus, it was determined that methylation of the v-mos gene was more inhibitory to transformation than methylation of the viral long terminal repeat.     

Developmental regulation of cytosine methylation in the nuclear ribosomal RNA genes of Pisum sativum

Prominent features of the cytosine methylation pattern of the Pisum sativum nuclear ribosomal RNA genes have been defined. Cytosine methylation within the C-C-G-G sequence was studied using the restriction enzymes HpaII and MspI and gel blot hybridizations of the restriction digests. The extent to which particular features of the methylation pattern change during seedling development has also been determined. Total cellular DNA, purified from defined sections of pea seedlings grown under different lighting conditions, was analyzed with DNA hybridization probes derived from different portions of a cloned member of the nuclear rRNA gene family. By use of an indirect end-labeling technique, a map of 23 cleavable HpaII and/or MspI sites in genomic rDNA was constructed. The map covers about 90% of the rDNA repeat including the entire non-transcribed spacer region and most of the rRNA coding sequences. One notable feature of the map is that the most prominent HpaII site, located about 800 base-pairs upstream from the 5' end of the mature 18 S rRNA, is cleaved only in one of the two most abundant rDNA length variants (the short variant). With a gel blot assay specific for cleavage at this site, we estimated the HpaII sensitivity of DNA preparations from several stages of pea seedling development. We find that, while methylation is generally low in young seedlings, DNA obtained from the apical buds of pea seedlings is highly methylated. Further, the methylation level of rDNA within the pea bud decreases as the buds are allowed to develop under continuous white light. Our data, taken together with published studies on pea seedling development, indicate that cytosine methylation levels may be related to the regulated expression of the nuclear rRNA genes in pea.