Oncotarget

Research Papers:

SIRT1 regulates Mxd1 during malignant melanoma progression

PDF |  HTML  |  Supplementary Files  |  How to cite  |  Order a Reprint

Oncotarget. 2017; 8:114540-114553. https://doi.org/10.18632/oncotarget.21457

Metrics: PDF 493 views  |   HTML 1667 views  |   ?  

Fabiana M. Meliso, Danilo Micali, Camila T. Silva, Thaís S. Sabedot, Simon G. Coetzee, Adrian Koch, Fabian B. Fahlbusch, Houtan Noushmehr, Regine Schneider-Stock and Miriam G. Jasiulionis _

Abstract

Fabiana M. Meliso1,*, Danilo Micali1,*, Camila T. Silva1, Thaís S. Sabedot2, Simon G. Coetzee2, Adrian Koch3, Fabian B. Fahlbusch4, Houtan Noushmehr2, Regine Schneider-Stock3 and Miriam G. Jasiulionis1

1Ontogeny and Epigenetics Laboratory, Pharmacology Department, Universidade Federal de São Paulo - UNIFESP, São Paulo, SP, Brazil

2Department of Genetics, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, SP, Brazil

3Experimental Tumorpathology, Institute of Pathology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany

4Department of Pediatrics and Adolescent Medicine, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany

*These authors share the first authorship

Correspondence to:

Miriam G. Jasiulionis, email: mjasiulionis@gmail.com

Keywords: sirtuin-1 (SIRT1), MYC/MAX/MXD1 network, epigenetic regulation, melanoma progression, DNMT3B

Received: July 06, 2016     Accepted: September 13, 2017     Published: October 03, 2017

ABSTRACT

In a murine melanoma model, malignant transformation promoted by a sustained stress condition was causally related to increased levels of reactive oxygen species resulting in DNA damage and massive epigenetic alterations. Since the chromatin modifier Sirtuin-1 (SIRT1) is a protein attracted to double-stranded DNA break (DSB) sites and can recruit other components of the epigenetic machinery, we aimed to define the role of SIRT1 in melanomagenesis through our melanoma model. The DNA damage marker, γH2AX was found increased in melanocytes after 24 hours of deadhesion, accompanied by increased SIRT1 expression and decreased levels of its target, H4K16ac. Moreover, SIRT1 started to be associated to DNMT3B during the stress condition, and this complex was maintained along malignant progression. Mxd1 was identified by ChIP-seq among the DNA sequences differentially associated with SIRT1 during deadhesion and was shown to be a common target of both, SIRT1 and DNMT3B. In addition, Mxd1 was found downregulated from pre-malignant melanocytes to metastatic melanoma cells. Treatment with DNMT inhibitor 5AzaCdR reversed the Mxd1 expression. Sirt1 stable silencing increased Mxd1 mRNA expression and led to down-regulation of MYC targets, such as Cdkn1a, Bcl2 and Psen2, whose upregulation is associated with human melanoma aggressiveness and poor prognosis. We demonstrated a novel role of the stress responsive protein SIRT1 in malignant transformation of melanocytes associated with deadhesion. Mxd1 was identified as a new SIRT1 target gene. SIRT1 promoted Mxd1 silencing, which led to increased activity of MYC oncogene contributing to melanoma progression.


SIRT1 regulates Mxd1 during malignant melanoma progression | Meliso | Oncotarget

INTRODUCTION

Sirtuin-1 (SIRT1) is an NAD-dependent class III histone deacetylase belonging to the protein family of sirtuins. It has usually been related to gene silencing through its deacetylase action on histones and non-histone targets. SIRT1 is well known to be responsive to stress and participates in the establishment and development of several tumors, including prostate, lung, breast, ovarian, hepatocellular carcinomas, and melanoma [1, 2, 3, 4, 5, 6]. Nevertheless, the role of SIRT1 in malignant transformation is controversially discussed and seems to be dependent on cell type and microenvironment (as reviewed by [7, 8]). It can promote tumor suppression [9, 10], but can also trigger malignant progression [3, 5, 11 12].

Carcinogenesis is associated with specific mechanisms of intracellular stress status, such as oxidative stress, and its consequences, including epigenetic modifications and DNA damage [13, 14, 15, 16]. The induction of a DNA double strand break in an exogenous E-cadherin promoter CpG island resulted in the relocation of SIRT1 to the site of damage, along with other components of the epigenetic machinery, such as DNA methyltransferase 1 (DNMT1), DNA methyltransferase 3B (DNMT3B), and Enhancer of zeste homolog 2 (EZH2) [17]. In this regard, it was reported that SIRT1 can relocate from regions non-rich in CGs to areas with high levels of CGs to increase their methylation [18]. Oberdoerffer and coworkers showed SIRT1 relocation in embryonic mesenchymal cells treated with hydrogen peroxide [19]. Here, SIRT1 dissociates from repetitive sequences and oncogenes and associates to DNA break regions, resulting in transcriptional changes similar to those observed during aging. Thus, SIRT1 relocation in response to oxidative stress might have a crucial role in the establishment of aberrant epigenetic patterns associated with the development of pathologies, including cancer.

Particularly for melanomas, it was shown that increased oxidative stress leads to several mutations, which are related to tumor initiation and progression (reviewed by [20, 21, 22]). Previously, we reported that the induction of stress in non-tumorigenic melanocyte lineage melan-a via several cycles of cell adhesion impediment resulted in epigenetic reprogramming and malignant transformation [23, 24, 25, 26]. Through this protocol, different cell lines were established, such as pre-malignant 4C melanocyte lineage, non-metastatic 4C11- and metastatic 4C11+ melanoma cell lines. In this in vitro melanocyte malignant transformation model, we demonstrated that after anchorage blockade, there was an increased level of reactive oxygen species, especially superoxide anion, which contributes to increased expression of DNMT1, DNMT3B and global content of methylated cytosine [27, 25]. These results suggested an important role of epigenetic modifications caused by environmental changes for the initiation of malignant transformation [23, 24]. Our previous studies have also shown increased SIRT1 expression during the initial phases of melanoma progression [25].

One of the major signaling pathways involved in cancer progression is the MYC/MAX/MAD network, which is related to cell growth, proliferation, differentiation, and apoptosis (reviewed by [29, 30]). The components MYC (Avian myelocytomatosis viral oncogene homolog) and MAD (or MXD, MAX dimerization protein 1) have antagonistic action, leading to transcriptional activation or silencing of several genes, respectively, by competing and interacting with MAX protein binding domain [31]. Myc acts as an oncogene and is frequently overexpressed in many types of cancer [32, 33, 34, 35]. In addition, Myc has been described as a candidate target for melanoma treatment [36]. It was shown that SIRT1 binds to the C-terminal region of MYC, which leads to MYC-MAX association and its subsequent activation [37].

