Report

Genomic Loss of microRNA-101 Leads to Overexpression of Histone Methyltransferase EZH2 in Cancer

See allHide authors and affiliations

Science  12 Dec 2008:
Vol. 322, Issue 5908, pp. 1695-1699
DOI: 10.1126/science.1165395

Abstract

Enhancer of zeste homolog 2 (EZH2) is a mammalian histone methyltransferase that contributes to the epigenetic silencing of target genes and regulates the survival and metastasis of cancer cells. EZH2 is overexpressed in aggressive solid tumors by mechanisms that remain unclear. Here we show that the expression and function of EZH2 in cancer cell lines are inhibited by microRNA-101 (miR-101). Analysis of human prostate tumors revealed that miR-101 expression decreases during cancer progression, paralleling an increase in EZH2 expression. One or both of the two genomic loci encoding miR-101 were somatically lost in 37.5% of clinically localized prostate cancer cells (6 of 16) and 66.7% of metastatic disease cells (22 of 33). We propose that the genomic loss of miR-101 in cancer leads to overexpression of EZH2 and concomitant dysregulation of epigenetic pathways, resulting in cancer progression.

Polycomb group proteins, including enhancer of zeste homolog 2 (EZH2), play a master regulatory role in controlling important cellular process such as maintaining stem cell pluripotency (13), cell proliferation (4, 5), early embryogenesis (6), and X chromosome inactivation (7). EZH2 functions in a multiprotein complex called polycomb repressive complex 2 (PRC2), which includes SUZ12 (suppressor of zeste 12) and EED (embryonic ectoderm development) (8, 9). The primary activity of the EZH2 protein complex is to trimethylate histone H3 lysine 27 (H3K27) at target gene promoters, leading to epigenetic silencing (10, 11). Mounting evidence suggests that EZH2 has properties consistent with those of an oncogene because overexpression promotes cell proliferation, colony formation, and increased invasion of benign cells in vitro (4, 5, 12) and induces xenograft tumor growth in vivo (13). Likewise, knockdown of EZH2 in cancer cells results in growth arrest (4, 13) as well as diminished tumor growth (10) and metastasis in vivo (14).

EZH2 was initially found to be elevated in a subset of aggressive clinically localized prostate cancers and almost all metastatic prostate cancers (4). Subsequently, EZH2 has also been found to be aberrantly overexpressed in breast cancer (12), melanoma (15), bladder cancer (16), gastric cancer (17), and other cancers (5). Thus, although EZH2 is broadly overexpressed in aggressive solid tumors and has properties of an oncogene, the genetic mechanism of EZH2 elevation in cancer is unclear.

Because microRNAs (miRNAs) have gained considerable attention as regulators of gene expression (18) and play important roles in cellular differentiation and embryonic stem cell development (19), we postulated that they may play a role in modulating EZH2 expression. To test whether miRNAs play a role in controlling EZH2 expression, we computationally nominated those that might contribute to EZH2 regulation. Because intersecting the results of multiple prediction algorithms can increase specificity at the cost of lower sensitivity (20), we integrated the results of the prediction software programs PicTar (21), TargetScan (22), miRanda (23), and miRInspector (24). Overall, only 29 miRNAs were found by any program to target EZH2, whereas only microRNA-101 (miR-101) and miR-217 were found by all four programs to be predicted to regulate EZH2 (Fig. 1A and table S1) (25). Furthermore, PicTar, miRanda, and TargetScan predicted two miR-101–binding sites within the EZH2 3′ untranslated region (3′UTR) (Fig. 1B), whereas PicTar and TargetScan predicted two miR-217 binding sites within the EZH2 3′UTR. Of the 34 miRNAs predicted to regulate EZH2 by at least one program (table S1), only miR-101 had a strong negative association with prostate cancer progression from benign to localized disease to metastasis.

Fig. 1.

miR-101 regulates EZH2 transcript and protein expression. (A) Venn diagram displaying miRNAs computationally predicted to target EZH2 by PicTar (blue), miRanda (red), TargetScan (green), and MicroInspector (orange). (B) Schematic of two predicted miR-101–binding sites in the EZH2 3′UTR. (C) miR-101 regulates EZH2 transcript expression. Quantitative RT-PCR of EZH2 in SKBr3 cells transfected with precursor miR-101 is shown. Control miR and other precursor miRNAs (miR-26a, miR-128a, and miR-217) were also used for transfection. (D) miR-101 regulates PRC2 protein expression. miR-101 down-regulates EZH2 protein as well as PRC2 members SUZ12 and EED in SKBr3 cells. Control miRs and EZH2-specific siRNA were also used for transfection. The experiment was performed three independent times and a representative result is displayed. GAPDH, glyceraldehyde-3-phosphate dehydrogenase.

