| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Department of Urology and Department of Biochemistry and Molecular Biology, University of Southern California, 1441 Eastlake Ave., Rm. 5420, Norris Cancer Center, Mail Stop 73, Los Angeles, California 90033
Received 7 December 2006/ Returned for modification 5 January 2007/ Accepted 11 January 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Although it has been demonstrated that DNA methylation affects transcription elongation more than initiation in fungi (1, 26), the impact of DNA methylation on transcription may be more complex in mammalian cells. Two studies using a recombinase-mediated transgenic system in mammalian cells proposed that DNA methylation may trigger a closed chromatin structure that blocks transcription initiation when the promoter and coding regions are both methylated (29); additionally, DNA methylation in the intragenic region may lead to histone hypoacetylation and decrease the efficiency of transcription elongation (19). Other studies using nonreplicating plasmids have indicated that the size of the methylated region (14) and methylation adjacent to the promoter (12, 13, 22) may inhibit transcription. While some important implications were made in these studies, the impact of the methylated DNA segments on transcription is difficult to interpret because the size and location of the methylated DNA region may affect transcription differently. While a methylated coding region has little impact on transcription initiation when it is about 1 kb downstream from the promoter, as observed by Lorincz et al. (19), it is conceivable that the same methylated sequence might have a much larger impact on transcription initiation if it were immediately downstream from the promoter.
In Saccharomyces cerevisiae and Drosophila melanogaster, the level of methylated histone H3 at lysine 4 was found to be enriched at the active gene sites (2, 30). Dimethylated and trimethylated histone H3 at lysine 4 (H3K4me2 and H3K4me3, respectively) appeared to be linked to active transcription in chicken erythrocytes and embryos (28, 31) and in human cells (17). H3K4me2 has been found to be present around the active promoters where TFIID and RNA polymerase were also present (15), as well as in the coding sequences of human cells (20). H3K4me2 and H3K4me3 were found at 46% and 66%, respectively, of the start sites of genes on human chromosomes 21 and 22 (3). The presence of H3K4me2 at some locations also appears to be cell type specific, indicating the association of H3K4me2 with markers of active chromosomal regions (3).
Studies of the effect of DNA methylation on local chromatin structure have been more limited. Levels of both H3K4me2 and H3K4me3 were found to be two- to threefold lower in the methylated coding region of the transgene where the depletion of RNA polymerase II (Pol II) was also observed (19). A lack of correlation among gene expression, histone modification, and DNA methylation has been described for some imprinted genes (36). DNA methylation and histone modification are two factors among many in the different pathways and mechanisms for gene suppression.
The present study was designed to determine whether DNA methylation dictates a closed chromatin structure that influences transcription initiation and elongation, depending on the location of the DNA methylation. Methylated DNA patches on the episomes are stably maintained for long intervals after transfection into human cells, as described previously (9, 11). We generated episomes with methylated DNA patches in the promoter, in the reporter coding region, in both the promoter and the reporter coding region, in the first 600 bp of the reporter coding region, and in the last 1.2 kb of the reporter coding region to study the impact of DNA methylation on chromatin structure and gene expression. The sizes of the DNA methylation patches in the promoter alone and in the first 600 bp of the coding region are similar, so we could distinguish whether the differences between methylation in the promoter and in the coding region were due to the size or the location of the methylation. The episome with DNA methylation in the first 600 bp of the coding region and the episome with DNA methylation in the last 1.2 kb of the coding region allowed us to discern whether the location of the DNA methylation patch in the coding region affected transcription differently. We analyzed the presence of H3K4me2, H3K4me3, and Pol II in various regions in correlation with the reporter gene expression of these episomes in human cells.
We found that the size of the DNA methylation patch alone does not determine the extent of transcriptional repression because methylated DNA patches of similar sizes impact transcriptional elongation more than transcriptional initiation. The presence of H3K4me2 correlates with the status of DNA methylation and clearly not with gene expression or Pol II presence. While the level of H3K4me3 also closely correlates with the status of DNA methylation, two exceptions were observed, indicating a contributory role of transcription to the presence of this modified histone. The findings here are consistent with previous findings that methylated H3K4 are enriched in actively transcribed regions. These findings further extend the understanding of how DNA methylation affects the promoter and the coding region. We conclude that DNA methylation in the promoter affects transcription initiation and that DNA methylation has a much larger impact on transcription elongation than on initiation. We also conclude that DNA methylation plays an essential role in the exclusion of methylated H3K4 and dictates a closed chromatin structure. Furthermore, the association of histone modification with DNA methylation has very distinct boundaries that do not spread into regions only 200 bp away. We propose that DNA methylation at the junction of transcriptional initiation and elongation is most critical in transcription suppression and that this effect is mechanistically mediated through chromatin structure.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Nomenclature of plasmids. Plasmids with methylation patches of 518 bp to 2.4 kb in length at different locations of the transcription unit were generated using the unique restriction sites available. Plasmids were designated with the addition of "Me" to the name of the patch that is methylated. These include pMeLTR, pMeLuc, pMeLTR/Luc, pMe5'Luc, and pMe3'Luc. The fully unmethylated plasmid is designated pCLH22, and the fully methylated plasmid is designated pMeCLH22.
In vitro DNA methylation. Plasmid or DNA fragments were methylated in vitro at all CpG sites by using SssI methylase (New England Biolabs) under the conditions recommended by the manufacturer. Following methylation, DNA was dialyzed with nitrocellulose filters (Millipore) after phenol-chloroform extraction. The methylation status of each DNA fragment was confirmed by HhaI restriction enzyme digestion.
In vitro ligation. The vector preparations were test ligated without the insert to identify any uncut or vector religation products. Only preparations with a background less than 0.01% were used in constructing the patch-methylated plasmids. Methylated and unmethylated portions of the plasmid were ligated in vitro to generate patch-methylated constructs. After ligation, the DNA was phenol-chloroform extracted and dialyzed with nitrocellulose filters. If the Rous sarcoma virus long terminal repeat (RSVLTR) promoter and the luciferase gene remained intact on the unligated vector, these unligated molecules could possibly recircularize in mammalian cells after transfection, thereby generating undesired luciferase expression. Therefore, exonuclease V (US Biochemical) treatment was carried out after ligation of those constructions to eliminate any unligated DNA. Ligated DNA was quantitated by transformation of Escherichia coli by using a known quantity of pCLH22 as a standard.
Cell line and transfection. The calcium phosphate transfection method (9, 35) was used to transfect 293/EBNA1 cells (9). All transfections were done at least in duplicate in each experiment. All experiments to determine luciferase expression and plasmid quantitation were performed at least twice, using products from independent ligations for confirmation. As described previously (11), each transfection was harvested at least three consecutive times over a 3-week interval, at least 1 week after transfection, to assess any possible changes occurring over time. For chromatin immunoprecipitation (ChIP) assays, transfected cells were selected with 200 µg/ml hygromycin, with the exception of cells transfected with pMeCLH22, which do not express the hygromycin resistance gene due to methylation.
Episome recovery and analysis. As described previously, 1.25% of the cells were harvested for the luciferase assay; 2.5% of the cells were replated on a 100-mm-diameter tissue culture plate; and the remaining cells were harvested by the Hirt method to recover episomal DNA each time transfected cells reached confluence (9). The same fraction (10%) of the DNA harvested from each transfection was digested with XbaI to linearize the plasmid for DNA quantitation and double digested with XbaI and HhaI to determine the methylation status using Southern blotting. Southern blots were exposed to a phosphorimaging screen, and the radioactivity was quantitated with a PhosphorImager (Bio-Rad GS525 or Bio-Rad FX). DNA quantitation was also confirmed by real-time quantitative PCR (Q-PCR) using at least two different TaqMan probe and primer sets.
Luciferase assay. An aliquot (1.25%) of transfected cells was harvested and lysed for the luciferase assay. Luciferase activities were analyzed using a Monolight 2020 luminometer (Analytical Luminescence) as described previously (9). After subtraction of the background reading, the luciferase reading was normalized for the amount of DNA from the same harvest, as determined by Southern blotting and/or Q-PCR as described previously (9-11). This ensured that the luciferase expression from the same quantity of plasmids with different methylation states from different transfections were being compared.
Chromatin immunoprecipitation. ChIP assays were performed by a modification of the protocol described by Braunstein et al. (4) and as described previously (11) with minor modifications. Exponentially growing cultures of 293/EBNA1 cells transfected with patch-methylated stable episomes were fixed with 1% formaldehyde for 10 min at room temperature. Approximately 5 x 106 fixed cells were resuspended in 1 ml radioimmunoprecipitation assay (RIPA) buffer (150 mM NaCl, 50 mM Tris-HCl [pH 8.0], 0.5% sodium deoxycholate, 1% Nonidet P-40, 0.1% sodium dodecyl sulfate [SDS]) that contained mammalian protease inhibitors (Roche). The cell suspension was sonicated and centrifuged as described previously (11). One 120-µl aliquot of soluble chromatin was used for each immunoprecipitation, and one aliquot was reserved as the total chromatin fraction (TCF) without immunoprecipitation steps for Q-PCR analysis. For immunoprecipitation, 1 µg each of anti-dimethyl-histone H3 (Lys4) and anti-trimethyl-histone H3 (Lys4) (Upstate Biotechnology and Abcam) were mixed with 120 µl of soluble chromatin, 6 µg of sheared salmon sperm DNA, and 20 µl of protein G-Sepharose slurry (in RIPA buffer) with rotary mixing at 4°C, overnight. Immunoprecipitation using antibodies against RNA polymerase II (clone 8WG16; Covance) was carried out as described above but mixed with 240 µl of soluble chromatin instead of 120 µl. An aliquot of soluble chromatin was processed without antibodies (no antibody [Ab]) in parallel as negative controls for the immunoprecipitation experiments. The Sepharose beads were collected by centrifugation and washed sequentially for 10 min each at 4°C, once with 0.5 ml RIPA buffer and twice with 0.5 ml phosphate-buffered saline. Immunocomplexes were eluted twice by incubation with 1% SDS in 0.1 M NaHCO3 for 30 min at room temperature, each time with rotary mixing. Eluates were incubated at 65°C for 5 h after the addition of NaCl to a 200 mM final concentration to reverse the formaldehyde cross-links. The TCF aliquot was included from the reverse cross-linking step and processed as the other immunoprecipitation samples. After cross-links were reversed, EDTA was added to a final concentration of 10 mM, Tris-HCl (pH 6.5) was added to a final concentration of 40 mM, and 4 µg proteinase K was added to each tube. After overnight incubation at 37°C, samples were phenol-chloroform extracted and ethanol precipitated with 4 µg of glycogen (Roche). The final DNA pellet was dissolved in 50 µl TE (Tris-EDTA; pH 8.0) for Q-PCR analysis. Five independent immunoprecipitations for each of the anti-H3K4me2, anti-H3K4me3 (four experiments with antibodies from Upstate and one experiment with antibodies from Abcam for these two antibodies), and the anti-RNA polymerase II were carried out for each patch-methylated episome.
Q-PCR and statistical analysis. Q-PCR was performed with a Bio-Rad iCycler using iQ Supermix (Bio-Rad) and 1 µl of each DNA sample. Fluorescence-labeled TaqMan probes for six regions of the plasmid pCLH22 were synthesized by Biosearch Technologies, and the primers were synthesized by Operon Technologies. Four of the primer and probe sequences were described previously (11). Primers 5'-ACAAGGAGAGAAAAAGCACCGT and 5'-GCCTTCCTAATAATTCACGATCGT and probe 5'-CATGCCGATTGGTGGAAGTAAGGTGG were used for the LTR3 region. Primers 5'-GTCACGCTCGTCGTTTGGT and 5'-TCATGTAACTCGCCTTGATCGT and probe 5'-TGGCTTCATTCAGCTCCGGTTCCC were used for the pBR backbone. Primers 5'-TGAATTCCCCAATGTCAAGCA and 5'-TCGGCACTTTGCATCGG and probe 5'-TTCCGGAATCGGGAGCGCG were used for the hygromycin 5 (Hyg5) region that is 121 bp downstream from the herpes simplex virus thymidine kinase promoter. All Q-PCR assays were carried out using the same two-step program: 95°C for 15 s and 60°C for 1 min for 40 cycles. Within each 96-well plate of Q-PCRs, titrations of known amounts of pCLH22 DNA were included as positive controls and for quantitation. DNA from the TCF and ChIP sample immunoprecipitated with no Ab was included in each 96-well plate of Q-PCRs with ChIP samples immunoprecipitated with antibodies. All Q-PCRs were performed in duplicate. The fraction of immunoprecipitated DNA (%IP) was calculated by dividing the amount of DNA from the sample immunoprecipitated with antibodies (after the background amplification from the no-Ab control was subtracted) by the amount of DNA in the corresponding TCF sample. The %IP of regions LTR1, LTR3, Luc1, Luc2, Hyg, and Hyg5 was then normalized by the %IP of the pBR region of each episome to derive the distribution of H3K4me2 and H3K4me3 on the episome. The %IP of regions LTR1, LTR3, Luc1, and Luc2 from each episome was normalized by the %IP of corresponding regions from pCLH22 to derive the distribution of RNA Pol II in the transcription unit on the episome. A t test was carried out to test the hypothesis that there is no significant difference between the %IP of a region from a patch-methylated episome and the %IP of the same region from pCLH22.
| RESULTS |
|---|
|
|
|---|
|
Compared with the pBR region of the episome, there was very little H3K4me2 associated with all regions studied on pMeCLH22, which is the episome with a high density of DNA methylation at CpG sites (Fig. 1B and data not shown). It is important to note that the percentage of DNA immunoprecipitated from the pBR region of pMeCLH22 by the anti-H3K4me2 antibody was more than fivefold less than that of pCLH22. Therefore, the level of H3K4me2 was more than 100-fold lower on pMeCLH22 than on pCLH22 in the RSVLTR region and the luciferase coding region. The luciferase gene activity of the fully methylated episome was more than 500-fold lower than that of the unmethylated episome, based on the luciferase assay (9-11). It is clear that H3K4me2 associates with actively transcribed DNA. This result is consistent with the findings by Schneider et al. (28) in which the level of H3K4me2 is increased in actively transcribed genes and that little or no H3K4me2 is found in most of the tissue-specific genes when they are in a silent state, in higher eukaryotes. This finding is also consistent with the observation that methylated H3K4 is associated with transcriptional activation in the chicken ß-globin locus (18). The finding that the level of H3K4me2 is substantially reduced on methylated DNA compared with that of unmethylated DNA is consistent with the finding of Lorincz et al. (19).
H3K4me3 distributes differently than H3K4me2 and is depleted from the transcribed region of the fully methylated episome. To explore the presence and distribution of H3K4me3 on unmethylated and fully methylated stable episomes pCLH22 and pMeCLH22 (in human cells), anti-H3K4me3 was used for ChIP assays, followed by analysis using the seven probe and primer sets described above. The controls, calculations, and normalizations were carried out as described above. It is not possible to compare the levels of H3K4me2 and H3K4me3 on the episome because the efficiencies of the two different antibodies in immunoprecipitation may be different. However, the distributions of H3K4me2 and H3K4me3 across the episome can be compared.
While the distribution of H3K4me3 showed a trend similar to that of H3K4me2 for most of the regions, Luc1 showed the lowest level of H3K4me3 among the seven regions studied (Fig. 1C and data not shown). The level of H3K4me3 was also considerably lower on pMeCLH22 than on pCLH22, except for the pBR region (Fig. 1C). It is uncertain why H3K4me3 was present at similar levels on pCLH22 and on pMeCLH22 in the pBR region, while none or very little of it was found in other regions of pMeCLH22. Otherwise, the distributions of H3K4me3 were mostly consistent with the findings described by Schneider et al. (28) and Lorincz et al. (19), as discussed above regarding H3K4me2. Therefore, there is a good correspondence between the stable replicating episome and the chromosome.
DNA methylation impacts the presence of H3K4me2 and H3K4me3 on the patch-methylated episomes. From the experiments described above, it is clear that H3K4me2 and H3K4me3 levels are much reduced on methylated episomes relative to those on unmethylated ones. We were interested in whether the respective distributions of H3K4me2 and H3K3me3 in the transcription unit are affected more by transcription or by DNA methylation. To test this, we carried out ChIP assays using the anti-H3K4me2 and anti-H3K4me3 antibodies on human cells harboring patch-methylated episomes pMeLTR, pMeLTR/Luc, and pMeLuc. The luciferase activity was approximately twofold lower for pMeLTR, which has CpG methylation only in the RSVLTR promoter, than for fully unmethylated pCLH22 (Fig. 2 and reference 11). The luciferase activity of pMeLTR/Luc, which has CpG methylation in the RSVLTR promoter and the luciferase coding region, was 50-fold lower than that of pCLH22 (Fig. 2 and reference 11). An eightfold reduction of luciferase activity was detected for pMeLuc, which has CpG methylation in the luciferase coding region only, compared with that of pCLH22.
|
On pMeLTR/Luc, the H3K4me2 level was much reduced in the methylated LTR1 region (more than 16-fold lower than in the pBR region), the LTR3 region (more than 6-fold lower than in the pBR region), the Luc1 region (20-fold lower than in the pBR region), and the Luc2 region (25-fold lower than in the pBR region). Similarly, the H3K4me2 level was much reduced in the methylated Luc1 and Luc2 regions, while the unmethylated LTR1 and LTR3 regions had higher levels of H3K4me2 on pMeLuc (Fig. 2). The levels of H3K4me2 remained high in the Hyg and Hyg5 regions that were unmethylated on all three patch-methylated episomes (data not shown). Patterns of H3K4me3 distribution on these three patch-methylated episomes were similar to those of H3K4me2 (Fig. 2). This finding is consistent with the observation of Lorincz et al. (19) that the methylated region is more than 500 bp downstream of the promoter.
The methylated luciferase gene region did not reduce the level of H3K4me2 or H3K4me3 in the unmethylated promoter of pMeLuc, even though the first methylated CpG was only 170 bp downstream from the LTR3 region. Likewise, promoter methylation on pMeLTR did not appear to reduce the H3K4Me2 or the H3K4Me3 level in the unmethylated coding region, despite the fact that the nearest CpG methylation site was only 115 bp upstream from the Luc1 region and that the luciferase gene activity was reduced twofold. These findings indicate that the chromatin structure associated with DNA methylation does not extend far upstream or downstream and that reduction of gene activity is not clearly correlated with the reduction of H3K4me2 and H3K4me3 levels. Importantly, the data here suggest a strong correlation between DNA methylation and the depletion of H3K4me2 and H3K4me3, even though a direct effect cannot be certain.
DNA methylation immediately downstream from the transcriptional start site has a major impact on gene activity. The experiments above did not allow us to clearly conclude whether the depletion of H3K4me2 and H3K4me3 is correlated with reporter gene activity, even though the findings suggest that DNA methylation in the coding region plays a major role in determining the presence of H3K4me2 and H3K4me3 locally. We generated two additional patch-methylated episomes, pMe5'Luc and pMe3'Luc, with different sizes and locations for the methylation patch in the luciferase coding region. The luciferase activity of pMe5'Luc, which has CpG methylation in the first 680 bp of the luciferase coding region, was reduced approximately fivefold from that of the fully unmethylated pCLH22 (Fig. 3). In contrast, the luciferase activity of pMe3'Luc, which has CpG methylation in the last 1,220 bp of the luciferase coding region, was very similar to that of the unmethylated pCLH22 (Fig. 3). This result is sufficient to conclude that the length of the patch of DNA methylation in the transcription unit is not the key determinant of transcription and that DNA methylation near the transcription start site has a major impact on gene activity. Compared with the luciferase activity of pMeLTR (Fig. 2), which had a methylated DNA patch of a size similar to that of pMe5'Luc, it is evident that the location of DNA methylation is critical in transcription suppression and that it is possible that DNA methylation impacts transcription elongation more than initiation.
|
Levels of H3K4me3 are reduced in both the methylated Luc1 region and the unmethylated Luc2 region of pMe5'Luc as well as in the methylated Luc2 region of pMe3'Luc (Fig. 3). The much reduced level of H3K4me3 found in the unmethylated Luc2 region of pMe5'Luc was surprising, but it was reproduced in all five experiments. This is one of the two instances where the presence of H3K4me2 and H3K4me3 is not concordant (Fig. 1, see the pBR region of pMeCLH22). These findings suggest that the gene activity may have some impact on the presence of H3K4me3 in addition to DNA methylation. The lack of H3K4me2 and H3K4me3 in the 3' portion of the coding region on pMe3'Luc does not affect gene expression significantly. Hence, H3K4me2 and H3K4me3 are not absolutely required for the passage of RNA Pol II, consistent with the implications of the study by Schneider et al. (28).
The presence of RNA polymerase II in the transcription unit is not correlated with the level of H3K4me2 or H3K4me3. To further define the relationships among DNA methylation, H3K4 modifications, and transcription, we carried out ChIP assays using antibodies against Pol II to determine the presence of Pol II on the patch-methylated episomes in human cells. Five independent ChIP assay experiments were carried out with cells harboring each of the six patch-methylated and fully unmethylated episomes. Q-PCR was carried out in duplicate with standards and controls for each run as described above. There was no appreciable signal above background level from the LTR1 and the pBR regions in all of the 35 ChIP assays (seven episomes, five experiments). This is consistent with the lack of transcription in the pBR region and the distances of the LTR1 region upstream from the TATA box (310 bp) and the transcription start site of the RSVLTR promoter (339 bp).
Unlike the assays for H3K4, no internal control region can be used to normalize the abundance of Pol II on each episome because no Pol II is in the regions without transcription. The ChIP assay results for Pol II on each patch-methylated episome were normalized relative to the results on fully unmethylated pCLH22, which had unaffected luciferase expression. Compared with unmethylated pCLH22, a more than fivefold reduction of Pol II was observed in the LTR3 region, and a two- to threefold reduction of Pol II was observed in the Luc1 and Luc2 regions of pMeLTR (Fig. 4). As described above, the levels of H3K4me2 and H3K4me3 remained high in the Luc1 and Luc2 regions of this stable episome. The reduction of Pol II in the unmethylated luciferase coding region is most likely due to the inhibition of transcription initiation by the DNA methylation and histone modification in the promoter.
|
There was a 25- to 100-fold reduction of Pol II in the Luc1 and Luc2 regions of pMe5'Luc, while only a small reduction of Pol II was observed in the LTR3 region (Fig. 4). This finding indicates that CpG methylation in the 600 bp of DNA immediately downstream from the promoter does not have a significant impact on the presence of Pol II in the promoter. Despite the nearly normal level of Pol II in the promoter of pMe5'Luc, very little Pol II was present in the coding region only 115 bp downstream from the transcriptional start site. This finding is consistent with the previous finding that DNA methylation in the coding region can reduce the efficiency of Pol II elongation (19). It is also important to note that the level of H3K4me2 was not reduced in the Luc2 region, even though a significant reduction of Pol II was observed in the same region of pMe5'Luc. This strongly indicates that transcriptional activity is not concordant with the presence of H3K4me2, and it is consistent with the observation for pMeLuc described above.
A nearly twofold reduction of Pol II was observed in the unmethylated LTR3 and Luc1 regions of pMe3'Luc, while a more than 10-fold reduction of Pol II was observed in the methylated Luc2 region of the same episome (Fig. 4). This finding adds further support to the interpretation here that DNA methylation can reduce Pol II elongation efficiency. As expected, there was nearly no Pol II found in all three regions of fully methylated pMeCLH22 (Fig. 4).
Taken together, the findings here strongly indicate that the presence of H3K4me2 correlates with DNA methylation and not the presence of Pol II. While the presence of H3K4me3 is not entirely concordant with DNA methylation, it also does not correlate with the presence of Pol II. These findings also indicate that DNA methylation immediately downstream from the promoter can reduce the Pol II presence at the promoter through mechanisms other than H3K4 modifications. This effect of DNA methylation is potentially size dependent. It is noteworthy that the luciferase gene activity correlates with the combined levels of Pol II in the Luc1 and Luc2 regions better than the combined levels of Pol II in the promoter and the coding region (LTR3, Luc1, and Luc2) (Fig. 4).
| DISCUSSION |
|---|
|
|
|---|
Methylation of the promoter (pMeLTR) and the 600 bp of the 5' region of the coding region (pMe5'Luc) have very different impacts on reporter gene activity, even though the DNA methylation patch sizes were similar. The level of Pol II was reduced about sevenfold in the promoter of pMeLTR, and that most likely leads to the two- to threefold reduction of Pol II in the coding region. In contrast, Pol II was nearly absent in the coding region of pMe5'Luc despite the relatively high level of Pol II detected in the promoter. These findings strongly suggest that DNA methylation impacts transcriptional initiation less than elongation.
DNA methylation at the 5' end of the coding region inhibits reporter gene expression much more than DNA methylation at the 3' portion of the coding region, even though the methylation patch at the 3' end was twice the size (Fig. 4, compare pMe5'Luc and pMe3'Luc). Also, DNA methylation at the 5' end of the coding region has an effect on reporter gene expression that is similar to DNA methylation in the entire coding region (Fig. 4, compare pMeLuc and pMe5'Luc). These findings further indicate that the location of the DNA methylation plays a more important role in transcriptional suppression than the size of the DNA methylation patch. It is clear that DNA methylation at the junction of transcriptional initiation and elongation may play a crucial role in gene suppression. DNA methylation in the region surrounding the transcriptional start site has been correlated with the lack of expression of genes in specific tissues or those on the inactivated X chromosome; these include the glutathione S-transferase gene (21), the hypoxanthine phosphoribosyltransferase gene (24), the Maspin gene (8), and the DEAD-box protein 4 (32). It has been proposed previously that methylation at five CpG sites surrounding the MAGE-A1 transcription start site might be required for the suppression of the gene (6). Our findings here may explain the role of DNA methylation in the suppression of these genes.
While DNA methylation in the promoter may have a small impact on transcriptional initiation, DNA methylation at the beginning of the transcriptional start site may inhibit transcription more dramatically. When DNA methylation is present in the region surrounding the transcriptional start site, which includes parts of the promoter and the coding region, transcription is effectively suppressed. The low reporter gene activity from the episome harboring DNA methylation in both the promoter and in the coding region (pMeLTR/Luc) supports this view, even though the size of the methylation patch cannot be completely ruled out as a contributing factor. We know that episomes harboring 120 bp of methylated DNA patch immediately downstream from the transcriptional start site showed similar reporter activity and no reduction of H3K4me2, H3K4me3, and Pol II levels compared with those of the fully unmethylated episome (data not shown). Therefore, a minimum size of DNA methylation (between 120 bp and 600 bp) at a critical region, such as the region of transition for transcriptional initiation and elongation, may efficiently suppress gene activity.
In the current study, DNA methylation of the entire 1.9-kb reporter coding region (pMeLuc) led to a more than eightfold reduction of Pol II, while the levels of H3K4me2 and H3K4me3 remained high in the promoter. However, 600 bp of DNA methylation at the 5' end of the coding region (pMe5'Luc) did not change the levels of Pol II, H3K4me2, and H3K4me3 in the promoter. Very little Pol II was found in the coding region of pMeLuc and pMe5'Luc, indicating the lack of transcription, particularly elongation. These findings indicate that the size of the DNA methylation patch may affect transcription initiation in the promoter immediately upstream through a mechanism independent of the presence of H3K4me2 and H3K4me3. This conclusion is not inconsistent with the finding of Lorincz et al. (19) that the methylated coding region, at about 1 kb downstream from the promoter, has little impact on transcriptional initiation. It is also noteworthy that an accumulation of Pol II in the promoter was not observed on either pMeLuc or pMe5'Luc, supporting the model in which the reduction of transcriptional elongation efficiency leads to premature Pol II dissociation, a conclusion also favored by Lorincz et al. (19). A much larger impact on reporter gene expression is observed when DNA methylation is present in the first 600 bp of the coding region (pMe5'Luc) than when the 1.2 kb of DNA at the 3' end of the coding region is methylated (pMe3'Luc). The level of Pol II is much lower in the methylated 3' region than the unmethylated 5' region on pMe3'Luc, indicating that DNA methylation affects transcription elongation. The transcription efficiency in the 3' coding region may play a less important role in overall gene activity, based on the observation that the reporter gene activity of pMe3'Luc is nearly normal (despite the absence of Pol II in the methylated 3' coding region).
It has been reported that methylated H3K4 is found in the actively transcribed gene region in higher eukaryotes (3, 15, 20, 28) and that the level of H3K4me2 is coupled with transcription (30). It has also been reported that the level of H3K4me2 is two- to threefold lower in the methylated transgene than in the unmethylated transgene (19). However, the connection between DNA methylation, H3K4 methylation, and transcription cannot be discerned from these previous studies. We observed a tight correlation between the depletion of H3K4me2 in the regions of DNA methylation, regardless of the reporter gene activity. The luciferase gene activity from pMe3'Luc is 20% reduced, but the level of H3K4me2 in the methylated 3' portion of the coding region is fivefold or more lower than those of episomes (pMeLTR, pMe5'Luc, and pCLH22) without DNA methylation in the same region. Conversely, the level of H3K4me2 remains high in the unmethylated DNA regions regardless of the presence of Pol II. For example, the reporter gene activity is more than fivefold suppressed and the Pol II level is more than 20-fold reduced, while the level of H3K4me2 remains high in the unmethylated 3' coding region of pMe5'Luc. Furthermore, DNA methylation approximately 200 bp downstream does not affect the presence of H3K4me2 in the unmethylated promoter on pMeLuc and pMe5'Luc. This distinct boundary effect strongly supports the view that DNA methylation affects the presence of H3K4me2 directly.
We propose here that DNA methylation dictates a closed chromatin structure that is devoid of H3K4me2 and inhibits transcription and that the presence of H3K4me2 marks an open chromatin structure that would permit transcription if all other conditions for active transcription were fulfilled (Fig. 5). Based on this model, H3K4me2 can be detected where active transcription is observed, and it can also be detected at other regions with open chromatin when no active transcription occurs in these regions. This model is consistent with the observations that H3K4me2 correlates with active chromosomal regions (3) and that H3K4me2 marks monoallelically expressed genes (25). This model would provide an explanation for the detection of H3K4me2 in the inactive ß-globin locus in the study by Schneider et al. (28), and it explains the high level of H3K4me2 observed in the unmethylated pBR region on the episomes that had no detectable transcriptional activity in the current study. This model is also consistent with the finding that methylated H3K4 is absent from the methylated DNA in the imprinted gene, even though DNA methylation is not concordant with gene expression (36). This model may provide further mechanistic understanding of the findings in a recent study describing the correlation between histone methylation and transgene expression (37).
|
Using antibodies from two independent sources, we consistently detected H3K9me3 on fully methylated pMeCLH22 in all regions examined (data not shown). However, the relationship among DNA methylation, transcription, and the presence of methylated H3K9 on the patch-methylated episomes appears to be much more complex. Methylated H3K9 has been reported to be associated with condensed heterochromatin, but H3K9 methylation has also been found to be enriched in actively transcribed regions of some genes (7, 33). It is possible that transcriptional elongation plays a significant role in the H3K9 modification. Further analysis of the interplay between DNA methylation and histone modifications is necessary.
| ACKNOWLEDGMENTS |
|---|
This work was supported by NIH grant RO1 GM54781.
| FOOTNOTES |
|---|
Published ahead of print on 22 January 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
2. Bernstein, B. E., E. L. Humphrey, R. L. Erlich, R. Schneider, P. Bouman, J. S. Liu, T. Kouzarides, and S. L. Schreiber. 2002. Methylation of histone H3 Lys 4 in coding regions of active genes. Proc. Natl. Acad. Sci. USA 99:8695-8700.
3. Bernstein, B. E., M. Kamal, K. Lindblad-Toh, S. Bekiranov, D. K. Bailey, D. J. Huebert, S. McMahon, E. K. Karlsson, E. J. Kulbokas III, T. R. Gingeras, S. L. Schreiber, and E. S. Lander. 2005. Genomic maps and comparative analysis of histone modifications in human and mouse. Cell 120:169-181.[CrossRef][Medline]
4. Braunstein, M., A. B. Rose, S. G. Holmes, C. D. Allis, and J. R. Broach. 1993. Transcriptional silencing in yeast is associated with reduced nucleosome acetylation. Genes Dev. 7:592-604.
5. Cedar, H. 1988. DNA methylation and gene activity. Cell 53:3-4.[CrossRef][Medline]
6. De Smet, C., C. Lurquin, B. Lethé, V. Martelange, and T. Boon. 1999. DNA methylation is the primary silencing mechanism for a set of germ line- and tumor-specific genes with a CpG-rich promoter. Mol. Cell. Biol. 19:7327-7335.
7. Eissenberg, J. C., and A. Shilatifard. 2006. Leaving a mark: the many footprints of the elongating RNA polymerase II. Curr. Opin. Genet. Dev. 16:184-190.[CrossRef][Medline]
8. Futscher, B. W., M. M. Oshiro, R. J. Wozniak, N. Holtan, C. L. Hanigan, H. Duan, and F. E. Domann. 2002. Role for DNA methylation in the control of cell type specific maspin expression. Nat. Genet. 31:175-179.[CrossRef][Medline]
9. Hsieh, C.-L. 1994. Dependence of transcriptional repression on CpG methylation density. Mol. Cell. Biol. 14:5487-5494.
10. Hsieh, C.-L. 1997. Stability of patch methylation and its impact in regions of transcriptional initiation and elongation. Mol. Cell. Biol. 17:5897-5904.[Abstract]
11. Irvine, R. A., I. G. Lin, and C.-L. Hsieh. 2002. DNA methylation has a local effect on transcription and histone acetylation. Mol. Cell. Biol. 22:6689-6696.
12. Jones, P. L., G. J. Veenstra, P. A. Wade, D. Vermaak, S. U. Kass, N. Landsberger, J. Strouboulis, and A. P. Wolffe. 1998. Methylated DNA and MeCP2 recruit histone deacetylase to repress transcription. Nat. Genet. 19:187-191.[CrossRef][Medline]
13. Kaludov, N. K., and A. P. Wolffe. 2000. MeCP2 driven transcriptional repression in vitro: selectivity for methylated DNA, action at a distance and contacts with the basal transcription machinery. Nucleic Acids Res. 28:1921-1928.
14. Kass, S. U., J. P. Goddard, and R. L. Adams. 1993. Inactive chromatin spreads from a focus of methylation. Mol. Cell. Biol. 13:7372-7379.
15. Kim, T. H., L. O. Barrera, M. Zheng, C. Qu, M. A. Singer, T. A. Richmond, Y. Wu, R. D. Green, and B. Ren. 2005. A high-resolution map of active promoters in the human genome. Nature 436:876-880.[CrossRef][Medline]
16. Klose, R. J., S. A. Sarraf, L. Schmiedeberg, S. M. McDermott, I. Stancheva, and A. P. Bird. 2005. DNA binding selectivity of MeCP2 due to a requirement for A/T sequences adjacent to methyl-CpG. Mol. Cell 19:667-678.[CrossRef][Medline]
17. Liang, G., J. C. Lin, V. Wei, C. Yoo, J. C. Cheng, C. T. Nguyen, D. J. Weisenberger, G. Egger, D. Takai, F. A. Gonzales, and P. A. Jones. 2004. Distinct localization of histone H3 acetylation and H3-K4 methylation to the transcription start sites in the human genome. Proc. Natl. Acad. Sci. USA 101:7357-7362.
18. Litt, M. D., M. Simpson, M. Gaszner, C. D. Allis, and G. Felsenfeld. 2001. Correlation between histone lysine methylation and developmental changes at the chicken beta-globin locus. Science 293:2453-2455.
19. Lorincz, M. C., D. R. Dickerson, M. Schmitt, and M. Groudine. 2004. Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells. Nat. Struct. Mol. Biol. 11:1068-1075.[CrossRef][Medline]
20. Miao, F., and R. Natarajan. 2005. Mapping global histone methylation patterns in the coding regions of human genes. Mol. Cell. Biol. 25:4650-4661.
21. Millar, D. S., C. L. Paul, P. L. Molloy, and S. J. Clark. 2000. A distinct sequence (ATAAA)n separates methylated and unmethylated domains at the 5'-end of the GSTP1 CpG island. J. Biol. Chem. 275:24893-24899.
22. Nan, X., H. H. Ng, C. A. Johnson, C. D. Laherty, B. M. Turner, R. N. Eisenman, and A. Bird. 1998. Transcriptional repression by the methyl-CpG-binding protein MeCP2 involves a histone deacetylase complex. Nature 393:386-389.[CrossRef][Medline]
23. Ng, H. H., F. Robert, R. A. Young, and K. Struhl. 2003. Targeted recruitment of Set1 histone methylase by elongating Pol II provides a localized mark and memory of recent transcriptional activity. Mol. Cell 11:709-719.[CrossRef][Medline]
24. Park, J.-G., and V. M. Chapman. 1994. CpG island promoter region methylation patterns of the inactive-X-chromosome hypoxanthine phosphoribosyltransferase (Hprt) gene. Mol. Cell. Biol. 14:7975-7983.
25. Rougeulle, C., P. Navarro, and P. Avner. 2003. Promoter-restricted H3 Lys 4 di-methylation is an epigenetic mark for monoallelic expression. Hum. Mol. Genet. 12:3343-3348.
26. Rountree, M. R., and E. U. Selker. 1997. DNA methylation inhibits elongation but not initiation of transcription in Neurospora crassa. Genes Dev. 11:2383-2395.
27. Santos-Rosa, H., R. Schneider, A. J. Bannister, J. Sherriff, B. E. Bernstein, N. C. Emre, S. L. Schreiber, J. Mellor, and T. Kouzarides. 2002. Active genes are tri-methylated at K4 of histone H3. Nature 419:407-411.[CrossRef][Medline]
28. Schneider, R., A. J. Bannister, F. A. Myers, A. W. Thorne, C. Crane-Robinson, and T. Kouzarides. 2004. Histone H3 lysine 4 methylation patterns in higher eukaryotic genes. Nat. Cell Biol. 6:73-77.[CrossRef][Medline]
29. Schübeler, D., M. C. Lorincz, D. M. Cimbora, A. Telling, Y.-Q. Feng, E. E. Bouhassira, and M. Groudine. 2000. Genomic targeting of methylated DNA: influence of methylation on transcription, replication, chromatin structure, and histone acetylation. Mol. Cell. Biol. 20:9103-9112.
30. Schübeler, D., D. M. MacAlpine, D. Scalzo, C. Wirbelauer, C. Kooperberg, F. van Leeuwen, D. E. Gottschling, L. P. O'Neill, B. M. Turner, J. Delrow, S. P. Bell, and M. Groudine. 2004. The histone modification pattern of active genes revealed through genome-wide chromatin analysis of a higher eukaryote. Genes Dev. 18:1263-1271.
31. Skowronska-Krawczyk, D., M. Ballivet, B. D. Dynlacht, and J. M. Matter. 2004. Highly specific interactions between bHLH transcription factors and chromatin during retina development. Development 131:4447-4454.
32. Song, F., J. F. Smith, M. T. Kimura, A. D. Morrow, T. Matsuyama, H. Nagase, and W. A. Held. 2005. Association of tissue-specific differentially methylated regions (TDMs) with differential gene expression. Proc. Natl. Acad. Sci. USA 102:3336-3341.
33. Vakoc, C. R., S. A. Mandat, B. A. Olenchock, and G. A. Blobel. 2005. Histone H3 lysine 9 methylation and HP1gamma are associated with transcription elongation through mammalian chromatin. Mol. Cell 19:381-391.[CrossRef][Medline]
34. Vire, E., C. Brenner, R. Deplus, L. Blanchon, M. Fraga, C. Didelot, L. Morey, A. Van Eynde, D. Bernard, J. M. Vanderwinden, M. Bollen, M. Esteller, L. Di Croce, Y. de Launoit, and F. Fuks. 2006. The Polycomb group protein EZH2 directly controls DNA methylation. Nature 439:871-874.[CrossRef][Medline]
35. Wigler, M., R. Sweet, G. K. Sim, B. Wold, A. Pellicer, E. Lacy, T. Maniatis, S. Silverstein, and R. Axel. 1992. Transformation of mammalian cells with genes from prokaryotes and eukaryotes. Biotechnology 24:444-452.[Medline]
36. Yamasaki-Ishizaki, Y., T. Kayashima, C. K. Mapendano, H. Soejima, T. Ohta, H. Masuzaki, A. Kinoshita, T. Urano, K.-I. Yoshiura, N. Matsumoto, T. Ishimaru, T. Mukai, N. Niikawa, and T. Kishino. 2006. Role of DNA methylation and histone H3 lysine 27 methylation in tissue-specific imprinting of mouse Grb10. Mol. Cell. Biol. 27:732-742.
37. Yan, C., and D. D. Boyd. 2006. Histone H3 acetylation and H3 K4 methylation define distinct chromatin regions permissive for transgene expression. Mol. Cell. Biol. 26:6357-6371.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | J. Virol. | Eukaryot. Cell |
|---|
| Microbiol. Mol. Biol. Rev. | Clin. Vaccine Immunol. | All ASM Journals |
|---|