Considering the fact that sustained stress of melanocytes leads to their malignant transformation and taking into account the potential role of SIRT1 in melanoma progression, this work aimed to evaluate the role of SIRT1 in our melanoma transformation model. For the first time, we show an interaction between SIRT1 and MYC/MAX/MAD network in melanoma progression, highlighting the SIRT1 role in the regulation of important pathways related to cancer.

RESULTS

SIRT1 shows increased expression and associates with DNMT3B during cellular stress caused by melanocyte adhesion blockade

In a first step, we evaluated DNA damage and SIRT1 expression in melan-a melanocyte lineage subjected to adhesion blockade for 24 hours - a condition which, when repeated, results in malignant transformation [27]. An increase in γH2AX protein expression (Figure 1A), a characteristic mark of DNA double-strand breaks, as well as in SIRT1 protein expression (Figure 1B), was observed during deadhesion. This expression pattern indicates that deadhesion could lead to DNA damage and activates stress response via SIRT1. Indeed, as shown in Figure 1C, a decrease in the acetylation of lysine 16 of histone 4 (H4K16ac), a target of SIRT1, was observed during anchorage blockade, suggesting the increased SIRT1 expression could be associated with an increase in its deacetylase activity. Moreover, SIRT1 becomes associated with DNMT3B during deadhesion (Figure 1D).

SIRT1

Figure 1: SIRT1 shows increased expression and associates with DNMT3B during stress condition. The expression of γH2AX (A), SIRT1 (B) and H4K16ac (C) were evaluated by western blot in adhered (ma) and deadhered (D24h) melan-a melanocytes. The expression of β-ACTIN or ERK was used as endogenous control. The interaction between SIRT1 and DNMT3B (D) was determined by co-immunoprecipitation in adherent (ma) and deadhered (D24h) melanocytes. Total: Total protein extract used in IP assay; IP: immunoprecipitated protein fraction.

The cellular stress caused by melanocyte adhesion blockade leads to relocation of SIRT1 to DNA sequences associated with cell survival, proliferation, and migration

The SIRT1 relocation was analyzed by ChIP-seq to determine the DNA sequences where SIRT1 is bound during deadhesion compared to adhered melanocytes. The results showed that SIRT1 binds to a greater number of DNA sequences during the stressful condition compared to normal adhesion conditions (Supplementary Table 1). Our in silico analysis showed that in normal condition, SIRT1 tends to increase its association with DNA sequences next to genes related to pathways such as immune system process, transmembrane response and extracellular matrix organization (Figure 2A and Supplementary Table 2). While in the stress situation, SIRT1 has presented higher concentration value for sequences next to genes related to pathways involved in the establishment of cancer phenotype, like those related to migration, proliferation, angiogenesis, stress response, adhesion, invasion, carbohydrate metabolism, tumor growth, inflammation, chemotaxis, cell death, and cell survival (Figure 2A2B, Supplementary Table 3). Among those pathways, our in silico analysis identified five SIRT1-associated genes (App, Ncam1, Plat, Slc9a1 and Vegfc) that were concomitantly involved in the cellular processes migration, survival and proliferation (Figure 2C), further emphasizing the role of SIRT1 in the acquisition of malignant characteristics during cellular stress. Interestingly, an intronic region of the Mxd1 gene (represented by the purple bar in the Figure 3B), an important modulator of MYC function [38], is among the sequences to which SIRT1 was found to be associated during stress condition (Figure 3A).

Under

Figure 2: Under cellular stress, SIRT1 relocates to DNA sequences containing genes related to migration, survival, and proliferation. (A) Semantic similarities of Gene Ontology (GO) terms for the genes related to normal adhesion condition of melanocytes. The scatterplot shows the representative GO clusters in two-dimensional spaces. Bubble color indicates the p-value and the size indicates the frequency of GO terms. More semantically similar GO terms are closer in the plot. The mainly pathways found in adhesion were extracellular matrix organization, transmembrane transport and immune system process. The main pathways found under stress condition (deadhesion) were regulation of cell migration, proliferation, angiogenesis, negative regulation of cell death, regulation of growth, among others. (B) Biological pathways related to genes in which SIRT1 is associated during deadhesion. (C) Venn diagram showing the number of genes related to migration, survival and proliferation in which SIRT1 has associated to during deadhesion.

Increased

