SIRT6 histone deacetylase functions as a potential oncogene in human melanoma

Copyright: Garcia-Peterson et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC-BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

ABSTRACT

Melanoma is an aggressive skin cancer that can rapidly metastasize to become fatal, if not diagnosed early. Despite recent therapeutic advances, management of melanoma remains difficult. Therefore, novel molecular targets and strategies are required to manage this neoplasm. This study was undertaken to determine the role of the sirtuin SIRT6 in melanoma. Employing a panel of human melanoma cells and normal human melanocytes, we found significant SIRT6 mRNA and protein upregulation in melanoma cells. Further, using a tissue microarray coupled with quantitative Vectra analysis, we demonstrated significant SIRT6 overexpression in human melanoma tissues. Lentiviral short hairpin RNA-mediated knockdown of SIRT6 in A375 and Hs 294T human melanoma cells significantly decreased cell growth, viability, and colony formation, induced G1-phase arrest and increased senescence-associated beta-galactosidase staining. As autophagy is important in melanoma and is associated with SIRT6, we used a qPCR array to study SIRT6 knockdown in A375 cells. We found significant modulation in several genes and/or proteins (decreases in AKT1, ATG12, ATG3, ATG7, BAK1, BCL2L1, CLN3, CTSB, CTSS, DRAM2, HSP90AA1, IRGM, NPC1, SQSTM1, TNF, and BECN1; increases in GAA, ATG10). Our data suggests that increased SIRT6 expression may contribute to melanoma development and/or progression, potentially via senescence- and autophagy-related pathways.

INTRODUCTION

Melanoma is one of the deadliest forms of skin cancer, and can rapidly metastasize to become lethal if not diagnosed early or left untreated. It is estimated that 87,110 new cases of melanoma will be diagnosed in the US alone in 2017 [1]. Current options for melanoma management include surgical resection (if caught early without any sign of distant metastasis), radiation therapy, immunotherapy, and specific pathway targeted drugs such as B-Raf proto-oncogene, serine/threonine kinase (BRAF) inhibitors Vemurafenib (Zelboraf) and dabrafenib (Tafinlar) [2]. However, significant adverse events, as well as tumor resistance and recurrence frequently prevent these treatments from being completely successful [3-6]. Further, even with these advanced treatments, 9,730 deaths are estimated to occur in 2017 in the US due to this neoplasm [1]. Therefore, novel molecular targets and treatments are required for effective management of this cancer. Melanoma provides an ideal neoplasm for new molecularly-targeted drug development, since genetic alterations contribute greatly in the pathogenesis of this disease. In this direction, our laboratory has been interested in identifying the role and functional significance of the sirtuin family of proteins in melanoma [7-13].

The sirtuins (SIRTs) are a family of seven (SIRTs 1-7) NAD-dependent mammalian histone and protein deacetylases that regulate multiple pathways involved in aging, metabolism and cancer. These important enzymes have been found to impact a wide variety of cellular and organismal processes, including, metabolism, transcriptional regulation, DNA repair, and aging [13-16]. Since many of these processes are involved in cancer initiation and progression, sirtuins are considered an attractive target in cancer research. Although research to date indicates that sirtuins are involved in cancer, their role is extremely complex and appears to depend on the cellular context, with many sirtuins showing both tumor suppressor as well as oncogenic properties [17]. Like other sirtuins, nuclear SIRT6 is emerging as a critical player in a variety of cellular functions, including regulation of pathways related to gene transcription, glucose homeostasis, DNA repair and telomere integrity [18-20]. Interestingly, investigations into the exact role of SIRT6 in cancer appears to result in mixed answers, with available evidence for both tumor promoter as well as tumor suppressor roles of this sirtuin [18, 21-25]. However, the role of SIRT6 in melanoma is unknown. This study was undertaken to determine the role and functional significance of SIRT6 in melanoma.

RESULTS

SIRT6 is overexpressed in human melanoma cell lines and clinical tissues

As the first step of our investigation into understanding the role of SIRT6 in melanoma, we determined the expression profile of SIRT6 in normal human epidermal melanocytes (NHEM) and a panel of human melanoma cell lines (A375, Hs 294T, G361, SK-MEL-2, SK-MEL-28, WM35, WM115 and 451Lu). We found that compared to NHEM, most of the melanoma cell lines exhibited significantly higher expression of SIRT6 protein and mRNA as assessed by immunoblot and RT-qPCR analyses, respectively (Figure 1A-B). Next, to determine the clinical significance of SIRT6 in melanoma, we assessed the expression pattern of SIRT6 in human tissues using a human melanoma tissue microarray (TMA), containing a total of 62 cases of malignant melanoma, 20 metastatic malignant melanoma and 18 nevus tissues. Individual tissue cores were analyzed with Vectra/inForm quantitation. As shown in Figure 1C-E (and detailed in Supplementary Figure 1 and Supplementary Table 1), compared to melanocytic nevi, we found a significantly higher expression of SIRT6 in melanoma. Interestingly, we found SIRT6 overexpression both in the cytoplasm, as well as in nucleus of the cells (Figure 1E). Moreover, metastatic samples were found to have more SIRT6 than primary malignant melanoma samples. This differential expression of SIRT6 in melanoma encouraged us to further explore its role in melanoma.

