Genetic and epigenetic loss of microRNA-31 leads to feed-forward expression of EZH2 in melanoma
PDF | HTML | Supplementary Files | How to cite
Metrics: PDF 4049 views | HTML 3404 views | ?
Irfan A. Asangani1,2,*,Paul W. Harms1,2,3,*, Lois Dodson1, Mithil Pandhi1, Lakshmi P. Kunju2, Christopher A. Maher1,2,4, Douglas R. Fullen2,3, Timothy M. Johnson3, Thomas J. Giordano2, Nallasivam Palanisamy1,2,5,8, Arul M. Chinnaiyan1,2,5,6,7,8
1 Michigan Center for Translational Pathology
2 Department of Pathology, University of Michigan
3 Department of Dermatology, University of Michigan
4 Center for Computational Medicine and Bioinformatics
5 Comprehensive Cancer Center, University of Michigan Medical School
6 Howard Hughes Medical Institute, University of Michigan Medical School
7 Department of Urology, University of Michigan
8 Shared senior authors
* These authors contributed equally to this work
Arul M. Chinnaiyan, email:
Keywords: : microRNA-31, melanoma, tumor suppressor, EZH2, DZNep
Received: August 15, 2012, Accepted: August 28, 2012, Published: August 31, 2012
MicroRNAs (miRs) play a key role in cancer etiology by coordinately repressing numerous target genes involved in cell proliferation, migration and invasion. The genomic region in chromosome 9p21 that encompasses miR-31 is frequently deleted in solid cancers including melanoma; however the expression and functional role of miR-31 has not been previously studied in melanoma. Here, we queried the expression status and performed functional characterization of miR-31 in melanoma tissues and cell lines. We found that down-regulation of miR-31 was a common event in melanoma tumors and cell lines and was associated with genomic loss in a subset of samples. Down-regulation of miR-31 gene expression was also a result of epigenetic silencing by DNA methylation, and via EZH2-mediated histone methylation. Ectopic overexpression of miR-31 in various melanoma cell lines inhibited cell migration and invasion. miR-31 targets include oncogenic kinases such as SRC, MET, NIK (MAP3K14) and the melanoma specific oncogene RAB27a. Furthermore, miR-31 overexpression resulted in down-regulation of EZH2 and a de-repression of its target gene rap1GAP; increased expression of EZH2 was associated with melanoma progression and overall patient survival. Taken together, our study supports a tumor suppressor role for miR-31 in melanoma and identifies novel therapeutic targets.
Melanoma is the fifth most common cancer in the United States, with nearly 70,000 new diagnoses and 9000 deaths annually . Localized disease is managed by simple excision however advanced metastatic disease often requires chemotherapy, immune therapy, and/or small molecule inhibitors of BRAF, none of which are curative in most patients. Therefore, prognosis for advanced melanoma remains very poor .
Melanomas are genetically complex malignancies characterized by dysregulation of multiple signaling and tumor suppressor/oncogene pathways, including BRAF, NRAS, PTEN, and/or CDKN2A [3-5]. microRNAs (miRs) play an important role in multiple cancers including melanoma, and may act either as oncomiRs by targeting tumor suppressors, or as tumor-suppressive miRs by targeting oncogenes .
miR-31 is emerging as a complex player in a number of cancers. Evidence suggests miR-31 can acts as either an oncomiR or a tumor-suppressive miR in a tumor type-specific manner. As an oncomiR, miR-31 is overexpressed and promotes proliferation and/or invasion in carcinoma of the lung, colon, and head and neck [7-13], whereas it plays a tumor-suppressive role in regulating proliferation, invasion, and metastasis in other tumor types [13-18]. The parent gene for miR-31, hsa-mir-31, is located on 9p21.3, a chromosomal region frequently deleted in melanoma and other malignancies [18-20]. However, investigations of miR-31 expression and function in melanoma have been limited and not much is known about miR-31 target genes in melanoma.
EZH2, the catalytic subunit of Polycomb Repressive Complex 2 (PRC2) is overexpressed in various solid cancers and is associated with metastatic progression and poor prognosis [21-23]. EZH2 is a transcriptional repressor that catalyzes histone H3K27 trimethylation  and its upregulation in cancer is due in part to MEK-ERK-mediated transcriptional activation and/or by loss of miR-101 that negatively regulates EZH2 expression by binding to its 3’-UTR [25-28]. However, the expression status and mechanisms of EZH2 overexpression are also poorly understood in melanoma.
In this study, we characterized the expression status and functional role of miR-31 and validated several of its important target genes in melanoma. We found recurrent genomic loss of miR-31 and abrogated expression by epigenetic silencing through DNA methylation and EZH2-mediated histone H3K27 trimethylation. Transient overexpression of miR-31 in various melanoma cell lines inhibited cell migration and invasion, supporting a role for miR-31 as a tumor-suppressive miR that targets multiple oncogenes such as SRC, MET, NIK and RAB27a. Finally, we found that miR-31 negatively regulates the expression of EZH2 that may in turn mediate epigenetic silencing of miR-31 resulting in a mutually antagonistic feedback loop.
Recurrent deletion of miR-31 locus in melanoma
In order to better understand the genomic aberrations of melanoma at the copy number level, we performed array–comparative genomic hybridization (aCGH) analysis of twenty-six melanoma tissue samples and fifteen melanoma cell lines along with primary melanocytes as control. We focused primarily on the genomic region of chromosome 9p21 since this region containing cdkn2a and cdkn2b is known to be deleted in multiple solid and hematologic malignancies . We found deletion in cdkn2a and cdkn2b loci in 50% of the samples analyzed (21/42) that included both focal and larger homozygous and heterozygous deletion. Furthermore, a subset of samples (9/42, 21%) displayed larger deletions of chromosome 9p that includes the region encoding miR-31 (Figure 1A). Out of fifteen melanoma cell lines analyzed by aCGH, Malme3M and Sk-Mel5 displayed complete loss of miR-31 locus, whereas primary melanocytes did not harbor any deletions.
Figure 1: miR-31 displays decreased expression in melanoma, and is deleted in a subset of melanoma tumors. A, genomic aberrations at 9p21 (indicated by the red vertical line on the chromosome schematic in top panel) were examined by array comparative genomic hybridization, with the green and red dots indicating genomic location of the probes used in the microarray. The hsa-miR-31 locus is highlighted by a red box. B, genomic status of miR-31 in melanoma cell lines by qPCR of genomic DNA. The breast carcinoma cell line MDA-MB-231 intended to be a control displays homozygous deletion of miR-31. Normal melanocytes and HEK-293 genomic DNA served as control. C, miR-31 expression in normal melanocytes and melanoma cell lines as determined by qRT-PCR. D, expression of the oncogenic microRNA miR-182 in melanoma cell lines relative to melanocytes. E, miR-31 expression in dysplastic nevi tested (light blue), primary melanoma (dark blue) and metastatic melanoma (red) relative to benign melanocytes (green). Average fold miR31 expression and standard error of the mean for each group are shown. P-values are relative to benign melanocytes. (Inset) miR-31 expression in melanocytes (green) relative to other benign human tissues (blue).
Loss of lineage-specific miR-31 expression in primary and metastatic melanoma
Following the aCGH analysis in melanoma tissues and cell lines, we further validated the allelic imbalance of miR-31 by genomic PCR. As expected, no PCR-amplified product was observed in Malme3M and Sk-Mel5 cell lines that harbor homozygous deletion of miR-31 locus as determined by aCGH (Figure 1B). Furthermore, interphase fluorescence in situ hybridization (FISH) analysis of Malme3M showed a clear loss of miR-31 genomic region (Supplementary Figure 1). Analysis of two other cell lines MM603 and Sk-Mel94 revealed a one copy loss of miR-31; however, five additional melanoma cell lines retained normal copy numbers. Furthermore, we serendipitously found that the MDA-MB-231 breast cancer cell line that was shown to lack miR-31 expression  has a homozygous deletion of miR-31 locus (Figure 1B), while HEK-293 cells and melanocyte controls had normal copies. Additionally, mature miR-31 expression levels were evaluated in a panel of 15 melanoma cell lines and two primary melanocytes with MDA-MB-231 cell line serving as control. As depicted in Figure 1C, all melanoma cell lines except A375 exhibited lower levels of miR-31 compared to primary melanocytes, with Sk-Mel5, Malme3M and MDA-MB-231 showing no signal for mature miR-31. As a control, miR-182, an oncogenic miR in melanoma , was significantly overexpressed in all the melanoma cell lines compared to normal melanocytes, confirming the robustness and sensitivity of our PCR-based miR detection (Figure 1D).
miR-31 was reported as both an oncogene and a tumor suppressor gene based on its expression status in various cancers. To determine the expression level of miR-31 in melanoma tissues, we analyzed dysplastic nevus (n=2), primary melanoma (n=9) and metastatic melanoma (n=71) samples for transcript levels of mature miR-31 by PCR quantification. miR-31 expression levels in dysplastic nevus, primary and metastatic melanoma samples were significantly decreased compared to primary melanocytes (Figure 1E). Importantly, one metastatic melanoma sample harboring heterozygous loss of miR-31 (Mel-57) and three of the samples with homozygous deletion (Mel-54, 58 and 66) of miR-31 (Figure 1A) displayed remarkable concordance with miR-31 transcript levels. The residual miR-31 transcripts detected in these samples is likely due to normal stromal contamination of the samples. In addition, among various normal tissues, melanocytes displayed the highest miR-31 expression (Figure 1E inset). These data indicate that melanocytes are generally enriched with miR-31, and during the progression towards melanoma, miR-31 expression is abrogated either by genomic loss or by epigenetic silencing.
Figure 2: miR-31 expression is downregulated by epigenetic silencing in melanoma. A, qRT-PCR for miR-31 using total RNA extracted from G361 treated with 5’-aza-dC melanoma cells. IGFBP3 (green) is included as a positive control for drug treatment. B, Expression of EZH2 (blue) and miR-31 (black) upon EZH2 knockdown in D14 melanoma cells. C, expression of miR-31 after treatment with global histone methylation inhibitor deazaneplanocin A (DZNep) in various melanoma cell lines.
DNA demethylation as well as EZH2 depletion by siRNA or DZNep treatment induces miR-31 expression
A detailed query of miR-31 map from the UCSC genome browser revealed that the gene is encoded in an intron of MIR31HG (miR-31 host gene), a 105kB non-coding gene containing four exons (Supplementary Figure 2) located on chromosome 9p21.3. The transcriptional regulatory region of MIR31HG was enriched for the H3K27Ac mark (ENCODE ChIP-seq data) that is associated with transcriptionally active genes. Furthermore, this regulatory region consists of 77 CpG islands surrounding the transcription start site (TSS). ENCODE ChIP-seq tracks of NT2-D1 pluripotent human testicular embryonal carcinoma cell line were enriched for H3K27me3 transcriptional silencing mark in its regulatory region (Supplementary Figure 2). This observation suggests that miR-31 transcription could be silenced by CpG island methylation or by histone H3K27me3. Indeed, transcriptional silencing by CpG island methylation and H3K27me3 mediated by EZH2 have been reported in breast cancer and ATL leukemia [32, 33]. In order to test whether miR-31 locus is epigenetically silenced CpG island methylation of its regulatory region in melanoma, we treated G361 cells that have normal miR-31 copy number (Figure 1B) with a DNA methylation inhibitor, 5’-aza-2’deoxyCytidine (5-Aza-dC). The IGFBP3 gene, known to be hypermethylated in melanoma and other cancers [31, 34], was used as a positive control and was highly expressed upon 5-aza-dC treatment (Figure 2A, right panel). More importantly, miR-31 expression was also induced upon 5-aza-dC treatment (Figure 2A left panel) indicating that miR-31 transcription is silenced by DNA methylation in melanoma cells. Moreover, miR-31 was shown to be repressed by EZH2 through H3K27me3 in other cancer types [33, 35].
We next determined whether EZH2 regulates miR-31 expression in melanoma cells. D14 cells were transfected with two independently validated EZH2 siRNA (Figure 2B right panel) and expression levels of miR-31 were measured. The levels of miR-31 were markedly increased upon EZH2 knockdown (Figure 2B left panel). Next, treatment of five different melanoma cell lines with deazaneplanocin A (DZNep), a compound that, which depletes EZH2 and other polycomb repressive complex 2 (PRC2) components and thus attenuates H3K27me3 , led to the induction of miR-31 in all cell lines (Figure 2C). These results support the notion that miR-31 is epigenetically silenced by DNA methylation and EZH2-mediated histone methylation in melanomas.
Figure 3: Effect of miR-31 overexpression on proliferation, migration, and invasion of melanoma cells. A, qRT-PCR for miR-31 using total RNA from various melanoma cell lines after transient transfection of precursor miR-31 (pre-miR-31) or control miR. U6 snRNA serves as an endogenous control. Raw Ct values are plotted for both miR-31 and U6 snRNA. B, effect of miR-31 overexpression on proliferation in melanoma cell lines. Cells were quantitated by Coulter Counter at the indicated time points to monitor growth. C, effect of miR-31 overexpression on migration through Transwell and invasion through Matrigel in melanoma cell lines. Migrated and invaded cells were stained with crystal violet, solubilized and quantitated by measuring absorbance at 560 nm.* indicates p<0.05, Student’s t-test.
miR-31 inhibits migration/invasion of melanoma cells
We next examined the functional role of miR-31 in context of oncogenic phenotypes such as cell proliferation, migration and invasion of melanoma cell lines. Transient overexpression of miR-31 in Malme3M and Sk-Mel5 cells (Figure 3A) harboring genomic deletion of miR-31 resulted in divergent phenotypes for cell proliferation; Malme3M cell growth was significantly inhibited whereas Sk-Mel5 cells were unaffected (Figure 3B). However, both cell lines displayed significant reduction in migration and invasion upon miR-31 overexpression as measured by modified Boyden chamber assays (Figure 3C). Likewise, MM603 and Sk-Mel28 cells that have at least one normal copy of miR-31 also exhibited differential effect in cell proliferation upon miR-31 overexpression, with MM603 cell growth significantly inhibited while Sk-Mel28 cells were unaffected. Again, both these cell lines showed marked inhibition in their migration and invasion potential upon miR-31 overexpression. A375 cells that express high levels of endogenous miR-31 (Figure 1C) displayed proliferation advantage in the serum-free condition upon miR-31 knockdown (Supplementary Figure 3). Overall, these data suggest that miR-31 is involved in cell migration and invasion in melanoma similar to its role in breast cancer .
Figure 4: Targets of miR-31 in melanoma. A, schematic of predicted miR-31 binding sites in 3’ UTR of putative miR-31 targets, with complementary pairs shown in black and mismatches shown in red. B, immunoblot analysis of predicted miR-31 targets SRC, NIK, RAB27a, and MET in melanoma cell lines and benign melanocytes. PIK3R3, an indirect target of miR-31, is also included. β-Actin was used as loading control. C, immunoblot analysis of predicted miR-31 target proteins SRC, NIK, RAB27A, and MET, upon miR-31 overexpression. PIK3R3 as an indirect target of miR-31 displays reduced levels. β-Actin served as a loading control. D, HEK-293 cells transfected with control antagomiR (black) or anti-miR-31 antagomiR (green) along with reporter constructs containing 3’-UTRs of the predicted miR-31 targets SRC, RAB27A, NIK, and MET. RhoA 3’UTR reporter construct was used as a positive control. E, same as (D) but with control miR (black), miR-31 (red) and reporter constructs containing 3’ UTRs of the predicted miR-31 targets SRC, RAB27A, NIK, MET, and RhoA.
miR-31 targets include melanoma oncogenes SRC, NIK, RAB27a and MET
In order to identify miR-31 target genes in melanoma, we took a multi-pronged approach utilizing, (a) microarray gene expression profiling of Malme3M cells transfected with miR-31, (b) review of full list of miR-31 targets predicted by TargetScan, PITA and miRanda, and (c) melanoma literature search (Supplementary Figure 4). We focused on MET, RAB27a, SRC, and NIK (also known as MAP3K14) as putative target genes as they contain cognate miR-31 binding sites in their 3’-UTR (Figure 4A). Moreover, NIK has recently been reported as a miR-31 target in ATL leukemia . Although the PIK3R3 gene does not harbor a miR-31 binding site in its 3’-UTR, we included it as a potential downstream indirect target of miR-31 since its expression was downregulated in our microarray profile analysis (Supplementary Figure 4C). We first performed immunoblot analysis of all four potential target protein expressions in a panel of nine melanoma cell lines and two normal primary melanocyte culture (Figure 4B). SRC protein was readily detectable in multiple melanoma cell lines and NIK was overexpressed in melanoma cell lines but was detected at low levels in normal melanocytes (Figure 4B). Likewise, MET expression was detected in many of the analyzed cells with the highest expression in SK-Mel5 cells. Interestingly, PIK3R3 was exclusively detected in melanoma cells, with low or undetectable expression in normal melanocytes (Figure 4B). Similarly RAB27a, a gene involved in the regulation of membrane trafficking and exosome formation [37, 38] was highly expressed in a number of melanoma cells.
Three cell lines, MM603, SK-Mel28 and SK-Mel5, were selected for validation since all four miR-31targets SRC, NIK, RAB27a and MET along with PIK3R3 were substantially expressed in these cells. Overexpression of miR-31 by transient transfection in these cells led to significant reduction in SRC, NIK, RAB27a and MET protein levels (Figure 4C). Moreover, PIK3R3, an indirect target of miR-31 also showed reduction in all three cell lines (Figure 4C). Depletion of EZH2 activity by DZNep treatment led to an increase in miR-31 levels (Figure 3C) accompanied by a concomitant reduction in SRC, NIK and RAB27a (Supplementary Figure 5). Reciprocal experiments with miR-31 knockdown were not possible since these cell lines already express very low endogenous miR-31 or in the case of SK-Mel5 cells, miR-31 copies are deleted. To determine whether miR-31 directly targets the 3’UTR of SRC, RAB27a, NIK and MET, we cloned the full length 3’UTR of these genes into a luciferase reporter. We used 3’UTR of RhoA as a positive control since RhoA 3’UTR was reported as a target of miR-31 . Co-transfection of the luciferase reporter 3’UTR plasmids with anti-miR-31 led to significant de-repression of luciferase activity in HEK-293 cells (Figure 4D). Reciprocally, overexpression of miR-31 led to a significant decrease in luciferase activity of these 3’UTR reporters (Figure 4E), demonstrating the inhibitory effect of miR-31 on SRC, RAB27a, NIK and MET translation by binding to their 3’UTR.
Figure 5: The histone methyltransferase EZH2 is negatively regulated by miR-31 in melanoma. A, (top panel) qRT-PCR analysis of EZH2 mRNA in melanoma cell lines relative to benign melanocytes. (bottom panel) immunoblot analysis of EZH2 and MMSET in benign melanocytes and melanoma cell lines. GAPDH was used as loading control. B, qRT-PCR analysis for EZH2 expression in miR-31 (red) or control miR transfected cells (gray). Expression was normalized with gapdh transcript levels. C, immunoblot analysis for EZH2, MMSET and Rap1GAP using cell lysates from miR-31 and control miR transfected melanoma cell lines. GAPDH served as loading control.
miR-31 downregulates EZH2 expression
Upon further analysis of Malme3M gene expression profiling data, we noted that EZH2 expression was downregulated by approximately 30% upon miR-31 overexpression (Supplementary Figure 4C). Given that EZH2 acts as a repressor of miR-31 expression, we sought to further characterize the mutually antagonistic regulatory axis between miR-31 and EZH2. We found that high miR-31 expressing primary melanocyte cultures expressed low levels of EZH2 mRNA and protein by qRT-PCR and immunoblotting, respectively. However, melanoma cell lines displayed overexpression of EZH2 mRNA and protein (Figure 5A) and notably the EZH2 mRNA and protein were highly correlated (Spearman r= 0.88, p= 0.0003, Supplementary Figure 6A). The histone methyltransferase MMSET, which is highly correlated with EZH2 in many other solid cancers (Asangani et al., manuscript in review), also demonstrated increased expression in melanoma cell lines relative to primary melanocytes (Figure 5A), and was correlated with EZH2 expression in melanoma cells (Supplementary Figure 6B). Importantly, overexpression of miR-31 resulted in downregulation of EZH2 mRNA and protein (Figures 5B and 5C), and its downstream target MMSET in all three melanoma lines tested. Furthermore, Rap1GAP, known to be repressed by EZH2 , was upregulated upon miR-31 overexpression in SK-Mel5 and MM603 cells (Figure 5C). Thus, overexpression of miR-31 resulted in downregulation of EZH2 and subsequently to the modulation of its targets MMSET and Rap1GAP; altered expression of these targets might have contributed to loss of cellular invasion and migration phenotypes as described above (Figure 3).
Figure 6: EZH2 is overexpressed in melanoma and correlates with patient survival. A, qRT-PCR for EZH2 transcript in a cohort of primary melanoma (blue), metastatic melanoma (red) and benign melanocytes (green). B, representative images depicting immunostaining pattern for EZH2 in benign nevi and melanocytic lesions. Inset magnification 40x showing nuclear positivity for EZH2. C, quantitation of EZH2 immunohistochemical staining, overall staining for EZH2 was measured by multiplying staining percentage (0-100%) by staining intensity on a numerical scale (none=1, weak=2, moderate=3, strong=4), resulting in an overall product score. D, Kaplan-Meier survival analysis showing association between strong EZH2 expression in locoregional metastases and decreased overall survival relative to metastases with moderate/low EZH2 expression (p=0.044, Mantel-Cox log-rank test).
Increased expression of EZH2 is associated with melanoma progression and survival
We next measured the expression levels of EZH2 mRNA and protein in a large cohort of melanoma samples. Quantitative RT-PCR analysis demonstrated that EZH2 mRNA was overexpressed in metastatic melanomas (n=71) relative to primary melanocyte cultures (p-value =0.002) (Figure 6A). Primary melanoma tumors (n=5) displayed a trend toward increased EZH2 mRNA expression relative to primary melanocytes that did not reach statistical significance due to small sample size (Figure 6A). EZH2 protein was determined in situ by immunohistochemistry utilizing melanoma tissue microarrays (TMA) and individual sections consisting of benign nevus (n=16), dysplastic nevus (n=13), primary melanoma (n=66) and metastatic melanoma (n=127) samples. We found the expression of EZH2 to be low or absent in benign nevi (Figure 6B and C), and EZH2 expression displayed progressive increase from benign nevi, dysplastic nevi, localized and metastatic melanoma samples (p-value =0.0015, =1.7E-16, =0.05 respectively) (Figure 6B and C). To assess whether EZH2 protein expression is associated with overall survival, Kaplan-Meier analysis was carried out for the samples with follow-up data. In univariate Kaplan-Meier analysis, high EZH2 protein expression in locoregional metastases was associated with decreased 5-year survival (p=0.044, Mantel-Cox log-rank test) (Figure 6D). Taken together, these data suggests EZH2 as a highly expressed histone methyltransferase which might be playing a role in melanoma progression leading to poor survival.
Figure 7: Model for mechanisms of interaction between miR-31 and EZH2 in melanoma cells. Targets of miR-31 activate cell signaling pathways that promote EZH2 expression and/or activity. Genomic loss or epigenetic silencing of miR-31 in melanoma leads to increased expression of miR-31 targets and EZH2 activity, which results in further suppression of miR-31 expression and subsequent cancer cell proliferation and metastasis.
In this study we have determined the expression status and regulatory mechanism of miR-31 in melanoma and performed functional characterization of its role in various oncogenic phenotypes and validated novel targets. Although miR-31 has been shown to act as either an oncomiR or tumor suppressor in various tumor types , in this study we clearly establish miR-31 as a tumor suppressor, demonstrate that it targets multiple oncogenes such as SRC, MET, NIK and RAB27a and regulates the expression of EZH2 in melanoma.
The gene encoding miR-31 is located within the intron of a long noncoding RNA located on chromosome 9p21. Deletions of chromosome 9p21 are common in melanoma and other tumors and have been shown to include the miR-31 locus in cancers including mesothelioma and urothelial carcinoma [18-20]. By array CGH we observed both heterozygous and homozygous deletions involving the miR-31 locus in 21% of melanomas. The expression of miR-31 was profoundly suppressed in most of the melanoma tissues and cell lines; the loss of expression is attributed to both genomic deletion and epigenetic silencing – a common mechanism of downregulation of tumor suppressors in cancer. Epigenetic modification of chromatin via methylation of cytosine residues in CpG islands or by post-translational modification of histones marks transcriptional availability of genes. Polycomb repressor complex 2, of which EZH2 is the catalytic component, mediates trimethylation of histone 3 at lysine K27 (H3K27) and subsequent repression of target genes  including miR-31 expression in adult T cell leukemia and prostate cancer [33, 35], and ENCODE ChipSeq data indicates an enrichment of trimethylated H3K27 at the miR-31 promoter region. In addition, miR-31 is silenced in triple-negative breast carcinoma by CpG island methylation . We observed a significant induction of miR-31 expression upon depletion of EZH2 either by RNA interference or by DZNep treatment as well as by treatment with the DNA methylation inhibitor, 5’aza-dC.
The inhibitory effect of miR-31 on cell proliferation has been reported for ATL and mesothelioma [18, 33], however cell growth of only a subset of melanoma cell lines were inhibited upon ectopic overexpression of miR-31 in our studies. Conversely, knockdown of miR-31 led to increased proliferation only in serum-free conditions. miR-31-mediated attenuation of local invasion and metastasis without affecting primary tumor growth was elegantly demonstrated in breast cancer [13, 30], and we observed a similar effect on migration and invasion upon miR-31 overexpression in all melanoma cell lines tested, supporting its role as a metastasis suppressor in melanomas.
In order to identify targets that are regulated by miR-31 to mediate tumor suppressive effects in melanoma, we employed an integrative approach utilizing gene expression microarrays, miR target prediction algorithms, and literature searches. We nominated SRC, MET, NIK and RAB27a as possible miR-31 targets in melanoma; these genes are known to play pro-tumorigenic roles in melanoma [37, 39-42]. The validation of these four melanoma oncogenes as a direct miR-31 target represents a significant finding and elucidates the molecular pathways affected by chromosome 9p21 deletion in melanoma. Moreover, NIK (NF-kB inducing kinase) which promotes cell proliferation and survival in melanoma has recently been validated as a miR-31 target in ATL .
Downregulation of EZH2 may represent an additional tumor suppressive role for miR-31 since EZH2 plays a pro-tumorigenic role in melanoma . It is unlikely that EZH2 expression is directly targeted by miR-31 as EZH2 3’UTR lacks a miR-31 binding site. However, EZH2 expression and/or function can be regulated through an alternate mechanism such as MEK-ERK-Elk-1 activation pathway [26, 27]. We observed the presence of multiple cis-regulatory elements for transcriptional factors Elk-1, Sp1, AP1 (c-Jun) and NF-kB in the proximal promoter of EZH2 (Supplementary Figure 7) that are known to be activated by miR-31 targets SRC, MET, and NIK proteins. Therefore, miR-31 may indirectly inhibit EZH2 via repression of these upstream activators. Taken together, our findings indicate that EZH2 and miR-31 engage in a mutually antagonistic regulatory circuit that can accelerate tumorigenic progression (Figure 7).
In conclusion, we found recurrent deletion and epigenetic silencing of miR-31 in melanoma. And miR-31 acts as a tumor-suppressive microRNA by exerting inhibitory effects on cell motility and invasion. Importantly, we identified SRC, MET, NIK and RAB27a to be direct targets of miR-31. Moreover, we found that miR-31 also engages in an antagonistic regulatory circuit with EZH2, and loss of miR-31 expression may contribute to EZH2 overexpression in melanoma and provides a potential therapeutic axis which can be targeted by small molecule inhibitor.
If not indicated otherwise, melanoma cell lines were obtained from the American Type Culture Collection (ATCC) and cultured in DMEM with 10% FBS. SK-Mel5, SK-Mel28, and MM603 melanoma cell lines were provided by M. Soengas, and cultured in DMEM with 10% FBS. Primary foreskin melanocyte cultures were provided by M. Verahaegen and cultured in Media 254CF (Cascade Biologics) with human melanocyte growth supplement (S-002-5, Cascade Biologics).
All human subjects studies were approved by the Institutional Review Board of the University of Michigan. All tumor tissue was procured from patients who underwent excision at the University of Michigan Health System between the years of 2007 and 2010. At time of collection, tumor tissue was flash-frozen in liquid nitrogen and stored at -80oC until RNA or DNA extraction. Formalin-fixed paraffin-embedded tissue for tissue microarray construction was obtained from archival tissue blocks. Human epidermal melanocyte (HEM) RNA was obtained from Cell Applications, or collected from cultured primary melanocytes.
array–Comparative Genomic Hybridization (aCGH)
aCGH of tumor and cell line samples was performed using gDNA on Agilent’s 244K aCGH microarrays (Human Genome CGH 244K Oligo Microarray) as previously described .
RNA and DNA isolation
For patient tumor samples, areas containing at least 70% tumor were targeted for RNA isolation using hematoxylin and eosin (H&E) stains obtained on frozen sections for each specimen. Representative areas were removed from the tissue block and homogenized in the presence of Trizol reagent (Life Technologies) and total cellular RNA was purified according to the manufacturer’s standard protocol. For cell lines, total RNA was isolated using RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. After purification, RNA quality was assessed by Agilent Bioanalyzer. DNA was isolated using the DNeasy Blood & Tissue kit (Qiagen).
Quantitative real-time PCR
All PCR reactions were performed on StepOne real-time PCR system (Applied Biosystems) using SYBR Green or Taqman detection. For analysis of relative genomic copy number, primers for pri-miR-31 (Applied Biosystems assay ID Hs03302684_pri) or the reference gene PGK1 (Roche 5046173001) were utilized in reactions with 50ng genomic DNA. Stem-loop RT-PCR analysis of miR-31 expression was performed using the hsa-miR31 microRNA Assay (Applied Biosystems assay ID 002279) as described previously . RNU6B expression for normalization was prepared and analyzed using the RNU6B microRNA Assay (Applied Biosystems assay ID 001093). Results were analyzed using the 2−DDCT-method. SYBR Green PCR primers used in this study are listed in Table S1.
RNA interference and microRNA transfections
Knockdown of EZH2 was accomplished by RNA interference using two different commercially available and well-characterized siRNA duplexes (Dharmacon, Lafayette, CO). Precursors of miR-31 (Ambion PM11465), negative control miR (Ambion 17110), anti-miR-31 (Ambion 11465) and negative control anti-miR (Ambion 17010) were purchased from Ambion (Austin, TX). Transfections were performed with OptiMEM (Invitrogen) and oligofectamine (Invitrogen) as previously described .
Antibodies and Immunoblot analyses
Antibodies used in the study are listed in Table S2 and were employed at dilutions recommended by the manufacturers. For miR-31 overexpression studies, cell lines SK-Mel5 (having miR-31 genomic deletion), SK-Mel28, and MM603 were transfected with miR-31 as described above, and protein was extracted 72 hours post-transfection. Cells were homogenized in RIPA buffer containing complete proteinase inhibitor (Roche, Indianapolis, IN, USA), and protein concentrations were determined by BCA kit. For immunoblot analysis, 200ug total protein extract was boiled in sample buffer and 10-15ug aliquots were separated by SDS-PAGE and transferred onto Polyvinylidene Difluoride membrane (GE Healthcare). The membrane was incubated for half an hour in blocking buffer [Tris-buffered saline, 0.1% Tween (TBS-T), 5% nonfat dry milk] followed by incubation overnight at 4ºC with the primary antibody. Following a wash with TBS-T, the blot was incubated with horseradish peroxidase-conjugated secondary antibody and signals were visualized by enhanced chemiluminescence system as per manufacturer’s protocol (GE Healthcare).
Cell Proliferation assay
Cells were transfected with precursor miR-31 or negative miR control. After 24 hours, cells were plated into 6-well plates at a density of 2x104 cells/well. Cells were trypsinized and cell counts were analyzed on a Coulter counter (Beckman Coulter, Fullerton, CA) at indicated time points.
Motility and basement membrane invasion assay
Cells were transfected with precursor miR31 or negative control miR. After 24 hours, cells were seeded on 8 µm inserts (BD Falcon) coated with Matrigel (for invasion assays) or uncoated (for motility assays) in a 24-well culture plates. Fetal bovine serum was added to the lower chamber as a chemoattractant. After 72 hours non-invading cells were gently removed by cotton swab. Invasive/migrated cells on the lower side of the chamber were stained with crystal violet and air-dried. Relative invasion/migration was quantitated by solubilization of crystal violet dye and measurement of absorbance at 560 nm.
5’ aza-deoxycytidine and DZNep treatment
Cells were seeded in 6-well plates and treated the following day with the indicated concentrations of 5-aza-2’-deoxycytidine (EMD Biosource) or an equivalent volume of DMSO vehicle. 5’ azacytidine was renewed daily for five days of treatment. Similarly, DZNep was added to the cells and cultured for four days. Total RNA and protein were extracted for PCR or immunoblot analysis.
Construction of 3’UTR-luciferase plasmids and reporter assays
The SRC 2.0kb full length 3′ UTR Luciferase reporter plasmid was purchased from GeneCopoeia (Rockville, MD). The full length 3’UTR of NIK, MET, RAB27a and RhoA was PCR amplified using cDNA and inserted into Topo-TA cloning vector (Invitrogen), and after sequence confirmation subcloned into pEZX-MT01 vector (GeneCopoeia). For 3’UTR luciferase reporter assays, HEK-293 cells were cultured in 24-well plates and co-transfected with the 1µg reporter construct and 50nM precursor miR or negative control miR in duplicates or triplicates using FuGENE® HD (Promega). 48 hours post-transfection luciferase activities were measured using Luc-Pair miR Luciferase Assay kit (GeneCopoeia, Rockville, MD) and normalized by Renilla luciferase activity from matched wells.
For immunohistochemical analysis three tissue microarrays (TMA) were used, containing a total of 4 benign nevi, 1 dysplastic nevus, 63 primary melanomas, and 127 metastatic melanomas, with each case represented by triplicate 0.6 mm cores. Immunohistochemistry was also performed on individual slides for 12 benign nevi, 12 dysplastic nevi, and 4 primary melanomas. Tumor content of each core or slide was verified by H&E stain. Immunohistochemistry was performed as previously described , using anti-EZH2 antibody at 1:100 dilution. Slides were scored independently by two pathologists (PH and LK) for percentage of positive cells and intensity of staining. All positive cases displayed nuclear staining. Photomicrographs were taken with a SPOT Insight Color camera (Diagnostic Instruments) on an Olympus BX41 microscope with Olympus UPlanFL 10x and 40x objectives using SPOT Basic software.
All statistical analysis was carried out using GraphPad Prism software program. For quantitative data, treatment groups were reported as mean ± SEM and compared using the unpaired Student’s t-test. For Kaplan-Meir survival analysis, EZH2 expression values were categorized into low/moderate (<300 product score) and high (>300 product score). Statistical significance was established at p ≤ 0.05 unless otherwise noted.
We thank Monique Verhaegan and Andrzej Dlugosz for providing primary melanocytes. RAB27A mouse monoclonal 4B12 was a kind gift from Chiara Recchi and Miguel Seabra (UK). We thank Kalpana Ramnarayanan, Xiaojun Jing, Xuhong Cao, Jyoti Athanikar and Karen Giles for technical assistance, critical reading, and submission of the manuscript.
A.M.C. is supported by the Doris Duke Charitable Foundation Clinical Scientist Award. A.M.C. is an American Cancer Society Research Professor and A. Alfred Taubman Scholar. N.P. is supported by Development award form Melanoma Research Alliance.
Conflicts of interest:
A.M.C. is a co-founder of Compendia Biosciences and served as a consultant to Ventana Medical Systems. None of these companies were involved in these studies or approved of its content.
1. Siegel R, Ward E, Brawley O and Jemal A. Cancer statistics, 2011: the impact of eliminating socioeconomic and racial disparities on premature cancer deaths. CA Cancer J Clin. 2011; 61(4):212-236.
2. Nikolaou VA, Stratigos AJ, Flaherty KT and Tsao H. Melanoma: new insights and new therapies. J Invest Dermatol. 2012; 132(3 Pt 2):854-863.
3. Sharpless E and Chin L. The INK4a/ARF locus and melanoma. Oncogene. 2003; 22(20):3092-3098.
4. Davies H, Bignell GR, Cox C, Stephens P, Edkins S, Clegg S, Teague J, Woffendin H, Garnett MJ, Bottomley W, Davis N, Dicks E, Ewing R, Floyd Y, Gray K, Hall S, et al. Mutations of the BRAF gene in human cancer. Nature. 2002; 417(6892):949-954.
5. Hodis E, Watson IR, Kryukov GV, Arold ST, Imielinski M, Theurillat JP, Nickerson E, Auclair D, Li L, Place C, Dicara D, Ramos AH, Lawrence MS, Cibulskis K, Sivachenko A, Voet D, et al. A landscape of driver mutations in melanoma. Cell. 2012; 150(2):251-263.
6. Perera RJ and Ray A. Epigenetic regulation of miRNA genes and their role in human melanomas. Epigenomics. 2012; 4(1):81-90.
7. Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M and Garcia-Foncillas J. Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues. Mol Cancer. 2006; 5:29.
8. Slaby O, Svoboda M, Fabian P, Smerdova T, Knoflickova D, Bednarikova M, Nenutil R and Vyzula R. Altered expression of miR-21, miR-31, miR-143 and miR-145 is related to clinicopathologic features of colorectal cancer. Oncology. 2007; 72(5-6):397-402.
9. Cottonham CL, Kaneko S and Xu L. miR-21 and miR-31 converge on TIAM1 to regulate migration and invasion of colon carcinoma cells. J Biol Chem. 2010; 285(46):35293-35302.
10. Gao W, Shen H, Liu L, Xu J and Shu Y. MiR-21 overexpression in human primary squamous cell lung carcinoma is associated with poor patient prognosis. J Cancer Res Clin Oncol. 2011; 137(4):557-566.
11. Liu X, Sempere LF, Ouyang H, Memoli VA, Andrew AS, Luo Y, Demidenko E, Korc M, Shi W, Preis M, Dragnev KH, Li H, Direnzo J, Bak M, Freemantle SJ, Kauppinen S, et al. MicroRNA-31 functions as an oncogenic microRNA in mouse and human lung cancer cells by repressing specific tumor suppressors. J Clin Invest. 2010; 120(4):1298-1309.
12. Wang CJ, Stratmann J, Zhou ZG and Sun XF. Suppression of microRNA-31 increases sensitivity to 5-FU at an early stage, and affects cell migration and invasion in HCT-116 colon cancer cells. BMC Cancer. 2010; 10:616.
13. Valastyan S and Weinberg RA. miR-31: a crucial overseer of tumor metastasis and other emerging roles. Cell Cycle. 2010; 9(11):2124-2129.
14. Guo J, Miao Y, Xiao B, Huan R, Jiang Z, Meng D and Wang Y. Differential expression of microRNA species in human gastric cancer versus non-tumorous tissues. J Gastroenterol Hepatol. 2009; 24(4):652-657.
15. Bhatnagar N, Li X, Padi SK, Zhang Q, Tang MS and Guo B. Downregulation of miR-205 and miR-31 confers resistance to chemotherapy-induced apoptosis in prostate cancer cells. Cell Death Dis. 2010; 1:e105.
16. Zhang Y, Guo J, Li D, Xiao B, Miao Y, Jiang Z and Zhuo H. Down-regulation of miR-31 expression in gastric cancer tissues and its clinical significance. Med Oncol. 2010; 27(3):685-689.
17. Schaefer A, Jung M, Mollenkopf HJ, Wagner I, Stephan C, Jentzmik F, Miller K, Lein M, Kristiansen G and Jung K. Diagnostic and prognostic implications of microRNA profiling in prostate carcinoma. Int J Cancer. 2010; 126(5):1166-1176.
18. Ivanov SV, Goparaju CM, Lopez P, Zavadil J, Toren-Haritan G, Rosenwald S, Hoshen M, Chajut A, Cohen D and Pass HI. Pro-tumorigenic effects of miR-31 loss in mesothelioma. J Biol Chem. 2010; 285(30):22809-22817.
19. Flores JF, Walker GJ, Glendening JM, Haluska FG, Castresana JS, Rubio MP, Pastorfide GC, Boyer LA, Kao WH, Bulyk ML, Barnhill RL, Hayward NK, Housman DE and Fountain JW. Loss of the p16INK4a and p15INK4b genes, as well as neighboring 9p21 markers, in sporadic melanoma. Cancer Res. 1996; 56(21):5023-5032.
20. Veerla S, Lindgren D, Kvist A, Frigyesi A, Staaf J, Persson H, Liedberg F, Chebil G, Gudjonsson S, Borg A, Mansson W, Rovira C and Hoglund M. MiRNA expression in urothelial carcinomas: important roles of miR-10a, miR-222, miR-125b, miR-7 and miR-452 for tumor stage and metastasis, and frequent homozygous losses of miR-31. Int J Cancer. 2009; 124(9):2236-2242.
21. Bachmann IM, Halvorsen OJ, Collett K, Stefansson IM, Straume O, Haukaas SA, Salvesen HB, Otte AP and Akslen LA. EZH2 expression is associated with high proliferation rate and aggressive tumor subgroups in cutaneous melanoma and cancers of the endometrium, prostate, and breast. J Clin Oncol. 2006; 24(2):268-273.
22. Kleer CG, Cao Q, Varambally S, Shen R, Ota I, Tomlins SA, Ghosh D, Sewalt RG, Otte AP, Hayes DF, Sabel MS, Livant D, Weiss SJ, Rubin MA and Chinnaiyan AM. EZH2 is a marker of aggressive breast cancer and promotes neoplastic transformation of breast epithelial cells. Proc Natl Acad Sci U S A. 2003; 100(20):11606-11611.
23. Varambally S, Dhanasekaran SM, Zhou M, Barrette TR, Kumar-Sinha C, Sanda MG, Ghosh D, Pienta KJ, Sewalt RG, Otte AP, Rubin MA and Chinnaiyan AM. The polycomb group protein EZH2 is involved in progression of prostate cancer. Nature. 2002; 419(6907):624-629.
24. Margueron R and Reinberg D. The Polycomb complex PRC2 and its mark in life. Nature. 2011; 469(7330):343-349.
25. Banerjee R, Mani RS, Russo N, Scanlon CS, Tsodikov A, Jing X, Cao Q, Palanisamy N, Metwally T, Inglehart RC, Tomlins S, Bradford C, Carey T, Wolf G, Kalyana-Sundaram S, Chinnaiyan AM, et al. The tumor suppressor gene rap1GAP is silenced by miR-101-mediated EZH2 overexpression in invasive squamous cell carcinoma. Oncogene. 2011; 30(42):4339-4349.
26. Fujii S, Fukamachi K, Tsuda H, Ito K, Ito Y and Ochiai A. RAS oncogenic signal upregulates EZH2 in pancreatic cancer. Biochem Biophys Res Commun. 2012; 417(3):1074-1079.
27. Fujii S, Tokita K, Wada N, Ito K, Yamauchi C, Ito Y and Ochiai A. MEK-ERK pathway regulates EZH2 overexpression in association with aggressive breast cancer subtypes. Oncogene. 2011; 30(39):4118-4128.
28. Varambally S, Cao Q, Mani RS, Shankar S, Wang X, Ateeq B, Laxman B, Cao X, Jing X, Ramnarayanan K, Brenner JC, Yu J, Kim JH, Han B, Tan P, Kumar-Sinha C, et al. Genomic loss of microRNA-101 leads to overexpression of histone methyltransferase EZH2 in cancer. Science. 2008; 322(5908):1695-1699.
29. Lindor NM, Jalal SM, Kumar S and Lteif AN. Multiple primary tumors associated with chromosome 9p deletion. Am J Med Genet A. 2007; 143(1):95-97.
30. Valastyan S, Reinhardt F, Benaich N, Calogrias D, Szasz AM, Wang ZC, Brock JE, Richardson AL and Weinberg RA. A pleiotropically acting microRNA, miR-31, inhibits breast cancer metastasis. Cell. 2009; 137(6):1032-1046.
31. Oy GF, Slipicevic A, Davidson B, Solberg Faye R, Maelandsmo GM and Florenes VA. Biological effects induced by insulin-like growth factor binding protein 3 (IGFBP-3) in malignant melanoma. Int J Cancer. 2010; 126(2):350-361.
32. Augoff K, McCue B, Plow EF and Sossey-Alaoui K. miR-31 and its host gene lncRNA LOC554202 are regulated by promoter hypermethylation in triple-negative breast cancer. Mol Cancer. 2012; 11:5.
33. Yamagishi M, Nakano K, Miyake A, Yamochi T, Kagami Y, Tsutsumi A, Matsuda Y, Sato-Otsubo A, Muto S, Utsunomiya A, Yamaguchi K, Uchimaru K, Ogawa S and Watanabe T. Polycomb-Mediated Loss of miR-31 Activates NIK-Dependent NF-kappa B Pathway in Adult T Cell Leukemia and Other Cancers. Cancer Cell. 2012; 21(1):121-135.
34. Tomii K, Tsukuda K, Toyooka S, Dote H, Hanafusa T, Asano H, Naitou M, Doihara H, Kisimoto T, Katayama H, Pass HI, Date H and Shimizu N. Aberrant promoter methylation of insulin-like growth factor binding protein-3 gene in human cancers. Int J Cancer. 2007; 120(3):566-573.
35. Cao Q, Mani RS, Ateeq B, Dhanasekaran SM, Asangani IA, Prensner JR, Kim JH, Brenner JC, Jing X, Cao X, Wang R, Li Y, Dahiya A, Wang L, Pandhi M, Lonigro RJ, et al. Coordinated regulation of polycomb group complexes through microRNAs in cancer. Cancer Cell. 2011; 20(2):187-199.
36. Tan J, Yang X, Zhuang L, Jiang X, Chen W, Lee PL, Karuturi RK, Tan PB, Liu ET and Yu Q. Pharmacologic disruption of Polycomb-repressive complex 2-mediated gene repression selectively induces apoptosis in cancer cells. Genes Dev. 2007; 21(9):1050-1063.
37. Akavia UD, Litvin O, Kim J, Sanchez-Garcia F, Kotliar D, Causton HC, Pochanard P, Mozes E, Garraway LA and Pe’er D. An integrated approach to uncover drivers of cancer. Cell. 2010; 143(6):1005-1017.
38. Peinado H, Aleckovic M, Lavotshkin S, Matei I, Costa-Silva B, Moreno-Bueno G, Hergueta-Redondo M, Williams C, Garcia-Santos G, Ghajar CM, Nitadori-Hoshino A, Hoffman C, Badal K, Garcia BA, Callahan MK, Yuan J, et al. Melanoma exosomes educate bone marrow progenitor cells toward a pro-metastatic phenotype through MET. Nat Med. 2012.
39. Dhawan P and Richmond A. A novel NF-kappa B-inducing kinase-MAPK signaling pathway up-regulates NF-kappa B activity in melanoma cells. J Biol Chem. 2002; 277(10):7920-7928.
40. Homsi J, Cubitt C and Daud A. The Src signaling pathway: a potential target in melanoma and other malignancies. Expert Opin Ther Targets. 2007; 11(1):91-100.
41. Otsuka T, Takayama H, Sharp R, Celli G, LaRochelle WJ, Bottaro DP, Ellmore N, Vieira W, Owens JW, Anver M and Merlino G. c-Met autocrine activation induces development of malignant melanoma and acquisition of the metastatic phenotype. Cancer Res. 1998; 58(22):5157-5167.
42. Thu YM, Su Y, Yang J, Splittgerber R, Na S, Boyd A, Mosse C, Simons C and Richmond A. NF-kappaB inducing kinase (NIK) modulates melanoma tumorigenesis by regulating expression of pro-survival factors through the beta-catenin pathway. Oncogene. 2012; 31(20):2580-2592.
43. Fan T, Jiang S, Chung N, Alikhan A, Ni C, Lee CC and Hornyak TJ. EZH2-dependent suppression of a cellular senescence phenotype in melanoma cells by inhibition of p21/CDKN1A expression. Mol Cancer Res. 2011; 9(4):418-429.
44. Lonigro RJ, Grasso CS, Robinson DR, Jing X, Wu YM, Cao X, Quist MJ, Tomlins SA, Pienta KJ and Chinnaiyan AM. Detection of somatic copy number alterations in cancer using targeted exome capture sequencing. Neoplasia. 2011; 13(11):1019-1025.
45. Asangani IA, Rasheed SA, Nikolova DA, Leupold JH, Colburn NH, Post S and Allgayer H. MicroRNA-21 (miR-21) post-transcriptionally downregulates tumor suppressor Pdcd4 and stimulates invasion, intravasation and metastasis in colorectal cancer. Oncogene. 2008; 27(15):2128-2136.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.