To examine whether miR-101 regulates the 3′UTR of EZH2, we generated luciferase reporters encoding the normal, antisense, and mutated versions of the EZH2 3′UTR. Overexpression of miR-101, but not miR-217 or control miRNA, decreased the activity of the luciferase reporter encoding the 3′UTR of EZH2 (fig. S1). Similarly, the antisense and mutant EZH2 3′UTR activities were not inhibited by miR-101. To explore whether the 3′UTR binding by miR-101 results in down-regulation of the EZH2 transcript, we transfected SKBr3 breast cancer cells (which express high levels of endogenous EZH2) with precursors of miR-101, miR-217, and a control miRNA, as well as several other unrelated miRNAs. Quantitative reverse transcription polymerase chain reaction (RT-PCR) demonstrated a reduction in EZH2 transcript levels by miR-101 (Fig. 1C) but not miR-217 or other control miRs.

To determine whether miR-101 represses EZH2 protein expression, we performed immunoblot analysis using an EZH2-specific antibody as well as antibodies to other PRC2 members, including EED and SUZ12 (Fig. 1D). In addition to miR-101, we included other miRNAs that were predicted to regulate EZH2, including miR-217 and miR-26a. Control miR-495 was predicted by TargetScan to target the PRC1 component BMI-1. Only miR-101 and EZH2 small interfering RNA (siRNA) attenuated EZH2 protein expression. miR-101 overexpression also leads to repression of EZH2's tight binding partners in the PRC2 complex: EED and, to a lesser extent, SUZ12. These proteins are thought to form a coregulated functional complex, and altering the expression of one affects the expression of the others (5, 26, 27). In this particular case, upon further inspection of the 3′UTRs of the PRC2 components, miR-101 binding sites were found in EED (fig. S2) but not in SUZ12. Because miRNAs are known to regulate multiple target genes, and in some cases hundreds of genes (18), we used the prediction algorithm TargetScan to nominate targets of miR-101. In addition to EZH2 and EED, we tested four predicted targets of miR-101 (table S2) that have been implicated in cancer pathways, including n-Myc, c-Fos, AT-rich interactive domain 1A (also called SWI-like and ARID1A), and fibrillin 2 (FBN2). None of these putative miR-101 targets were affected by overexpression of miR-101 (Fig. 1D). To support the findings from our miR-101–overexpression experiments, we employed antagomiR technology (28) to specifically inhibit miR-101 expression in benign immortalized breast epithelial cells (fig. S3). Two independent antagomiRs to miR-101 (i and ii) induced expression of EZH2 protein in benign breast epithelial cells.

To determine whether miR-101 affects EZH2 and PRC2 function, we evaluated cellular proliferation, a property known to be regulated by EZH2 (4, 5). miR-101 overexpression in SKBr3 and DU145 cells markedly attenuated cell proliferation (Fig. 2A and fig. S4). Overexpression of EZH2 (without an endogenous 3′UTR) rescued the inhibition of cell growth by miR-101, which suggests target specificity.

Fig. 2.

The role of miR-101 in regulating cell proliferation, invasion, and tumor growth. (A) miR-101 overexpression reduces cell proliferation. A cell growth assay of SKBr3 cells treated with either precursor miR-101 or siRNA targeting EZH2 is shown. Cell growth relative to the control miRNA and control siRNA duplex was measured. Rescue experiments were performed by overexpressing EZH2 (minus its endogenous 3′UTR) in miR-101–treated cells. (B) miR-101 expression decreases cell invasion of DU145 prostate carcinoma cells. We transfected cells with miR-101, EZH2-specific siRNA, control miR, and nontargeting siRNA. miR-101 was also overexpressed in those cells that overexpressed EZH2 by andenoviral infection. All cells were subjected to a matrigel invasion assay. (C) AntagomiRs to miR-101 induce the invasiveness of benign immortalized H16N2 breast epithelial cells. Representative fields of invaded and stained cells are shown in the inset. P values were calculated between control antagomiR, antagomiR-101i, and antagomiR-101ii. (D) Overexpression of miR-101 attenuates prostate tumor growth. Overexpression of miR-101 reduces DU145 tumor growth in a mouse xenograft model. Plot of mean tumor-volume trajectories over time for the mice inoculated with (red) miR-101– and (green) vector-expressing DU145 cells. Error bars represent the SE of the mean at each time point. The inset displays the decrease of EZH2 protein levels in miR-101–expressing cell lines.