SIRT6 knockdown inhibits growth and clonogenic survival of A375 and Hs 294T human melanoma cell lines

As an initial assessment to understand the functional significance of SIRT6 in melanoma, we employed lentiviral-mediated transient transfection of short hairpin RNA (shRNA) to knockdown the SIRT6 gene expression in A375 and Hs 294T cells. To determine the optimal shRNA for knockdown, we tested five different shRNA constructs and chose the one with the greatest knockdown as compared to nonsense shRNA (shNS) (Supplementary Figure S2) for subsequent studies. As shown in Figure 2A-B, we achieved significant knockdown of SIRT6 in both A375 and Hs 294T cells, both at the protein and mRNA levels. Next, we assessed the effects of SIRT6 knockdown on the proliferative potential of melanoma cells using the Promega RealTime-Glo MT Cell Viability Assay (Figure 2C) and the trypan blue exclusion assay (Figure 2D). Our data demonstrated that SIRT6 knockdown significantly decreased cell growth and viability in both of the melanoma cell lines tested. We further assessed the long-term clonogenic survival of the cells following SIRT6 knockdown. As seen in Figure 2E, SIRT6 knockdown resulted in a marked decrease in colony formation abilities of A375 and Hs 294T melanoma cells, suggesting that knockdown of SIRT6 in melanoma cells results in markedly diminished long-term proliferative potential.

SIRT6 knockdown induces G1-phase arrest and senescence-like phenotypes in human melanoma cells

To determine if the observed anti-proliferative response of SIRT6 knockdown in melanoma cells was associated with a dysregulation in cell cycle, we performed a cell cycle distribution analysis by flow cytometry (Figure 3A-B). We found that SIRT6 knockdown resulted in an enhanced accumulation of cells in G0/G1 phase in both A375 and Hs 294T human melanoma cell lines. This increase in G1 phase was accompanied by a concomitant decrease of cells in S phase and G2/M phase. The finding of increased G1 arrest, paired with observation of a multinucleated phenotype in our cell cultures after knockdown (data not shown), suggested the possible induction of senescence in the cells. To confirm this, we determined the effects of SIRT6 knockdown on senescence-associated beta-galactosidase (SA-β-Gal) staining in A375 and Hs 294T. As shown in Figure 3C, we found that SIRT6 knockdown resulted in an appreciable increase in SA-β-Gal positive cells, suggesting senescence induction.

SIRT6 knockdown causes modulation of autophagy-related genes in melanoma cells

Because autophagy and senescence are believed to be functionally intertwined processes [26], we next wanted to assess the effects of SIRT6 knockdown on autophagy-related pathways in A375 melanoma cells. To achieve this, we employed a commercially available autophagy PCR array containing 84 genes that are related to autophagy. A number of genes in the array were found to be modulated following SIRT6 knockdown in A375 cells. Detailed information and fold change in genes following SIRT6 knockdown are provided in Supplementary Table 4. A heat map generated from the data obtained is shown in Figure 4A. We found that 17 of the 84 genes tested were significantly modulated by two-fold or more (Figure 4B). Of these 17 genes, 2 were upregulated (ATG10 and GAA) and 15 were downregulated (AKT1, ATG12, ATG3, ATG7, BAK1, BCL2L1, CLN3, CTSB, CTSS, DRAM2, HSP90AA1, IRGM, NPC1, SQSTM1 and TNF). Out of the 17 genes, only CTSS and TNF were found to have more than 4-fold changes.

Modulations in autophagy-related genes are associated with the antiproliferative response of SIRT6 inhibition

In order to understand and connect the SIRT6-mediated modulations in autophagy-related genes, the 17 genes that were ≥2-fold modulated were analyzed using IPA and a network pathway was generated that connected 12 of the 17 genes that were found to be modulated in PCR array (Figure 5A). Interestingly, many of the network genes were found to have links with cancer (shown with a pink boundary). IPA was also able to identify links to genes which were not part of the PCR array (denoted as uncolored nodes), suggesting that they may be associated with SIRT6 signaling. IPA was further used to explore the predicted functions of the genes affected by SIRT6 knockdown. It was found that the SIRT6 knockdown-affected genes are associated with cell transformation, tumor invasion and movement of tumor cells (Figure 5B). The highest modulated gene, TNF, appears to affect all three functions predicted by IPA. However, the finding that BAK1 is downregulated (depicted with yellow lines), appears to be inconsistent with the predicted effect of SIRT6 inhibition.

To confirm the findings of our PCR array data, we performed RT-qPCR validation of the significantly regulated autophagy-related genes identified above. As shown in Figure 6A, the changes in these 17 genes due to SIRT6 knockdown were comparable to the array data. Further, to determine the modulations in autophagy markers at the protein level, we performed immunoblot analyses using specific antibodies against selected autophagy marker proteins. As shown in Figure 6B, we found significant modulation in various autophagy-related proteins in SIRT6 knockdown cells compared to nonsense control cells. SIRT6 knockdown was found to cause a marked decrease in the protein levels of LC3 II (autophagy related form of LC3), ATG3, ATG7, SQSTM1 and BECN in melanoma cells. However, GAA and ATG10 protein levels were found to be increased, following SIRT6 knockdown in melanoma cells. These results further confirmed that SIRT6 knockdown-mediated anti-proliferative response is associated with modulations in autophagy in melanoma cells.

