Abstract
The H3K9me2 histone methyltransferases G9a and GLP repress Mage-a class cancer germ-line (CG) antigen gene expression in murine embryonic stem (ES) cells, but the role of these enzymes in CG antigen gene regulation in human cancer cells is unknown. Here we show that whereas independent or dual knockdown of G9a and GLP in human cancer cells leads to reduced global and CG antigen promoter–associated H3K9me2 levels, it does not activate CG antigen gene expression. Moreover, CG antigen gene repression is maintained following pharmacologic targeting of G9a or treatment of G9a knockdown cells with the histone deacetylase inhibitor trichostatin A. However, G9a knockdown cells display increased sensitivity to CG antigen gene activation mediated by the DNA methyltransferase inhibitor decitabine. To account for these findings, we examined DNA methylation at CG antigen gene promoters in both cell types. We found robust DNA hypomethylation in G9a/GLP targeted murine ES cells but a lack of DNA methylation changes in G9a/GLP targeted human cancer cells; intriguingly, this distinction also extended to markers of global DNA methylation. These data reveal that G9a/GLP is required for DNA methylation of CG antigen genes and genomic DNA in murine ES cells, but not human cancer cells, and implicate DNA methylation status as the key epigenetic mechanism involved in CG antigen gene repression. (Mol Cancer Res 2009;7(6):851–62)
Introduction
Cancer germ-line (CG) antigens (also known as cancer-testis antigens) are targets of tumor vaccines currently undergoing clinical evaluation worldwide (1). The rationale for vaccine-based targeting of CG antigens in cancer includes both the inherent immunogenicity of these antigens and their restricted expression in normal human tissues as well as widespread expression in human tumors (1, 2). CG antigen genes seem to be principally regulated at the transcriptional level by epigenetic signals. Lines of evidence for this are as follows: (a) Promoter DNA hypomethylation correlates with CG antigen gene expression in human tumor and normal tissues (3, 4); (b) treatment of nonexpressing cancer cells with epigenetic modulatory drugs, including DNA methyltransferase (DNMT) and histone deacetylase (HDAC) inhibitors, induces CG antigen gene expression (5-7); and (c) specific histone H3 modification patterns are associated with CG antigen gene expression status in human cancer cells (8-10). More detailed elucidation of the epigenetic mechanisms regulating CG antigen gene expression in human cancer cells and tumors should help facilitate the clinical development of CG antigen-targeted therapeutic approaches (2).
The enzymes that mediate histone H3 lysine 9 (H3K9) methylation have been identified, and these include class Suv39 enzymes, which target heterochromatic loci and catalyze H3K9 trimethylation (H3K9me3), and G9a and GLP (also known as EuHMTaseI), which target euchromatic loci and catalyze H3K9 dimethylation (H3K9me2; refs. 11-14). Similar to DNA methyltransferase I (DNMT1) knockout, G9a or GLP knockout causes embryonic lethality in mice, implicating a crucial role for epigenetic repression in early mammalian development (12, 13, 15). Gene expression profiling of G9a knockout murine embryonic stem (ES) cell lines identified Mage-a class CG antigen genes as targets of epigenetic repression by G9a, and subsequent experiments showed a similar role for GLP (12, 13). H3K9me2 levels are reduced both genome-wide and at Mage-a promoters in G9a or GLP knockout ES cells (12, 13). These data suggest that H3K9me2 may play a primary role in mediating CG antigen gene repression in murine ES cells, which seems to conflict with the prevailing view that DNA methylation is the primary mediator of CG antigen gene silencing in human cancer (1, 4, 8).
We previously established the human colon adenocarcinoma cell lines HCT116 and RKO as useful models to study the mechanisms of CG antigen gene repression in human cancer (5, 8). CG antigen genes are transcriptionally silenced in HCT116 and RKO cells by DNA hypermethylation and are activated by pharmacologic or genetic targeting of DNMT enzymes in these cells (5, 8). Moreover, CG antigen gene induction in DNMT-deficient HCT116 cells coincides with both DNA hypomethylation and remodeling of histone code patterns at CG antigen gene promoters (8). The fact that CG antigen genes are infrequently expressed in human colorectal cancer also supports the relevance of this experimental model (1).
In the current study, we examined the role of G9a and GLP in regulating CG antigen expression in human cancer cells, using HCT116 and RKO colorectal cancer cells as models. We hypothesized that G9a/GLP may repress CG antigen genes in nonexpressing human cancer cell lines because H3K9me2 and DNA methylation signals are often interdependent and G9a propagates epigenetic repression at the mammalian replication fork via its interaction with DNMT1 (16). Our data reveal that in human cancer cells, unlike mouse ES cells, G9a and GLP are dispensable for CG antigen gene repression, despite playing a role in the maintenance of H3K9 methylation in both cell types. Most importantly, we find that CG antigen genes are DNA hypomethylated in G9a or GLP null murine ES cells, but their DNA methylation status is unchanged in human cancer cells following G9a and/or GLP targeting. These data reinforce the notion that epigenetic control processes are regulated in distinct ways in ES cells and cells of somatic origin, including cancer cells. Furthermore, they implicate DNA methylation as the primary repression mechanism for CG antigen gene expression.
Results
CG Antigen Gene Expression and Histone Modifications in G9a Knockdown Human Cancer Cells
To assess the role of G9a in CG antigen gene repression in human cancer cells, we initially developed an effective transient siRNA G9a targeting strategy in RKO cells (Fig. 1A). Low siRNA transfection efficiency prevented implementation of this strategy in HCT116 cells (data not shown). Reverse transcription-PCR (RT-PCR) analysis revealed that efficient G9a knockdown did not induce expression of three representative CG antigen genes, MAGE-A1, XAGE-1, and NY-ESO-1 (8), in RKO cells (Fig. 1B). Because sustained knockdown might be required to affect CG antigen gene expression, and to validate these findings in another cell type, we established clonal G9a stable knockdown cell lines in both HCT116 and RKO cells (Fig. 1C). Interestingly, during derivation of the these cell lines, we observed that cancer cells expressing G9a shRNA showed substantially reduced cell clonogenicity, particularly in RKO cells (Fig. 1D; Table 1). We also observed reduced cell viability following siRNA-mediated G9a knockdown in RKO (data not shown). These data suggest that G9a contributes to cell growth and/or survival in colorectal cancer cells, consistent with previous reports using mice and prostate cancer cells (12, 17). The role of G9a in cell survival may explain why it was difficult to achieve complete G9a knockdown in human cancer cells (Fig. 1A and C).
Transient and stable G9a knockdown in human cancer cells. A. Transient siRNA knockdown of G9a in RKO cells. Control and G9a-specific siRNAs were administered at either 5 or 50 nmol/L; cells were harvested at day 5; and G9a expression was measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of B. B. End-point RT-PCR analysis of MAGE-A1, NY-ESO-1, and XAGE-1 for the experiment shown in A. GAPDH expression was analyzed to control for cDNA input, and amplification of cDNA from RKO cells treated with 5 μmol/L 5-aza-2′-deoxycytidine (decitabine; DAC) served as a positive control for CG antigen gene expression. CG antigen gene amplification was done for 35 cycles, whereas GAPDH was done for 25 cycles. The sample key is shown at the bottom. LF2000, Lipofectamine-only transfection control. C. Western blot analysis of G9a protein expression in control shRNA (shCon) and G9a shRNA (shG9a) stable RKO and HCT116 cell lines. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. D. Cells stably expressing shG9a have reduced clonogenicity. Representative images of crystal violet-stained colonies of control and G9a knockdown stable RKO cells.
Transient and stable G9a knockdown in human cancer cells. A. Transient siRNA knockdown of G9a in RKO cells. Control and G9a-specific siRNAs were administered at either 5 or 50 nmol/L; cells were harvested at day 5; and G9a expression was measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of B. B. End-point RT-PCR analysis of MAGE-A1, NY-ESO-1, and XAGE-1 for the experiment shown in A. GAPDH expression was analyzed to control for cDNA input, and amplification of cDNA from RKO cells treated with 5 μmol/L 5-aza-2′-deoxycytidine (decitabine; DAC) served as a positive control for CG antigen gene expression. CG antigen gene amplification was done for 35 cycles, whereas GAPDH was done for 25 cycles. The sample key is shown at the bottom. LF2000, Lipofectamine-only transfection control. C. Western blot analysis of G9a protein expression in control shRNA (shCon) and G9a shRNA (shG9a) stable RKO and HCT116 cell lines. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. D. Cells stably expressing shG9a have reduced clonogenicity. Representative images of crystal violet-stained colonies of control and G9a knockdown stable RKO cells.
Quantification of Colony Number in Control and G9a Knockdown Stable RKO and HCT116 Cell Lines
Cell Line . | Average No. of Colonies . | %Change . | |
---|---|---|---|
shCon (n = 2) . | shG9a (n = 2) . | ||
RKO | 1,521 | 894 | ↓41.2 |
HCT116 | 643 | 598 | ↓7.0 |
Cell Line . | Average No. of Colonies . | %Change . | |
---|---|---|---|
shCon (n = 2) . | shG9a (n = 2) . | ||
RKO | 1,521 | 894 | ↓41.2 |
HCT116 | 643 | 598 | ↓7.0 |
To explore the phenotype of stable G9a knockdown human cancer cells, we used Western blot analysis to measure global H3K9me2 levels and H3K27me2 levels because these two modifications have been previously linked to G9a function (18, 19). Notably, H3K9me2 levels were reduced in both cell types whereas H3K27me2 levels were not significantly altered (Fig. 2A). A lack of activity of G9a toward H3K27me2 is consistent with the results of recent studies of G9a function in different model systems (20-23). To assess whether G9a knockdown altered histone modifications at CG antigen gene promoters, we used quantitative chromatin immunoprecipitation (ChIP)-PCR to determine H3K9me2 levels at the MAGE-A1 and XAGE-1 promoters as described previously (8). Similar to the effect on global levels, H3K9me2 levels are reduced at both MAGE-A1 and XAGE-1 in G9a stable knockdown HCT116 and RKO cells (Fig. 2B). In contrast, H3K27me2 levels are unchanged at these promoters (Fig. 2C). Whereas reduced H3K9me2 suggested that gene expression might be altered, CG antigen genes remained repressed in both G9a knockdown cell lines (Fig. 2D). This observation was repeatable as independently derived stable G9a knockdown RKO and HCT116 cell clones also did not express CG-X genes (data not shown). Stable G9a knockdown cells also showed slight increases in the global level of the active histone mark H3K9ac (Fig. 3A). However, consistent with a lack of CG antigen expression in the stable G9a knockdown cells, ChIP analyses revealed that the levels of H3K9ac and H3K4me2 at CG gene promoters, which are associated with CG gene expression (8), were not increased (Fig. 3B-E).
Characterization of stable G9a knockdown human cancer cells. A. Western blot analysis of H3K9me2 and H3K27me2 levels in control shRNA– and G9a shRNA–expressing stable cell lines. Coomassie staining confirmed equivalent protein input, and band densitometry was done as described in Materials and Methods. B and C. Quantitative ChIP-PCR analysis of H3K9me2 (B) and H3K27me2 (C) levels at the MAGE-A1 and XAGE-1 5′ CpG island regions. Bars, 1 SD. D. RT-PCR analysis of MAGE-A1, XAGE-1, and NY-ESO-1 expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom.
Characterization of stable G9a knockdown human cancer cells. A. Western blot analysis of H3K9me2 and H3K27me2 levels in control shRNA– and G9a shRNA–expressing stable cell lines. Coomassie staining confirmed equivalent protein input, and band densitometry was done as described in Materials and Methods. B and C. Quantitative ChIP-PCR analysis of H3K9me2 (B) and H3K27me2 (C) levels at the MAGE-A1 and XAGE-1 5′ CpG island regions. Bars, 1 SD. D. RT-PCR analysis of MAGE-A1, XAGE-1, and NY-ESO-1 expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom.
H3K9ac and H3K4me2 in G9a knockdown human cancer cells. A. Western blot analysis of H3K9ac and H3K4me2 levels in control shRNA– and G9a shRNA–expressing stable RKO and HCT116 cell lines. Coomassie staining confirmed equivalent protein input, and band densitometry was done as described in Materials and Methods. B and C. Quantitative ChIP-PCR analysis of H3K9ac levels at the MAGE-A1 and XAGE-1 5′ CpG island regions in G9a targeted RKO (B) and HCT116 (C) cell lines. D and E. Quantitative ChIP-PCR analysis of H3K4me2 levels at the MAGE-A1 and XAGE-1 5′ CpG island regions in G9a targeted RKO (D) and HCT116 (E) cell lines. RKO cells treated with decitabine (DAC) served as a positive control for both histone modifications in RKO cells, whereas DNMT1−/−, DNMT3b−/− HCT116 cells (DKO) served as a positive control for both histone modifications in HCT116. Bars, 1 SD.
H3K9ac and H3K4me2 in G9a knockdown human cancer cells. A. Western blot analysis of H3K9ac and H3K4me2 levels in control shRNA– and G9a shRNA–expressing stable RKO and HCT116 cell lines. Coomassie staining confirmed equivalent protein input, and band densitometry was done as described in Materials and Methods. B and C. Quantitative ChIP-PCR analysis of H3K9ac levels at the MAGE-A1 and XAGE-1 5′ CpG island regions in G9a targeted RKO (B) and HCT116 (C) cell lines. D and E. Quantitative ChIP-PCR analysis of H3K4me2 levels at the MAGE-A1 and XAGE-1 5′ CpG island regions in G9a targeted RKO (D) and HCT116 (E) cell lines. RKO cells treated with decitabine (DAC) served as a positive control for both histone modifications in RKO cells, whereas DNMT1−/−, DNMT3b−/− HCT116 cells (DKO) served as a positive control for both histone modifications in HCT116. Bars, 1 SD.
Because G9a knockdown was robust but not complete in G9a shRNA knockdown clones, we further targeted G9a in stable G9a knockdown RKO cells using transient siRNA transfection (Fig. 4A). Whereas G9a siRNA transfection resulted in further diminishment of G9a protein, CG-X genes remained repressed (Fig. 4A and B).
CG-X antigen genes are repressed in RKO cells following dual shRNA/siRNA G9a knockdown or dual G9a/GLP knockdown. A. RKO cells stably expressing control shRNA or G9a shRNA were mock transfected or transfected with a G9a-specific siRNA. Cells were harvested 5 d after transfection and G9a protein levels were measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of B. B. RT-PCR analysis of CG antigen gene expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom. C. RKO cells stably expressing control shRNA or G9a shRNA were mock transfected or transfected with a GLP-specific siRNA. Five days after transfection, G9a and GLP protein levels were measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of D. D. RT-PCR analysis of CG antigen gene expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom.
CG-X antigen genes are repressed in RKO cells following dual shRNA/siRNA G9a knockdown or dual G9a/GLP knockdown. A. RKO cells stably expressing control shRNA or G9a shRNA were mock transfected or transfected with a G9a-specific siRNA. Cells were harvested 5 d after transfection and G9a protein levels were measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of B. B. RT-PCR analysis of CG antigen gene expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom. C. RKO cells stably expressing control shRNA or G9a shRNA were mock transfected or transfected with a GLP-specific siRNA. Five days after transfection, G9a and GLP protein levels were measured by Western blot. α-Tubulin expression was measured as control for protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of D. D. RT-PCR analysis of CG antigen gene expression. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown at the bottom.
Dual Targeting of G9a and EuHMTaseI/GLP in Human Cancer Cells
GLP (EuHMTaseI) has similar substrate specificity as G9a and has a nonredundant and essential role in mammalian development (13). Despite this nonredundancy in mouse development, it remained plausible that GLP could compensate for G9a loss in human cancer cells with regard to CG antigen gene repression. To test this possibility, we used a siRNA approach. Treatment of RKO cells with GLP siRNA resulted in effective GLP knockdown and, notably, also reduced G9a protein levels (Fig. 4C). This may be explained by the fact that GLP and G9a form an intracellular complex that stabilizes the G9a protein (13). Consistent with this idea, GLP knockdown did not affect G9a mRNA expression (data not shown). Despite effective reduction of both G9a and GLP in RKO cells, CG antigen genes remained repressed (Fig. 4D). Furthermore, we did not observe an additional decrease of H3K9me2 levels at CG antigen gene promoters following GLP knockdown in stable G9a knockdown cells (data not shown). These data suggest that redundancy between G9a and GLP does not account for the repression of CG antigen genes in cancer cells sustaining G9a knockdown.
Pharmacologic Targeting of G9a in Human Cancer Cells
Recently, a novel small-molecule inhibitor of G9a was identified by high-throughput screening (23). This agent, BIX-01294, was shown to selectively target G9a and, to a lesser extent, GLP, and to reduce H3K9me2 levels in mammalian cells at low micromolar concentrations (23). Notably, BIX-01294 treatment was also shown to activate Mage-A2 expression in wild-type mouse ES cells (23). Here, we sought to complement the genetic knockdown approach described above by examining the effect of BIX-01294 treatment on CG antigen gene expression in human cancer cells. BIX-01294 treatment reduced global H3K9me2 levels in both control and G9a stable knockdown HCT116 and RKO cells, validating the activity of this agent in human cancer cells (Fig. 5A and B). However, BIX-01294 treatment failed to induce CG-X gene expression in either cell type (Fig. 5C and D). These data provide additional evidence that G9a is not required for CG antigen gene repression in human cancer cells.
Combined genetic and pharmacologic targeting of G9a in human cancer cells. A and B. RKO (A) and HCT116 (B) cells were treated with DMSO (vehicle) or BIX-01294; cells were harvested 48 h later; and H3K9me2 levels were determined by Western blot analyses. Coomassie staining confirmed equal protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of the figure. C and D. RT-PCR analysis of CG antigen gene expression in the experiment described in A and B. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown below the figure.
Combined genetic and pharmacologic targeting of G9a in human cancer cells. A and B. RKO (A) and HCT116 (B) cells were treated with DMSO (vehicle) or BIX-01294; cells were harvested 48 h later; and H3K9me2 levels were determined by Western blot analyses. Coomassie staining confirmed equal protein input, and band densitometry was done as described in Materials and Methods. The sample key is shown at the bottom of the figure. C and D. RT-PCR analysis of CG antigen gene expression in the experiment described in A and B. PCR conditions and controls are the same as described in Fig. 1B, and the sample key is shown below the figure.
Effect of Combined G9a Knockdown and HDAC Inhibitor Treatment on CG Antigen Gene Expression in Human Cancer Cells
The observation that G9a knockdown results in reduced H3K9me2 at CG antigen promoters without altering H3K9ac raised the question of whether forced induction of H3K9ac in G9a knockdown cells could elicit CG antigen expression. To test this, we used the classic HDAC inhibitor trichostatin A (TSA) to treat control and stable G9a knockdown cells. As anticipated, TSA treatment markedly increased global H3 acetylation and H3K9 acetylation levels in both control and stable G9a knockdown cells (Fig. 6A). Moreover, TSA induced H3K9ac levels at CG antigen promoters (Fig. 6B and C). Despite this effect, MAGE-A1 remained silent in both RKO and HCT116 cells, whereas XAGE-1 remained silent in RKO (Fig. 6D). For XAGE-1 expression in HCT116, TSA treatment led to a moderate level of induction, suggesting that histone acetylation regulates this gene in certain cellular contexts (Fig. 6D). Taken together, however, these data suggest that remodeling the histone code at CG antigen gene promoters is generally insufficient to provoke CG antigen gene expression. A possible explanation for this phenomenon is that CG antigen gene promoters are densely DNA methylated in these cell types (8).
Effect of combined G9a knockdown and HDAC inhibitor treatment on CG antigen gene expression. RKO and HCT116 stable control shRNA and stable G9a knockdown cells were treated with the HDAC inhibitor TSA at 600 nmol/L concentration for 24 h before harvesting. A. Histone protein extracts were prepared, as described in Materials and Methods, and used for Western blot analysis of acetylated histone H3 or acetylated histone H3K9 levels. Coomassie staining confirmed equivalent protein loading. B and C. Acetylated H3K9 levels at the MAGE-A1 (B) and XAGE-1 (C) promoters were measured by quantitative ChIP-PCR. D. MAGE-A1 and XAGE-1 mRNA expression was measured by end-point RT-PCR. n and p, negative (no template) and positive (decitabine-treated cells) PCR controls, respectively.
Effect of combined G9a knockdown and HDAC inhibitor treatment on CG antigen gene expression. RKO and HCT116 stable control shRNA and stable G9a knockdown cells were treated with the HDAC inhibitor TSA at 600 nmol/L concentration for 24 h before harvesting. A. Histone protein extracts were prepared, as described in Materials and Methods, and used for Western blot analysis of acetylated histone H3 or acetylated histone H3K9 levels. Coomassie staining confirmed equivalent protein loading. B and C. Acetylated H3K9 levels at the MAGE-A1 (B) and XAGE-1 (C) promoters were measured by quantitative ChIP-PCR. D. MAGE-A1 and XAGE-1 mRNA expression was measured by end-point RT-PCR. n and p, negative (no template) and positive (decitabine-treated cells) PCR controls, respectively.
G9a Knockdown Sensitizes Human Cancer Cells to 5-Aza-2′-Deoxycytidine (Decitabine)–Mediated CG Antigen Gene Activation
Whereas the results above indicate that inhibition of G9a/GLP is not sufficient to provoke CG antigen gene expression in human cancer cells, either alone or in combination with TSA, it remained plausible that G9a targeting could potentiate CG antigen gene activation by DNMT inhibitors. A similar relationship exists between HDAC inhibitors and DNMT inhibitors in that HDAC inhibition can facilitate the activation of methylation silenced genes by the DNMT inhibitor 5-aza-2′-deoxycytidine (24). We thus tested the response of control and G9a knockdown human cancer cells to treatment with decitabine. Notably, stable G9a knockdown HCT116 cells treated with decitabine showed significantly increased MAGE-A1 and XAGE-1 gene expression and moderately increased NY-ESO-1 expression as compared with control cells treated with decitabine alone (Fig. 7A-C). In contrast to HCT116, an enhanced CG antigen expression effect was only marginal in decitabine-treated G9a stable knockdown RKO cells (data not shown). These data show that G9a knockdown can synergize with DNMT inhibition to promote CG antigen gene expression, but that this effect is, in part, cell context dependent.
Stable G9a knockdown HCT116 cells show increased sensitivity to decitabine-mediated CG antigen gene activation. HCT116 cell lines stably expressing control or G9a-specific shRNAs were treated with the indicated concentrations of decitabine (5-aza-CdR) for 48 h; RNA extracts were harvested; and MAGE-A1 (A), XAGE-1 (B), and NY-ESO-1 (C) expression was measured using quantitative RT-PCR. Expression data from all three genes were normalized to GAPDH as a cDNA input control. Columns, mean expression values of triplicate data points; bars, 1 SD.
Stable G9a knockdown HCT116 cells show increased sensitivity to decitabine-mediated CG antigen gene activation. HCT116 cell lines stably expressing control or G9a-specific shRNAs were treated with the indicated concentrations of decitabine (5-aza-CdR) for 48 h; RNA extracts were harvested; and MAGE-A1 (A), XAGE-1 (B), and NY-ESO-1 (C) expression was measured using quantitative RT-PCR. Expression data from all three genes were normalized to GAPDH as a cDNA input control. Columns, mean expression values of triplicate data points; bars, 1 SD.
CG Antigen Gene Regulation and DNA Methylation in G9a and GLP Knockout Mouse ES Cells
In contrast to our data using human cancer cells, G9a and GLP knockout mouse ES cells show robust expression of Mage-a class CG antigen genes (12, 13). Because loss of H3K9me2 globally and at CG antigen gene promoters is observed in both mouse ES cells and human cancer cells responding to G9a and/or GLP targeting (refs. 12, 13, 17; Fig. 2), reduction of this repressive epigenetic histone mark does not seem to explain the differential response. Of note, previous studies have observed altered DNA methylation in G9a knockout mouse ES cells (18, 25, 26). We hypothesized that differences in promoter DNA methylation status at CG antigen genes might explain the distinct response of mouse ES and human cancer cells to G9a/GLP targeting. To test this, we initially used G9a and GLP knockout murine ES cell lines (Fig. 8A). Using pan-Mage-a primers that amplify murine Mage-a2, Mage-a5, Mage-a6, and Mage-a8, we observed that G9a and GLP knockout ES cells express high levels of Mage-a genes (Fig. 8B). Of the murine Mage-a gene family members, only Mage-a2 and Mage-a8 contain classic 5′ CpG islands (27), and thus we focused further attention on these two genes. Both Mage-a2 and Mage-a8 were induced in G9a and GLP knockout ES cells (Fig. 8B). We next used methylation-specific PCR to determine the methylation status of the 5′ CpG region of Mage-a2 in wild-type, G9a knockout, and GLP knockout ES cells. As shown in Fig. 8C, Mage-a2 was hypermethylated in wild-type ES cells and becomes dramatically hypomethylated in GLP and G9a knockout cells. Moreover, reintroduction of either GLP or G9a into the corresponding knockout cell line led to partial remethylation of Mage-a2 (Fig. 8C). This remethylation correlated with repression of Mage-a2 (12, 13). Similar to the Mage-a2 result, bisulfite pyrosequencing analysis indicated that Mage-a8 was dramatically hypomethylated in G9a and GLP null ES cells and remethylated following reintroduction of the relevant protein (Fig. 8D). To determine if the effect on DNA methylation in G9a and GLP knockout ES cells extended to markers of global DNA hypomethylation, we used two different measures of global DNA methylation status: (a) 5-methyl-2′-deoxycytidine (5mdC) levels and (b) bisulfite pyrosequencing of the murine repetitive element B1 (28). Both knockout cell types display sizable reduction of 5mdC levels, which was partially restored following reintroduction of the corresponding proteins (Fig. 8E). A similar phenotype was observed for methylation of the repetitive B1 element (Fig. 8F). These data reveal that G9a and GLP are required for DNA methylation of CG antigen genes as well as other regions of the genome in murine ES cells, and provide a plausible explanation for the activation of CG antigen genes in this cell type.
CG antigen gene expression and DNA methylation in G9a and GLP knockout mouse ES cells. A. Western blot analysis of G9a and GLP protein expression. α-Tubulin expression was measured as a control for protein input. B. RT-PCR analysis of Mage-a, Mage-a2, and Mage-a8 expression. Gapdh was amplified as a control for cDNA input. All reactions used 30 PCR cycles. Specific amplification of Mage-a2 and Mage-a8 was confirmed by DNA sequencing of the PCR amplicons. C. Methylation-specific PCR analysis of Mage-a2 5′ region methylation. Results are shown for control DNAs, wild-type ES cells, and G9a and GLP knockout ES cells before and after reintroduction of the knocked out protein. U and M lanes correspond to PCR amplifications specific for unmethylated and methylated DNA, respectively. U and M control DNAs are described in Materials and Methods. D. Bisulfite pyrosequencing analysis of Mage-a8 5′ region methylation. Samples are as described in C. E. Liquid chromatography mass spectrometry analysis of total genomic 5mdC. Samples are same as described in C. F. Bisulfite pyrosequencing analysis of B1 repetitive element methylation. Samples are as described in C. Statistical analysis shown in D to F used one-tailed t test (GraphPad Prism) and the relevant P values are indicated in the figure.
CG antigen gene expression and DNA methylation in G9a and GLP knockout mouse ES cells. A. Western blot analysis of G9a and GLP protein expression. α-Tubulin expression was measured as a control for protein input. B. RT-PCR analysis of Mage-a, Mage-a2, and Mage-a8 expression. Gapdh was amplified as a control for cDNA input. All reactions used 30 PCR cycles. Specific amplification of Mage-a2 and Mage-a8 was confirmed by DNA sequencing of the PCR amplicons. C. Methylation-specific PCR analysis of Mage-a2 5′ region methylation. Results are shown for control DNAs, wild-type ES cells, and G9a and GLP knockout ES cells before and after reintroduction of the knocked out protein. U and M lanes correspond to PCR amplifications specific for unmethylated and methylated DNA, respectively. U and M control DNAs are described in Materials and Methods. D. Bisulfite pyrosequencing analysis of Mage-a8 5′ region methylation. Samples are as described in C. E. Liquid chromatography mass spectrometry analysis of total genomic 5mdC. Samples are same as described in C. F. Bisulfite pyrosequencing analysis of B1 repetitive element methylation. Samples are as described in C. Statistical analysis shown in D to F used one-tailed t test (GraphPad Prism) and the relevant P values are indicated in the figure.
CG Antigen Gene and Global Genomic DNA Methylation in G9a Knockdown Human Cancer Cells
We next examined DNA methylation in stable G9a and GLP knockdown human cancer cells. We used quantitative bisulfite pyrosequencing to determine the methylation level of the 5′ CpG island promoters of MAGE-A1, NY-ESO-1, and XAGE-1. As expected, in both the RKO and HCT116 cell lines, stable G9a knockdown did not alter CG antigen promoter methylation levels (Fig. 9A and B). Standard sodium bisulfite sequencing and methylation-specific PCR analyses confirmed the pyrosequencing data (data not shown). In addition, stable G9a knockdown in HCT116 and RKO cells did not lead to altered global DNA methylation, as assessed either by 5mdC levels or by pyrosequencing of the human LINE-1 repetitive DNA element (Fig. 9C and D). Consistent with the findings using stable G9a knockdown cells, transient knockdown of G9a or GLP in the stable G9a knockdown RKO cell line or combined stable G9a knockdown and TSA treatment in either HCT116 or RKO also did not alter CG antigen or global DNA methylation (data not shown).
DNA methylation in stable G9a knockdown human cancer cells. A and B. Bisulfite pyrosequencing analysis of MAGE-A1, NY-ESO-1, and XAGE-1 5′ promoter CpG island region methylation in RKO (A) and HCT116 (B) control and stable G9a knockdown cell lines. C and D. Global DNA methylation analyses of RKO (C) and HCT116 (D) control and stable G9a knockdown cell lines. Two measurements of global DNA methylation were done: bisulfite pyrosequencing analysis of LINE-1 repetitive element methylation and liquid chromatography mass spectrometry determination of total 5mdC levels.
DNA methylation in stable G9a knockdown human cancer cells. A and B. Bisulfite pyrosequencing analysis of MAGE-A1, NY-ESO-1, and XAGE-1 5′ promoter CpG island region methylation in RKO (A) and HCT116 (B) control and stable G9a knockdown cell lines. C and D. Global DNA methylation analyses of RKO (C) and HCT116 (D) control and stable G9a knockdown cell lines. Two measurements of global DNA methylation were done: bisulfite pyrosequencing analysis of LINE-1 repetitive element methylation and liquid chromatography mass spectrometry determination of total 5mdC levels.
Discussion
We previously reported that disruption of DNMT enzymes in human cancer cells results in reduced H3K9me2 levels at CG antigen gene promoters, correlating with increased gene expression (8). In addition, it was previously shown that murine ES cells deficient for G9a or GLP express Mage-a class CG antigen genes, coincident with reduced H3K9me2 levels at their promoters (12, 13). These data suggested that H3K9me2 and/or G9a/GLP may play a primary role in CG antigen gene repression in human cancer cells. However, our current data argue against this model. First, either transient or stable G9a knockdown in human cancer cell lines does not induce CG antigen gene expression, despite loss of H3K9me2 both globally throughout the genome and specifically at CG antigen promoter regions. Second, simultaneous knockdown of G9a and GLP and combined genetic and pharmacologic inhibition of G9a also failed to induce CG gene expression. Third, combined G9a knockdown and HDAC inhibitor treatment generally failed to induce CG gene expression. Moreover, we previously found that unlike double DNMT1/3b knockout cells, single DNMT1−/− or DNMT3b−/− HCT116 cells show reduced H3K9me2 at CG-X gene promoters in the absence of DNA methylation changes or expression of CG genes (8). Collectively, these data suggest that G9a/GLP knockdown and/or altered histone modification states are not sufficient alone to reactivate the expression of CG antigen genes in human cancer cells. Moreover, our data argue that promoter DNA hypomethylation is the chief epigenetic signal leading to CG antigen expression in human cancer cells. In agreement with our findings for CG antigen genes, it was recently reported that G9a knockdown in human prostate cancer cells leads to the up-regulation of very few genes in a genome-wide microarray screen and does not cause DNA hypomethylation of methylation-silenced tumor suppressor genes (17).
Based on our data from human cancer cells, we reexamined the mechanism of CG antigen gene activation in G9a- or GLP-deficient murine ES cells (12, 13). We hypothesized that G9a or GLP loss could lead to CG gene promoter DNA hypomethylation in murine ES cells, but not in human cancer cells, thus explaining the differential outcomes in the two systems. Our data validated this hypothesis. In both G9a and GLP null ES cells, the Mage-a2 and Mage-a8 promoters were severely hypomethylated coincident with gene activation. That this effect was specific was validated by experiments that showed that reintroduction of wild-type G9a or GLP led to remethylation and silencing of these genes. Interestingly, G9a and GLP null ES cells also displayed both a global loss of methylation at the non–long terminal repeat (LTR) short interspersed repetitive element (SINE) B1 and a reduction of total genomic 5mdC levels. In stark contrast, in human cancer cells in which G9a and/or GLP was down-regulated by RNA interference, both CG promoter and global DNA methylation levels (either of the LINE-1 element or 5mdC levels) were unaltered. This was true despite significant reduction of both global and promoter-specific H3K9me2 levels in these cells. These data emphasize a primary regulatory role for DNA methylation in CG antigen gene regulation. The apparently tight dependence of CG gene silencing on DNA methylation may be accounted for by direct inhibition of Ets transcription factor binding, the recruitment of methylated DNA binding proteins at CG antigen gene promoter regions, or other as yet undefined mechanisms (5, 29, 30).
Our finding of DNA hypomethylation in G9a/GLP null murine ES cells is consistent with other reports (18, 25, 26). DNA hypomethylation in G9a null ES cells was initially described for the Prader-Willi imprinting center, and it was later shown that methylation of the Oct3/4 gene during differentiation is blocked in G9a null ES cells (25, 26). In addition, a recent study reported that 32 of 1,300 (∼2.5%) genomic loci are hypomethylated in G9a−/− ES cells, implying a more widespread effect (18). Here we show that G9a or GLP null ES cells display significant reductions in 5mdC levels and hypomethylation of the murine B1 SINE repetitive DNA element. Similarly, Dong et al. (31) recently reported that murine LTR and non-LTR LINE elements are DNA hypomethylated in G9a null ES cells. Together, these data suggest that G9a or GLP loss disrupts maintenance and/or de novo DNA methylation. Loss of DNA methylation in G9a and GLP knockout ES cells does not seem to result from decreased expression of Dnmt genes or proteins (31, 32). Instead, the mechanistic basis for the DNA hypomethylation effect may reflect either reduced recruitment of Dnmt1 to target sites during DNA replication or reduced Dnmt3a and/or Dnmt3b targeting to nonreplicating chromatin (16, 31, 32). Intriguingly, recent data suggest that DNA hypomethylation in G9a or GLP null murine ES cells may result from a loss of a function unrelated to their histone methyltransferase activity (31, 32).
It is important to note that although H3K9me2 may not be sufficient for CG antigen gene repression in human cancer cells, it seems to play a secondary role. Notably, we found that certain G9a knockdown human colorectal cancer cells show an increased sensitivity to decitabine-mediated CG antigen gene induction. Additionally, it was recently shown that G9a and DNMT1 knockdown cooperatively facilitate induction of the DNA methylation–silenced tumor suppressor gene MASPIN in breast cancer cells (33). In ES cells, unlike human cancer cells (or other somatic cells), it seems that either H3K9 methylation or DNA methylation is sufficient for transcriptional repression (17, 32). This may reflect a “dual silencing” mechanism to maintain silent states during cellular differentiation, which is uniquely important in ES cells (17, 32).
We show here that murine ES cells are more susceptible to the loss of DNA methylation at CG antigen genes following G9a or GLP loss than are human cancer cells. One possible explanation for the differential effect in the two cell types is based on the structure of CG antigen genes in the two species. In humans, CG antigen genes have promoter regions of relatively high CpG density, whereas murine Mage-a class CG genes generally have CpG poor promoter regions (ref. 8 and data not shown). To mitigate this issue, our study focused on the two murine Mage-a genes (Mage-a2 and Mage-a8) that have CpG dense promoter regions. Therefore, it is unlikely that the CpG density of the studied CG gene promoters accounts for the different results observed in the two cell types. Of note, the orthologous human genes MAGE-A2 and MAGE-A8 were not activated in human cancer cells sustaining G9a knockdown (data not shown). Moreover, the effects of G9a or GLP loss on global DNA methylation (e.g., 5mdC levels) in the two cell types are also distinct, despite widespread conservation of the mouse and human genomes. A second possible explanation for the differential effect relates to the two cell types under study. We speculate that the dramatic loss of DNA methylation in G9a or GLP null ES cells reflects a greater level of “epigenetic plasticity” in this cell type than is seen in somatic cells (including cancer cells). ES cells may require higher epigenetic plasticity because the ultimate fate of these cells is epigenetic reprogramming along specific cell lineages (34). One example of the epigenetic plasticity of ES cells is the bivalent chromatin domains consisting of H3K4me3 and H3K27me3 often found at the promoters of developmentally regulated genes in ES cells (35). Also consistent with this explanation, G9a knockout embryos (somatic cells) do not display DNA hypomethylation at target loci, despite loss of H3K9 methylation at these regions (32). Finally, studies of the effect of DNA methylation on chromatin structure and histone modifications suggest that there are important distinctions between epigenetic regulation in ES cells and somatic cells (36).
In summary, our data emphasize a primary regulatory role for DNA methylation in CG antigen gene regulation in human cancer. There has been increasing interest in the utilization of epigenetic modulatory drugs to promote CG antigen gene expression to augment CG antigen vaccine efficacy in human cancer patients (2, 37). Our data suggest that these approaches may need to include direct inhibition of DNMT enzymes as a required step to promote CG antigen gene expression. However, it remains plausible that agents that modify histone modification status, including HDAC or H3K9 methyltransferase inhibitors, may prove useful to enhance the effect of DNMT inhibitors in this context.
Materials and Methods
Cell Culture and Drug Treatments
Colon adenocarcinoma cell lines RKO and HCT116 (American Type Culture Collection) were cultured as described previously (8). Mouse ES cells [wild-type (TT2), G9a−/− (22-10), GLP−/− (CD-10)] and G9a and GLP reintroduction cell lines were described previously (12, 13). Control or stable G9a knockdown HCT116 or RKO cell lines (described below) were treated with 5-aza-2′-deoxycytidine (decitabine; Sigma) dissolved in PBS and were harvested 48 h posttreatment. In other experiments, cells were treated for 24 h with 600 nmol/L TSA (Sigma) solubilized in DMSO. The small-molecule G9a inhibitor BIX-01294 (23) was solubilized in DMSO and used to treat control or stable G9a knockdown HCT116 or RKO cell lines for 48 h.
Transient and Stable RNA Interference Knockdown Experiments
For transient knockdown of G9a or GLP (EuHMTaseI), RKO cells were transfected with 5 to 50 nmol/L siRNA (Ambion) using Lipofectamine 2000 (Invitrogen). As a negative control, RKO cells were transfected with 5 to 50 nmol/L of siControl Non-Targeting siRNA #2 (Dharmacon). siRNA transfections occurred on day 0 and day 2, and cells were harvested on day 5. For stable knockdown, a shRNA specific for human G9a was designed using the shRNA insert design tool (Ambion). The shRNA duplex was ligated into the pSilencer 3.1-H1 neo vector (Ambion). After sequence verification, the G9a shRNA plasmid or a control shRNA vector (Ambion) was transfected into RKO and HCT116 cells using Lipofectamine 2000 (Invitrogen). Antibiotic selection was done using 1.5 mg/mL (RKO) or 400 μg/mL (HCT116) G418 (Mediatech) and was continued until all mock-transfected control cells were killed. G9a shRNA– and control shRNA–expressing cell clones were picked and expanded while being maintained in normal media supplemented with G418 at one half of the antibiotic kill concentration. In combination shRNA/siRNA experiments, control or stable G9a knockdown RKO cells were transfected with various siRNAs using Lipofectamine 2000 on day 0 and day 2, and cells were harvested on day 5. siRNA and shRNA sequences are shown in Supplementary Table S1.
Colony Formation Assays
Equivalent numbers of RKO and HCT116 cells were transfected with control or G9a shRNA plasmids, and antibiotic selection was done on cells grown in 15-cm tissue culture plates. For colony counts, cells were washed in PBS, fixed in 10% formaldehyde, and stained with 0.1% crystal violet in methanol. Colonies were counted manually.
RT-PCR and Quantitative RT-PCR
Total RNA was extracted using Trizol (Invitrogen), and cDNA was generated using the First Strand cDNA Synthesis Kit (Fermentas). RT-PCR and quantitative RT-PCR amplifications of human MAGE-A1, NY-ESO-1, XAGE-1, and GAPDH were done as described previously (8). Primer sequences for murine Mage-a, Mage-a2, Mage-a8, and Gapdh were designed using Primer 3 (38) and are shown in Supplementary Table S1. Murine pan Mage-a primers were adapted from a previous study (39) and used to amplify mouse Mage-a2, Mage-a5, Mage-a6, and Mage-a8. Primers were obtained from IDT.
Western Blot Analyses of G9a and GLP
Western blots of whole-cell protein extracts were done as described previously (8). G9a was detected using the D141-3 primary antibody (MBL International Corp.) at 1:100 dilution followed by protein A-horseradish peroxidase (Amersham Biosciences) at 1:1,000 dilution. Alternatively, G9a was detected using a polyclonal antibody (kind gift from Dr. Sriharsa Pradhan, New England Biolabs, Ipswich, MA) at 1:1,000 dilution followed by antirabbit IgG (Santa Cruz Biotechnology) at 1:10,000 dilution. GLP (EuHMTaseI) was detected using primary antibody D220-3 (MBL International) at 1 μg/mL, in combination with protein A-horseradish peroxidase at 1:1,000 dilution. α-Tubulin was detected as described previously (5). Proteins were detected using Western Lightning Chemiluminescence Reagent (Perkin-Elmer Life Sciences) and the Versadoc imaging system (Bio-Rad).
Measurement of Histone Modifications
Global levels of histone H3 lysine 9 dimethylation (H3K9me2), lysine 27 dimethylation (H3K27me2), and lysine 9 acetylation (H3K9ac) were measured by Western blotting of acid-soluble nuclear protein extracts. Extracts were prepared according to the Upstate Acid Extraction protocol (Upstate Biotechnology) and dialysis was accomplished using extra-strength 10,000 MWCO Slide-A-Lyzer Dialysis Cassettes (Pierce). Extracts were fractionated by SDS-PAGE, and Coomassie staining was used to confirm equivalent histone protein loading. Primary antibodies for Western blot and ChIP analysis specific for H3K9me2, H3K27me2, H3K9ac, and H3K4me2 were obtained from Upstate Biotechnology, and secondary antibodies were obtained from Amersham Biosciences. Western blot detection was done as described above. Quantitative ChIP assays were used to detect histone code modifications at the MAGE-A1 and XAGE-1 promoters as described previously (8).
Densitometry Analysis of Western Blots
Protein densitometry was used to quantify the extent of siRNA and shRNA knockdowns and histone modification levels on Western blots. Chemiluminescent imaging of Western blots was done with VersaDoc (Bio-Rad) and signal quantification was done with Quantity One software (Bio-Rad), using background signal correction. G9a and GLP levels in total protein extracts were normalized by the level of α-tubulin in sequentially probed blots. Histone modification levels in acid-soluble protein extracts were not normalized; however, Coomassie staining was used to visually confirm equivalent protein input. All Western blots were loaded with equivalent amounts of protein per lane.
DNA Methylation Analysis of CG Antigen Genes
DNA methylation levels of human MAGE-A1, NY-ESO-1, and XAGE-1 were determined using quantitative sodium bisulfite pyrosequencing (40). NY-ESO-1 pyrosequencing primers were described previously (41). Pyrosequencing assays specific for the 5′ CpG island or CpG-rich regions of human MAGE-A1 and XAGE-1 were designed using the Pyrosequencing Assay Design Software (Biotage) and the primers are shown in Supplementary Table S1. All primers were obtained from IDT. Genomic DNA samples were isolated using the Puregene kit (Gentra Systems), and sodium bisulfite conversion and quantitative pyrosequencing were accomplished as described previously (10, 41). Of the murine Mage-a genes, only Mage-a2 (region: −179 to +54, 234 bp, 51% GC, obs/exp = 0.6) and Mage-a8 (region: −127 to +92, 220 bp, 52% GC, obs/exp = 0.6) contain 5′ CpG islands, meeting the classic definition (ref. 27; data not shown), and thus were selected for further analysis. Mage-a2 methylation was determined by methylation-specific PCR (42) using the primers shown in Supplementary Table S1. Mage-a8 methylation was determined by bisulfite pyrosequencing as described above, using the primers shown in Supplementary Table S1. Controls included DNA methylated in vitro with SssI (New England Biolabs), and unmethylated DNA was generated as described previously (43).
Global Genomic DNA Methylation Analyses
Global genomic DNA methylation was measured by determination of 5mdC levels using a quantitative liquid chromatography mass spectrometry method (44) or by sodium bisulfite pyrosequencing of the murine SINE B1 or human LINE-1 repetitive elements as described previously (28, 41). Primers were obtained from IDT.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Acknowledgments
We thank Dr. E. Michael August of Boehringer Ingelheim (Ridgefield, CT) for kindly providing BIX-01294.
References
Competing Interests
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.