We previously showed that upon overexpression, EZH2 can induce cell invasion in matrigel-coated basement membrane invasion assays (12). Here we show that miR-101 overexpression markedly inhibits the in vitro invasive potential of DU145 prostate-cancer cells (Fig. 2B) and SKBr3 breast cancer cells (fig. S5). Similarly, stable expression of miR-101 in DU145 cells showed a reduction in EZH2 expression and reduced invasion (fig. S6, A and B). Overexpression of EZH2 rescued the inhibition that was mediated by miR-101. Another in vitro readout for tumorigenic potential, increased cell migration, was also inhibited by miR-101 (fig. S7). Because overexpression of miR-101 attenuates cancer invasion, inhibition of miR-101 should enhance this neoplastic phenotype. Two independent antagomiRs targeting miR-101 (i and ii) induced an invasive phenotype when transfected into benign immortalized breast epithelial cell lines H16N2 or HME (Fig. 2C and fig. S8).

To assess whether miR-101 inhibits anchorage-independent growth, we used a soft-agar assay. DU145 prostate cancer cells stably overexpressing miR-101 exhibited markedly reduced colony formation relative to the parental cells or vector controls (fig. S9). Furthermore, in vivo, DU145 cells expressing miR-101 grew significantly slower than the vector control xenografts (P = 0.0001) (Fig. 2D), demonstrating that miR-101 has properties consistent with that of a tumor suppressor in these particular assays.

Because EZH2 and PRC2 regulate gene expression by trimethylating H3K27, we hypothesized that miR-101 overexpression would result in decreased overall H3K27 trimethylation in cancer cells. SKBr3 breast cancer and DU145 prostate cancer cells transfected with miR-101 or EZH2 siRNA for 7 days displayed a global decrease in trimethyl H3K27 levels (fig. S10A). The effect of miR-101 on H3K27 methylation was negated by overexpression of EZH2 (fig. S10B).

To test the level of promoter occupancy of the H3K27 histone mark, we performed chromatin immunoprecipitation (ChIP) assays in cancer cells overexpressing miR-101. We found significant reduction in the trimethyl H3K27 histone mark at the promoter of known PRC2 target genes such as ADRB2, DAB2IP, CIITA, and WNT1 in miR-101–overexpressing SKBr3 cells and EZH2 siRNA–treated cells (Fig. 3A and fig. S11). To determine whether the reduced promoter occupancy by H3K27 results in concomitant reduction of gene expression, we performed quantitative RT-PCR on the PRC2 targets tested by ChIP assay. As expected, there was a significant increase in target gene expression in both miR-101– and EZH2 siRNA–treated cells (Fig. 3B). To further explore miR-101 regulation of EZH2 and its possible similarity with EZH2-specific RNA interference (RNAi), we examined whether miR-101 overexpression and EZH2 knockdown generated similar gene expression profiles. To assess this, we conducted gene-expression array analysis of SKBr3 cells transfected with either miR-101 or EZH2 siRNA duplexes. Genes that were overexpressed at the twofold threshold were significantly overlapping in both the miR-101– and EZH2 siRNA–transfected cells (P = 6.08 × 10–17) (fig. S12). Similarly, those genes that were repressed also had significant overlap (P = 3.24 × 10–27).

Fig. 3.

miR-101 regulation of the cancer epigenome through EZH2 and H3K27 trimethylation. (A) ChIP assay of the trimethyl H3K27 histone mark when miR-101 is overexpressed. Known PRC2 repression targets were examined in SKBr3 cells. ChIP was performed to test H3K27 trimethylation at the promoters of ADRB2, DAB2IP, CIITA, RUNX3, CDH1, and WNT1. GAPDH, KIAA0066, and NUP214 gene promoters served as controls. (B) Quantitative RT-PCR of EZH2 target genes was performed with SKBr3 cells transfected with miR-101. The EZH2 transcript and its known targets, including ADRB2, DAB2IP, CIITA, RUNX3, and E-cadherin (CDH1) were measured.

We next investigated whether miR-101 expression inversely correlates with EZH2 levels in human tumors. A meta-analysis of a majority of the publicly available miRNA expression data sets suggested that miR-101 is significantly underexpressed in prostate, breast, ovarian, lung, and colon cancers (table S3). Because EZH2 was initially found to be overexpressed in a subset of aggressive clinically localized prostate cancers and almost universally elevated in metastatic disease (4), we examined miR-101 in a similar context of prostate cancer progression by doing quantitative PCR analysis (Fig. 4A and fig. S13). As expected, metastatic prostate cancers expressed significantly higher levels of EZH2 as compared with those of clinically localized disease or benign adjacent prostate tissue (P < 0.0001). Consistent with a functional connection between miR-101 and EZH2, miR-101 expression was significantly decreased in metastatic prostate cancer relative to that in clinically localized disease or benign adjacent prostate tissue (P < 0.0001). miR-217, which like miR-101 was predicted to regulate EZH2, did not exhibit significant differences between metastatic disease and clinically localized prostate cancer or benign prostate tissue (P = 0.35 and 0.13, respectively).