DISCUSSION

The objective of this study was to determine the role of the nuclear sirtuin SIRT6 in melanoma. Based on limited available information to date, the role of SIRT6 in cancer is somewhat controversial. While some studies have suggested a tumor suppressor function of SIRT6 in liver and intestinal cancers [18, 21, 22], others have shown that it has oncogenic properties in prostate cancer and non-melanoma skin cancer [23-25]. Thus, the role of SIRT6 in cancer appears to be quite complex and possibly context- and tissue-specific. Interestingly, SIRT6 has been suggested to be an oncogene in keratinocyte-derived skin cancers. Ming et al demonstrated that skin-specific SIRT6 deletion inhibited tumorigenesis in mice [25]. Mechanistically, it was found that SIRT6 promoted COX-2 activation via repressing AMPK signaling, thereby increasing cell proliferation and survival in the epidermis [25]. Another finding supporting the pro-proliferative role of SIRT6 in keratinocytes stems from the proposed inverse association between SIRT6 and the tumor suppressor microRNA34a in keratinocyte differentiation [23]. These findings suggested that SIRT6 could be oncogenic in skin, especially in the keratinocyte-derived cells. Our study focused on melanocytic cells and we found strong evidences supporting the pro-proliferative function of SIRT6 in melanoma.

Our data demonstrated that SIRT6 is upregulated, both at mRNA and protein levels, in melanoma cell lines as well as clinical tissue samples of human melanoma. Thus, our data appears to be consistent with the pro-proliferative function of SIRT6 as seen in prostate and non-melanoma skin cancer [23-25]. In addition, we found that RNA interference-mediated knockdown of SIRT6 reduces cellular growth, viability, and colony formation in A375 and Hs 294T melanoma cells, further underlining the potential pro-proliferative function of SIRT6 in melanoma. Interestingly, our data demonstrated that SIRT6 knockdown induces senescence in melanoma cells as evident by morphological changes (characterized by multinucleated phenotype), as well as increased SA-β-Gal staining [27]. Our data also demonstrated that SIRT6 knockdown produced a marked accumulation of cells in G1 phase at the expense of G2/M and S phases. Senescent cells are generally believed to arrest in G1 or G0 stage of the cell cycle, thus indicating the possibility that SIRT6 knockdown causing the cells to become senescent. This is quite similar to the results we have reported earlier for SIRT3 [7].

In an effort to uncover the mechanism behind the anti-proliferative effects shown by the SIRT6 knockdown in melanoma we performed an autophagy PCR array. Autophagy is an extremely orchestrated process involving activities of various autophagy-related proteins known as ATGs. Recent studies have implicated SIRT6 in senescence [28-30] as well as autophagy [31, 32]. Indeed, autophagy and senescence are believed to have multiple and often overlapping and complementary functions in cancer, both in terms of influencing tumor development and in modulating the response to chemotherapy and/or radiation [33]. Autophagy regulatory proteins work by activating autophagy and beclin 1 regulator 1 (AMBRA1), ATG8/microtubule-associated light chain 3 (LC3), and Sequestosome1 (SQSTM1), followed by the inactivation of mammalian target of rapamycin (mTOR) to induce activation [34]. Interestingly, autophagy has been shown to suppress tumorigenesis by degrading oncogenic and toxic proteins [35, 36] or promote tumor development by reducing cancer cell vulnerability to stress and maintaining tumor cells homeostasis [37]. Our results demonstrate that SIRT6 knockdown significantly alters genes and pathways related to autophagy in melanoma cells both at mRNA and protein levels. Upon SIRT6 knockdown, modulations in various autophagy markers including BECN1, SQSTM1, ATG3, ATG7, ATG10 and GAA were confirmed by PCR array, RT-qPCR and/or immunoblotting. Moreover, SIRT6 knockdown reduced the conversion of LC3 protein from its free form LC3-I to phosphatidylethanolamine-conjugated form LC3-II, which represents an inhibition of a crucial step in autophagosome formation [38]. These data suggest that inhibition of SIRT6 by shRNA inhibits autophagy in melanoma cells. Indeed, SIRT6 has been previously suggested to affect the autophagy in human bronchial epithelial cells and esophageal cancer cells [31, 39, 40]. Based on our data, it appears that SIRT6 positively regulates autophagy in melanoma cells. This is an interesting observation because melanoma researchers seem to have come to a consensus that autophagy is tumor-suppressive in the early stages of cancer and tumor-promoting in established tumors (discussed in [41]). Indeed, future studies are required to carefully study the role of SIRT6 and autophagy in early versus late melanomas.

Overall, our study provides strong evidence suggesting a tumor promoter function of SIRT6 in melanoma, warranting future detailed investigations into the role, mechanism and therapeutic targetability of SIRT6 in melanoma. In vivo studies in relevant melanoma models are needed to further validate our in vitro findings.

MATERIALS AND METHODS

Cell culture

