Previous Article | Next Article ![]()
Eukaryotic Cell, January 2006, p. 18-25, Vol. 5, No. 1
1535-9778/06/$08.00+0 doi:10.1128/EC.5.1.18-25.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Molecular and Human Genetics,1 Graduate Program in Structural Computational Biology and Molecular Biophysics,2 Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas,3 Faculty of Computer and Information Science, University of Ljubljana, Ljubljana, Slovenia4
Received 13 September 2005/ Accepted 25 October 2005
|
|
|---|
|
|
|---|
In the past 15 years, it has been accepted that DNA methylation does not occur in Dictyostelium. Smith and Ratner (40) reported that methylation was absent from CCGG sites in repetitive DNA and in the actin multigene family. Drosophila melanogaster was also thought to be an exception for a long time, but recent evidence demonstrated a functional DNA methylation system in Drosophila (14, 25). A small amount of 5mC, consisting of 0.1 to 0.2% of the total cytosine residues, has been detected by methylcytosine antibodies and by high-performance liquid chromatography (HPLC) analysis.
DNA methylation in Drosophila is mediated by the DNA methyltransferase Dnmt2 (32). The Dnmt2 methyltransferase family is highly conserved from yeast to humans, but its trans-methyl activity is rather weak. The biological function of this family is unknown because disruption of the genes caused no readily observable phenotypes in the mouse, yeast, or fly (20, 28, 45). Nevertheless, the evolutionary conservation of this gene family suggests an ancestral origin and an essential function in eukaryotes.
The Dictyostelium genome-sequencing consortium reported an unusual distribution of G+C-rich regions throughout the genome and an underrepresentation of CpG dinucleotides relative to the isomer GpC (9). Such a bias is believed to reflect methylation of cytosine in CpGs, probably because methylated cytosine promotes the mutagenic transition from CpG to TpG. In addition, the Dictyostelium genome sequence revealed the existence of a DNA methyltransferase for 5mC (30). The Dictyostelium DnmA (dictyBase identification no. DDB0231095) is highly similar to other members of the Dnmt2 subfamily. These observations suggested that methylation of cytosine may occur in Dictyostelium and that it may serve as a useful model system for the study of Dnmt2 transmethylases.
We show here that the Dictyostelium genome does contain 5mC, albeit at very low levels. We also show that DNA methylation is developmentally regulated and that deletion of the dnmA gene results in reduced methylation and in developmental defects. We found that CpG dinucleotides have a unique distribution in the Dictyostelium genome and that 5mC residues are found around some of the DIRS transposable elements and in the guaB gene.
|
|
|---|
Purification of genomic DNA and dot blot analysis. Genomic DNA was purified using three methods. The CTAB method (47) was used with minor modifications. Nuclei were lysed in 100 mM EDTA and 5% sodium lauryl sarcosyl at 55°C for 20 min. Genomic DNA was incubated in CTAB solution (1% CTAB, 0.7 M NaCl, 10 mM EDTA, 50 mM Tris-HCl [pH 8.0], 0.5% polyvinylpyrrolidone) at 65°C for 5 min and extracted with chloroform and phenol-chloroform and purified by ethanol precipitation. For the sodium dodecyl sulfate (SDS)-proteinase K method, nuclei were lysed in 1% SDS and 0.2 mg/ml proteinase K for 60 min at 60°C. Genomic DNA was extracted with phenol and phenol-chloroform and purified by ethanol precipitation. For extraction with plant DNAzol (Invitrogen), nuclei from vegetative cells or whole cells from developing cultures were treated according to the manufacturer's recommended protocol. Genomic DNA from all three methods was treated with 100 µg/ml RNase A (Sigma) for 1 h at 37°C. The different methods produced essentially identical results, supporting the notion that the signals observed were not an artifact of the purification method. PCR fragments of the thymidine kinase gene thyB (dictyBase identification no. DDB0191436) were used as positive or negative controls for dot blots. To generate a fully methylated positive control fragment, 5-methyl-dCTP (Roche Applied Science) was used instead of dCTP in the PCR. These DNA fragments were used as standards to quantify the amount of 5mC in the genomic DNA samples (see Fig. S1 in the supplemental material).
Detection and quantification of 5mC in genomic DNA were performed by dot blot analysis with an antibody against 5mC. This sensitive immunological method has been used by others to detect 5mC in Entamoeba histolytica (10) and in liver tumors (42) and is especially suitable for detection and quantification of small amounts and small changes in 5mC levels. Genomic DNA (1 to 25 µg) and PCR fragments were denatured with 0.4 N NaOH at 100°C for 10 min, neutralized with ammonium acetate (pH 7.0), and blotted on a nitrocellulose membrane (Schleicher & Schuell). Dot blots were incubated with a 1:750 dilution of anti-5-methylcytosine antibody (Megabase Research), followed by incubation with a 1:5,000 dilution of alkaline phosphatase-conjugated goat anti-rabbit immunoglobulin G antibody (Sigma), or with a 1:1,000 dilution of anti-5-methylcytosine antibody (Abcam, Inc.), followed by incubation with a 1:8,000 dilution of alkaline phosphatase-conjugated donkey anti-sheep immunoglobulin G (heavy plus light) antibody (Jackson ImmunoResearch). Signals were detected by chemiluminescence with the Western-Star system (Applied Biosystems) according to the manufacturer's recommended protocol. Antibodies were then removed by incubation with 62.5 mM Tris-HCl (pH 6.8), 2% SDS, and 100 mM 2-mercaptoethanol at 50°C for 30 min, and the blots were hybridized with a radiolabeled DNA probe for DIRS-1 (37) and the control probe thyB.
Analysis of DNA methylation by restriction endonuclease digests and quantitative PCR. Dictyostelium genomic DNA was digested with methylation-sensitive (HpaII and/or AvaI) or -insensitive endonucleases (MspI and/or BsoBI) and used as templates for PCRs. PCR primer pairs were designed to amplify regions flanking the restriction sites (HpaII/MspI and AvaI/BsoBI) in the guaB gene locus. Quantitative PCR was performed using the SYBR green dye in an MJ Research Opticon 2 instrument according to the manufacturer's recommended protocol. Standard curves were obtained using a control primer set based on serial dilutions of undigested genomic DNA. End products were resolved on a 1% agarose gel to verify the results. For Southern blot analyses, the genomic DNA was also digested with the XbaI restriction endonuclease and the DNA was analyzed as described previously (17) with the DIRS-1 probe or with the guaB probe.
Computational analysis of CpG dinucleotides in the Dictyostelium genome. We used the "newcpgreport" tool from the EMBOSS package (33) to calculate the observed/expected ratio and to find CpG islands. The definition of CpG islands we used is similar to the one used in mammalian systems with a correction for low G+C content (13). We required the G+C content to be greater than 40% and the observed/expected ratio of CpG to be greater than 0.6. The CpG observed/expected ratios were computed for 100-bp windows, sliding at 1-base increments. The minimal length of a CpG island is 200 bp. We also mapped the codistribution of CpG islands and 124 retrotransposable elements (80 DIRS, 33 SKIPPER, 11 TDD3).
For the computational analysis and corresponding visualizations, we used three different metrics and computed them from the genome sequence of all the genes (13,629) and a subset of 3,745 genes for which expression data are available (43). First, we counted the number of CpG dinucleotides at each position relative to the translational START or STOP sites in the given group of genes, from 10 kb 5' to 10 kb 3' of the site and displayed the total number as a histogram, where the width of the bin is 250 bp. Second, the G+C content (%) at a position in a gene was computed by centering a window of 100 bp and then counting the number of G's and C's in the window. This was done for all the genes and then averaged into a single value that was plotted at the corresponding position in a line graph. The values were computed every 30 bp. Third, the observed/expected ratio was computed in each window of 100 bp. We obtained the observed number by counting the number of CpGs in the window. The following formula as implemented in the "newcpgreport" tool was used to compute the expected number of CpGs: (number of G's x number of C's)/window length. The averages of the values were computed every 30 bp and plotted as a line graph.
|
|
|---|
![]() View larger version (58K): [in a new window] |
FIG. 1. DNA methylation is developmentally regulated. Genomic DNA was prepared from cells at different developmental stages and dot blotted on nitrocellulose membranes in the indicated amounts. (A) DNA methylation was detected by reacting the membranes with an antibody directed against 5-methylcytosine ( 5mC), and the total amount of DNA was estimated by hybridization with a radioactive probe against DIRS-1 (DIRS). (B) The antibody and hybridization signals were quantified, and methylation levels were normalized to the amounts of DNA and plotted as a function of developmental time. The plot indicates that DNA methylation is increased during Dictyostelium development. Results are the averages and standard deviations from 3 replications.
|
Unusual distribution of CpG islands in the Dictyostelium genome. DNA methylation often occurs on cytosines in the dinucleotide CpG. One of the reasons to suspect that Dictyostelium DNA was indeed methylated was the observation that the dinucleotide CpG was underrepresented in the genome (9). In other organisms, methylation of CpG has a profound regulatory role that is associated with the CpG islands (7). CpG islands have been reported around RTEs in other genomes, and the function of methylation at RTEs is a defense mechanism against the expression of these parasitic elements (5, 46). Therefore, we tested whether they were also overrepresented near RTEs in the Dictyostelium genome. Figure 2A shows the distribution of RTEs and CpG islands along the 6 chromosomes and several unassigned contigs in the Dictyostelium genome. We observed that the CpG islands are not distributed evenly throughout the genome and that they are more frequently found around RTEs (Fig. 2A). The most striking correlation is found at the centromeric end of chromosomes 1 and 6 and at the opposite end of chromosome 2 (Dictyostelium chromosomes are telocentric and the map places the centromeres on the left). We also searched for CpG islands using a model that accounts for the bimodal distribution of C+G in coding and noncoding regions, and the results were essentially identical to the ones shown (data not shown).
![]() View larger version (48K): [in a new window] |
FIG. 2. Uneven distribution of CpG dinucleotides in the Dictyostelium genome. CpG islands were identified in the Dictyostelium genome. (A) CpG islands (blue vertical lines) and RTEs (green vertical lines) were plotted along the 6 chromosomes, indicating a high degree of overlap. All of the 13,629 predicted genes in the genome were aligned by their ORFs such that the translational start site (B) or the translational termination site (C) was aligned. CpG dinucleotide frequency was computed in a 20-kb region centered on the respective border of the ORF (open bars). The ratio between the observed and expected frequencies of CpG dinucleotides is shown as a green line. The G+C content (%) is indicated by a red line.
|
Due to the unique nucleotide composition of the genome, Dictyostelium ORFs have significantly higher G+C contents than their flanking and intergenic regions (9). Plotting the average G+C content around the ORF borders revealed the expected drop 5' of the start site, followed by a sharp increase 3' of the start site (Fig. 2B). The end of the ORF exhibited the expected opposite trend with a mild increase 5', a sharp decline immediately 3', and a gradual increase back to the average level following the ORF end (Fig. 2C). We therefore tested the ratio between the observed number of CpG dinucleotides and the number expected from the G+C content. We found that the observed/expected ratio increased sharply immediately 5' of the ORF start site, sharply decreased immediately 3', and gradually returned to average further 3' of the start site (Fig. 2B). At the end of the ORF, the CpG ratio gradually decreased 5' to the end, sharply increased at the end, and then sharply decreased back to the average level 3' to the end of the ORF (Fig. 2C). We tested the distribution of all other dinucleotides but found none that behaved like CpG (data not shown). Finally, we performed the same analysis on about 3,700 genes that were known to be expressed from microarray data (43). The results were essentially identical to those shown for the entire predicted ORFeome (data not shown). These results support the notion that CpG dinucleotides have a special role around the borders of open reading frames.
Evidence for selective DNA methylation in the Dictyostelium genome. The correlation between RTEs and CpG islands (Fig. 2A) prompted us to test whether DNA methylation could be observed specifically in RTEs. We digested genomic DNA from growing and developing cells with restriction endonucleases that are sensitive or resistant to the presence of 5mC and separated the restriction fragments by gel electrophoresis. Examination of the banding pattern after ethidium bromide staining revealed that the methylation-sensitive and the methylation-insensitive enzyme reactions were carried out with the same efficiency (Fig. 3A). We therefore conclude that the reactions were carried out to completion. The DNA was then subjected to Southern blot analysis with a probe against one of the most abundant RTEs, DIRS-1 (6, 9). The data in Fig. 3B show the expected multitude of bands but also show that one of the bands was methylated. Comparing the pattern obtained with the 5mC-sensitive endonucleases in lane 3 with the pattern obtained with the 5mC-resistant endonucleases in lane 4, it is clear that at least one type of DIRS RTEs was protected from digestion by the methylation-sensitive endonuclease. This finding was best observed at 16 h of development (Fig. 3B). It is also likely that most of the DIRS RTEs were not methylated on CpG, since their Southern blot patterns were indistinguishable by this method.
![]() View larger version (116K): [in a new window] |
FIG. 3. Selective DNA methylation in the DIRS retrotransposon and the guaB gene. Genomic DNA was prepared from vegetative (0 h, lanes 1 and 2) and developing (16 h, lanes 3 and 4) cells and digested with the XbaI restriction endonuclease and with restriction endonucleases that are sensitive (both HpaII and AvaI, lanes 1 and 3) or insensitive (both MspI and BsoBI, lanes 3 and 4) to 5mC. Digested DNA was separated by electrophoresis, along with the 1-kb DNA ladder as a size marker, and stained with ethidium bromide to visualize the banding pattern (A). The DNA was analyzed by Southern blotting with the DIRS probe (B) or with the guaB probe (C). Open arrows indicate the full-size fragments, and filled arrows indicate the resulting fragments.
|
To further test the notion that Dictyostelium DNA methylation is rather sparse, we performed quantitative analysis of DNA methylation at the guaB locus. The guaB locus contains 3 putative methylation sites that can be detected by differential restriction endonuclease digestion (Fig. 4A). We designed PCR primer pairs that flank these sites as well as one pair to amplify a control fragment near the 5' border of the ORF. Genomic DNA was digested with the appropriate restriction endonucleases and purified, and quantitative PCRs were performed and normalized (Fig. 4B). We found that site 1 was not methylated, since both types of endonuclease digestion eliminated the ability to generate a PCR product. Both site 2 and site 3 were protected from digestion with 5mC-sensitive endonucleases but not protected from the 5mC-insensitive endonuclease. We therefore conclude that the sites within the guaB ORF are methylated, whereas the site 5' of the ORF is not. The data also support the observation that DNA methylation is developmentally regulated because the phenomenon was more pronounced in DNA from developing cells than in DNA from vegetative cells (Fig. 4B).
![]() View larger version (40K): [in a new window] |
FIG. 4. Quantitative analysis of DNA methylation in the guaB locus. The guaB gene contains several putative methylation sites. (A) A physical map of the guaB locus with the 5' border on the left. The ORF is indicated by a thick gray arrow. Quantitative PCR was done with four primer sets, indicated by small adjacent black arrows pointing at each other. Sites 1 and 3 are restriction sites for the methylation-sensitive endonuclease HpaII and the methylation-insensitive endonuclease MspI. Site 2 is a restriction site for the methylation-sensitive endonuclease AvaI and the methylation-insensitive endonuclease BsoBI. The control pair of primers (cont) was used for normalization. (B) Genomic DNA from vegetative (0 h) and terminally developed (24 h) cells was digested with methylation-sensitive (hatched bars) or -insensitive (white bars) endonucleases. The fraction (% of total) of DNA that remained undigested was measured by quantitative PCR across the respective sites as indicated in panel A and calculated from the average of two replications. Significant and developmentally regulated DNA methylation can be seen in sites 2 and 3.
|
![]() View larger version (27K): [in a new window] |
FIG. 5. Physical maps of the dnmA locus. (A) Restriction maps of the dnmA locus in the wild-type strain (AX4) and in the dnmA knockout strain (KO). Probe 1 consists of the dnmA gene, and probe 2 consists of the plasmid DNA used to replace the gene in the knockout mutant. Open boxes and arrowhead indicate the dnmA exons. Restriction sites are indicated as follows: EV, EcoRV; EI, EcoRI; Cl, ClaI. (B) Genomic DNA from the wild-type strain (AX4) and from the dnmA knockout strain (KO) was digested as indicated and subjected to Southern blot analysis with probes 1 and 2.
|
![]() View larger version (33K): [in a new window] |
FIG. 6. Reduced global DNA methylation in the dnmA knockout strain. Genomic DNA was extracted from vegetative (0 h) and developed (24 h) cells of the wild-type (AX4) and dnmA knockout (KO) strains. DNA methylation was detected by dot blotting with the anti-5mC antibody as described in the legend to Fig. 1. (A) DNA methylation signal intensity during growth (black bar) and development (white bar) are the averages of results from two independent replications of each of two independent knockout strains (KO) and the averages of results from three independent replications of the wild-type strain, showing a reduction in the developmentally regulated methylation. (B) Representative dot blot data comparing the dnmA knockout strain (KO) and the wild-type strain (AX4) signals in growth and development.
|
![]() View larger version (54K): [in a new window] |
FIG. 7. dnmA is responsible for selective DNA methylation at the guaB locus. Genomic DNA was prepared from developed wild-type (AX4) and dnmA knockout (KO) cells. The DNA was digested with a methylation sensitive (S, AvaI) or a methylation-insensitive (IS, BsoBI) endonuclease and analyzed by PCR using the primer set indicated in the legend to Fig. 4. PCR products were resolved on a 1% agarose gel containing ethidium bromide and visualized by fluorescence under UV light. Photographs are shown, and control PCR amplification products are shown below the experimental lanes.
|
dnmA activity is essential for proper culmination. Gross morphological aberrations are a good measure of developmental function in Dictyostelium. We therefore developed the dnmA-knockout mutant side by side with wild-type cells and compared their morphologies. The mutant and wild-type morphological and temporal progression through development were essentially identical for the first 20 h of development (data not shown). During culmination however, the dnmA knockout mutant exhibited subtle but reproducible defects. Many of the mutant sori were fragmented along the stalk, and they were clearer than the wild-type sori (Fig. 8). Clear sori sometimes indicate inefficient sporulation or precocious germination, and we indeed found that the mutant formed fewer spores than the wild type (data not shown). However, the sporulation phenotype was less penetrant and less reproducible than the sorus fragmentation phenotype. Overall, it is clear that dnmA is essential for proper terminal development, consistent with the temporal increase in DNA methylation during development.
|
View larger version (20K): [in a new window] |
FIG. 8. dnmA is essential for proper development. Wild-type (AX4) and dnmA knockout (KO) cells were starved and developed on black nitrocellulose filters for 24 h. In the mutant cells, sori were found at various points along the stalk, whereas the wild type has one sorus at the top of the stalk. Mutant sori were also more transparent than wild-type sori.
|
|
|
|---|
The function of the Dnmt2 family of DNA methyltransferases is somewhat enigmatic. In most organisms, these enzymes have little or no activity, and mutations in the respective genes have almost no consequences (20, 28, 45). DNA methylation in Dictyostelium is fairly rare, in that only 0.03 to 0.14% of the total cytosines are methylated. The fraction of methylated cytosines is higher in other organisms with about 0.1% in Drosophila melanogaster, 2 to 10% in mammals, and more than 30% in some plants (2, 14). In addition, the G+C content in the Dictyostelium genome is very low (22.43%), so there is just one methylated cytosine per 7 to 30 kb in the genome. This low level of DNA methylation is consistent with the weak activity of Dnmt2 enzymes. Nevertheless, mutating dnmA resulted in a subtle yet obvious developmental defect such as fragmentation of the sori midway along the stalk. The positioning of the sorus to the top of the stalk is likely to have a selective advantage in spore dispersal. Therefore, a mutant lacking dnmA would probably have a competitive disadvantage against wild-type cells. These findings provide a demonstration of function for a Dnmt2 DNA methyltransferase. Another recent example is the finding that Drosophila Dnmt2 might be involved in longevity and aging (22).
Dictyostelium DnmA contains all 10 conserved motifs, including the methyl donor (S-adenosyl-L-methionine) binding motif and the active site, and an invariant polypeptide in the target-recognizing domain, TRD, which is the putative DNA recognition site (21). The latter domain is not found in other Dnmt family members, so it is thought that Dnmt2 recognizes a specific kind of target through this TRD domain (8). In this paper, we found that the Dictyostelium DnmA functions as one of the DNA methyltransferase and methylates its targets in a selective manner. In D. melanogaster, most methylation was seen at CpT, and only a small fraction of the 5mC was detected in CpG dinucleotides (25). Overexpression of Dnmt2 in D. melanogaster increased the methylation levels of cytosine in the nonsymmetrical CpT and CpA dinucleotides (20). It is therefore conceivable that the cytosines in CpT and CpA are also recognized by the Dictyostelium DnmA.
The low level of DNA methylation and the developmental regulation explain why previous studies have failed to find it (40). A random choice of probes for Southern blot analysis would have to be very lucky to reveal a methylated gene, and physical methods like HPLC analysis are limited to fairly high proportions of 5mC. We have tried to analyze Dictyostelium DNA by HPLC, followed by mass spectrometry, but found that the proportion of 5mC in the genome was below the limit of detection of the method (data not shown). We were fortunate that the antibody we used is highly sensitive and highly specific to 5mC in DNA. We also had the advantage of knowing the genome sequence as a guide in searching for good probes for Southern blot analysis.
DNA methylation in Dictyostelium is developmentally regulated, and the highest degree of methylation is observed in cells at the end of development. It is tempting to speculate that the timing of highest methylation and the late morphological phenotype of the dnmA mutant are causally correlated. It is possible that the methylation of a small subset of genes is required for regulation of late gene expression and that disruption of the process leads to developmental defects. Another possibility is that the correlation is indirect. Terminal morphology is the result of proper execution of all developmental pathways. Therefore, we cannot rule out the possibility that the dnmA knockout phenotype is the result of subtle defects in processes that are executed early in development and are manifested only later on.
The unusual distribution of CpG islands in the Dictyostelium genome is also quite interesting. The association of CpG islands with RTEs in other organisms is usually attributed to the silencing effects of DNA methylation (5, 46). It is reasonable to hypothesize that silencing RTE transcription by DNA methylation may have a selective advantage in maintaining genome stability and reducing the metabolic burden associated with expressing the RTE genes. This assertion may have to be reexamined in light of our findings. In Dictyostelium, CpG islands are also associated with RTEs, but most of them are probably not methylated because we could not detect methylation around most DIRS1 RTEs. It is therefore possible that the CpG islands have a role that is independent of their methylation in addition to their proven, methylation-dependent role in the regulation of chromatin structure and gene expression. In that context, CpG dinucleotides have an unusual distribution around the borders of ORFs. It is conceivable that most of them are not methylated in Dictyostelium, so we propose that CpGs and CpG islands have a role that is independent of methylation. What that role might be is currently a matter of speculation, but the fact is that the unusual distribution is also found in other organisms (3, 38).
Finally, we did not find evidence for the role of DNA methylation in the regulation of gene expression in Dictyostelium, but we only examined one gene. A more comprehensive study might reveal whether DNA methylation is important for gene regulation and which genes and pathways are regulated by that mechanism.
In an independent study, Nellen and coworkers have detected 5mC DNA methylation of CpA and CpT dinucleotides. In their hands, deletion of dnmA did not result in morphological defects, but their parental strain (AX2) was different from the one we used (AX4). These researchers found that DNA methylation was necessary for transcriptional silencing of the Skipper RTE and showed evidence for the correlation of DNA methylation and RNA interference (19a). Together, our data establish that Dictyostelium cells utilize DNA methylation for regulation of gene expression and of development.
Work at BCM was supported by a grant from the National Institute of Child Health and Human Development (PO1 HD39691). M.K. is a recipient of a research fellowship from the Uehara Memorial Foundation. Work at University of Ljubljana was supported by a Program Grant (P2-0209) and a Scientific Research Collaboration Grant between the United States and the Republic of Slovenia, both from the Slovenian Research Agency.
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»