Fig. 4.

Genomic loss of the miR-101 locus may explain overexpression of EZH2 in solid tumors. (A) miR-101 transcript levels are inversely correlated with EZH2 expression in prostate cancer progression. We performed quantitative PCR for miR-101 and miR-217 by using total RNA from benign adjacent prostate, prostate cancer (PCA), and metastatic (MET) prostate cancer tissue. EZH2 expression was analyzed from the same RNA samples. (B) Genomic PCR of miR-101–1 and miR-101–2 in prostate cancer progression. Vertical axes represent log (base 2) relative quantification values; dashed lines are shown at the deletion threshold of log2(0.7) ≈ –0.51. For clarity, points have been horizontally displaced within each sample class. (C) Heat-map representation of matched normal, tumor, and metastatic samples (from right to left) in which miR-101 transcript, EZH2 transcript, and both miR-101–1 and miR-101–2 relative copy number were assessed. miR-101 and EZH2 expression is represented by a color scale highlighting down-regulation (green), no alteration (black), and up-regulation (red) of transcripts. miR-101–1 and miR-101–2 relative quantitation (RQ) of copy number are represented as homozygous loss (<0.3; bright green), single-copy loss (<0.7; light green), no copy number change (≥0.7 and ≤1.3; black), single-copy gain (>1.3; light red), and double-copy gain (>1.7; bright red). (D) Evidence that the miR-101–1 locus is somatically lost in tumors samples relative to matched normal samples. Nine metastatic prostate cancers were chosen that had copy number loss in the miR-101–1 locus, and matched normal tissue was analyzed for comparison. Bar heights represent differences in log2(RQ) values between metastatic and matched normal tissues.

To investigate the mechanism for miR-101 transcript loss in prostate cancer progression, we performed quantitative genomic PCR for miR-101. miR-101 has two genomic loci that are on chromosome 1 (miR-101–1) and chromosome 9 (miR-101–2) (fig. S14, A and B). Based on genomic PCR, 2 of 16 clinically localized prostate cancer s and 17 of 33 metastatic prostate cancers exhibited loss of the miR-101–1 locus (Fig. 4B). Similarly, 4 of 16 clinically localized prostate cancers and 8 of 33 metastatic prostate cancers displayed loss of miR-101–2 (Fig. 4B). Figure 4C displays a heat-map representation of matched samples in which miR-101 transcript, EZH2 transcript, miR-101–1 genomic loci, and miR-101–2 genomic loci were monitored. EZH2 transcript levels were inversely associated with miR-101 transcript levels across prostate cancer progression to metastasis (P < 0.0001). EZH2 tended to be uniformly elevated in samples in which the miR-101–1 or miR-101–2 genomic loci exhibited a loss in copy number (P = 0.004, permutation test).

To formally demonstrate that genomic loss of miR-101 loci was somatic in nature, we identified nine metastatic prostate cancers that exhibited loss of miR-101–1 and obtained DNA from matched normal tissue. As expected, eight of nine cases exhibited a marked decrease in relative levels of miR-101–1 copy number in the cancer as compared with that in matched normal tissue (Fig. 4D). We also explored miR-101 genomic loss in other tumor types. Using a number of experimental platforms, we demonstrated focal loss (∼20 kB) of miR-101–1 in a subset of breast, gastric, and prostate cancers (figs. S15 to S17). Furthermore, we explored public-domain high-density array comparative genomic hybridization and single-nucleotide polymorphism array data sets and observed a genomic loss of one or both miR-101 loci in a subset of glioblastoma multiforme, lung adenocarcinoma, and acute lymphocytic leukemia (fig. S18) (25).

miR-101, by virtue of its regulation of EZH2, may have profound control over the epigenetic pathways that are active not only in cancer cells but in normal pluripotent embryonic stem cells. Overexpression of miR-101 may configure the histone code of cancer cells to that associated with a more benign cellular phenotype. Because the loss of miR-101 and concomitant elevation of EZH2 are most pronounced in metastatic cancer, we postulate that miR-101 loss may represent a progressive molecular lesion in the development of more aggressive disease. Approaches to reintroduce miR-101 into tumors may have therapeutic benefit by reverting the epigenetic program of tumor cells to a more normal state.

Supporting Online Material

www.sciencemag.org/cgi/content/full/1165395/DC1

Materials and Methods

Figs. S1 to S18

Tables S1 to S10

References

References and Notes

View Abstract

Navigate This Article