All cell lines are routinely (upon thawing, or if a cell line is in culture for more than 3 months) authenticated using STR testing through our core facility, the University of Wisconsin Translational Research Initiatives in Pathology (TRIP) lab. For our experiments, human melanoma cell lines (A375, Hs 294T, G361, SK-MEL-2, SK-MEL-28, WM35, WM115, 451Lu) and normal adult human melanocyte (NHEM) were purchased from ATCC, and HEK293T Lenti-X cells were purchased from TaKaRa Bio. SK-MEL-28, WM35 and SK-MEL-2 were maintained in EMEM supplemented with 1 mM sodium pyruvate and non-essential amino acids, G361 was maintained in McCoy’s 5a medium (Corning) and HEK293T Lenti-X, Mel-ST, A375, Hs 294T were maintained in DMEM (Corning), all with 10% fetal bovine serum (Sigma-Aldrich) and NHEMs were cultured in dermal cell medium supplemented with adult growth supplement (ATCC). All the cells were maintained at standard cell culture conditions (37°C, 5% CO₂ in humidified incubator).

Reverse transcription-quantitative real-time PCR (RT-qPCR)

For RT-qPCR analysis, total RNA was isolated from cells using the RNeasy Mini Kit (Qiagen), and first strand cDNA was transcribed using Oligo(dT) primers, dNTPs and M-MLV reverse transcriptase (Promega). RT-qPCR was perform using the BioRad CFX96 Touch and the BioRad iTaq Universal SYBR Green Supermix with first strand cDNA, forward and reverse gene-specific primer pairs. Melt curve analysis was performed to ensure the specificity of target amplicon. Relative mRNA levels were analyzed using GAPDH as endogenous control and ΔΔCT algorithm using the BioRad CFX Manager software.

Immunohistochemical (IHC) staining

