|
|
|
RESEARCH PAPER
1 McArdle Laboratory for Cancer Research, University of Wisconsin Medical School, Madison, Wisconsin 53706, USA; 2 Howard Hughes Medical Institute, Division of Nucleic Acids Enzymology, Department of Biochemistry, Robert Wood Johnson Medical School, Piscataway, New Jersey 08854, USA; 3 NimbleGen Systems Inc., Madison, Wisconsin 53711, USA
| Abstract |
|---|
|
|
|---|
[Keywords: Suz12; histone methylation; polycomb; Eed; RNA interference; chromatin immunoprecipitation]
Received March 3, 2004; revised version accepted April 27, 2004.
The PcG complexes have been proposed to control gene activity via transcriptional repression. Recently, in vitro assays have provided insight concerning the mechanism underlying PcG-mediated transcriptional repression. Characterization of the PRC2/3 complexes indicates that each complex contains an intrinsic histone lysine methyltransferase (HKMT) activity that is mediated by the SET [Su-(var)39;E(z);Trithorax] domain of Ezh2 (Kuzmichev et al. 2002
). In vitro, the Ezh2 protein within the PRC2/3 complexes can methylate Lys 9 (H3-K9) and Lys 27 (H3-K27) of histone H3 and Lys 26 of histone H1 (H1-K26), depending on whether the oligonucleosomes contain histone H1 (Kuzmichev et al. 2004
). The efficiency with which the different histones (H1 or H3) are methylated depends on the specific form of Eed present in the complex (Kuzmichev et al. 2004
). The PRC2 complex contains the longest form of Eed (Eed1) and methylates both H1-K26 and H3-K27. However, PRC2 preferentially methylates H1-K26 when nucleosome arrays contain histone H1. The PRC3 complex, which contains the two shortest forms of Eed (Eed3/4), methylates H3-K27. The intermediate form of Eed (Eed2) is also present in yet another distinct PRC complex (A. Kirmizis, R. Margueron, A. Kuzmichev, P. Farnham, and D. Reinberg, unpubl.). Based on these and previous discoveries, a model for the mechanism of PcG-mediated transcriptional silencing has been proposed. The PRC2 or PRC3 complex is thought to first catalyze the addition of methyl groups to H3-K27 that serve as signals for the recruitment of PRC1. Binding of PRC1 is then proposed to either block the recruitment of transcriptional activating factors, such as SWI/SNF (a Trithorax complex), facilitating the establishment of a stable, repressive chromatin structure or to prevent transcription initiation by prebound factors (Simon 2003
; Dellino et al. 2004
).
Although an attractive model, it has not yet been demonstrated that the recruitment of HKMT-containing PRC complexes results in histone modification in mammalian cells owing to the fact that no target genes for these complexes have been identified. In contrast, target genes for PcG complexes have been identified in Drosophila. In the fruit fly, PcG target genes are regulated by specific DNA regions called Polycomb response elements (PREs; Simon et al. 1993
). PREs are
23 kb long and are composed of multiple different DNA motifs (Bloyer et al. 2003
). Usually PREs contain binding sites for GAGA and PHO, two site-specific DNA-binding proteins that interact with PcG proteins (Poux et al. 2001
) and function to recruit Drosophila PcG complexes to DNA (for review, see Pirrotta et al. 2003
). Mammalian PREs have not yet been identified, and therefore, it remains unclear how mammalian PcG complexes are recruited to chromatin to regulate expression of specific target genes. Previous reports have demonstrated an interaction between PcG proteins and members of the general transcriptional machinery, suggesting that binding of PcG complexes to promoters may be stabilized by these interactions (Breiling et al. 2001
; Saurin et al. 2001
). However, it is likely that a HKMT-containing PRC complex must first be recruited by a sequence-specific DNA-binding protein.
Because none of the HKMT-containing PRC complexes contain a site-specific DNA-binding protein, target genes cannot be identified using a bioinformatics approach to search for consensus binding sites in the genome. Therefore, we have used two global screening methods to identify genes that are regulated by mammalian PRC2/3 complexes. First, using oligonucleotide cDNA arrays, we identified genes whose expression changes upon siRNA-mediated depletion of Suz12. Second, we used a combination of chromatin immunoprecipitation and genomic microarrays to identify promoters bound by Suz12. Interestingly, we show that multiple components of the PRC2/3 complexes colocalize on these newly identified target genes and that binding of these proteins to target genes correlates with methylation of H3-K27.
| Results |
|---|
|
|
|---|
Although Drosophila PcG proteins have been mostly implicated in transcriptional control during development, human PcG proteins have been linked to cancer development. For example, previous studies have demonstrated that Ezh2 is overexpressed in various human cancers (Varambally et al. 2002
; Bracken et al. 2003
; Kleer et al. 2003
). Additionally, in endometrial stromal tumors, a fragment of the Suz12 gene is frequently fused to a second gene encoding a zinc finger protein (Koontz et al. 2001
). In our previous studies, we have observed up-regulation of Suz12 mRNA in tumors of the colon, breast, and liver (Kirmizis et al. 2003
). Specifically, we found that eight of nine colon cancer patients tested had high levels of Suz12 mRNA in the tumor tissue as compared with the normal tissue. Furthermore, we showed that the up-regulation of Suz12 mRNA is mediated by the TCF/
-catenin transcription complex (Kirmizis et al. 2003
), whose increased activity is responsible for the development of the majority of human colon cancers (for review, see Bienz and Clevers 2000
). The fact that Suz12 mRNA is present at high levels in colon tumors suggests that colon cancer provides a good model system to identify PRC target genes. However, it was first necessary to determine if the increased Suz12 mRNA expression in colon tumors correlates with increased Suz12 protein levels. We developed an antibody against the Suz12 protein and demonstrated that it recognizes the correct protein, using cells transfected with a plasmid expressing an HA-tagged Suz12. Western analysis showed that the HA antibody recognized a 93-kDa protein in cells transfected with HA-Suz12 but not in cells transfected with the empty vector (Fig. 1A, left panel). The purified Suz12 antibody recognized the same protein in both transfected cell populations (Fig. 1A, right panel). This result is expected because the endogenous protein runs the same size as the HA-tagged protein. We next analyzed whole-tissue extracts from five different colon cancer patients for Suz12 protein levels. As shown in Figure 1B, four out of five samples showed increased levels of Suz12 in the tumor tissue as compared with adjacent normal tissue. One patient had high Suz12 levels in both the normal and tumor tissues, and, therefore, it is possible that the normal sample was contaminated with tumor tissue. These findings are in agreement with our previous results showing up-regulation of Suz12 mRNA in colon tumors. In three of our extract preparations, the antibody recognized the 93-kDa band plus a lower band of
85 kDa. This lower band has also been detected by others using a different Suz12 antibody (Tie et al. 2003
). It is possible that the lower band represents a Suz12 degradation product that is produced during extraction of the tumor tissue because we have not detected it in nuclear extracts prepared from cell lines (data not shown).
|
Recent studies have shown that Suz12 is a component of the PRC2/3 complexes (Cao et al. 2002
; Czermin et al. 2002
; Kuzmichev et al. 2002
, 2004
; Muller et al. 2002
), suggesting that up-regulation of Suz12 in colon cancers may modulate the tumor phenotype by modifying histones and thus changing gene expression profiles. To examine this hypothesis, we used RNA interference (RNAi) to decrease Suz12 expression in colon cancer cells. We designed small interfering RNAs (siRNAs) to target the Suz12 transcript and tested them using two human colon cancer cell lines (Fig. 2; data not shown). We observed a decrease in both Suz12 mRNA and protein levels by 24 h after transfecting cells with the Suz12 siRNAs (Fig. 2A,B, cf. lanes 1 and 4). The down-regulation of Suz12 was more prominent after 72 h of transfection (Fig. 2A,B, cf. lanes 3 and 6). Lamin A/C and green fluorescent protein (GFP) siRNA duplexes did not affect the amounts of Suz12 verifying the specificity of the RNAi approach. We next examined the phenotype of the colon cancer cells having reduced Suz12 expression. We did not detect any changes in the morphology or the cell cycle profile of the cells after treatment with Suz12 siRNAs for 72 h (data not shown). However, RNAi does not eliminate Suz12 expression entirely from the cells, and the remaining Suz12 may be sufficient to mask any obvious phenotypes, especially in a short-term assay. It is also possible that certain functions of Suz12 in these cells are redundant with another protein, and, therefore, removing only Suz12 is not enough to cause changes in cell proliferation.
|
|
The experiments described above show that down-regulation of Suz12 alters the expression level of multiple genes. However, these findings do not demonstrate if Suz12 controls the expression of these genes directly by binding to their promoters. It is possible that removal of Suz12 from the cells changes the expression of these genes through an indirect mechanism. Therefore, we used chromatin immunoprecipitation (ChIP) to investigate whether Suz12 is recruited to the promoters of the deregulated genes. Because Suz12 is a component of transcriptional repressive complexes, we focused our initial investigation on the genes whose expression increased upon depletion of Suz12 from the cells. We performed ChIP experiments in SW480 cells using antibodies to Suz12, RNA polymerase II (Pol II), and a nonspecific IgG control. Using PCR primers that spanned the transcriptional start site, we monitored binding of Suz12 and Pol II to six promoters that drive expression of mRNAs that were up-regulated upon loss of Suz12. Two of the genes, Myelin transcription factor 1 (MYT1) and N-ethylmaleimide-sensitive factor (NSF), are expressed from two different transcriptional start sites separated by several kilobases. Therefore, we monitored binding of Suz12 to both the proximal and distal promoter regions of these two genes. Suz12 showed robust binding to the distal promoter region of MYT1 but not to the proximal promoter region (Fig. 4, lane 3). We did not detect strong binding of Suz12 to the distal or the proximal promoter region of the NSF gene or at the other four tested promoters. Studying specific protein/DNA interactions through ChIP assays requires the use of a highly specific antibody. The specificity of the Suz12 antibody was first illustrated by the fact that it recognized a single protein band by immunoblot analysis in nuclear extract preparations of cultured cells (data not shown) and by the fact that the protein detected on the Western blot is diminished after introduction of Suz12 siRNAs to the cells (Fig. 2B). To further demonstrate the specificity of the Suz12 antibody, we performed a ChIP assay in which the antibody was preincubated with the Suz12 immunogen (a fragment of Suz12 protein used to develop the antibody) at a 10-fold excess by weight. The immunogen abrogated the ability of the Suz12 antibody to enrich the MYT1 distal promoter (Fig. 4, lane 1). Collectively, these results demonstrate that the Suz12 antibody used in the ChIP experiments specifically enriches selected chromatin fragments. Therefore, the ChIP analysis demonstrated that MYT1 expression is directly regulated by Suz12, whereas expression of the other five tested genes may be indirectly regulated by Suz12.
|
500 bp from the Suz12- and IgG-immunoprecipitated chromatin samples as well as from an input chromatin reference sample (see Fig. 7C, right panel [below]). The amplicons were labeled with a fluorescent dye and hybridized onto separate custom promoter arrays. To identify the genes that are directly regulated by Suz12, we compared the hybridization intensities between the Suz12 array and the input array and calculated a Suz12/Input ratio for each oligonucleotide probe. Regions with a Suz12/Input ratio higher than 2 were considered positive for Suz12 binding. To eliminate from consideration promoter regions that are enriched nonspecifically during the ChIP procedure, we also calculated an IgG/Input ratio for each oligonucleotide probe. We performed three independent analyses and found eight promoters that were repeatedly enriched by the Suz12 antibody suggesting that those genes are directly regulated by Suz12 (Table 1). Profiles from a representative experiment of four of the eight positive promoters are shown in Figure 5. Use of the high-density oligonucleotide arrays allowed for more precise localization of Suz12 recruitment to its target loci than did the standard ChIP experiments. The small size of the amplicons in combination with the high density of the custom oligonucleotide arrays mapped the binding site of Suz12 within 500-bp regions. The results from these ChIP-chip experiments were validated in independent ChIP experiments using PCR primers that spanned the regions showing the highest enrichment by the Suz12 antibody (Fig. 5, inserts). As expected, analysis of the MYT1 promoter showed specific recruitment of Suz12. Interestingly, Suz12 binding to MYT1, as well as to some of the other promoters, occurred at two distinct promoter regions; one near the transcriptional start site and one farther upstream (e.g., MYT1 and SYBL1). The two occupied regions on these promoters were separated by several hundred bases, suggesting that Suz12 may associate with the transcriptional start site via a DNA looping mechanism mediated by proteinprotein interactions that is initiated by specific binding of Suz12 to an upstream region. Such a looping mechanism has been proposed to mediate the interaction of Drosophila PcG proteins with their target promoters (Orlando et al. 1998
|
|
|
Although the experiments described above identified a small set of direct target genes of Suz12, they did not identify numerous high-affinity binding sites. As noted above, this could be caused by the fact that most of the genes identified in the siRNA experiments are indirectly regulated by Suz12. Therefore, to identify additional Suz12 targets, we used a ChIP-CpG-island microarray technique that has previously been used in our laboratory to identify genomic loci bound by E2F family members (Weinmann et al. 2002
; Oberley et al. 2003
; Wells et al. 2003
). This technique will identify direct Suz12 targets that are within several kilobases of the CpG islands present on the microarray. Using the Suz12 antibody and preimmune (PI) serum as a control, we performed a standard ChIP assay in SW480 cells. We next prepared amplicons from the Suz12- and PI-immunoprecipitated chromatin samples as well as from an input chromatin reference sample. The amplicons from the immunoprecipitated samples were labeled with one fluorescent dye, whereas the reference sample (input) was labeled with a different fluorescent dye. To identify Suz12-bound loci, we applied the Suz12 and input amplicons to an array containing
12,000 CpG islands. We selected all the CpG clones that had intensity signals from the Suz12 chromatin at least threefold higher than the input control (Fig. 6A). To control for CpG clones that were enriched nonspecifically by the ChIP procedure, we hybridized a second array with the preimmune versus input amplicons. CpG islands that were enriched by the preimmune serum were discarded from further analysis. This procedure was repeated twice, and 52 CpG clones were selected as positives. Next, we used the University of California-Santa Cruz Genome Browser (http://genome.ucsc.edu) to determine the genomic location of all the CpG clones. Out of 52 clones, 41 contained sequences that matched to locations near promoters or 3'-ends of known and predicted genes (Supplementary Table S2). Because there is some degree of redundancy of the CpG islands printed on the arrays, three gene promoters were identified twice as unique clones. Also, three clones contained sequences that matched to two different gene promoters. Probably during the construction of the CpG-island library, DNA fragments representing different CpG islands were introduced in tandem within one clone. Additionally, two CpG clones did not match to any genomic loci, three matched to sequences that were not near annotated genes, and six matched to repeat elements; these 11 clones are not listed in Supplementary Table S2.
|
PRC-mediated repression of Suz12 target genes is associated with methylation of H3-K27
Recent studies have shown that Suz12 is a component of two different PRC complexes, each of which contains Ezh2, an HKMT that is associated with transcriptional silencing, and at least one form of Eed (Kuzmichev et al. 2004
). Therefore, we examined binding of Suz12, Ezh2, and Eed along a 5-kb region around the distal MYT1 promoter (Fig. 7A). For these experiments, we generated amplicons from ChIP samples by LMPCR and used them in conventional PCR using promoter-specific primers. This type of analysis allows precise localization of the binding of a factor to a promoter region because, as shown in Figure 7C, amplicons have a smaller range in size (200500 bp) as compared with immunoprecipitated DNA fragments (13 kb) isolated from a standard ChIP assay. Using this analysis, we showed that Suz12 binds to two regions in the 5-kb sequence; one near the distal transcriptional start site and one several kilobases upstream (Fig. 7B). These results are in complete agreement with our findings from the experiments using the custom oligonucleotide arrays (cf. Figs. 7B and 5). Using this high-resolution PCR analysis, we also found that Ezh2 and Eed colocalize with Suz12 on the MYT1 promoter, confirming that MYT1 is a target gene of the PRC complexes (Fig. 7B). To determine if the recruitment of the PRC complexes correlates with the HKMT activity of Ezh2, we next examined the methylation of H3-K27 at the 5-kb region of the MYT1 promoter. We show that binding of Suz12, Ezh2, and Eed to the MYT1 promoter correlates with, but does not precisely overlap, trimethylation of H3-K27. This result suggests that the PRC complexes can methylate histones located at the PRC-binding sites and in the surrounding region. Interestingly, no methylation of H3-K9 was detected at the repressed MYT1 locus, indicating that in vivo PRC-mediated silencing is associated with trimethylation of Lys 27 but not Lys 9 of histone H3 (Fig. 7B), as previously suggested (Kuzmichev et al. 2002
). The results of the siRNA experiments (demonstrating increased MYT1 mRNA in the absence of Suz12) and ChIP experiments (demonstrating methylation of H3-K27) suggest that the MYT1 promoter is kept in a repressed state by its association with the PRC2/3 complexes. Interestingly, we could not detect binding of Pol II to any location within the 5-kb sequence surrounding the MYT1 distal promoter. This lack of Pol II suggests a model for the PRC-mediated regulation of the MYT1 promoter; first, recruitment of PRC2/3 causes methylation of H3-K27, which then results in the formation of a chromatin structure that is incompatible with binding of Pol II, and thus the formation of a transcription initiation complex is prevented.
The robust Suz12 and weak Pol II signals seen on four of the promoters identified on the CpG-island arrays (Fig. 6C; data not shown) suggest that they may be regulated by a PRC complex in a manner similar to that of the MYT1 promoter. To test this hypothesis, we examined binding of Ezh2 to the AK056349, Wnt1, KCNA1, and CNR1 promoters. As shown in Figure 7D, Ezh2 is recruited to all of these Suz12 target promoters. Interestingly, binding of Suz12 and Ezh2 correlates with methylation of H3-K27, but not of H3-K9, at these four promoters. To demonstrate the specificity of our antibodies, we analyzed the methylation status of two other promoters, one active (the largest subunit of Pol II) and one inactive (the pericentromeric gene HTF6), in SW480 cells. As expected, the Pol II promoter was bound by Pol II but not by Suz12 or Ezh2 and showed no methylation of H3-K27 or H3-K9. In contrast, the silenced pericentromeric HTF6 promoter was characterized by methylation of H3-K27 and H3-K9. Methylation of H3-K27 at the HTF6 promoter may be mediated by a methyltransferase activity other than that of the PRC complexes because Suz12 and Ezh2 were not bound to the promoter (Fig. 7D). Alternatively, the kinetics of binding of an Ezh2-containing complex to this promoter may be different from the other promoters studied.
Our RNAi experiments described above showed that removal of Suz12 from colon cancer cells can result in increased expression of some target genes. Because of the correlation between Suz12 binding and H3-K27 methylation, we hypothesized that the increased expression is caused by the loss of the histone methyltransferase activity after depletion of Suz12 from cells. To test this hypothesis, we examined the binding pattern of Ezh2 and the status of H3-K27 methylation on five target promoters after siRNA-mediated removal of Suz12 from SW480 cells. We incubated cells with Suz12 siRNAs for 6 d, and then we performed ChIP assays using antibodies to Suz12, Ezh2, trimethyl H3-K27, and IgG as a control. A 6-d treatment was performed to allow for a greater loss of the methylated histones from the target promoters. As expected, we found that depletion of Suz12 from cells decreases the binding of Suz12 to the tested promoters (Fig. 8, cf. lanes 1 and 7). Interestingly, depletion of Suz12 from cells also resulted in a dramatic decrease of Ezh2 binding and H3-K27 methylation at the Suz12 target promoters. (Fig. 8, cf. lanes 2 and 8, and lanes 3 and 9). These results demonstrate that repression of the identified Suz12 target promoters is associated with methylation of H3-K27. Also, these data show that the decreased expression of Suz12 in cultured cells disrupts the function of the PRC2/3 complexes by preventing their recruitment to target promoters.
|
| Discussion |
|---|
|
|
|---|
The Suz12 target gene MYT1 was originally cloned from a human brain cDNA library on the basis of its ability to bind cis-regulatory elements of the glia-specific myelin proteolipid protein (PLP) gene and is suggested to be the prototype of the C2HC-type zinc finger protein family (Kim and Hudson 1992
). More recently, the Xenopus ortholog of MYT1 (X-MYT1) was identified as a transcriptional activator because it could induce expression of an N-tubulin promoter reporter construct in transient transfection assays (Bellefroid et al. 1996
). In the same study, the authors demonstrated that dominant-negative forms of X-MYT1 inhibited normal neurogenesis, suggesting that X-MYT1 is essential for inducing neuronal differentiation. Intriguingly, a recent report showed that the Xenopus ortholog of Ezh2 (XEZ) is expressed exclusively in the anterior neural plate during early Xenopus embryogenesis and postulated that XEZ might be involved in delaying anterior neuronal differentiation (Barnett et al. 2001
). Based on our findings, it is possible that Ezh2 delays neuronal differentiation, via the PRC2/3 complexes, by repressing the activity of the MYT1 gene. Besides MYT1, we have identified four additional promoters as being robustly bound by components of the PRC2/3 complexes; each of these promoters is also characterized by high levels of H3-K27. Although a link between components of the PRC2/3 complexes and Wnt1, the cannabinoid receptor (CNR1), or the potassium channel KCNA1 have not been previously reported, it is intriguing to note that these mRNAs are expressed at very low levels in most human tissues (see http://expression.gnf.org), suggesting that they may be generally silenced by the PRC complexes. In support of this hypothesis, we have shown that some of these target genes are bound by PRC2/3 components in other cell lines, such as the human MCF7 and mouse F9 (data not shown).
Previous studies have used mRNA expression profiling to identify Ezh2 target genes. In one study, gene expression was monitored in cells that lacked Ezh2 (Bracken et al. 2003
), whereas in the other study (Varambally et al. 2002
) the analysis was performed in cells that overexpressed Ezh2. The genes identified in these previous studies do not overlap with the PRC2/3 target genes that we have now identified. It is possible that many, if not all, of the previously identified genes are indirectly regulated by Ezh2. It is also possible that the distinct sets of identified genes might be the result of using different cell types and/or different DNA microarrays. Therefore, our future studies will aim at the identification of a large set of mammalian PRC target genes, using various different normal and cancer cell types.
Although we have described the first known direct target genes of a human PRC complex, identifying a larger set of targets is essential to better understand the mechanisms that mediate the functions of Polycomb complexes. For example, it is unclear if all Suz12 and/or Ezh2 target genes will be regulated by the PRC2/3 complexes or if different Eed isoforms (and thus different PRC complexes) will predominate in different cell or tumor types. Therefore, it remains possible that methylation of histone H1, but not histone H3, will occur at certain PRC target genes. Also, it is possible that binding of components of the PRC2/3 complexes may sometimes lead to activation, rather than repression of gene expression. A PcG protein that has both activation and repression functions is now classified in the Enhancer of Trithorax and Polycomb (ETP) group (for review, see Brock and van Lohuizen 2001
). Previous studies have shown that loss of expression of Ezh2 in Drosophila results in both increased and decreased gene expression (LaJeunesse and Shearn 1996
). Similar findings were presented in another study in which siRNA-mediated depletion of Ezh2 or Eed in mammalian cells resulted in decreased expression of several genes (Bracken et al. 2003
). Unfortunately, the previous studies did not distinguish between direct and indirect regulation by mammalian Ezh2 or Eed. However, our findings provide support for the hypothesis that components of the PRC2/3 complex can directly function as activators. Using ChIP, we showed that Suz12 and Ezh2 are recruited to the promoters of some genes (i.e., SYBL1) that are down-regulated upon Suz12 depletion from cells (Fig. 6; data not shown), suggesting that the PRC complexes normally activate the expression of these genes. Future studies should determine the mechanism of transcriptional activation mediated by the PRC complexes.
Finally, it is still unclear how the PRC complexes are directed to their target loci. Because no site-specific DNA-binding factor has been purified in the PRC complexes, it is not yet possible to identify a PRC target promoter by sequence inspection. In fact, even identifying a gene as being regulated by a PRC complex does not provide information as to exactly where the complex is recruited because the PRC-specific element (i.e., the PRE) might be located at a distance from the start of transcription. For example, we found that Suz12 was bound at -4000 in some promoters and at +1 in others. Interestingly, we also showed that many Suz12 target genes had multiple binding sites. In summary, this current study, which has identified the first set of human PRC2/3 target genes and demonstrated the utility of custom oligonucleotide arrays in presenting a detailed profiling of the binding of PRC components to large regions of target promoters, provides a significant step toward the goal of defining a mammalian PRE. Once we have identified a large number of direct target genes, we will use the information derived from the analysis of the custom arrays to search for common sequence elements in small, defined regions of the PRC target promoters.
| Materials and methods |
|---|
|
|
|---|
Using PCR, a 366-bp fragment from the N terminus and a 318-bp fragment from the C terminus of Suz12 ORF were generated and cloned into the Pet15b expression vector (catalog no. 69661, Novagen). The sequence of the primers can be found in Supplementary Table S3. The KIAA0160 (alternative name for Suz12) human cDNA clone, a gift from Takahiro Nagase (Kazusa DNA Research Institute, Osaka, Japan) was used as a template. Protein expression and purification were performed using the Bugbuster His-bind Purification Kit (catalog no. 70793, Novagen) according to the manufacturer's instructions. To prepare the Suz12 immunogen for rabbit injections, the N- and C-terminal protein fragments were combined in equal amounts to a final concentration of 1 µg/µL. Rabbit immunizations and antibody purifications were performed at Covance Inc. The Suz12 antibody is now available at Abcam, Inc. (catalog no. ab12201).
Plasmid constructs and transient transfections
To prepare a construct expressing HA-tagged Suz12, primers containing EcoRI sites on their 5'-ends were used in PCR to generate full-length cDNA of Suz12 (2219 bp). The KIAA0160 human cDNA clone was used as a template. The sequence of the primers can be found in Supplementary Table S3. The PCR product was digested with EcoRI and ligated into the EcoRI site of the pCMV-HA vector (Clontech). Bacterial colonies were screened for inserts and sequenced at the University of Wisconsin Biotechnology Sequencing Facility to ensure proper sequence and orientation. For transient transfections, SW480 human colon cancer cells were cultured in McCoy's 5A medium (Invitrogen) supplemented with 10% fetal bovine serum (Nova-Tech) and 1% Penicillin/Streptomycin (Invitrogen), and incubated at 37°C in a humidified 5% CO2 incubator. One day before the transfection, SW480 cells were plated at 2 x 105 cells per 60-mm dish. The cells were transfected with 2 µg of the HA-Suz12-expressing construct (pCMV-HA-Suz12) or 2 µg of the control construct (pCMV-HA-empty) using the FuGENE 6 transfection reagent (Roche) according to the manufacturer's protocol. Cells were harvested 48 h after transfection and were subjected to immunoblot analysis.
Immunoblot analysis
Normal and diseased human colon tissues obtained from the Cooperative Human Tissue Network (CHTN) were minced into small pieces using a razor blade. About 300 mg of minced tissue was resuspended in 3 mL of 1x PBS buffer and disaggregated using a Medimachine (catalog no. 340587, Becton Dickinson) to achieve single cells. The cells were then diluted in a 2x volume of whole-cell lysis buffer (0.1 M HEPES, 0.5 M KCl, 5 mM MgCl2, 0.5 mM EDTA, 35% glycerol, 5 mM NaF, 1 mM DTT, 1 mM PMSF) and were alternately frozen and thawed three times in liquid nitrogen to rupture the cell membranes. The samples were incubated for 30 min on ice to lyse the nuclei and then centrifuged at 4000 x g for 5 min to pellet the cell membranes. The protein content of the resulting supernatant was quantitated using a standard Bradford assay. Equal amounts of protein (2040 µg) from the tumor and normal tissue samples of each patient were subjected to standard Western analysis. In the case of the RNAi experiments (see below) and the transient transfections (see above), 5 x 105 cells were harvested from each treatment and examined by a standard Western analysis. Antibodies used for probing included anti-Suz12 (described above), anti-Lamin A/C (sc-7292, Santa Cruz Biotechnology), anti-HA (catalog no. MMS-101R, Covance Research Products), and antiactin (sc-1615, Santa Cruz Biotechnology).
RNA interference
Two different 21-nt duplex siRNAs for Suz12 (target sequences of 5'-CCCGGAAATTTCCCGTCCC-3' and 5'-GAGATGACCT GCATTGCCC-3', custom orders), Suz12 smart-pool (M0069570050) and control siRNAs for Lamin A/C (target sequence 5'-CTG GACTTCCAGAAGAACA-3', D0010000105) or GFP (target sequence 5'-GGCTACGTCCAGGAGCGCACC-3', D0013000120) were synthesized by Dharmacon. For the immunoblot analysis (see above) and RNA extraction (see below), SW480 and HT29 human colon cancer cells were plated at 2 x 105 cells per well in a 12-well plate. Twenty-four hours after plating, the cells were transfected either with both Suz12 siRNA duplexes together (100 nmoles each) or with each control siRNA (200 nmoles each) using oligofectamine (Invitrogen) or trans-IT TKO (Mirus) according to the manufacturers' instructions. At various time points after transfection, the cells were treated with trypsin and harvested for immunoblot analysis, or RNA preparation. For the RNAi-ChIP experiments, SW480 cells were plated at a density of 1.2 x 106 cells per 100-mm dish, and after 24 h they were transfected with 100 nm of Suz12 smart-pool using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. The cells were incubated with the Suz12 smart-pool for 72 h and then were harvested and replated for another 72 h of transfection with 100 nm of Suz12 smart-pool. At the end of the second 72-h incubation, the cells were harvested and used in ChIP assays as described below.
RTPCR analysis
Total RNA was prepared from 1 x 106 SW480 cells of each indicated treatment using the RNeasy kit (QIAGEN) as described by the manufacturer. For each RTPCR reaction, 100 ng of RNA was used, and the amplifications were performed as described previously (Graveel et al. 2001
). A list of all the primer sequences used in RTPCR is available as Supplementary Table S3. All the primers used in this study were obtained from the University of Wisconsin Biotechnology Center.
RNA microarray analysis
Total RNA was extracted from cells transfected with either Suz12 siRNAs or control GFP siRNA as described above. For each treatment, RNA was prepared from three independent experiments and was combined into one sample to reduce the experimental variation. Targets for microarray hybridization were generated from the RNA according to the supplier's instructions (Affymetrix). The human U133A gene chip (Affymetrix) was used for gene expression profiling, which represents
33,000 transcripts. Four U133A arrays were used for this analysis; two arrays were hybridized with cRNA obtained from cells treated with siRNAs to Suz12, and two arrays were hybridized with cRNA obtained from cells treated with siRNA to GFP. Hybridization, washing, scanning, and analysis of gene chips were performed at the University of Wisconsin Gene Expression Center according to the manufacturer's instructions. Expression levels were analyzed by the statistical algorithm in the Microarray Analysis Suite (MAS) 5.0 software (Affymetrix) using the default parameters. The data from the GFP siRNA treatment were used as a baseline expression for comparison with the Suz12 siRNA-treated sample.
ChIPs
Formaldehyde cross-linking and ChIP in SW480 cells were performed as described before (Weinmann et al. 2001
) with the following exceptions. Precleared chromatin from 2 x 106 cells was used for each ChIP sample. After reversing the cross-links, each individual ChIP sample was purified using the Qiaquick PCR purification kit (QIAGEN), and samples were eluted with 30 µL of elution buffer. The antibodies used in the ChIP assays included purified Suz12 antibody (see above), Ezh2, Eed-M26, trimethyl H3-K27, trimethyl H3-K9 (ab8898, Abcam Inc.), RNA Polymerase II (sc-899, Santa Cruz Biotechnology), histone H3 (ab1791, Abcam Inc.), rabbit anti-goat (catalog no. 55335, ICN Pharmaceuticals; used as an IgG control), and rabbit preimmune serum (Covance; used as a control in CpG-island array experiments). All antibodies were used at 1 µg per IP except the Suz12 antibody, which was used at 3 µg per IP. For the immunogen IP control (experiments shown in Fig. 4), Suz12 antibody was mixed with a 10-fold excess by weight of immunogen (bacterially expressed Suz12 protein fragments) and incubated with rotation at 4°C overnight prior to use in the ChIP assay. A detailed protocol of the ChIP procedure and the sequences of all the primers used in PCR can be found in Supplementary Table S3.
LMPCR
The generation of amplicons from individual ChIP samples was adapted from Ren et al. (2000
). A detailed protocol can be found in Oberley et al. (2004
) or at http://mcardle.oncology.wisc.edu/farnham.
Oligonucleotide promoter array analysis
Oligonucleotide promoter array analysis was performed by NimbleGen Systems as part of a Chromatin Immunoprecipitation (ChIP) Array Service. The custom arrays, each having 362,890 60-nt oligonucleotides, contained probes that tiled through 5 kb of the promoter regions of the selected genes. The probes represented both the forward and reverse strands, were spaced every 15 bases, and were printed at five (reverse strand) or six (forward strand) random locations on the array. For labeling, hybridization, and data analysis, refer to the Supplemental Material.
CpG-island microarray analysis
For detailed protocols on these steps, see Oberley et al. (2004
) or http://mcardle.oncology.wisc.edu/farnham. The Cy5 dye was coupled to the Suz12 (or the preimmune) amplicon, and Cy3 dye was coupled to the input reference amplicon, using standard methods. Labeled chromatin was hybridized with human CpG-island microarrays obtained from the Ontario University Health Network (http://www.microarrays.ca). Hybridizations were performed as described previously (Oberley et al. 2003
). The hybridized microarrays were analyzed using the Genepix Pro 4.1 (Axon Instruments) software package. This provided a set of raw values for each feature on each array. To identify clones that are selectively enriched during ChIP relative to the starting population, the Cy5 and Cy3 channels were normalized across the entire array, by taking the ratio of the medians for all features and normalizing them to unity. After normalization, features of poor intensity (<500) and those that had obvious blemishes were manually flagged and removed from the putative positive list. For all quality features, a ratio was generated for each feature that was the intensity in the Cy5 channel minus background divided by the intensity in the Cy3 channel minus background. Features with Suz12/Input intensity ratios above 3 and PI/Input intensity ratios below 1.3 were selected for further analysis.
| Acknowledgments |
|---|
|
|
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| Footnotes |
|---|
Article and publication are at http://www.genesdev.org/cgi/doi/10.1101/gad.1200204.
4 Corresponding author. E-MAIL farnham{at}oncology.wisc.edu; FAX (608) 262-2824. ![]()
| References |
|---|
|
|
|---|
Bellefroid, E.J., Bourguignon, C., Hollemann, T., Ma, Q., Anderson, D.J., Kintner, C., and Pieler, T. 1996. X-MyT1, a Xenopus C2HC-type zinc finger protein with a regulatory function in neuronal differentiation. Cell 87: 1191-1202.[CrossRef][Medline]
Bienz, M. and Clevers, H. 2000. Linking colorectal cancer to Wnt signaling. Cell 103: 311-320.[CrossRef][Medline]
Bloyer, S., Cavalli, G., Brock, H.W., and Dura, J.M. 2003. Identification and characterization of polyhomeotic PREs and TREs. Dev. Biol. 261: 426-442.[CrossRef][Medline]
Bracken, A.P., Pasini, D., Capra, M., Prosperini, E., Colli, E., and Helin, K. 2003. EZH2 is downstream of the pRB-E2F pathway, essential for proliferation and amplified in cancer. EMBO J 22: 5323-5335.[CrossRef][Medline]
Breiling, A., Turner, B.M., Bianchi, M.E., and Orlando, V. 2001. General transcription factors bind promoters repressed by Polycomb group proteins. Nature 412: 651-655.[CrossRef][Medline]
Brock, H.W. and van Lohuizen, M. 2001. The Polycomb groupNo longer an exclusive club? Curr. Opin. Genet. Dev. 11: 175-181.[CrossRef][Medline]
Cao, R., Wang, L., Xia, L., Erdjument-Bromage, H., Tempst, P., Jones, R.S., and Zhang, Y. 2002. Role of histone H3 lysine 27 methylation in Polycomb-group silencing. Science 298: 1039-1043.
Czermin, B., Melfi, R., McCabe, D., Seitz, V., Imhof, A., and Pirrotta, V. 2002. Drosophila enhancer of Zeste/ESC complexes have a histone H3 methyltransferase activity that marks chromosomal Polycomb sites. Cell 111: 185-196.[CrossRef][Medline]
Dellino, G.I., Schwartz, Y.B., Farkas, G., McCabe, D., Elgin, S.C.R., and Pirrotta, V. 2004. Polycomb silencing blocks transcription initiation. Mol. Cell. 13: 887-893.[CrossRef][Medline]
Graveel, C.R., Jatkoe, T., Madore, S.J., Holt, A.L., and Farnham, P.J. 2001. Expression profiling and identification of novel genes in hapatocellular carcinomas. Oncogene 20: 2704-2712.[CrossRef][Medline]
Jacobs, J.J. and van Lohuizen, M. 1999. Cellular memory of transcriptional states by Polycomb-group proteins. Semin. Cell Dev. Biol. 10: 227-235.[CrossRef][Medline]
____. 2002. Polycomb repression: From cellular memory to cellular proliferation and cancer. Biochim. Biophys. Acta 1602: 151-161.[Medline]
Jurgens, G. 1985. A group of genes controlling the spatial expression of the bithorax complex in Drosophila. Nature 316: 153-155.[CrossRef]
Kim, J.G. and Hudson, L.D. 1992. Novel member of the zinc finger superfamily: A C2-HC finger that recognizes a gliaspecific gene. Mol. Cell. Biol. 12: 5632-5639.
Kirmizis, A., Bartley, S.M., and Farnham, P.J. 2003. Identification of the Polycomb group protein SU(Z)12 as a potential molecular target for human cancer therapy. Mol. Cancer Ther. 2: 113-121.
Kleer, C.G., Cao, Q., Varambally, S., Shen, R., Ota, I., Tomlins, S.A., Ghosh, D., Sewalt, R.G., Otte, A.P., Hayes, D.F., et al. 2003. EZH2 is a marker of aggressive breast cancer and promotes neoplastic transformation of breast epithelial cells. Proc. Natl. Acad. Sci. 100: 11606-11611.
Koontz, J.I., Soreng, A.L., Nucci, M., Kuo, F.C., Pauwels, P., van den Berghe, H., Dal Cin, P., Fletcher, J.A., and Sklar, J. 2001. Frequent fusion of the JAZF1 and JJAZ1 genes in endometrial stromal tumors. Proc. Natl. Acad. Sci. 98: 6348-6353.
Kuzmichev, A., Nishioka, K., Erdjument-Bromage, H., Tempst, P., and Reinberg, D. 2002. Histone methyltransferase activity associated with a human multiprotein complex containing the Enhancer of zeste protein. Genes & Dev. 16: 2893-2905.
Kuzmichev, A., Jenuwein, T., Tempst, P., and Reinberg, D. 2004. Different Ezh2-containing complexes target methylation of histone H1 or nucleosomal histone H3. Mol. Cell 14: 183-193.[CrossRef][Medline]
LaJeunesse, D. and Shearn, A. 1996. E(z): A polycomb group gene or a trithorax group gene? Development 122: 2189-2197.[Abstract]
Muller, J., Hart, C.M., Francis, N.J., Vargas, M.L., Sengupta, A., Wild, B., Miller, E.L., O'Connor, M.B., Kingston, R.E., and Simon, J.A. 2002. Histone methyltransferase activity of a Drosophila Polycomb group repressor complex. Cell 111: 197-208.[CrossRef][Medline]
Oberley, M.J., Inman, D., and Farnham, P.J. 2003. E2F6 negatively regulates BRCA1 in human cancer cells without methylation of histone H3 on lysine 9. J. Biol. Chem. 278: 42466-42476.
Oberley, M.J., Tsao, J., Yau, P., and Farnham, P.J. 2004. High throughput screening of chromatin immunoprecipitates using CpG island microarrays. Methods Enzymol. 376: 316-335.
Orlando, V. 2003. Polycomb, epigenomes, and control of cell identity. Cell 112: 599-606.[CrossRef][Medline]
Orlando, V., Jane, E.P., Chinwalla, V., Harte, P.J., and Paro, R. 1998. Binding of Trithorax and Polycomb proteins to the bithorax complex: Dynamic changes during early Drosophila embryogenesis. EMBO J. 17: 5141-5150.[CrossRef][Medline]
Pasini, D., Bracken, A.P., and Helin, K. 2004. Polycomb group proteins in cell cycle progression and cancer. Cell Cycle 3: 22-26.[Medline]
Pirrotta, V., Poux, S., Melfi, R., and Pilyugin, M. 2003. Assembly of Polycomb complexes and silencing mechanisms. Genetica 117: 191-197.[CrossRef][Medline]
Poux, S., Melfi, R., and Pirrotta, V. 2001. Establishment of Polycomb silencing requires a transient interaction between PC and ESC. Genes & Dev. 15: 2509-2514.
Ren, B., Robert, F., Wyrick, J.J., Aparicio, O., Jennings, E.G., Simon, I., Zeitlinger, J., Schreiber, J., Hannett, N., Kanin, E., et al. 2000. Genome-wide location and function of DNA binding proteins. Science 290: 2306-2309.
Saurin, A.J., Shao, Z., Erdjument-Bromage, H., Tempst, P., and Kingston, R.E. 2001. A Drosophila Polycomb group co