Figure 3: Increased SIRT1 association to Mxd1 during stress condition. (A) Heatmap showing the intensity in which SIRT1 was found to bind to Mxd1 DNA sequence in adhered (AD) and deadhered (D24h) melan-a melanocytes. (B) Mouse gene Mxd1 sequence. The purple rectangle represents the intronic region of the Mxd1 sequence where SIRT1 was differentially associated to before and during stress condition. The oligonucleotide sequence of the peak is represented below, where the orange rectangle represents the regions over which the primers forward and reverse were designed to confirm SIRT1 and DNMT3B association to Mxd1 sequence, by Real time-qPCR. The green rectangle represents a CpG island where Mxd1 transcription start site lies on. The blue schedule represents complete Mxd1 gene, located in the chromosome 6. Figure modified from UCSC Genome browser (https://genome.ucsc.edu/).

SIRT1 regulates Mxd1 during melanoma progression

Since SIRT1 participates on stress response [18, 19], has elevated expression, increased Mxd1 association in melanocytes submitted to stress condition, and is usually involved in gene silencing, as expected for histone deacetylase, we suggest that SIRT1 could participate in Mxd1 repression along melanocyte malignant transformation. To test this hypothesis, we evaluated Mxd1 gene expression by RT-qPCR. As shown in Figure 4A, the Mxd1 expression decreased along melanoma progression, being significantly reduced in pre-malignant 4C melanocytes and both non-metastatic 4C11- and metastatic 4C11+ tumor cell lines compared with parental melan-a melanocytes. By ChIP assay, using primers targeting the same Mxd1 sequence found to be enriched in ChIP-seq (Figure 3B, primer sequences highlighted in orange), we showed that SIRT1 was more associated with Mxd1 in melanoma cell lines (non-metastatic 4C11- and metastatic 4C11+ cells) compared to non-tumorigenic ones (melan-a melanocytes and pre-malignant 4C melanocytes) (Figure 4B). To further evaluate the effect of SIRT1 on Mxd1 expression, the melanoma cells 4C11- and 4C11+ were stable silenced for SIRT1 (Figure 4C and 4D). Sirt1 silencing increased expression of Mxd1 in both cell lines compared to non-silenced cells (Figure 4E and 4F), pointing to a Mxd1 repression by SIRT1 in melanoma cells.

SIRT1

Figure 4: SIRT1 regulates Mxd1 during melanoma progression. (A) Mxd1 expression in cell lines representing different phases of melanoma progression determined by RT-qPCR. (B) Association of SIRT1 to Mxd1 gene along malignant transformation was evaluated by ChIP. SIRT1 protein expression in non-metastatic 4C11- (C) and metastatic 4C11+ (D) cell lines after stable SIRT1 silencing using shRNA analyzed by western blot. Mxd1 expression in 4C11- (E) and 4C11+ (F) melanoma cell lines silenced for SIRT1 evaluated by RT-qPCR. ma: melan-a melanocyte lineage; 4C: pre-malignant melanocyte lineage; 4C11-: non-metastatic melanoma cell line; 4C11+: metastatic melanoma cell line; WT: wild type; Non-target: non-target shRNA control; sh#: SIRT1 silenced clones. *p < 0.01, **p < 0.001, ****p < 0.00001.

Demethylating agent treatment reverses Mxd1 repression

As shown above, SIRT1 formed a new complex with DNMT3B during melanocyte deadhesion. To evaluate if this SIRT1-DNMT3B complex is maintained along melanoma progression, we performed protein co-immunoprecipitation. In fact, protein-protein interaction between SIRT1 and DNMT3B was observed both in 4C pre-malignant melanocytes and in 4C11- and 4C11+ tumor cells (Figure 5A). Interestingly, DNMT3B was found to be associated with Mxd1 gene at the same site where SIRT1 interacts, since we have used the same oligonucleotides designed for the sequence of Mxd1 peak where SIRT1 has associated during stress condition (Figure 3B, primer sequences highlighted in orange). Similar to SIRT1, DNMT3B was stronger associated with Mxd1 in the tumor cell lines (Figure 5B). Moreover, Mxd1 repression was reversed after treating tumor cell lines with the demethylating agent 5-Aza-2’deoxycytidine (Figure 5C5E). Analyses by pyrosequencing of Mxd1 promoter showed no significant difference in DNA methylation status, at least in this evaluated DNA sequence (Data not shown).

Demethylating

Figure 5: Demethylating agent treatment reverses Mxd1 repression by SIRT1. (A) The interaction between SIRT1 and DNMT3B determined by protein co-immunoprecipitation assay. (B) Association of DNMT3B to Mxd1 gene along the malignant transformation was evaluated by ChIP. Mxd1 expression after treatment for 48h with 10 μM 5AzaCdR analyzed by RT-qPCR in 4C (C), 4C11- (D) and 4C11+ (E) cell lines. ma: non-tumorigenic melan-a melanocytes, 4C: pre-malignant melanocytes, 4C11-: non-metastatic melanoma cell line, 4C11+: metastatic melanoma cell line. Total: total protein extract used in IP assay; IP: immunoprecipitated protein fraction. *p < 0.01.

Down regulation of Mxd1 favors MYC activation

It is known that MXD1 competes with MYC oncoprotein for binding to MAX (MYC associated factor X). Since Mxd1 is downregulated by SIRT1 in 4C11- and 4C11+ tumor cell lines (Figure 4E and 4F), we have postulated that MAX could become more available to MYC binding in these melanoma cells. If a larger amount of MYC binds to MAX, we would expect an increased activation of MYC target genes. Among MYC targets, there are genes involved in tumorigenic processes, such as Cdkn1 (Cyclin Dependent Kinase Inhibitor 1A), Bcl2 (B-cell lymphoma 2) and Psen2 (Presenilin-2). In our melanoma model, the expression of these genes is increased in the tumor cell lines (Figure 6A) and, as hypothesized, there was a significant reduction of Cdkn1a, Bcl2 and Psen2 expression in SIRT1-silenced 4C11+ melanoma cell line (Figure 6B6D, respectively) compared to non-silenced cell lines. Together, these results show that SIRT1 is triggering MYC target gene expression by inhibiting MXD1 function, which may contribute to cancer phenotype acquisition. In line with this finding, in silico analysis of gene interactions between Sirt1, Myc and Max revealed several genes linked to tumor incidence (Figure 6E).

Down-regulation

Figure 6: Down-regulation of SIRT1 leads to an increase of Mxd1 expression and subsequent downregulation of MYC target genes. (A) Heatmap showing increased expression of MYC targets in 4C11+ melanoma cells compared to melan-a melanocytes obtained by microarray assays (color key chart with respective values is given). Arrows indicate Cdkn1a, Bcl2 and Psen2 as some of MYC targets in 4C11+ melanoma cells. Expression of Cdkn1a (B), Bcl2 (C) and Psen2 (D) was determined by RT-qPCR in non-metastatic 4C11- and metastatic 4C11+ melanoma cell lines, control and silenced for SIRT1. (E) Illustration of the interrelation among MYC and SIRT1 target genes altered in our melanoma model by Ingenuity functional pathway analysis (IPA®) graphical database of networks of interacting genes (Ingenuity Knowledge Base, IKB®). Only genes with interactions are displayed and were further subjected to the built-in “Grow to Diseases & Functions” feature to detect genes related with “tumor incidence”. WT: wild type; non-target: non-target shRNA control; sh#: SIRT1 silenced clones. ***p < 0.0001, ****p < 0.00001.

Decreased levels of Mxd1 and increased levels of Myc-target genes in human melanomas are related to poor prognosis

From the Oncomine microarray public database analysis for human tumors [39], corroborating the results described above, Mxd1 seems to be reduced in melanoma cells and in metastatic cells, while Myc target genes are increased in cutaneous melanoma and seem to be associated with a poor prognosis (Figure 7).

Decreased

Figure 7: Decreased levels of Mxd1 and increased levels of Myc-target genes in human melanomas are related to poor prognosis. Data obtained from Oncomine database [40, 41, 42, 43] are shown across multiple independently published microarray studies as indicated. Boxes represent the interquartile range (25th–75th percentile). The bars denote the median. P <-1E4

DISCUSSION

Besides many advantages of existing mouse melanoma models, there are specific shortcomings regarding their recapitulation of natural tumor progression, from increasing proliferation to invasion and metastasis [44]. Our in vitro melanocyte malignant transformation model based on repeated cycles of anchorage blockade allows studying different stages of melanoma progression from normal melanocytes to non-metastatic and metastatic melanoma cells [27]. We do not use any carcinogenic agent or genetic engineering as stimulus, as the natural ability to survive anchorage blockade is a characteristic for disseminating tumor cells and finally leads to tumor transformation [45].

Previously, we have identified massive epigenetic alterations along tumor progression in our model [25]. Particularly, DNA methyltransferases were deregulated after anchorage blockade as a result of increased levels of superoxide anion [27, 28]. Particularly, we observed a SIRT1 up-regulation in our melanoma progression model [25]. Similarly, Ohanna and coworkers reported SIRT1 overexpression in resistant BRAFV600E-mutated melanoma cells [6]. The role of the deacetylase SIRT1 in cancer is not clearly understood and seems to be dependent on cell type and stimulus. To unravel the role of SIRT1 over-expression after anchorage blockade, we aimed to identify novel SIRT1 targets in different stages of melanoma development. For the first time, we unraveled that SIRT1 recruits DNMT3B and relocates to Mxd1 gene under stress. MXD1 is a MYC inhibitor, competing with MYC for the binding of MAX co-activator. MYC/MAX/MXD1 network has been shown to play a major role in the development of neuroblastoma and melanoma [36, 46]. A high c-MYC and SIRT1 protein co-expression has been demonstrated to be associated with malignant transformation of specific colorectal cancer subtypes [47]. Inversely, Li and coworkers reported a positive regulation of SIRT1 by the c-MYC pathway in acute myeloid leukemia [48]. Since in our model SIRT1 represses Mxd1 gene expression, we suggest that MYC is activated in melanoma progression. As MYC is closely associated with tumor aggressiveness such as self-renewal capacity [49], proliferation [6] and epithelial to mesenchymal transition [50], we hypothesize that this newly identified interaction might lead to early initiation and finally to tumor progression in our model. Furthermore, our data suggest a new mechanism explaining how SIRT1 inactivates gene expression through recruitment of DNMT, i.e. by recruitment of DNMT3B to the Mxd1 gene promoter. DNMT3B is a de novo DNA methyltransferase, that is responsible for establishing new methylation patterns and is involved in gene silencing [51, 52]. Our data indicate that SIRT1 and DNMT3B co-participate in Mxd1 epigenetic silencing in melanoma.

Our ChIP-seq data showed a SIRT1/Mxd1 interaction already after 24 hours of anchorage blockade. Surprisingly, the ChIP assay along melanoma progression stages did not confirm this interaction in pre-malignant 4C melanocytes, but SIRT1/Mxd1 interaction was remarkably evident at later stages. This genomic interaction might reflect high methylation dynamics along melanoma transformation and progression. The concept of high plasticity in methylation/demethylation in cancer has been highlighted very recently by Bell and coworkers [53] and Vizoso & Esteller [54]. Cells that are exposed to stressful conditions, like ROS induction by anchorage blockade in our model, need to adapt to these environmental alterations, which is mainly achieved by epigenetic reprogramming [55]. Indeed, SIRT1/DNMT3B complex was verified after 24h of anchorage blockade and at late stages of melanoma progression, but not at the pre-malignant 4C cells, a result that we aforementioned for SIRT1/Mxd1 interaction. Corresponding to this, treatment of cells with DNMT inhibitor led to a continuous increase in Mxd1 expression along melanoma progression, suggesting sequential acquisition of aberrant epigenetic marks. Moreover, Menssen and collaborators [56] found a positive feedback loop between SIRT1 and MYC activation in colorectal cancer. They described that SIRT1 increases transcription activity of c-MYC, whereas c-MYC induces SIRT1 deacetylase activity. As a limitation of our study, we have to take into account that ChIP-seq showed 500 bp long DNA binding region for SIRT1 at Mxd1 gene, whereas the ChIP assay verified only a smaller 100 bp DNA sequence in this region and we cannot exclude sub-sequential stronger or transient binding of SIRT1 along the different stages of tumor progression. Depending on the ROS-induced DNA damage, the SIRT1 relocation at Mxd1 sequence and stress-induced SIRT1 recruited repressive protein complexes might vary along the different steps of melanoma progression. Altogether, we give experimental evidences that in our melanocyte malignant transformation model, based on sustained anchorage blockade SIRT1, is a major player for the acquisition of an aggressive phenotype. We propose a new mechanism of MYC oncogene activation by SIRT1-dependent epigenetic silencing of Mxd1 gene.

Considering the reversibility of epigenetic changes, the participation of SIRT1 and other epigenetic components in reducing Mxd1 expression and in increasing MYC oncogenic activity might have prognostic and therapeutic potential in melanoma.

MATERIALS AND METHODS

In vitro melanoma progression model

The studies were performed using an in vitro malignant melanoma progression model, previously established by our research group [27]. In this model, the melanocyte malignant transformation occurs without using physical or chemical carcinogenic agents and without genetic manipulation. The non-tumorigenic melanocyte lineage melan-a, derived spontaneously from epidermal melanoblasts of mouse embryos [57], was submitted to a sustained stress situation, represented by cycles of successive adhesion blockade. The melan-a cells, that in normal condition grow adherent to the plate ground, were maintained in suspension for 96 hours. As expected for immortalized non-tumorigenic cell line, the vast majority of cells died by anoikis after 96 hours in suspension. After this time, the few surviving cells were plated again in adherent condition. Those cells were called 1C (subjected to one cycle of anchorage blockade). When the plate achieved 80% of cell confluence, the cells were trypsinized and replated in deadhesion condition for 96 hours and the surviving cells were plated again in adherent condition and called 2C. These steps were repeated a third and fourth time, giving rise to 3C and 4C cell lines, submitted to three and four cycles of deadhesion, respectively. Subsequently, 4C cells were subjected to a fifth cycle of anchorage blockade for 96 hours and the surviving spheroids were submitted to a limiting dilution. After this, all randomly selected clones (among them 4C11-) were shown to be tumorigenic when injected into syngeneic mice and some of them were able to form metastatic colonies in the lung when injected in the caudal vein (as 4C11+). Besides forming tumoral colonies in the lung after intravenous inoculation, metastatic foci in axillary lymph nodes, confirmed by histological analysis, were found after subcutaneous inoculation of 4C11+ cells, demonstrating its metastatic phenotype. All the lineages derived from melan-a have shown stable phenotype, since they have maintained their morphology, expression of a panel of genes and proteins, and the same behavior in functional assays (such as proliferation, migration, invasion, anoikis resisitance, etc) along time.

Cell culture

The cell lines used in this study were the non-tumorigenic melan-a melanocytes, the pre-malignant 4C melanocytes, the non-metastatic 4C11- melanoma cell line, and the metastatic 4C11+ melanoma cell line. The cells were cultured in RPMI pH 6.9 (Gibco, Carlsbad, CA), with 5% fetal bovine serum (FBS) (Gibco) and, just for melan-a, 200 nM phorbol 12-myristate 13-acetate – PMA (Amresco) was added. The PMA was added to the medium only for melan-a because it is necessary for proliferation and survival of cultured non-tumoral melanoblasts and melanocytes. Melan-a cells do not form tumors in syngeneic or athymic mice even after a long time being cultured in the presence of PMA (Bennett, Cooper, and Hart 1987). On the other hand, as described for melanoma cells, the cell lines obtained after four 96-hours cycles of adhesion blockade were capable of proliferating and avoiding senescence in the absence of PMA. All cell lines were grown at 37oC, under humidified atmosphere with 5% of CO2 until utilized in the experiments. Unless mentioned, all cell lines were cultivated in adherent conditions.

Western blot

Protein extracts, prepared with RIPA buffer (50 mM Tris HCl (pH 7,4); 50 mM NaCl; 1% NP40 10%; 0.5% sodium deoxycolate 10%; 0.1% SDS 1%) containing protease inhibitors (PMSF 1 mM, Na3VO4 1 mM, pepstatin 10 mg/ml, leupeptin 10 mg/ml, aprotinin 10 mg/ml), were separated by electrophoresis in SDS-PAGE polyacrylamide gel and transferred to a PVDF membrane. Then, each membrane was blocked for 40 minutes in a solution with 5% of non-fat milk. Next, primary antibody was added to the membranes and incubated according to manufacturer’s instructions. The membranes were washed and the secondary antibody was incubated for one hour, according to manufacturer’s instructions. Primary antibodies used were: anti-phospho-H2AX (Ser139) (Merck Millipore, 07–164), anti-SIRT1 (Abcam, #Ab12193), anti-DNMT1 (Imgenex, IMG-261A), anti-DNMT3B (Active Motif, #39207), anti-H4K16ac (Millipore, Temecula, CA, #07–329). After secondary antibodies incubation (KPL), the membranes were washed with TBS-T, detected with chemiluminescent SuperSignal (ThermoScientific) and visualized by luminescence reading device UVITEC (Cambridge, www.uvitec.co.uk).

Protein co-immunoprecipitation

The association between SIRT1 and DNMT3B proteins was verified by protein co-immunoprecipitation assay. For this purpose, Protein G – agarose beads (10007D, Life Technologies) were used and the assay was performed according to the manufacturer’s instructions. Immunoprecipitation was performed with anti-SIRT1 (Abcam, #ab12–193), anti-DNMT1 (Imgenex, IMG-261A) or anti-DNMT3B (Active Motif, #39207), following datasheet specifications. The samples were used for subsequent analyses following western blot procedure.

Chromatin immunoprecipitation (ChIP)

ChIP assay was performed with Active Motif kit (ChIP-IT® High Sensitivity, catalog no. 53040), according to the manufacturer’s instructions. The ChIP assay was performed in independent duplicates.

Chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq)

ChIP-seq assay was performed with Active Motif kit (ChIP-IT® High Sensitivity, catalog no. 53040), according to the manufacturer’s instructions. ChIP-seq assay was performed in independent duplicates. DNA from ChIP was prepared from melan-a melanocytes and from melan-a melanocytes submitted to anchorage blockade during 24 hours, as described by ChIP-IT® High Sensitivity (Active Motif – 53040) protocol. After ChIP, enrichment for known targets was verified with INPUT (total DNA extract) and IgG samples by qPCR assay before sequencing. For pPCR reactions, we have used 1/10 input dilution. The qPCR results were normalized with the percentage of input amplification (sample signal value from IP assay was divided by input sample signal value). The DNA samples immunoprecipitated with anti-SIRT1 have been sequenced on a SOLiD4 machine from Life Technologies (Life Technologies, Carlsbad, CA, USA). Libraries for ChIP-seq were prepared following protocols recommended by SOLiD4. For alignment and mapping, the LifeScope mapper was used to map the reads across reference GRCm38/mm10 mouse genome. Then, Samtool package was used to sort and to remove the duplicates reads to improve specificity [58]. The reads quality was verified by FASTQC and the good mapping quality scores were considered Phred values ≥ 20 [59]. From BAM files alignment, Peak Caller MACs (Model basic analysis of ChIP-seq) was used to find regions in the genome with a significant number of mapped reads [60]. After normalizing and eliminating background, the redundant peaks were removed. Next, we used the open source software package Diffbind (Differential binding analysis of ChIP-seq peak data) from Bioconductor to count, adjust the contrast and to analyze the DNA sequences where SIRT1 was differentially bound [61]. DSEQ2, a Bioconductor package [62], was used to identify ChIP peaks associated with SIRT1. Finally, the packages Genomic Ranges and org.Mm were used to identify the gene closest to each of those peaks, considering GRCm38/mm10 mouse as the reference genome [63].

Gene Ontology (GO) analysis

The GO analysis was performed to evaluate the mainly pathways that the closest genes to the peaks where SIRT1 has associated in normal and suspension condition were involved. For this, firstly we have used DAVID tool to obtain GO terms and p-values for each gene, according the procedure recommended [64, 65]. Then, we have used the obtained GO terms and p-values (less or equal to 0.05) for REVIGO analysis. Based on algorithms that reduce redundance, REVIGO finds a subset of representative and non-redundant GO terms as described previously [66]. The results were represented by two scatterplots, one for adhesion and other for deadhesion condition, showing the representatives clusterings remaining after the redundance reduction. More semantically similar GO terms are closer in the plot. The bubble size indicates the frequency of the GO term, while the color represents the p-value according to legend.

In silico functional pathway analysis

For in silico functional pathway analysis of pathways related to SIRT1-associated genes, as well as interaction analysis of SIRT1/MYC/MAX we used Ingenuity Pathway Analysis (IPA) software (Qiagen, Redwood City, CA, USA), as previously described [67].

5-AzaCdR and TSA treatment

The compound 5-AZA-2’deoxycytidine (5-AzaCdR) is an inhibitor of DNA methylation, while Trichostatin A (TSA) inhibits HDACs of class I, II and IV. To evaluate the effect of methylation and histone modification on gene expression, the melan-a, 4C, 4C11- and 4C11+ cell lines were treated with 5-AzaCdR and TSA, respectively, or with 5-AzaCdR and TSA concomitantly to verify if the effect of the treatment with both drugs would be amplified. Treatment was started at 40% confluence. The compound 5-AzaCdR was added to the plate at 10 μM final concentration in 10 ml of RPMI medium enriched with 5% FBS and the cells were cultivated for 24 hours in a humidified incubator at 37oC with 5% CO2. Every 24 hours, the medium was changed with the addition of 5-AzaCdR at the same conditions. For the treatment with TSA, the compound was added to the plate at the final concentration of 40 nM and the cells were maintained at the same conditions described before. For the treatment with 5-AzaCdR and TSA, the cells were primarily treated with 5-AzaCdR and, after 96 hours the medium was removed and the TSA was added, as described before. Subsequently, total RNA was extracted from the cells following TRIzol® recommended protocol.

Quantitative Real Time PCR

Quantitative analysis of mRNA expression was performed using the Corbett Rotor-Gene 6000 system detection using Rotor-Gene kit Fast Syber Green PCR Master Mix® (Qiagen, Dusseldorf, Germany).

Pyrosequencing

Pyrosequencing assay was performed with PyroMark Gold Q24 Qiagen kit (catalog number 1055272), according to the manufacturer’s instructions.

Sirt1 stable silencing

The stable silencing of Sirt1 was performed by MISSION® shRNA Lentiviral Transduction Particles from Sigma-Aldrich (SHCLNV-NM_019812). Lentiviral particles were used to establish cell lines with knockdown of Sirt1 or expressing the non-target shRNA sequence as a control. All transductions were performed using a MOI (multiplicity of infection) of 0.5, in the presence of 8 μg/mL polybrene. After transduction, puromycin-resistant cells were selected. The target silencing was checked by RT-qPCR.

Author contributions

Fabiana M. Meliso wrote the paper, performed all experiments and analyzed the experimental data. She also analyzed ChIPseq data supported by Thais S. Sabedot, Simon G. Coetzee and Houtan Noushmehr. Danilo Micali Pereira was responsible for SIRT1 silencing assay and contributed with Real Time PCR reactions and immunoprecipitation assays. Camila T. Silva assisted in primer design and pyrosequencing assay. Fabian B. Fahlbusch conducted the in silico IPA Analysis. Thaís S. Sabedot and Simon Coetzee collaborated with bioinformatics analyzes. Houtan Noushmehr supervised all ChIPseq data analysis by bioinformatics. Regine S. Stock supervised ChIP and ChIPseq assays. Miriam G. Jasiulionis designed this study and oversaw all this study and paper.

ACKNOWLEDGEMENTS

We thank Camila Ferreira de Souza for analyzing gene expression data by microarray with cells representing our melanoma progression model. We also thank Professor Alexandre de Castro Keller, who provided the equipment UVITEC for our use and the technician Ingrid Mons, who supported us in Pyrosequencing assays. We thank Dr. A. Ekici and Dr. S. Uebe, from the Institute of Human Genetics, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Germany, for performing the ChIP-Sequencing data collection.

CONFLICTS OF INTEREST

This work does not present any conflicts of interest.

GRANT SUPPORT

This work was supported by FAPESP (2011/01663-8, 2011/12306-1, 2014/13663-0, 2015/07925-5, and 2016/06488-3), DAAD (PKZ A/12/79134) and FAPESP/BAYLAT (2012/51300-7).

REFERENCES

1. Cui Y, Li J, Zheng F, Ouyang Y, Chen X, Zhang L, Chen Y, Wang L, Mu S, Zhang H. Effect of SIRT1 Gene on Epithelial-Mesenchymal Transition of Human Prostate Cancer PC-3 Cells. Med. Sci Monit. 2016; 22:380–386.

2. Mvunta DH, Miyamoto T, Asaka R, Yamada Y, Ando H, Higuchi S, Ida K, Kashima H, Shiozawa T. Overexpression of SIRT1 is Associated With Poor Outcomes in Patients With Ovarian Carcinoma. Appl Immunohistochem Mol Morphol. 2016; epub ahead of print.

3. Chung SY, Jung YY, Park IA, Kim H, Chung YR, Kim JY, Park SY, Im SA, Lee KH, Moon HG, Noh DY, Han W, Lee C, et al. Oncogenic role of SIRT1 associated with tumor invasion, lymph node metastasis, and poor disease-free survival in triple negative breast cancer. Clin. Exp. Metastasis 2016; 33:179–185.

4. Wang Y, Cui R, Zhang X, Qiao Y, Liu X, Chang Y, Yu Y, Sun F, Wang J. SIRT1 increases YAP- and MKK3-dependent p38 phosphorylation in mouse liver and human hepatocellular carcinoma. Oncotarget. 2016; 7:11284–98. https://doi.org/10.18632/oncotarget.7022 https://www.ncbi.nlm.nih.gov/pubmed/26824501.

5. Lu J, Zhang M, Huang Z, Sun S, Zhang Y, Zhang L, Peng L, Ma A, Ji P, Dai J, Cui T, Liu H, Gao J. SIRT1 in B[a]P-induced lung tumorigenesis. Oncotarget. 2015; 6:27113–29. https://doi.org/10.18632/oncotarget.4729.

6. Ohanna M, Bonet C, Bille K, Allegra M, Davidson I, Bahadoran P, Lacour JP, Ballotti R, Bertolotto C. SIRT1 promotes proliferation and inhibits the senescence-like phenotype in human melanoma cells. Oncotarget. 2014; 5:2085–95. https://doi.org/10.18632/oncotarget.1791.

7. Yang H, Bi Y, Xue L, Wang J, Lu Y, Zhang Z, Chen X, Chu Y, Yang R, Wang R, Liu G. Multifaceted Modulation of SIRT1 in Cancer and Inflammation. Crit. Rev. Oncog. 2015; 20:49–64.

8. Lin Z, Fang D. The Roles of SIRT1 in Cancer. Genes Cancer. 2013; 4:97–104. https://doi.org/10.1177/1947601912475079.

9. Kikuchi K, Noguchi A, Kasajima R, Miyagi Y, Hoshino D, Koshikawa N, Kubota A, Yokose T, Takano Y. Association of SIRT1 and tumor suppressor gene TAp63 expression in head and neck squamous cell carcinoma. Tumour Biol. 2015; 36:7865–72.

10. Di Sante G, Pestell TG, Casimiro MC, Bisetto S, Powell MJ, Lisanti MP, Cordon-Cardo C, Castillo-Martin M, Bonal DM, Debattisti V, Chen K, Wang L, He X, et al. Loss of Sirt1 promotes prostatic intraepithelial neoplasia, reduces mitophagy, and delays PARK2 translocation to mitochondria. Am. J. Pathol. 2015; 185:266–79.

11. Zhang S, Zhang D, Yi C, Wang Y, Wang H, Wang J. MicroRNA-22 functions as a tumor suppressor by targeting SIRT1 in renal cell carcinoma. Oncol. Rep. 2016; 35:559–67.

12. Grbesa I, Pajares MJ, Martínez-Terroba E, Agorreta J, Mikecin AM, Larráyoz M, Idoate MA, Gall-Troselj K, Pio R, Montuenga LM. Expression of sirtuin 1 and 2 is associated with poor prognosis in non-small cell lung cancer patients. PLoS One 2015; 10:e0124670.

13. Nassour J, Martien S, Martin N, Deruy E, Tomellini E, Malaquin N, Bouali F, Sabatier L, Wernert N, Pinte S, Gilson E, Pourtier A, Pluquet O, et al. Defective DNA single-strand break repair is responsible for senescence and neoplastic escape of epithelial cells. Nat. Commun. 2016; 7:10399.

14. Mileo AM, Miccadei S. Polyphenols as Modulator of Oxidative Stress in Cancer Disease: New Therapeutic Strategies. Oxid Med Cell Longev. 2016; 2016:6475624. https://doi.org/10.1155/2016/6475624.

15. Matt S, Hofmann TG. The DNA damage-induced cell death response: a roadmap to kill cancer cells. Cell Mol Life Sci. 2016; 73:2829–50.

16. Roos WP, Thomas AD, Kaina B. DNA damage and the balance between survival and death in cancer biology. Nat. Rev. Cancer 2016; 16:20–33.

17. O’Hagan HM, Mohammad HP, Baylin SB. Double strand breaks can initiate gene silencing and SIRT1-dependent onset of DNA methylation in an exogenous promoter CpG island. PLoS Genet. 2008; 4:e1000155.

18. O’Hagan HM, Wang W, Sen S, Destefano Shields C, Lee SS, Zhang YW, Clements EG, Cai Y, Van Neste L, Easwaran H, Casero RA, Sears CL, Baylin SB. Oxidative damage targets complexes containing DNA methyltransferases, SIRT1, and polycomb members to promoter CpG Islands. Cancer Cell. 2011; 20:606–19.

19. Oberdoerffer P, Michan S, Mcvay M, Mostoslavsky R, Vann J, Park SK, Hartlerode A, Stegmuller J, Hafner A, Loerch P, Wright SM, Mills KD, Bonni A, et al. SIRT1 redistribution on chromatin promotes genomic stability but alters gene expression during aging. Cell 2008; 135:907–18.

20. Meierjohann S. Oxidative stress in melanocyte senescence and melanoma transformation. Eur J Cell Biol. 2014; 93:36–41.

21. Denat L, Kadekaro AL, Marrot L, Leachman SA, Abdel-Malek ZA. Melanocytes as instigators and victims of oxidative stress. J Invest Dermatol. 2014; 134:1512–18. https://doi.org/10.1038/jid.2014.65.

22. Sander CS, Hamm F, Elsner P, Thiele JJ. Oxidative stress in malignant melanoma and non-melanoma skin cancer. Br J Dermatol. 2003; 148:913–22.

23. Oba-Shinjo SM, Correa M, Ricca TI, Molognoni F, Pinhal MA, Neves IA, Marie SK, Sampaio LO, Nader HB, Chammas R, Jasiulionis MG. Melanocyte transformation associated with substrate adhesion impediment. Neoplasia. 2006; 8:231–41.

24. Campos ACE, Molognoni F, Melo FHM, Galdieri LC, Carneiro CRW, D’Almeida V, Correa M, Jasiulionis MG. Oxidative stress modulates DNA methylation during melanocyte Anchorage blockade associated with malignant transformation. Neoplasia 2007; 9:1111–1121.

25. Molognoni F, Cruz AT, Meliso FM, Morais AS, Souza CF, Xander P, Bischof JF, Costa FF, Soares MB, Liang G, Jones PA, Jasiulionis MG. Epigenetic reprogramming as a key contributor to melanocyte malignant transformation. Epigenetics 2011; 6:450–464.

26. de Souza CF, Xander P, Monteiro AC, Silva AG, da Silva DC, Mai S, Bernardo V, Lopes JD, Jasiulionis MG. Mining gene expression signature for the detection of pre-malignant melanocytes and early melanomas with risk for metastasis. PLoS One 2012; 7:e44800.

27. Melo FHM, Molognoni F, Morais AS, Toricelli M, Mouroc MG, Higac EMS, Lopes JD, Jasiulionis MG. Endothelial nitric oxide synthase uncoupling as a key mediator of melanocyte malignant transformation associated with sustained stress conditions. Free Radic Biol Med. 2011; 50:1263–1273.

28. Molognoni F, De Melo FH, Da Silva CT, Jasiulionis MG. Ras and Rac1, frequently mutated in melanomas, are activated by superoxide anion, modulate Dnmt1 level and are causally related to melanocyte malignant transformation. PLoS One 2013; 8:e81937.

29. Grandori C, Cowley SM, James LP, Eisenman RN. The Myc/Max/Mad network and the transcriptional control of cell behavior. Annu. Rev. Cell Dev. Biol. 2000; 16:653–99.

30. Linka JM, Hurlin PJ. The activities of MYC, MNT and the MAX-interactome in lymphocyte proliferation and oncogenesis. Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms 2015; 1849:554–562.

31. Muhle-Goll C, Nilges M, Pastore A. The leucine zippers of the HLH-LZ proteins Max and c-Myc preferentially form heterodimers. Biochemistry 1995; 34:13554–64.

32. Adams JM, Harris AW, Pinkert CA, Corcoran LM, Alexander WS, Cory S, Palmiter RD, Brinster RL. The c-myc oncogene driven by immunoglobulin enhancers induces lymphoid malignancy in transgenic mice. Nature 1985; 318:533–538.

33. Schwab M, Alitalo K, Karl-Heinz K, Varmus HE, Bishop JM, Gilbert F, Brodeur G, Goldstein M, Trent J. Amplified DNA with limited homology to myc cellular oncogene is shared by human neuroblastoma cell lines and a neuroblastoma tumour. Nature 1983; 305:245 - 248.

34. Chakraborty AA, Scuoppo C, Dey S, Thomas LR, Lorey SL, Lowe SW, Tansey WP. A common functional consequence of tumor-derived mutations within c-MYC. Oncogene 2015; 34:2406–2409.

35. Shroff EH, Eberlin LS, Dang VM, Gouw AM, Gabay M, Adama SJ, Bellovin DI, Tran PT, Philbrick WM, Garcia-Ocana A, Casey SC, Lia Y, Dang CV, et al. MYC oncogene overexpression drives renal cell carcinoma in a mouse model through glutamine metabolism. Proc. Natl. Acad. Sci. USA 2015; 112:6539–6544.

36. Pan D, Kim B, Hu G, Gupta DS, Senpan A, Yang X, Schmieder A, Swain C, Wickline SA, Tomasson MA, Lanza GM. A strategy for combating melanoma with oncogenic c-Myc inhibitors and targeted nanotherapy. Nanomedicine (Lond). 2015; 10:241–251.

37. Mao B, Zhao G, Lv X, Chen HZ, Xue Z, Yang B, Liu DP, Liang CC. Sirt1 deacetylates c-Myc and promotes c-Myc/Max association. Int J Biochem Cell Biol. 2011; 43:1573–81.

38. Diolaiti D, McFerrin L, Carroll PA, Eisenman RN. Functional interactions among members of the MAX and MLX transcriptional network during oncogenesis. Biochim Biophys Acta. 2015; 1849:484–500.

39. Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia 2004, 6:1–6.

40. Talantov D, Mazumder A, Yu JX, Briggs T, Jiang Y, Backus J, Atkins D, Wang Y. Novel genes associated with malignant melanoma but not benign melanocytic lesions. Clin. Cancer Res. 2005; 11:7234–42.

41. Haqq C, Nosrati M, Sudilovsky D, Crothers J, Khodabakhsh D, Pulliam BL, Federman S, Miller JR 3rd, Allen RE, Singer MI, Leong SP, Ljung BM, Sagebiel RW, et al. The gene expression signatures of melanoma progression. Proc. Natl. Acad. Sci. USA. 2005; 102:6092–7.

42. Xu L, Shen SS, Hoshida Y, Subramanian A, Ross K, Brunet JP, Wagner SN, Ramaswamy S, Mesirov JP, Hynes RO. Gene expression changes in an animal melanoma model correlate with aggressiveness of human melanoma metastases. Mol. Cancer Res. 2008; 6:760–9.

43. Bogunovic D, O’Neill DW, Belitskaya-Levy I, Vacic V, Yu YL, Adams S, Darvishian F, Berman R, Shapiro R, Pavlick AC, Lonardi S, Zavadil J, Osman I, et al. Immune profile and mitotic index of metastatic melanoma lesions enhance clinical staging in predicting patient survival. Proc. Natl. Acad. Sci. USA 2009; 106:20429–34.

44. Kuzu OF, Nguyen FD, Noory MA, Sharma A. Current State of Animal (Mouse) Modeling in Melanoma Research. Cancer Growth Metastasis 2015; 8:81–94

45. Paoli P, Giannoni E, Chiarugi P. Anoikis molecular pathways and its role in cancer progression. Biochim. Biophys. Acta 2013; 1833:3481–98.

46. Cascón A, Robledo M. MAX and MYC: a heritable breakup. Cancer Res. 2012; 72:3119–24.

47. Kriegl L, Vieth M, Kirchner T, Menssen A. Up-regulation of c-MYC and SIRT1 expression correlates with malignant transformation in the serrated route to colorectal cancer. Oncotarget. 2012; 3:1182–93. https://doi.org/10.18632/oncotarget.628.

48. Li L, Osdal T, Ho Y, Chun S, McDonald T, Agarwal P, Lin A, Chu S, Qi J, Li L, Hsieh YT, Dos Santos C, Yuan H, et al. SIRT1 activation by a c-MYC oncogenic network promotes the maintenance and drug resistance of human FLT3-ITD acute myeloid leukemia stem cells. Cell Stem Cell. 2014; 15:431–46. https://doi.org/10.1016/j.stem.2014.08.001.

49. Hishida T, Nozaki Y, Nakachi Y, Mizuno Y, Iseki H, Katano M, Kamon M, Hirasaki M, Nishimoto M, Okazaki Y, Okuda A. Sirt1, p53, and p38(MAPK) are crucial regulators of detrimental phenotypes of embryonic stem cells with Max expression ablation. Stem Cells 2012; 30:1634–44.

50. Aprelikova O, Chen K, El Touny LH, Brignatz-Guittard C, Han J, Qiu T, Yang HH, Lee MP, Zhu M, Green JE. The epigenetic modifier JMJD6 is amplified in mammary tumors and cooperates with c-Myc to enhance cellular transformation, tumor progression, and metastasis. Clin. Epigenetics 2016; 14:8–38.

51. Okano M, Bell DW, Haber DA, Li E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999; 99:247–57.

52. Rhee I, Bachman KE, Park BH, Jair KW, Yen RWC, Schuebel KE, Cui H, Feinberg AP, Lengauer C, Kinzler KW, Baylin SB, Vogelstein B. DNMT1 and DNMT3b cooperate to silence genes in human cancer cells. Nature 2002; 416:552–556.

53. Bell RE, Golan T, Sheinboim D, Malcov H, Amar D, Salamon A, Liron T, Gelfman S, Gabet Y, Shamir R, Levy C. Enhancer methylation dynamics contribute to cancer plasticity and patient mortality. Genome Res. 2016; 26:601–11.

54. Vizoso M, Esteller M. DNA methylation plasticity contributes to the natural history of metastasis. Cell Cycle 2015; 14:2863–4.

55. Goding CR, Pei D, Lu X. Cancer: pathological nuclear reprogramming? Nat. Rev. Cancer 2014; 14:568–73.

56. Menssen A, Hydbring P, Kapelle K, Vervoorts J, Diebold J, Lüscher B, Larsson LG, Hermeking H. The c-MYC oncoprotein, the NAMPT enzyme, the SIRT1-inhibitor DBC1, and the SIRT1 deacetylase form a positive feedback loop. Proc Natl Acad Sci USA. 2012; 24:E187–96.

57. Bennett DC, Cooper PJ, Hart IR. A line of non-tumorigenic mouse melanocytes, syngeneic with the B16 melanoma and requiring a tumour promoter for growth. Int J Cancer. 1987; 39:414–18.

58. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, and 1000 Genome Project Data Processing Subgroup. The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics. 2009; 25:2078–79.

59. Cock PJ, Fields CJ, Goto N, Heuer ML, Rice PM. The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids Res. 2010; 38:1767–71. https://doi.org/10.1093/nar/gkp1137.

60. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS. Model-based Analysis of ChIP-Seq (MACS). Genome Biol. 2008; 9:137.

61. Stark R, Brown G. DiffBind: differential binding analysis of ChIP-Seq peak. 2011. data. http://bioconductor.org/packages/release/bioc/vignettes/DiffBind/inst/doc/DiffBind.pdf.

62. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15:550.

63. Lawrence M, Huber W, Pagès H, Aboyoun P, Carlson M, Gentleman R, Morgan MT, Carey VJ. Software for computing and annotating genomic ranges. PLOS Comput Biol. 2013; 9:e1003118.

64. Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009a; 4:44–57.

65. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009b; 37:1–13.

66. Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011; 6:e21800.

67. Beinder L, Faehrmann N, Wachtveitl R, Winterfeld I, Hartner A, Menendez-Castro C, Rauh M, Ruebner M, Huebner H, Noegel SC, Doerr HG, Rascher W, Fahlbusch FB. Detection of expressional changes induced by intrauterine growth restriction in the developing rat mammary gland via exploratory pathways analysis. PLoS One 2014; 9:e100504.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.
PII: 21457