Human melanoma TMA was purchased from US Biomax (Cat #ME1004c), and was stained as described previously [7]. Briefly, slides were deparaffinized in xylene, and rehydrated through a standard graded ethanol series. Antigen retrieval was performed using IHC-Tek epitope retrieval solution and steamer set (IHC World, LLC.). The slides were then immersed in 3% H₂ O₂ for endogenous peroxidase quenching. After blocking with goat serum, Vectastain ABC-AP Kit and Vector Red Alkaline Phosphatase Substrate Kit I (Vector Labs) were used for tissue staining as per manufacturer protocol. SIRT6 primary antibody (Cell Signaling, Cat#12486) was used at 1:150 dilution. The slides were then counterstained with Harris modified hematoxylin (ThermoFisher Scientific Inc.), dehydrated with graded ethanol and xylenes, and mounted with Permount before undergoing image analysis with the Vectra system, as described below.

TMA image acquisition and analysis

For quantitative staining analysis, the Vectra automated quantitative tissue imaging system, including an automated slide scanner and state-of-the-art software (Nuance and inForm; PerkinElmer) was used, as described previously [8]. This system takes advantage of the unique spectral characteristics (curve) of each chromogen via creation of a spectral library from single-chromogen slides in order to determine the relative contribution of each chromogen to the stained area on a slide. The stained TMA was imaged with the Vectra slide scanner, using a scanning protocol that was created based on the core size and layout, as well as acquisition of the spectral library. For each slide, an 8-bit image cube from each of the TMA tissue cores was acquired and inForm software was used to segment tissues (melanoma versus others) and cells (nucleus versus cytoplasm) to analyze the protein levels. Then, the target SIRT6 signals (Vector Red staining) were quantitated within the tissue and subcellular compartment(s) after unmixing the spectral curves with inForm software. Using these systems allowed for elimination of signal noises and cross-talk. Continuous signal intensity data (mean optical density per pixel) was generated for each TMA core or cellular compartment.

Lentiviral production and transduction

For viral creation, Lenti-X 293T cells were passaged to attain 30–40% confluency the next day and then transfected using the CaPO₄ method, as described previously [42]. Five SIRT6-targeting shRNA constructs were purchased and the clone with the best knockdown was chosen for further study. Briefly, 8 μg shRNA plasmid DNA (nonsense or SIRT6-targeting shRNA) (Sigma Aldrich, sequence information in Supplementary Table S2), 0.5 μg of VSV-G, and 0.6 μg of Δ8.2 plasmids were mixed with sterile ddH₂ O and CaCl_{2. The DNA mix was} bubbled, and an equal volume of 2x HEPES buffered saline (pH 7.05) was added dropwise. The transfection mixture was added to the appropriate dishes and incubated overnight. After 24 h, the transfection media was removed and fresh media was added. After an additional 24 and 48 h, cell media containing shRNA lentivirus was collected and filtered before use. For target cell transduction, 1.5-2.0 x 10⁴ melanoma cells (A375 and Hs 294T) were plated in each well of a 6-well plate. Viral media was added to triplicate sets of cells along with 8 μg/ml of polybrene a total of four times over 2 days. After 48 hours of transduction, the viral media was removed and cells were allowed to recover for 48 hours before being collected for further experiments.

Trypan blue exclusion and RealTime-Glo MT cell viability assays

For trypan blue exclusion assay, following treatments and recovery time as outlined above, the cells were trypsinized, an aliquot was stained with trypan blue dye (ThermoFisher Scientific), and the viable, dead, and total cells were determined using the Countess II FL Automated Cell Counter (ThermoFisher Scientific). To further confirm the effects of SIRT6 knockdown, we used RealTime-Glo MT cell viability assay (Promega), according to the manufacturer’s protocol. For this purpose, following treatments, 500 cells were plated per well in 96 well plates and transduced as above. Twenty-four hours after transduction, RealTime-Glo reagent was added to the wells, and luminescence was measured at 48 h using the BioTek Synergy H1 plate reader.

Clonogenic cell survival assay

Forty eight hours after transductions, 200 viable cells from each treatment group were plated in separate wells of a 6-well plate. Cells were maintained under standard tissue culture conditions, changing media every 3 days for approximately 10 days, until the cells reached optimal density. Cells were then stained at room temp for 30 m with 0.5% crystal violet solution in 1:1 methanol:water, rinsed 2 times with phosphate buffered saline (pH 7.4), and air dried followed by digital photography.

Senescence-associated β-galactosidase (SA-β-Gal) staining

Senescence-associated beta-galactosidase catalyzes the hydrolysis of β-galactosides into monosaccharides in senescent cells. To determine the induction of senescence, 1000 cells were plated per well in 6 well plates and transduced as above. Forty eight hours after transduction, cells were stained using the Senescence β-Galactosidase Staining Kit (Cell Signaling) per manufacturer’s instructions. The next day, the cells were washed with PBS and counterstained with Nuclear Fast Red (Vector Laboratories) for 30 sec. The images of SA-β-gal positive cells (with blue color sedimentation in cytoplasm) were captured with a Nikon DS-Fi1 microscope (Nikon Instruments Inc.) using NIS Elements AR 3.1 software.

Cell cycle analysis

To determine the impact of SIRT6 knockdown on cell cycle, 15000-18000 cells were plated per well in 6-well plates and transduced as above. Forty-eight hours after transduction, cells were trypsinized, washed twice with PBS and fixed in ice-cold 70% ethanol overnight. Cells were then pelleted, washed twice with PBS, treated with RNase, and stained with propidium iodide (Invitrogen). Cell cycle distribution was analyzed on a FACScan benchtop cytometer at the University of Wisconsin Paul P. Carbone Comprehensive Cancer Center (UWCCC) Flow cytometry facility. Cell cycle distribution was analyzed using CellQuest software (BD Biosciences). Cell cycle analysis was performed using ModFit software (Verity Software House)

Immunoblot analysis

Protein lysates were prepared, quantified and subjected to SDS-PAGE followed by transferring to nitrocellulose membranes. The proteins were probed using primary antibodies (SIRT6, #12486; β-tubulin, #2128; ATG3, #3415; ATG7, #8558; Beclin1, #3495; β-Actin, #4970 from Cell Signaling Technologies, LC3B, #NB600-1384 from Novus Biologicals and ATG10, #A9356; GAA, #G2048 from Sigma-Aldrich) and appropriate secondary antibodies (anti-rabbit, #7404; anti-mouse, #7076, Cell Signaling Technologies) followed by chemiluminescence detection with Pierce ECL Western Blotting Substrate (Thermo Scientific). Blots were imaged using the ImageQuant LAS 4000 (GE Healthcare).

Autophagy RT2 profiler PCR array analysis

RT² Profiler™ PCR Arrays for Human Autophagy were purchased from Qiagen (Cat. #PAHS-084Z) and were used and analyzed according to the manufacturer’s instructions, employing the TaKaRa SYBR Premix ExTaq II polymerase. Five reference genes (B2M, HPRT1, RPLP0, GAPDH, and ACTB) were used for normalizing data between shNS-A375 cells and shSIRT6-A375 cells. Data for PCR array were analyzed using SABiosciences RT² Profiler PCR Data Analysis software and were considered significant at ≥2 fold change and P<0.05. Further analyses of differentially expressed genes were performed using Ingenuity Pathway Analysis (IPA). The predicted gene-gene interaction networks and functions network were generated and analyzed using inputs of gene identifiers and fold-changes (shNS versus shSIRT6). Further, the PCR array data were validated using RT-qPCR analyses as detailed in section 2.2. Primer pairs for AKT1, ATG10, ATG12, ATG3, ATG7, BAK1, BCL2L1, CLN3, CTSB, CTSS, DRAM2, GAA, HSP90AA1, IRGM, NPC1, SQSTM1 and TNF were retrieved from PrimerBank [43] and are detailed in Supplementary Table 3.

Statistical Analysis

All quantitative data were analyzed using GraphPad PRISM 5.0 software (GraphPad Software, Inc., La Jolla, CA). For TMA, statistical significance was determined by one-way ANOVA with Tukey’s post hoc test or by two-tailed Student’s t-test. For all other quantitative data, statistical analyses were performed with two-tailed unpaired Student’s t-test between two experimental groups and one-way analysis of variance for more than two experimental groups followed by Sidak post hoc test. Data were shown as mean ±SE unless otherwise specified. A p-value < 0.05 was considered to be statistically significant.

ABBREVIATIONS

AKT1: AKT serine/threonine kinase 1 AMPK: Activated Protein kinase ATG10 : autophagy related 10 ATG12: Autophagy related 12 ATG3: Autophagy related 3 ATG7: Autophagy related 7 BAK1: BCL2 antagonist/killer 1 BCL2L1: BCL2 like 1 BRAF: B-Raf proto-oncogene, serine/threonine kinase

CLN3: CLN3, battenin COX-2: Cyclooxygenase-endoperoxide synthase 2 CTSS: Cathepsin S CTSB: Cathepsin B DRAM2: DNA damage regulated autophagy modulator 2

GAA: Glucosidase alpha, acid GAPDH: Glyceraldehyde-3-phosphate dehydrogenase

HSP90AA1: Heat shock protein 90 alpha family class A member 1

IRGM: Immunity related GTPase M LC3: Microtubule-associated protein 1A/1B-light chain 3

NHEM: Normal human epidermal melanocytes NPC1: NPC intracellular cholesterol transporter 1 PI: Propidium iodide RT-qPCR: Reverse transcription, real-time quantitative polymerase chain reaction

SA-β-Gal : Senescence-associated beta- galactosidase

shRNA: Short hairpin RNA SIRT: Sirtuin SQSTM1: Sequestosome 1 TMA: Tissue microarray TNF: Tumor necrosis factor

AUTHORS’ CONTRIBUTIONS

LG, MN, CKS, GC, and NA were involved in the conception and design of the study. LG, MN, CKS, GC, and WH were involved in data collection. LG, MN, CKS, GC, WH, and NA were involved in the data analysis and writing.

ACKNOWLEDGEMENTS

The authors would like to thank Ms. Grace Souksavath for her technical help.

CONFLICTS OF INTEREST

The authors have no conflicts of interest to declare.

FUNDING

This work was partially supported by funding from the NIH (R01AR059130 and R01CA176748 to NA), and the Department of Veterans Affairs (VA Merit Review Awards I01BX001008 and I01CX001441; and a Research Career Scientist Award IK6BX003780 to NA). This work was also supported by the core facilities in the Skin Diseases Research Center (SDRC) Core Grant P30AR066524 from NIH/NIAMS, as well as the University of Wisconsin Carbone Cancer Center Support Grant P30CA014520.

• 1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. 2017; 67: 7-30. [PubMed]
• 2. Maverakis E, Cornelius LA, Bowen GM, Phan T, Patel FB, Fitzmaurice S, He Y, Burrall B, Duong C, Kloxin AM, Sultani H, Wilken R, Martinez SR, et al. Metastatic melanoma - a review of current and future treatment options. Acta Derm Venereol. 2015; 95: 516-24. [PubMed]
• 3. Mavropoulos JC, Wang TS. Managing the skin toxicities from new melanoma drugs. Curr Treat Options Oncol. 2014; 15: 281-301. [PubMed]
• 4. Peuvrel L, Quereux G, Saint-Jean M, Brocard A, Nguyen JM, Khammari A, Knol AC, Varey E, Dreno B. Profile of vemurafenib-induced severe skin toxicities. J Eur Acad Dermatol Venereol. 2016; 30: 250-7. [PubMed]
• 5. McArthur GA, Chapman PB, Robert C, Larkin J, Haanen JB, Dummer R, Ribas A, Hogg D, Hamid O, Ascierto PA, Garbe C, Testori A, Maio M, et al. Safety and efficacy of vemurafenib in BRAF(V600E) and BRAF(V600K) mutation-positive melanoma (BRIM-3): extended follow-up of a phase 3, randomised, open-label study. Lancet Oncol. 2014; 15: 323-32. [PubMed] https://doi.org/10.1016/S1470-2045(14)70012-9.
• 6. Spagnolo F, Ghiorzo P, Queirolo P. Overcoming resistance to BRAF inhibition in BRAF-mutated metastatic melanoma. Oncotarget. 2014; 5: 10206-21. [PubMed] https://doi.org/10.18632/oncotarget.2602.
• 7. George J, Nihal M, Singh CK, Zhong W, Liu X, Ahmad N. Pro-Proliferative Function of Mitochondrial Sirtuin Deacetylase SIRT3 in Human Melanoma. J Invest Dermatol. 2016; 136: 809-18. [PubMed] https://doi.org/10.1016/j.jid.2015.12.026.
• 8. Wilking-Busch MJ, Ndiaye MA, Huang W, Ahmad N. Expression profile of SIRT2 in human melanoma and implications for sirtuin-based chemotherapy. Cell Cycle. 2017; 16: 574-7. [PubMed] https://doi.org/10.1080/15384101.2017.1288323.
• 9. Garcia-Peterson LM, Wilking-Busch MJ, Ndiaye MA, Philippe CGA, Setaluri V, Ahmad N. Sirtuins in Skin and Skin Cancers. Skin Pharmacol Physiol. 2017; 30: 216-24. [PubMed] https://doi.org/10.1159/000477417.
• 10. Wilking MJ, Singh C, Nihal M, Zhong W, Ahmad N. SIRT1 deacetylase is overexpressed in human melanoma and its small molecule inhibition imparts anti-proliferative response via p53 activation. Arch Biochem Biophys. 2014; 563: 94-100. [PubMed] https://doi.org/10.1016/j.abb.2014.04.001.
• 11. Singh CK, George J, Nihal M, Sabat G, Kumar R, Ahmad N. Novel downstream molecular targets of SIRT1 in melanoma: a quantitative proteomics approach. Oncotarget. 2014; 5: 1987-99. [PubMed] https://doi.org/10.18632/oncotarget.1898.
• 12. Wilking MJ, Singh CK, Nihal M, Ndiaye MA, Ahmad N. Sirtuin deacetylases: a new target for melanoma management. Cell Cycle. 2014; 13: 2821-6. [PubMed] https://doi.org/10.4161/15384101.2014.949085.
• 13. Singh CK, Chhabra G, Ndiaye M, Garcia-Peterson LM, Mack NJ, Ahmad N. The Role of Sirtuins in Antioxidant and Redox Signaling. Antioxid Redox Signal. 2017: [Epub ahead of print]. [PubMed] https://doi.org/10.1089/ars.2017.7290.
• 14. Osborne B, Bentley NL, Montgomery MK, Turner N. The role of mitochondrial sirtuins in health and disease. Free Radic Biol Med. 2016; 100: 164-74. [PubMed]
• 15. Dai Y, Faller DV. Transcription Regulation by Class III Histone Deacetylases (HDACs)-Sirtuins. Transl Oncogenomics. 2008; 3: 53-65. [PubMed]
• 16. Kraemer KH, Levy DD, Parris CN, Gozukara EM, Moriwaki S, Adelberg S, Seidman MM. Xeroderma pigmentosum and related disorders: examining the linkage between defective DNA repair and cancer. J Invest Dermatol. 1994; 103: 96S-101S. [PubMed]
• 17. Bosch-Presegue L, Vaquero A. The dual role of sirtuins in cancer. Genes Cancer. 2011; 2: 648-62. [PubMed] https://doi.org/10.1177/1947601911417862.
• 18. Sebastian C, Zwaans BM, Silberman DM, Gymrek M, Goren A, Zhong L, Ram O, Truelove J, Guimaraes AR, Toiber D, Cosentino C, Greenson JK, MacDonald AI, et al. The histone deacetylase SIRT6 is a tumor suppressor that controls cancer metabolism. Cell. 2012; 151: 1185-99. [PubMed] https://doi.org/10.1016/j.cell.2012.10.047.
• 19. Tennen RI, Bua DJ, Wright WE, Chua KF. SIRT6 is required for maintenance of telomere position effect in human cells. Nat Commun. 2011; 2: 433. [PubMed] https://doi.org/10.1038/ncomms1443.
• 20. Lombard DB, Schwer B, Alt FW, Mostoslavsky R. SIRT6 in DNA repair, metabolism and ageing. J Intern Med. 2008; 263: 128-41. [PubMed] https://doi.org/10.1111/j.1365-2796.2007.01902.x.
• 21. Marquardt JU, Fischer K, Baus K, Kashyap A, Ma S, Krupp M, Linke M, Teufel A, Zechner U, Strand D, Thorgeirsson SS, Galle PR, Strand S. Sirtuin-6-dependent genetic and epigenetic alterations are associated with poor clinical outcome in hepatocellular carcinoma patients. Hepatology. 2013; 58: 1054-64. [PubMed] https://doi.org/10.1002/hep.26413.
• 22. Min L, Ji Y, Bakiri L, Qiu Z, Cen J, Chen X, Chen L, Scheuch H, Zheng H, Qin L, Zatloukal K, Hui L, Wagner EF. Liver cancer initiation is controlled by AP-1 through SIRT6-dependent inhibition of survivin. Nat Cell Biol. 2012; 14: 1203-11. [PubMed]
• 23. Lefort K, Brooks Y, Ostano P, Cario-Andre M, Calpini V, Guinea-Viniegra J, Albinger-Hegyi A, Hoetzenecker W, Kolfschoten I, Wagner EF, Werner S, Dotto GP. A miR-34aSIRT6 axis in the squamous cell differentiation network. EMBO J. 2013; 32: 2248-63. [PubMed] https://doi.org/10.1038/emboj.2013.156.
• 24. Liu Y, Xie QR, Wang B, Shao J, Zhang T, Liu T, Huang G, Xia W. Inhibition of SIRT6 in prostate cancer reduces cell viability and increases sensitivity to chemotherapeutics. Protein Cell. 2013; 4: 702-10. [PubMed] https://doi.org/10.1007/s13238-013-3054-5.
• 25. Ming M, Han W, Zhao B, Sundaresan NR, Deng CX, Gupta MP, He YY. SIRT6 promotes COX-2 expression and acts as an oncogene in skin cancer. Cancer Res. 2014; 74: 5925-33. [PubMed] https://doi.org/10.1158/0008-5472.CAN-14-1308.
• 26. White E, Lowe SW. Eating to exit: autophagy-enabled senescence revealed. Genes Dev. 2009; 23: 784-7. [PubMed]
• 27. Dimri GP, Lee X, Basile G, Acosta M, Scott G, Roskelley C, Medrano EE, Linskens M, Rubelj I, Pereira-Smith O, et al. A biomarker that identifies senescent human cells in culture and in aging skin in vivo. Proc Natl Acad Sci U S A. 1995; 92: 9363-7. [PubMed] https://doi.org/10.1073/pnas.92.20.9363.
• 28. Tasselli L, Xi Y, Zheng W, Tennen RI, Odrowaz Z, Simeoni F, Li W, Chua KF, Lerrer B, Gertler AA, Cohen HY. SIRT6 deacetylates H3K18ac at pericentric chromatin to prevent mitotic errors and cellular senescence. Nat Struct Mol Biol. 2016; 23: 434-40. [PubMed] https://doi.org/10.1038/nsmb.3202.
• 29. Tasselli L, Zheng W, Chua KF. SIRT6: Novel Mechanisms and Links to Aging and Disease. Trends Endocrinol Metab. 2017; 28: 168-85. [PubMed] https://doi.org/10.1016/j.tem.2016.10.002.
• 30. Vitiello M, Zullo A, Servillo L, Mancini FP, Borriello A, Giovane A, Della Ragione F, D’Onofrio N, Balestrieri ML, Lee N, Ryu HG, Kwon JH, Kim DK, et al. Multiple pathways of SIRT6 at the crossroads in the control of longevity, cancer, and cardiovascular diseases. Ageing Res Rev. 2017; 35: 301-11. [PubMed]
• 31. Huang N, Liu Z, Zhu J, Cui Z, Li Y, Yu Y, Sun F, Pan Q, Yang Q, Shao J, Yang X, Liu T, Zhang T, et al. Sirtuin 6 plays an oncogenic role and induces cell autophagy in esophageal cancer cells. Tumour Biol. 2017; 39: 1010428317708532. [PubMed]
• 32. Shun CT, Lin SK, Hong CY, Lin CF, Liu CM, Takasaka N, Araya J, Hara H, Ito S, Kobayashi K, Kurita Y, Wakui H, Yoshii Y, et al. Sirtuin 6 modulates hypoxia-induced autophagy in nasal polyp fibroblasts via inhibition of glycolysis. Am J Rhinol Allergy. 2016; 30: 179-85. [PubMed]
• 33. Gewirtz DA. Autophagy and senescence in cancer therapy. J Cell Physiol. 2014; 229: 6-9. [PubMed]
• 34. Meijer AJ, Codogno P. Regulation and role of autophagy in mammalian cells. Int J Biochem Cell Biol. 2004; 36: 244562. [PubMed]
• 35. Gomes LC, Odedra D, Dikic I, Pohl C. Autophagy and modular restructuring of metabolism control germline tumor differentiation and proliferation in C. elegans. Autophagy. 2016; 12: 529-46. [PubMed] https://doi.org/10.1080/15548627.2015.1136771.
• 36. Gozuacik D, Kimchi A. Autophagy as a cell death and tumor suppressor mechanism. Oncogene. 2004; 23: 2891906. [PubMed]
• 37. Degenhardt K, Mathew R, Beaudoin B, Bray K, Anderson D, Chen G, Mukherjee C, Shi Y, Gelinas C, Fan Y, Nelson DA, Jin S, White E. Autophagy promotes tumor cell survival and restricts necrosis, inflammation, and tumorigenesis. Cancer Cell. 2006; 10: 51-64. [PubMed] https://doi.org/10.1016/j.ccr.2006.06.001.
• 38. Tanida I, Ueno T, Kominami E. LC3 and Autophagy. Methods Mol Biol. 2008; 445: 77-88. [PubMed]
• 39. Lu J, Sun D, Liu Z, Li M, Hong H, Liu C, Gao S, Li H, Cai Y, Chen S, Li Z, Ye J, Liu P. SIRT6 suppresses isoproterenol-induced cardiac hypertrophy through activation of autophagy. Transl Res. 2016; 172: 96-112 e6. [PubMed]
• 40. Takasaka N, Araya J, Hara H, Ito S, Kobayashi K, Kurita Y, Wakui H, Yoshii Y, Yumino Y, Fujii S, Minagawa S, Tsurushige C, Kojima J, et al. Autophagy induction by SIRT6 through attenuation of insulin-like growth factor signaling is involved in the regulation of human bronchial epithelial cell senescence. J Immunol. 2014; 192: 958-68. [PubMed]
• 41. Ndoye A, Weeraratna AT. Autophagy- An emerging target for melanoma therapy. F1000Res. 2016; 5. [PubMed] https://doi.org/10.12688/f1000research.8347.1.
• 42. Schmit TL, Nihal M, Ndiaye M, Setaluri V, Spiegelman VS, Ahmad N. Numb regulates stability and localization of the mitotic kinase PLK1 and is required for transit through mitosis. Cancer Res. 2012; 72: 3864-72. [PubMed] https://doi.org/10.1158/0008-5472.CAN-12-0714.
• 43. Wang X, Spandidos A, Wang H, Seed B. PrimerBank: a PCR primer database for quantitative gene expression analysis, 2012 update. Nucleic Acids Res. 2012; 40: D11449. [PubMed] https://doi.org/10.1093/nar/gkr1013.