Methylation-dependent SOX9 expression mediates invasion in human melanoma cells and is a negative prognostic factor in advanced melanoma
© Cheng et al.; licensee BioMed Central. 2015
Received: 14 October 2014
Accepted: 23 January 2015
Published: 22 February 2015
Melanoma is the most fatal skin cancer displaying a high degree of molecular heterogeneity. Phenotype switching is a mechanism that contributes to melanoma heterogeneity by altering transcription profiles for the transition between states of proliferation/differentiation and invasion/stemness. As phenotype switching is reversible, epigenetic mechanisms, like DNA methylation, could contribute to the changes in gene expression.
Integrative analysis of methylation and gene expression datasets of five proliferative and five invasion melanoma cell cultures reveal two distinct clusters. SOX9 is methylated and lowly expressed in the highly proliferative group. SOX9 overexpression results in decreased proliferation but increased invasion in vitro. In a B16 mouse model, sox9 overexpression increases the number of lung metastases. Transcriptional analysis of SOX9-overexpressing melanoma cells reveals enrichment in epithelial to mesenchymal transition (EMT) pathways. Survival analysis of The Cancer Genome Atlas melanoma dataset shows that metastatic patients with high expression levels of SOX9 have significantly worse survival rates. Additional survival analysis on the targets of SOX9 reveals that most SOX9 downregulated genes have survival benefit for metastatic patients.
Our genome-wide DNA methylation and gene expression study of 10 early passage melanoma cell cultures reveals two phenotypically distinct groups. One of the genes regulated by DNA methylation between the two groups is SOX9. SOX9 induces melanoma cell invasion and metastasis and decreases patient survival. A number of genes downregulated by SOX9 have a negative impact on patient survival. In conclusion, SOX9 is an important gene involved in melanoma invasion and negatively impacts melanoma patient survival.
Melanoma is an aggressive skin cancer that originates from melanocytes, that is, pigment cells that reside in the basal layer of the epidermis and are derived from the neural crest during early development . It is the most life-threatening neoplasm of the skin and is considered a major health problem due to rising incidence and mortality rates [2,3]. Melanoma is a tumor with a high degree of heterogeneity and this phenotypic heterogeneity is reversible [4-7]. In addition to being a challenge for basic research, melanoma plasticity is a major hurdle for successful treatment . Investigating the molecular basis of phenotypic heterogeneity is crucial to better understand melanoma progression and should provide useful insights for the development of more effective therapies.
In an effort to elucidate the molecular mechanisms of melanoma progression, significant differences have been detected between melanoma cells from the same lesion [4,6,9]. We and others have found that melanoma cells generally express two distinct gene expression signatures, that these signatures correlate with in vitro characteristics and these phenotypes are reversible depending on their cellular microenvironments [10-12]. One signature is characterized by the upregulation of several melanocytic genes like MITF, TYR, and DCT. These melanoma cells are highly proliferative and weakly invasive in vitro and so are named the proliferative phenotype. The other signature is characterized by the upregulation of many mesenchymal genes such as WNT5A, TGFβ, and FGF2. In contrast to the proliferative cells, these cells are highly invasive but have a low proliferative capacity in vitro and are thus named the invasive phenotype. Meta-analysis of all available melanoma microarray datasets available on the NCBI GEO database confirmed these two gene signatures in 86% of the 536 melanomas . Immunohistochemical analyses of MITF and WNT5A, markers of the proliferative and invasive phenotype, respectively, of human primary and metastatic melanomas displayed an anti-correlative staining pattern confirming that these phenotypes exist in vivo . Together these findings culminated in the phenotype switching model for melanoma progression, in which melanoma cells respond to changing micro-environmental signals, such as hypoxia, by reprogramming their gene expression patterns to switch between states of proliferation and invasion [9,15]. Thus, phenotype switching has important implications in melanoma progression. Invasive phenotype cells characterized by low MITF expression, have stem-like properties , including the ability to initiate tumors with high efficiency . Consequently, tumors comprise a mix of MITF positive and negative melanoma cells .
DNA methylation provides a stable and heritable gene regulatory mechanism for which melanoma cells could alter the expression of many genes . Aberrant DNA methylation is a mechanism known to cause tumorigenesis . Tumor suppressor genes become silenced by hypermethylation of their promoter region, thus promoting tumorigenesis. Global hypomethylation has been observed in many cancers, including melanoma, to decrease with progression of the disease [21-23]. DNMT3a and DNMT3b, the de novo DNA methyltransferases, were shown to have increased expression in metastatic melanomas compared to primary melanomas . Another group showed that DNMT3a is required for melanoma development and metastasis in a melanoma mouse model . Several signaling pathways have been shown to be deregulated as a result of aberrant DNMT-dependent methylation in melanoma, which include MAPK, WNT, PI3K, pRB, and pathways in cell cycle, apoptosis, invasion, and metastasis . Progressive global DNA hypomethylation has been observed in malignant melanocyte transformation, and surprisingly transformation was blocked in the presence of 5-Aza-2-deoxycytidine (decitabine, Aza), a DNMT inhibitor . It would suggest that targeted hypomethylation is required for malignant transformation and not overall global hypomethylation caused by Aza treatment. This is supported by our observation that treating proliferative melanoma cells with Aza had no measureable effect on their invasive abilities (data not shown). 5-Aza-2-deoxycytidine treatment of various melanoma cell lines was shown to increase SOX9 expression and induce expression of p27 and p21 . SOX9 is a transcription factor involved in neural crest specification  and SOX9 overexpression in melanoma cell lines have been shown to induce cell cycle arrest in a p21 dependent manner . Taken together, it would suggest that DNA methylation has a crucial role in malignant transformation and progression by altering the landscape of the methylome to promote tumor progression, and SOX9 is one of the targets of DNA methylation that induces cell cycle arrest.
In this study, we examine the expression of the DNMTs between the proliferative and invasive melanoma cell cultures and describe the differential melanoma methylome by MeDIP-chip. We confirm that SOX9 expression is regulated by DNA methylation and has a role in cell cycle regulation, invasion in vitro and in vivo and could be a prognostic marker for overall survival in metastatic melanoma patients.
Proliferative melanoma cells have higher levels of global DNA-methylation
A 73-gene signature is significantly differentially methylated and expressed in proliferative melanoma cells
The difference in global methylation levels and protein expression of DNMT1 and DNMT3b prompted us to investigate the methylation profiles of the proliferative and invasive phenotype melanoma cells. To determine which CpG islands were differentially methylated between the proliferative and invasive phenotypes, we immunoprecipitated methylated DNA from five proliferative and five invasive melanoma cell cultures by MeDIP  followed by hybridization to Nimblegen Human DNA Methylation 3x720K CpG Island Plus RefSeq Promoter Arrays. This array contains 720,000 probes for 22,532 promoter regions and 27,728 CpG islands. We calculated the differential methylation levels between the five proliferative and five invasive melanoma cells with a sliding window ANOVA test with the R package DMR supplied from Nimblegen. We found 406 gene promoters to be significantly and differentially methylated between the proliferative and invasive phenotypes (Figure 1C). A total of 320 promoter regions were hypermethylated in the proliferative phenotype and 86 promoter regions were hypermethylated in the invasive phenotype. The greater number of hypermethylated regions in the proliferative phenotype would be consistent with the global methylation data.
Gene expression data for the 10 melanoma cell cultures previously generated by us  were reanalyzed for differential gene expression between the proliferative and invasive phenotypes using the R package limma . A total of 1,750 genes were differentially expressed between the proliferative and invasive phenotype (fold change >2, FDR corrected P <0.05) (Figure 1D). We then analyzed the relationship between the promoter methylation status and mRNA expression levels for all genes in both datasets. Genes were filtered for a peak score >2 for methylation, fold change >2 for gene expression and an FDR-corrected P value <0.05. A total of 73 genes showed both significant differential DNA methylation and significant differential expression between the proliferative and invasive phenotype (Additional file 2: Table S1). Sixty-two genes from the proliferative phenotype had hypermethylated promoters and low RNA expression and 11 genes in the invasive phenotype had hypermethylated promoters and low RNA expression as compared to the proliferative phenotype. This suggests that methylation has a role in regulating a portion of the genes differentially expressed between the proliferative and invasive phenotype. We hypothesized that the 73 genes with both differential DNA methylation and mRNA expression between the proliferative and invasive melanoma cells were likely to be true targets of epigenetic regulation in melanoma. To determine which groups of genes were functionally important, we performed pathway analysis of the 73 genes on MetaCore. We looked for enrichment of pathways under GO processes, process networks, and Pathway maps (Additional file 3: Table S2). Interestingly, we observed significant enrichment in pathways involved in EMT, melanoma, and cell differentiation. We decided to focus on SOX9 for validation due to its known function in melanocyte differentiation and melanoma progression [29,33,34].
Sox9 expression is silenced in proliferative melanoma cells through promoter DNA methylation
To confirm that promoter DNA hypermethylation correlated with transcriptional silencing of SOX9, we assessed mRNA levels using real-time RT-PCR in the 10 melanoma cell cultures. SOX9 mRNA was expressed robustly in the invasive phenotype melanoma cell cultures compared to the proliferative phenotype melanoma cell cultures (P <0.05) (Figure 2B). Protein expression of SOX9 was detected in all invasive phenotype melanoma cell lysates, but little to no expression of SOX9 was seen in the proliferative phenotype melanoma cell lysates (Figure 2C). To validate that SOX9 is indeed regulated by DNA methylation, we treated the five proliferative phenotype melanoma cell cultures with 5 μM 5-Aza-2′-deoxycytidine (a DNMT inhibitor) for 72 h. Re-expression of SOX9 was detected by western blot (Figure 2D). Thus, SOX9 expression is regulated by DNA methylation between the proliferative and invasive phenotype.
SOX9 mediates proliferation and invasion in melanoma cell cultures
Consistent with this observation, we overexpressed SOX9 in proliferative phenotype melanoma cell cultures (that is, M010817 and M980513) by lentiviral transfection (Figure 3D). Overexpression of SOX9 was previously shown to drive melanoma cells into cell cycle arrest . We measured proliferation and cell cycle progression by EdU and PI staining, respectively, and observed the cells transfected with vector have 50.5% in G1 phase, 37% in S phase, and 12.5% in G2/M phase, whereas the cells overexpressing SOX9 have 64.8% in G1 phase, 17.7% in S phase, and 17.5% in G2/M phase (P <0.05) (Figure 3E and F). The invasive capacity of M010817 and M980513 were significantly increased from 0.95% to 9.0% (P <0.05) and from 0.67% to 4.8% (P <0.05) from SOX9 overexpression (Figure 3G and H). In concordance with previously published data on SOX9 in melanoma, we also see G1/G0 arrest when SOX9 is overexpressed along with increased invasion.
SOX9 induces a partial invasive phenotype in proliferative melanoma cells
In vivo function of SOX9 overexpression
Validation with TCGA melanoma dataset
Clinical parameters of metastatic melanoma patients
Regional cutaneous or subcutaneous metastasis
Multivariate cox regression on metastatic melanoma patients
SOX9 (SOX9 low = 0, SOX9 high =1)
Gender (female = 0, male =1)
T0 (used as reference)
N0 (used as reference)
M0 (used as reference)
Lymph Node (used as reference)
Regional Cutaneous or Subcutaneous Metastasis
We tested all 31 overlapping genes for survival benefit in the metastatic dataset (Additional file 8: Table S7). Twelve of the genes of this set were upregulated and 19 were downregulated when SOX9 expression was high. Most of the genes, 92% (11/12), that were upregulated when SOX9 expression is high had no influence on patient survival, only one gene 7% (1/12) was associated with improved survival (Figure 7D). This suggests that the genes upregulated by SOX9 are not direct factors for patient survival. Interestingly, 58% (11/19) of the genes downregulated when SOX9 expression is high were associated with improved survival and the other 42% (8/19) had no influence on patient survival (Figure 7D). This suggests that SOX9 represses a group of genes important for patient survival. Taken together, high SOX9 expression leads to poor survival possibly due to the downregulation of several genes that influence patient survival.
SOX9 binds to the promoter regions of its target genes
Using the SOX9 binding motif to screen for potential binding sites on the promoter regions from the 31 overlapping targets of the SOX9 microarray and the TCGA melanoma dataset, we found 19 of the 31 genes had a potential SOX9 binding site within a 3 kb region upstream from the TSS. To determine whether SOX9 directly binds to the promoter regions of these genes, we performed ChIP analysis using SOX9 antibodies on M010817-SOX9 cells and measured SOX9 occupancy at promoter regions of TMEM158, TBX3, and FYB, for which we could design specific primers for qPCR (Figure 7E). The specificity of this assay was demonstrated by the enrichment of three known SOX9 target sequences, COL2A1 intron 1 , p21 , and SOX10 ), as compared to a non-target gene (IP10). We observed a specific association of SOX9 with TMEM158, TBX3, and FYB, suggesting that TMEM158, TBX3, and FYB are direct targets of SOX9 in melanoma.
Phenotypic, genetic, and epigenetic heterogeneity is a common feature in human melanomas [4,5,14,38]. Tumor subpopulations can be transient and have been seen to switch between phenotypic states in vivo [4,6,9,12,39]. We have previously described two subpopulations in melanoma, the proliferative phenotype and the invasive phenotype, which are defined by specific gene signatures, in vitro characteristics, and response to drug treatment [11,13,30,35,40]. Briefly, proliferative phenotype melanoma cells are distinguished by a high proliferative capacity and low invasive capacity and the invasive phenotype melanoma cells are distinguished by a low proliferative capacity and high invasive capacity. In this study, we found specific DNA methylation signatures for the proliferative and invasive melanoma phenotypes. We observed the invasive phenotype melanoma cell cultures had modest decrease of 5% in global methylation compared to the proliferative phenotype melanoma cell cultures. This may be due to decreased DNMT1 protein expression in the invasive phenotype melanoma cells. Global methylation levels have been observed to decrease as a cancerous lesion progresses from a benign tumor to metastasis [22,41], and we observed our invasive phenotype melanoma cell cultures had decreased DNA methylation levels and were more invasive than the proliferative phenotype melanoma cell cultures, suggesting the invasive cell cultures have progressed further in malignancy. Differential expression of the de novo DNA methyltransferase DNMT3b was also seen between the proliferative and invasive phenotype. These data are consistent with a model in which DNMT1 and DNMT3b have phenotype specificity and contribute to transcriptional heterogeneity by altering the methylation landscape of a melanoma cell in the context of melanoma phenotype switching. Pathway analysis of the 73 gene signature from the DNA methylation and gene expression array lead to the discovery of many transcription networks involved in development. These transcription factors were found to be hypomethylated and highly expressed in the invasive phenotype, which would suggest the invasive melanoma cell cultures may revert to a dedifferentiated state.
A number of other studies that have looked at genome-wide DNA methylation in melanoma have indicated that several tumor suppressors are silenced by DNA methylation compared to normal melanocytes [42-44] and compared to benign nevi . Also, a recent study investigating 5-hydroxymethylation (5-hmC) in melanoma found a global decrease of 5-hmC was necessary for melanoma formation . The results from these studies indicate aberrant DNA methylation is an important process in melanoma development and progression. In our work, we looked at the differences in DNA methylation landscape between 10 primary melanoma cell cultures and uncovered two distinct populations, as previously demonstrated by gene expression microarray analysis from our group . Surprisingly, the targets we found to be differentially methylated between the two phenotypes do not overlap with the targets found to be differentially methylated between normal melanocytes and melanoma, and benign nevi and melanoma. We did not detect any differential methylation in validated methylation gene sets such as COL1A2, NPM2, HSPB6, DDIT4L, and MTIG from Koga et al.  or UCHL1, COL1A2, THBS1, and TNFRSF10D from Bonazzi et al. . As those studies were comparing the methylation state of normal melanocytes to melanoma and in this study we compare within melanoma phenotypes, this might indicate that a different set of pathways are activated or silenced by DNA methylation in melanoma progression compared to melanoma initiation. In either case, it is clear that epigenetic modifications such as DNA methylation play an important role in melanoma initiation as well as progression, and embryonic developmental program reactivation may be one of the critical outcomes of this modulatory activity.
In our study, a subset of our melanoma cell cultures had lower SOX9 expression due to a hypermethylated promoter and the other subset with high SOX9 expression had a hypomethylated promoter. We confirmed that SOX9 is regulated by DNA methylation by treating low SOX9 expressing cells with 5-aza-2-deoxycytidine treatment and saw re-expression of SOX9. To determine if the regulation of SOX9 by DNA methylation is a common mechanism in melanoma or just seen within our melanoma cell cultures, we interrogated the melanoma TCGA dataset for SOX9 and found that SOX9 gene expression and DNA methylation are anti-correlated at three consecutive methylation probes in 293 samples. This provides strong evidence that specific DNA methylation is the molecular mechanism that regulates SOX9 expression in melanoma. Alcazar et al. demonstrated that after decitabine treatment of A375 and B16 melanoma cells, the promoter of SOX9 becomes hypomethylated and SOX9 is re-expressed with induction of p27 and p21 for cell cycle arrest . Passeron et al. also observed that SOX9 was downregulated in some melanoma cell lines and induction of SOX9 expression in these melanoma cell lines resulted in lower proliferation due to upregulation of p21 . We also overexpressed SOX9 in low SOX9 melanoma cell cultures and observed G1/S cell arrest, which is consistent with the study from Passeron et al. Although the proliferation rate is reduced, the invasive capacity of these SOX9 overexpressing cells is increased, which phenocopies the endogenous SOX9 expressing cells. Conversely, knockdown of SOX9 in the invasive phenotype melanoma cells reduced the invasive capacity of the cells. Microarray analysis of SOX9 overexpression revealed an EMT-like transcriptional signature and had 10% overlap with invasive phenotype gene signature which supports the notion that SOX9 is a factor that contributes to the invasive phenotype. In vivo, sox9 expression in B16F1 cells increases their metastatic potential causing more tumor lung nodules in the tail vein injection assay. Taken together, SOX9 is a gene that is regulated by DNA methylation and functionally, SOX9 mediates cell cycle progression, invasion, and metastasis in melanoma.
TCGA is a great resource for clinical and next-generation sequencing data on human tumors. We took advantage of the melanoma dataset and demonstrated that SOX9 expression levels have a significant impact on survival of metastatic melanoma patients but SOX9 did not have a significant impact on survival in patients with primary melanomas. This could suggest that SOX9 is required for progression of primary melanoma into metastasis and metastatic tumors with high SOX9 are more aggressive to the patient. There were no clinical metrics that could distinguish SOX9 high or low in primary melanoma. Only T1 and N0 stage in metastatic melanomas were significant between SOX9 high and low. Survival analysis of metastatic patients with high SOX9 expression versus low SOX9 expression revealed a significant difference in the overall 10-year survival rates. Patients with high-SOX9 expressing tumors had a 2.3 times increased risk of death compared to patients with low SOX9 expressing tumors. Based on these findings and the invasive properties of high SOX9 expressing melanomas, it would suggest that SOX9 expression in melanomas could push the tumor toward more aggressive metastasis. Thus, SOX9 could potentially be a prognostic marker for metastatic melanoma.
We performed differential gene expression analysis on the RNAseq dataset where we defined the SOX9 high group as having a minimum of three-fold greater expression than the SOX9 low group. We only saw an overlap of 31 genes between both datasets; however, the overlap was significant as determined by hypergeometric distribution. The contribution of heterogeneity in the melanoma TCGA patient population would be one of the largest factors for the difference in gene signatures between our SOX9 microarray and TCGA RNAseq data. Nonetheless, the significant overlap of genes narrow down the potential targets of SOX9. To confirm that the targets of SOX9 have prognostic value for the patients, we performed survival analysis on all 31 genes. Surprisingly, 58% (11/19) of the genes downregulated by SOX9 were associated with improved survival, which strongly suggests SOX9 represses a set of genes that decrease tumor malignancy. Genes that were upregulated by SOX9 expression had little impact on patient survival which implies that SOX9 expression alone is sufficient to drive disease progression. From this list of 31 genes, 19 of them had a potential SOX9 binding site in its promoter. We could validate TMEM158, TBX3, and FYB as direct targets of SOX9 binding by chromatin immunoprecipitation. FYB is downregulated when SOX9 levels are high suggesting a repressive effect of SOX9 on this gene. FYB is required for inflammatory cytokine production  but no known link has been established with melanoma. TMEM158 is upregulated by SOX9 but no clear role has been established for the gene in melanoma. TBX3 is also upregulated by SOX9 and TBX3 is known to cause increased invasiveness in melanoma [48,49], suggesting TBX3 could be an effector gene that drives the invasive phenotype we see in SOX9 high cells and in patients.
In conclusion, we found SOX9 to be regulated by DNA methylation, and high SOX9 expression leads to poor survival in melanoma patients due to the activation of EMT-like genes and the downregulation of potential tumor suppressor genes in melanoma cells. This was confirmed in vivo, and new direct targets of SOX9 that may mediate its function in tumor progression were identified by transcriptional profiling and chromatin-immunoprecipitation. Future therapies targeting SOX9 could be beneficial for patients to prevent progression and especially when combined with therapies targeting cells of the proliferative phenotype. Further investigation would be required to determine if SOX9 would have early prognostic value for tumor malignancy.
Materials and methods
Melanoma cell cultures were established from surplus material from primary cutaneous melanoma and melanoma metastases removed by surgery . Written informed consent was approved by the local IRB (EK647 and EK800). Clinical diagnosis was confirmed by histology and immunohistochemistry. Melanoma cells were released from tissue biopsies and grown as previously described . Melanoma cell cultures were maintained in RPMI (Invitrogen, Carlsbad, CA, USA) supplemented with 5 mM glutamine, 1 mM sodium pyruvate, and 10% heat-inactivated fetal calf serum, and cultured at 37°C and 5% CO2. As RNA was extract previously from these cell cultures for gene expression array analysis, all cell cultures used for experiments in this paper were within five passages of the RNA isolation time point.
5-methylcytosine relative content analyses
Global DNA methylation level was evaluated by MethylFlash Methylated DNA Quantification Kit (Epigentek, Farmingdale, NY, USA) as per manufacturer’s instructions.
MeDIP assay and analysis
The MeDIP assay was performed as described . A monoclonal antibody to 5-Methylcytidine (BI-MECY-100, Eurogentec, Belgium) was used for immunoprecipitation. The immunoprecipated DNA and sonicated input DNA were differentially labeled with fluorescent dyes (Cy3 and Cy5, respectively) and hybridized to Human DNA Methylation 3x720K CpG Island Plus RefSeq Promoter Arrays (Roche Nimblegen, Madison, WI, USA). Acquisition and analysis was performed using Nimblescan 2.5 and R package DMR provided by Nimblegen. All data have been deposited into NCBI GEO GSE57971.
Gene expression analysis
Gene expression datasets were obtained from NCBI GEO GSE33728 , and analysis was performed by R using the limma package. P values were adjusted by FDR multiple hypothesis test correction.
Primers for bisulphite sequencing
Treatment with decitabine
Decitabine (5-Aza-2′-deoxycytidine, Sigma Chemical (Aza)) was dissolved in DMSO as a 10 mM stock solution, aliquoted, and kept at −20°C. Primary melanoma cell cultures were seeded in Petri dishes (approximately 5,000 cells/cm2) in RPMI untreated or treated with Aza (5 μM) for 72 h, with fresh drug-supplemented medium every 24 h.
mRNA expression analysis
Total RNA was isolated using Trizol according to manufacturer’s instructions (Invitrogen). In total, 1 μg aliquots of RNA were reverse transcribed with Reverse Transcription System (Promega) according to the manufacturer’s instructions. Data collection and analysis were performed by ABI Viia7 Fast Real-Time PCR Systems (Applied Biosystems). Gene expression values of averaged triplicate reactions were normalized to RPL28 expression levels. RPL28 primers are as follows: 5′-GCAATTGGTTCCGCTACAAC-3′ and 5′-TGTTCTTGCGGATCATGTGT-3′. The primers for RT-PCR were purchased from QIAGEN: SOX9 (Hs_SOX9_1_SG).
Cells were washed twice with cold phosphate-buffered saline (PBS) and lysed at 4°C in lysis buffer containing 20 mM Tris–HCl (pH 7.5), 1% Triton X-100 (Sigma-Aldrich, St Louis, MO, USA), 137 mM NaCl, 10% glycerol, and protease and phosphatase inhibitors (Roche, Basel, Switzerland). Proteins were separated by SDS-PAGE using the NuPAGE SDS-PAGE Gel System (Invitrogen) under reducing conditions. A total of 15 μg of protein was mixed with 9 μL of NuPage LDS sample buffer (4×) (Invitrogen, NP0007), 3.6 μL of NUPAGE Sample Reducing (Invitrogen, NP0009) and filled up to 36 μL with RIPA buffer. This mixture was incubated at 85°C for 10 min while shaking at 900 rpm. Samples were loaded on NuPage precast gels (Invitrogen). Membranes were probed with the following antibodies: SOX9 (GTX109661, GeneTex, Hsinchu City, Taiwan); DNMT1 (ab13537, Abcam, Cambridge, UK); DNMT3a (ab2850, Abcam, Cambridge, UK); DNMT3b (ab16049, Abcam, Cambridge, UK); GAPDH (ab9483, Abcam, Cambridge, UK);
Silencing RNA (siRNA) transfection of melanoma cells was carried out using INTERFERin transfection solution according to the manufacturer’s protocol (Polyplus-transfection, Illkirch, France). Cells were transfected with 5 nM of siRNA (Qiagen) for 72 h before RNA or protein was extracted. As control siRNA, the All-Star negative siRNA sequence (Qiagen) was used, and gene-specific siRNAs targeting siSOX9 (SI00007595, SI00007609) were obtained from Qiagen.
SOX9 lentiviral transfection
Lentiviral particles containing plasmids expressing full-length SOX9 cDNA or eGFP were transfected into melanoma cells for 48 h. Media supplemented with 4 ng/mL blasticidin was used for selection. After 1 week of selection, protein lysate was extracted and analyzed for SOX9 expression. Plasmids for eGFP and SOX9 were a kind gift from Dr. Thierry Passeron .
SOX9 transient transfection
Empty vector or vector containing murine sox9 (a kind gift from Prof. Lukas Sommer) was transfected into B16F1 cells with jetPEI (Polyplus, France) according to manufacturer’s protocol. Protein was isolated at 24 h, 72 h, 144 h, and 288 h after transfection and analyzed for sox9 expression by western blot.
Microarray and pathway analysis
Gene expression of eGFP and SOX9 transfected cells were analyzed using the Affymetrix Human Genome U133 Plus 2.0 Array at the Functional Genomics Center Zurich (FGCZ). Differential gene expression was determined by R package limma . Pathway analysis was performed using MetaCore (GeneGo Inc., New York, USA).
Proliferation and cell cycle analysis
For cell cycle analysis, the Click-iT EdU Alexa Fluor 647 Flow Cytometry Assay Kit (Invitrogen) was used. Cells were labelled with PI according to the manufacturer’s protocol and the DNA content was measured using a BD FACSCanto II flow cytometer (BD Biosciences) and BD FACSDiva software (BD Biosciences). Cell cycle analysis was performed in triplicate.
Boyden chamber invasion assay
Cells were seeded on FluoroBlok 24-multiwell Insert System (351157, BD Biosciences) and Biocoat Tumor Invasion System (354165, BD Biosciences). The invasion assay was performed as per manufacturer’s instructions. Migrated and invaded cells were labeled with Calcein AM fluorescent dye (354216, BD Biosciences) and fluorescence was measured with Tecan GENios (Tecan, Männendorf, Switzerland) using 485 nm excitation and 535 nm emission. Relative invasion was calculated as the ratio of the fluorescence of invading cells of the Biocoat Tumor Invasion System divided by the fluorescence of migrating cells of the FluoroBlok 24-multiwell Insert System. Boyden chamber assays were performed in triplicate.
Cells were seeded in 24-well microplates at a density of 2 × 104 cells, and cell growth was determined with a standard colorimetric assay measuring 3-(4,5-dimethylthiazol-2-yl)-2,5 diphenyltetrazolium bromide (MTT) (Sigma-Aldrich, St Louis, MO, USA) reactivity after 72 h. MTT assays were performed in triplicate.
Primers for ChIP
In vivo metastasis
B16F1 cells were transfected with empty vector or vector containing murine sox9. 2 × 105 cells were injected intravenously into C57BL/6 J mice, five mice per group (Harlan Laboratories). After 12 days mice were sacrificed and lungs were examined for metastasis. Statistics were performed using the Mann–Whitney U test. All animal experiments have been approved by the veterinary authorities of Canton of Zurich, Switzerland, and were performed in accordance with Swiss law.
The SKCM DNA methylation, RNA-seq, and clinical dataset were downloaded on 28 July 2014 for analysis. Normalized reads from the level 3 RNA-seq data were used for analysis. The dataset was segregated into primary tumors and metastatic tumors for analysis. Chi-squared test was performed on the clinical parameters between the SOX9 high and low groups. Differential expression was analyzed with voom from the limma package . Log rank test and Cox proportional hazard ratio were analyzed by the survival R package . DNA methylation β-values were calculated by minfi . Correlation was calculated by Spearman’s rank correlation coefficient.
This work was financially supported by the Georg und Bertha Schwyzer-Winiker Stiftung (RD), the Forschungskredit FK-13-042 (PFC), the Swiss Cancer League and the Swiss National Science Foundation (LS, MPL), the Promedica Foundation (OS), and the Verein für Hautkrebsforschung (MPL). The University Research Priority Program (URPP) in translational cancer research at the University of Zürich provided cell lines used in this work. We thank the Cancer Biology PhD Program of the University/ETH Zürich for support and advice. We would also like to thank The Cancer Genome Atlas (TCGA) for providing important clinical and genomic data. Finally, all research on human materials was done according to Swiss law and following the Declaration of Helsinki on medical research. Consenting patients contributed surplus melanoma biopsy material according to IRB approval numbers (Ek.647/800), and cell lines were generated following IRB KEK-ZH-Nr.2014-0425.
- Dupin E, Le Douarin NM. Development of melanocyte precursors from the vertebrate neural crest. Oncogene. 2003;22:3016–23.View ArticlePubMedGoogle Scholar
- Garbe C, Leiter U. Melanoma epidemiology and trends. Clin Dermatol. 2009;27:3–9.View ArticlePubMedGoogle Scholar
- MacKie RM, Hauschild A, Eggermont AM. Epidemiology of invasive cutaneous melanoma. Ann Oncol. 2009;20:vi1–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Quintana E, Shackleton M, Foster HR, Fullen DR, Sabel MS, Johnson TM, et al. Phenotypic heterogeneity among tumorigenic melanoma cells from patients that is reversible and not hierarchically organized. Cancer Cell. 2010;18:510–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Perego M, Tortoreto M, Tragni G, Mariani L, Deho P, Carbone A, et al. Heterogeneous phenotype of human melanoma cells with in vitro and in vivo features of tumor-initiating cells. J Invest Dermatol. 2010;130:1877–86.View ArticlePubMedGoogle Scholar
- Roesch A, Fukunaga-Kalabis M, Schmidt EC, Zabierowski SE, Brafford PA, Vultur A, et al. A temporarily distinct subpopulation of slow-cycling melanoma cells is required for continuous tumor growth. Cell. 2010;141:583–94.View ArticlePubMed CentralPubMedGoogle Scholar
- Landsberg J, Kohlmeyer J, Renn M, Bald T, Rogava M, Cron M, et al. Melanomas resist T-cell therapy through inflammation-induced reversible dedifferentiation. Nature. 2012;490:412–6.View ArticlePubMedGoogle Scholar
- Romano E, Schwartz GK, Chapman PB, Wolchock JD, Carvajal RD. Treatment implications of the emerging molecular classification system for melanoma. Lancet Oncol. 2011;12:913–22.View ArticlePubMedGoogle Scholar
- Hoek KS, Goding CR. Cancer stem cells versus phenotype-switching in melanoma. Pigment Cell Melanoma Res. 2010;23:746–59.View ArticlePubMedGoogle Scholar
- Freedman JA, Tyler DS, Nevins JR, Augustine CK. Use of gene expression and pathway signatures to characterize the complexity of human melanoma. Am J Pathol. 2011;178:2513–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Hoek KS, Schlegel NC, Brafford P, Sucker A, Ugurel S, Kumar R, et al. Metastatic potential of melanomas defined by specific gene expression profiles with no BRAF signature. Pigment Cell Res. 2006;19:290–302.View ArticlePubMedGoogle Scholar
- Hendrix MJ, Seftor EA, Seftor RE, Kasemeier-Kulesa J, Kulesa PM, Postovit LM. Reprogramming metastatic tumour cells with embryonic microenvironments. Nat Rev Cancer. 2007;7:246–55.View ArticlePubMedGoogle Scholar
- Widmer DS, Cheng PF, Eichhoff OM, Belloni BC, Zipser MC, Schlegel NC, et al. Systematic classification of melanoma cells by phenotype-specific gene expression mapping. Pigment Cell Melanoma Res. 2012;25:343–53.View ArticlePubMedGoogle Scholar
- Eichhoff OM, Zipser MC, Xu M, Weeraratna AT, Mihic D, Dummer R, et al. The immunohistochemistry of invasive and proliferative phenotype switching in melanoma: a case report. Melanoma Res. 2010;20:349–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Widmer DS, Hoek KS, Cheng PF, Eichhoff OM, Biedermann T, Raaijmakers MI, et al. Hypoxia contributes to melanoma heterogeneity by triggering HIF1alpha-dependent phenotype switching. J Invest Dermatol. 2013;133:2436–43.View ArticlePubMedGoogle Scholar
- Carreira S, Goodall J, Denat L, Rodriguez M, Nuciforo P, Hoek KS, et al. Mitf regulation of Dia1 controls melanoma proliferation and invasiveness. Genes Dev. 2006;20:3426–39.View ArticlePubMed CentralPubMedGoogle Scholar
- Cheli Y, Giuliano S, Botton T, Rocchi S, Hofman V, Hofman P, et al. Mitf is the key molecular switch between mouse or human melanoma initiating cells and their differentiated progeny. Oncogene. 2011;30:2307–18.View ArticlePubMedGoogle Scholar
- Goodall J, Carreira S, Denat L, Kobi D, Davidson I, Nuciforo P, et al. Brn-2 represses microphthalmia-associated transcription factor expression and marks a distinct subpopulation of microphthalmia-associated transcription factor-negative melanoma cells. Cancer Res. 2008;68:7788–94.View ArticlePubMedGoogle Scholar
- Plass C. Cancer epigenomics. Hum Mol Genet. 2002;11:2479–88.View ArticlePubMedGoogle Scholar
- Baylin SB, Jones PA. A decade of exploring the cancer epigenome - biological and translational implications. Nat Rev Cancer. 2011;11:726–34.View ArticlePubMed CentralGoogle Scholar
- Hansen KD, Timp W, Bravo HC, Sabunciyan S, Langmead B, McDonald OG, et al. Increased methylation variation in epigenetic domains across cancer types. Nat Genet. 2011;43:768–75.View ArticlePubMed CentralPubMedGoogle Scholar
- Molognoni F, Cruz AT, Meliso FM, Morais AS, Souza CF, Xander P, et al. Epigenetic reprogramming as a key contributor to melanocyte malignant transformation. Epigenetics. 2011;6:450–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Sigalotti L, Coral S, Nardi G, Spessotto A, Cortini E, Cattarossi I, et al. Promoter methylation controls the expression of MAGE2, 3 and 4 genes in human cutaneous melanoma. J Immunother. 2002;25:16–26.View ArticlePubMedGoogle Scholar
- Nguyen T, Kuo C, Nicholl MB, Sim MS, Turner RR, Morton DL, et al. Downregulation of microRNA-29c is associated with hypermethylation of tumor-related genes and disease outcome in cutaneous melanoma. Epigenetics. 2011;6:388–94.View ArticlePubMed CentralPubMedGoogle Scholar
- Deng T, Kuang Y, Wang L, Li J, Wang Z, Fei J. An essential role for DNA methyltransferase 3a in melanoma tumorigenesis. Biochem Biophys Res Commun. 2009;387:611–6.View ArticlePubMedGoogle Scholar
- Schinke C, Mo Y, Yu Y, Amiri K, Sosman J, Greally J, et al. Aberrant DNA methylation in malignant melanoma. Melanoma Res. 2010;20:253–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Alcazar O, Achberger S, Aldrich W, Hu Z, Negrotto S, Saunthararajah Y, et al. Epigenetic regulation by decitabine of melanoma differentiation in vitro and in vivo. Int J Cancer. 2012;131:18–29.View ArticlePubMed CentralPubMedGoogle Scholar
- Cheung M, Briscoe J. Neural crest development is regulated by the transcription factor Sox9. Development. 2003;130:5681–93.View ArticlePubMedGoogle Scholar
- Passeron T, Valencia JC, Namiki T, Vieira WD, Passeron H, Miyamura Y, et al. Upregulation of SOX9 inhibits the growth of human and mouse melanomas and restores their sensitivity to retinoic acid. J Clin Invest. 2009;119:954–63.PubMed CentralPubMedGoogle Scholar
- Hoek KS, Eichhoff OM, Schlegel NC, Dobbeling U, Kobert N, Schaerer L, et al. In vivo switching of human melanoma cells between proliferative and invasive states. Cancer Res. 2008;68:650–6.View ArticlePubMedGoogle Scholar
- Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, Lam WL, et al. Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005;37:853–62.View ArticlePubMedGoogle Scholar
- Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.PubMedGoogle Scholar
- Passeron T, Valencia JC, Bertolotto C, Hoashi T, Le Pape E, Takahashi K, et al. SOX9 is a key player in ultraviolet B-induced melanocyte differentiation and pigmentation. Proc Natl Acad Sci U S A. 2007;104:13984–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Cook AL, Smith AG, Smit DJ, Leonard JH, Sturm RA. Co-expression of SOX9 and SOX10 during melanocytic differentiation in vitro. Exp Cell Res. 2005;308:222–35.View ArticlePubMedGoogle Scholar
- Eichhoff OM, Weeraratna A, Zipser MC, Denat L, Widmer DS, Xu M, et al. Differential LEF1 and TCF4 expression is involved in melanoma cell phenotype switching. Pigment Cell Melanoma Res. 2011;24:631–42.View ArticlePubMedGoogle Scholar
- Bell DM, Leung KK, Wheatley SC, Ng LJ, Zhou S, Ling KW, et al. SOX9 directly regulates the type-II collagen gene. Nat Genet. 1997;16:174–8.View ArticlePubMedGoogle Scholar
- Shakhova O, Cheng P, Mishra PJ, Zingg D, Schaefer SM, Debbache J, et al. Antagonistic cross-regulation between Sox9 and Sox10 controls an anti-tumorigenic program in Melanoma. PLoS Genet. 2015;11:e1004877.View ArticlePubMed CentralPubMedGoogle Scholar
- Mihic-Probst D, Ikenberg K, Tinguely M, Schraml P, Behnke S, Seifert B, et al. Tumor cell plasticity and angiogenesis in human melanomas. PLoS One. 2012;7:e33571.View ArticlePubMed CentralPubMedGoogle Scholar
- Hendrix MJ, Seftor EA, Hess AR, Seftor RE. Molecular plasticity of human melanoma cells. Oncogene. 2003;22:3070–5.View ArticlePubMedGoogle Scholar
- Zipser MC, Eichhoff OM, Widmer DS, Schlegel NC, Schoenewolf NL, Stuart D, et al. A proliferative melanoma cell phenotype is responsive to RAF/MEK inhibition independent of BRAF mutation status. Pigment Cell Melanoma Res. 2011;24:326–33.View ArticlePubMedGoogle Scholar
- Fraga MF, Herranz M, Espada J, Ballestar E, Paz MF, Ropero S, et al. A mouse skin multistage carcinogenesis model reflects the aberrant DNA methylation patterns of human tumors. Cancer Res. 2004;64:5527–34.View ArticlePubMedGoogle Scholar
- Bonazzi VF, Nancarrow DJ, Stark MS, Moser RJ, Boyle GM, Aoude LG, et al. Cross-platform array screening identifies COL1A2, THBS1, TNFRSF10D and UCHL1 as genes frequently silenced by methylation in melanoma. PLoS One. 2011;6:e26121.View ArticlePubMed CentralPubMedGoogle Scholar
- Koga Y, Pelizzola M, Cheng E, Krauthammer M, Sznol M, Ariyan S, et al. Genome-wide screen of promoter methylation identifies novel markers in melanoma. Genome Res. 2009;19:1462–70.View ArticlePubMed CentralPubMedGoogle Scholar
- Muthusamy V, Duraisamy S, Bradbury CM, Hobbs C, Curley DP, Nelson B, et al. Epigenetic silencing of novel tumor suppressors in malignant melanoma. Cancer Res. 2006;66:11187–93.View ArticlePubMedGoogle Scholar
- Conway K, Edmiston SN, Khondker ZS, Groben PA, Zhou X, Chu H, et al. DNA-methylation profiling distinguishes malignant melanomas from benign nevi. Pigment Cell Melanoma Res. 2011;24:352–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Lian CG, Xu Y, Ceol C, Wu F, Larson A, Dresser K, et al. Loss of 5-hydroxymethylcytosine is an epigenetic hallmark of melanoma. Cell. 2012;150:1135–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Rajasekaran K, Kumar P, Schuldt KM, Peterson EJ, Vanhaesebroeck B, Dixit V, et al. Signaling by Fyn-ADAP via the Carma1-Bcl-10-MAP3K7 signalosome exclusively regulates inflammatory cytokine production in NK cells. Nat Immunol. 2013;14:1127–36.View ArticlePubMedGoogle Scholar
- Peres J, Prince S. The T-box transcription factor, TBX3, is sufficient to promote melanoma formation and invasion. Mol Cancer. 2013;12:117.View ArticlePubMed CentralPubMedGoogle Scholar
- Rodriguez M, Aladowicz E, Lanfrancone L, Goding CR. Tbx3 represses E-cadherin expression and enhances melanoma invasiveness. Cancer Res. 2008;68:7872–81.View ArticlePubMedGoogle Scholar
- Raaijmakers MIG, Widmer DS, Maudrich M, Koch T, Langer A, Flace A, et al. A new live cell biobank workflow efficiently recovers heterogeneous melanoma cells from native biopsies. Exp Dermatol. 2015.Google Scholar
- Geertsen RC, Hofbauer GF, Yue FY, Manolio S, Burg G, Dummer R. Higher frequency of selective losses of HLA-A and -B allospecificities in metastasis than in primary melanoma lesions. J Invest Dermatol. 1998;111:497–502.View ArticlePubMedGoogle Scholar
- Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T. BiQ Analyzer: visualization and quality control for DNA methylation data from bisulfite sequencing. Bioinformatics. 2005;21:4067–8.View ArticlePubMedGoogle Scholar
- Smyth GK. Limma: linear models for microarray data. New York: Springer; 2005.Google Scholar
- Santoro R. Analysis of chromatin composition of repetitive sequences: the ChIP-Chop assay. Methods Mol Biol. 2014;1094:319–28.View ArticlePubMedGoogle Scholar
- Kadaja M, Keyes BE, Lin M, Pasolli HA, Genander M, Polak L, et al. SOX9: a stem cell transcriptional regulator of secreted niche signaling factors. Genes Dev. 2014;28:328–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Mathelier A, Zhao X, Zhang AW, Parcy F, Worsley-Hunt R, Arenillas DJ, et al. JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles. Nucleic Acids Res. 2014;42:D142–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Therneau TM: A Package for Survival Analysis in S. In R package version 237–4; 2013.Google Scholar
- Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.