Zebrafish modeling reveals that SPINT1 regulates the aggressiveness of skin cutaneous melanoma and its crosstalk with tumor immune microenvironment
Journal of Experimental & Clinical Cancer Research volume 38, Article number: 405 (2019)
Skin cutaneous melanoma (SKCM) is the most lethal form of skin cancer and while incidence rates are declining for most cancers, they have been steadily rising for SKCM. Serine protease inhibitor, kunitz-type, 1 (SPINT1) is a type II transmembrane serine protease inhibitor that has been shown to be involved in the development of several types of cancer, such as squamous cell carcinoma and colorectal cancer.
We used the unique advantages of the zebrafish to model the impact of Spint1a deficiency in early transformation, progression and metastatic invasion of SKCM together with in silico analysis of the occurrence and relevance of SPINT1 genetic alterations of the SKCM TCGA cohort.
We report here a high prevalence of SPINT1 genetic alterations in SKCM patients and their association with altered tumor immune microenvironment and poor patient survival. The zebrafish model reveals that Spint1a deficiency facilitates oncogenic transformation, regulates the tumor immune microenvironment crosstalk, accelerates the onset of SKCM and promotes metastatic invasion. Notably, Spint1a deficiency is required at both cell autonomous and non-autonomous levels to enhance invasiveness of SKCM.
These results reveal a novel therapeutic target for SKCM.
Skin cutaneous melanoma (SKCM) originates from melanocytes, neural-crest derived pigment-producing cells located in the epidermis, where their major function is to protect keratinocytes from UV-induced DNA damage . The malignant transformation of melanocytes generates this fatal form of skin cancer with a complex multigenic etiology that becomes extremely difficult to treat once it has metastasized. SKCM is the deadliest form of skin cancer (75% of deaths related to skin cancer) and it is common in the Western world. Indeed, its global incidence is 15–25 per 100,000 individuals . While incidence rates are declining for most cancers, they have been steadily rising for SKCM worldwide . Early detection is fundamental, since localized, early stage SKCM can be surgically excised with little chance of recurrence with a 98.2% of patient survival rate after 5-year survival as reported by The Surveillance, Epidemiology, and End Results (SEER) . Metastatic SKCM, however, is still an often fatal disease with a 5-year survival rate of 15–20% .
SKCM is one of the most recurrent types of cancer and its genetic heterogeneity has led in recent years to join forces to determine SKCM causes and develop effective therapies. Transformation of melanocytes into primary and then metastatic SKCM requires a complex interplay of exogenous and endogenous events . More than 50% of the tumors originate from normal skin rather than from dysplastic nevi, suggesting that SKCM not only appears to be due to the transformation of mature melanocytes . In this way, the identification of SKCM initiating cells is really important to devising methods for early detection and eradication of SKCM [7, 8]. Moreover, SKCM stem cell populations have been characterized and associated with tumor progression, immunoevasion, drug resistance and metastasis .
Inflammation can play a key role in cancer, from initiation of the transformed phenotype to metastatic spread. Nevertheless, inflammation and cancer have a profound yet ambiguous relationship. Inflammation (especially chronic inflammation) has protumorigenic effects, but inflammatory cells also mediate an immune response against the tumor and immunosuppression is known to increase the risk for certain tumors . Nowadays, skin cancers are also attributed to chronically injured or non-healing wounds and scars or ulcers that occur at sites of previous burns, sinuses, trauma, osteomyelitis, prolonged heat and chronic friction. The incidence of malignancy in scar tissues is 0.1–2.5% and it is estimated that underlying infections and inflammatory responses are linked to 15–20% of all deaths from cancer worldwide . Furthermore, chronic inflammation contributes to about 20% of all human cancers .
Serine protease inhibitor, kunitz-type, 1 (SPINT1), also known as hepatocyte growth factor activator inhibitor 1 (HAI1), is a type II transmembrane serine protease inhibitor that plays a crucial role in the regulation of the proteolytic activity of both suppression of tumorigenicity 14 (ST14), also known as Matriptase-1 [13,14,15], and Hepatocyte growth factor activator (HGFA) . The functional linkage between ST14 and SPINT1 has important implications for the development of cancer. ST14 activity, which is only partially opposed by endogenous SPINT1, causes increased sensitivity to carcinogens and produces spontaneous tumorigenesis in the skin of keratin-5-matriptase transgenic mice, while increased epidermal SPINT1 expression fully counteracts the oncogenic effect of ST14 . Furthermore, the expression of ST14 has been demonstrated to be up-regulated in various human cancer histotypes such as breast, cervix, ovaries, prostate, esophagus and liver cancers .
The close functional relationship between ST14 and SPINT1 was also observed in a zebrafish model of skin inflammation, carrying a hypomorphic mutation of spint1a. Indeed epidermal hyperproliferation and neutrophil infiltration observed in mutant zebrafish larvae are both rescued by st14a gene knock-down, suggesting a novel role for the SPINT1-ST14 axis in regulating inflammation [19, 20]. Given the unique advantages of the zebrafish model for tumor cell transplantation, thanks to the transparency of embryos and even adults that allows the tracking of donor cells and the strong correlation between alterations of Spint1a-St14a levels with tumor progression, the spint1a mutant zebrafish represents an attractive model to study the role of SPINT1 and chronic inflammation in SKCM.
Our results show that genetic alterations of SPINT1 correlated with a poor prognosis of SKCM patients and provide evidence that SPINT1 expression positively correlated with tumor macrophage infiltration, but not neutrophils. In line with these clinical data, we show that Spint1a deficiency enhances at both cell autonomous and non-autonomous levels cell dissemination of SKCM in zebrafish models by promoting tumor dedifferentiation and altered immune surveillance.
Wild-type zebrafish (Danio rerio H. Cypriniformes, Cyprinidae) were obtained from the Zebrafish International Resource Center (ZIRC, Oregon, USA) and mated, staged, raised and processed as described in the zebrafish handbook . Zebrafish fertilized eggs were obtained from natural spawning of wild type and transgenic fish held at our facilities following standard husbandry practices. Animals were maintained in a 12 h light/dark cycle at 28 °C. Tg(kita:GalTA4,UAS:mCherry)hzm1 zebrafish were crossed with Tg(UAS:eGFP-H-RAS_G12V)io6 line  to express oncogenic human HRAS_G12 V driven by the melanocyte cell-specific promoter kita. Transgene expression was confirmed in 3 day post-fertilization larvae (dpf) by fluorescence microscopy and by analyzing melanocyte hyperproliferation . The hi2217 line, which carries a hypomorphic spint1a mutant allele that promotes skin inflammation , and transparent roya9/a9; nacrew2/w2 (casper)  of 4–8 month old were previously described.
Zebrafish larvae were anesthetized by a solution of 0.16 mg/ml buffered tricaine (Sigma-Aldrich) in embryo medium. Adult zebrafish were anesthetized by a dual anesthetic protocol to minimize over-exposure to tricaine, in long-term studies (up to 40 min) . Briefly, the anesthesia was firstly induced by tricaine and then fish were transferred to tricaine/isoflurane solution (forane in ethanol, 1:9).
Tumor sampling, disaggregation and cell sorting
Primary melanoma tumors were excised from adult zebrafish once they had reached between 3 and 5 mm in diameter. Some individuals were euthanized according the European Union Council and IUAC protocol and others were monitored and maintained still alive after the tumor biopsy treated with conditioners to reduce fish stress and heal damaged tissue and wounds (STRESS COAT, API), as well as to protect from bacterial (MELAFIX, API) and fungal infections (PIMAFIX, API).
The tumor was excised with a clean scalpel and razor blade, placed in 2 ml of dissection media, composed by DMEM/F12 (Life Technologies), 100 UI/ml penicillin, 100 μg/ml streptomycin, 0.075 mg/ml Liberase (Roche). After manually disaggregation with a clean razor blade and incubation at room temperature for 30 min, 5 ml of wash media, composed by DMEM/F12 (Life Technologies), penicillin-streptomycin (Life Technologies), and 15% heat-inactivated FBS (Life Technologies), were added to the tumor slurry and manually disaggregated one last time. Next, the tumor cell suspensions were passed through a 40 μm filter (BD) into a clean 50 ml tube. An additional 5 ml of wash media was added to the initial tumor slurry and additionally filtered. This procedure was repeated twice. Cells were counted with a hemocytometer and the tubes of resuspended cells were centrifuged at 800 g for 5 min at 4 °C. The pellets of tumor cells were resuspended in PBS containing 5% FBS and kept on ice prior to transplantation .
The resulting cell suspension from zebrafish melanoma tumors was passed through a 40 μm cell strainer and propidium iodide (PI) was used as a vital dye to exclude dead cells. The Cell Sorting was performed on a “Cell Sorter” SONY SH800Z in which eGFP positive cells were sorted from the negative ones of the same cell tumor suspension.
Larval allotransplantion assays
Melanomas were disaggregated, then labelled with 1,1′-di-octa-decyl-3,3,3′,3′-tetra-methyl-indo-carbo-cya-nine perchlorate (DiI, ThermoFisher) and finally resuspended in a buffer containing 5% FBS in PBS. Between 25 to 50 cells/embryo were then injected in the yolk sac of 60–90 Casper or spint1 mutants zebrafish larvae of 48 h post-fertilization (hpf) and after 5 days at 28 °C, larvae were scored manually in blind samples for zebrafish melanoma cells dissemination by fluorescence microscopy . Melanoma cell invasion score was calculated as the percentage of zebrafish melanoma cell-invaded larvae over the total number of larvae analyzed taking into account also the number of tumor foci per larvae. Three tumor foci were established to score a larva as positive for invasion. Furthermore, larvae positive for invasion were also distinguished in three groups considering the number of positive foci per larvae: 3–5 foci per larvae, 5–15 foci per larvae and > 15 foci per larvae.
Adult allotransplantion assays
Adult zebrafish used as transplant recipients were immunosuppressed to prevent rejection of the donor material. Thus, the recipients were anesthetized, using the dual anesthetic protocol described above, and treated with 30 Gy of split dose sub-lethal X-irradiation (YXLON SMART 200E, 200 kV, 4.5 mA) 2 days before the transplantation. Then the immunosuppressed fish were maintained in fresh fish water treated with conditioners preventing infections onset and the consequent recipient deaths.
Anesthetized fish (10–20 per tumor) were placed dorsal side up on a damp sponge and stabilized with one hand. Using the other hand, the needle of a 10 μl beveled, 26S-guaged syringe (Hamilton) was positioned midline and ahead to the dorsal fin. 30,000, 100,000, 300,000 and 500,000 cells resuspended in PBS were injected into the dorsal subcutaneous cavity. The syringe was washed in 70% ethanol and rinsed with PBS between uses. Following transplantation, fish were placed into a recovery tank of fresh fish water and kept off-flow with daily water changes for 7 days. Large and pigmented tumors engrafted and were observed to expand by 10 days post-transplantation.
SKCM imaging in adult zebrafish
Adult zebrafish were scored weekly for melanoma formation starting at the first appearance of raised lesions. Tumor scoring was blinded and experiments were independently repeated at least 3 times. Zebrafish were anesthetized, placed in a dish of fish water and photographed using a mounted camera (Nikon D3100 with a Nikon AF-S Micro Lens). The pigmented tumor size was represented as the number of pigmented pixels (Adobe Photoshop CS5).
Analysis of gene expression
Once zebrafish tumors reached between 3 and 5 mm of diameter, they were excised and total RNA was extracted with TRIzol reagent (Invitrogen), following the manufacturer’s instructions, and then treated with DNase I, amplification grade (1 U/μg RNA; Invitrogen). SuperScript III RNase H− Reverse Transcriptase (Invitrogen) was used to synthesize first-strand cDNA with oligo (dT)18 primer from 1 μg of total RNA at 50 °C for 50 min. Real-time PCR was performed with an ABI PRISM 7500 instrument (Applied Biosystems) using SYBR Green PCR Core Reagents (Applied Biosystems). Reaction mixtures were incubated for 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C, 1 min at 60 °C, and finally 15 s at 95 °C, 1 min 60 °C and 15 s at 95 °C. For each mRNA, gene expression was normalized to the ribosomal protein S11 gene (rps11) content in each sample Pfaffl method . The primers used are shown in Additional file 1: Table S1. In all cases, each PCR was performed with triplicate samples and repeated at least in two independent samples.
Human SKCM dataset analysis
Normalized gene expression, patient survival data, genetic alterations and neutrophil/macrophage infiltration were downloaded on 14th December 2017 from SKCM repository of The Cancer Genome Atlas (TCGA, Provisional) from cBioPortal database (https://www.cbioportal.org/). This cohort includes 479 patient samples. No cleaning filter was applied and all patient data available for each study were used: 469 samples for expression analysis, 461 samples for survival analysis and 468 for immune cell infiltration. Transcript levels of SPINT1 in human samples from normal skin, benign and malignant melanoma was collected from Gene Expression Omnibus (GDS1375 dataset and 202826_at probe). Gene expression plots and regression curves for correlation studies were obtained using GraphPad Prism 5.03 (GraphPad Software).
Data are shown as mean ± SEM and they were analyzed by analysis of variance (ANOVA) and a Tukey multiple range test to determine differences between groups. The survival curves were analysed using the log-rank (Mantel-Cox) test. All the experiments were performed at least three times, unless otherwise indicated. The sample size for each treatment is indicated in the graph and/or in the figure legend. Statistical significance was defined as p < 0.05.
SPINT1 genetic alterations are associated with poor prognosis of SKCM patients and altered tumor immune microenvironment
To study the impact of SPINT1 in promoting SKCM progression and aggressiveness, an in silico analysis of human SKCM samples of the TCGA cohort was performed. This analysis revealed that genetic alterations occurred in 10% of SKCM patients; a relevant percentage comparing with major SKCM driven oncogenes and tumor suppressors (Fig. 1a). Among these genetic alterations, an increased mRNA level was the most prevalent alteration (7%), while 1.9% missense mutations of unknown significance and 1.9% deep deletions were also observed. Notably, these genetic alterations of SPINT1 significantly correlated with poor SKCM patient prognosis (Fig. 1b) and SPINT1 expression was significantly inhibited in human SKCM comparing with nevus and normal skin (Fig. 1c). In addition, the analysis of patient survival with increased SPINT1 transcript levels also revealed their poor prognosis but, unfortunately, the analysis of deep deletions and missense mutations separately gave no statistical significance because the number of patients for each condition was very low (data not shown). We next performed a GO enrichment analysis of biological process (Fig. 1d), analyzing the differentially expressed genes in SKCM samples of the TCGA cohort with missense mutations or copy-number alterations of SPINT1. The results showed that regulation of immune system, inflammatory response, cell cycle, cell adhesion, and extracellular matrix organization represent key pathways significantly affected in human SKCM with these SPINT1 genetic alterations. Collectively, these results suggest that both high and low levels of SPINT1 result in an unbalanced crosstalk between tumor cells and their microenvironment promoting higher aggressiveness.
The tumor microenvironment contains diverse leukocyte populations, including neutrophils, eosinophils, dendritic cells, macrophages, mast cells and lymphocytes . It is known that tumor-associated macrophages (TAM) are able to interact with tumor cells and can promote cancer progression [27,28,29]. As shown in Fig. 1e, the number of TAM in human SKCM samples correlated with the mRNA levels of SPINT1 in metastatic SKCM. However, the number of tumor-associated neutrophils (TAN) was independent of SPINT1 levels in both primary and metastatic SKCM. These data further confirmed the role of SPINT1 in the regulation of the crosstalk between tumor and inflammatory cells in human SKCM.
The expression of SPINT1 positively correlates with both inflammation and macrophage markers in human SKCM biopsies
In order to further understand the role of SPINT1 in SKCM, the RNA Seq database of the large TCGA cohort of SKCM was analyzed in terms of the expression of SOX10, TYR and DCT genes, that have been shown to be important in melanocyte development [30, 31]. In addition, SOX10 is a recognized biomarker for the diagnosis of SKCM . It was found that SPINT1 expression positively correlated with those of SOX10 and TYR, while a negative correlation was found between the expression of SPINT1 and DCT (Fig. 2a). The expression of the epithelial to mesenchymal transition (EMT) markers ZEB1, ZEB2 and TWIST1, but not TWIST2, negatively correlated with that of SPINT1 in SKCM (Fig. 2b).
SKCM cells release several cytokines and chemokines that recruit and polarize macrophages . Therefore, several inflammation markers were analyzed and only the expression of the genes encoding the receptor of the pro-inflammatory cytokine TNFα (TNFR1) and the receptor of the pro-inflammatory chemokine interleukin 8 (CXCR2), positively correlated with SPINT1 levels (Fig. 2c). Notably, the macrophage marker MFAP4 also positively correlated with SPINT1 expression (Fig. 2d). However, the M2 polarization marker CD163 (Fig. 2d) and several interferon-stimulated genes (ISGs) (Fig. 2e) were all unaffected by SPINT1 levels. Collectively, these results further suggest that SPINT1 regulates SKCM differentiation and aggressiveness, and macrophages infiltration.
Inflammation accelerates the onset of SKCM in zebrafish
Given the strong correlation between alterations of SPINT1 levels with the progression of SKCM and the crosstalk with the tumor immune microenvironment, we crossed the zebrafish line kita:Gal4;HRAS-G12V, which expresses the human oncogene HRAS-G12V in melanocytes and spontaneously develops SKCM , with the zebrafish mutant line spint1ahi2217Tg/hi2217Tg , which presents chronic skin inflammation (Fig. 3a). Firstly, we quantified by fluorescence microscopy the number of early oncogenically transformed goblet cells, which also expressed the kita promoter , in spint1a-deficient larvae and their wild type siblings (Fig. 3b). The results showed that spint1a deficiency resulted in increased number of HRAS-G12V+ cells (Fig. 3c). To determine if the enhanced Spint1a deficiency-driven oncogenic transformation was also able to promote SKCM aggressiveness, SKCM development in spint1ahi2217Tg/hi2217Tg fish were compared with wild type (spint1a+/hi2217Tg) from the end of metamorphosis stage (between 28 and 30 dpf) to 120 dpf (adult stage) (Fig. 3d-f). The resulting Kaplan-Meier curve showed a significant increased incidence of melanoma in the Spint1a-deficient fish, which developed SKCM in more than 50% of cases at 50–60 dpf compared with their wild type siblings which reached only 30% at this age (Fig. 3f). These data suggest that Spint1a deficiency increases oncogenic transformation and accelerates SKCM onset in vivo.
Spint1a deficiency is required at cell autonomous and non-autonomous levels to enhance SKCM cell dissemination in a zebrafish larval allotrasplantation model
To assess the in vivo role of Spint1a deficiency in SKCM invasiveness, SKCM tumors from spint1ahi2217Tg/hi2217Tg; kita:Gal4;HRAS-G12V and spint1a+/hi2217Tg; kita:Gal4;HRAS-G12V were disaggregated, after staining the cells were transplanted into the yolk sac of 2 dpf casper zebrafish larvae (Fig. 4a). The results showed that Spint1a deficiency in SKCM cells enhanced the dissemination of SKCM, assayed as the percentage of invaded larvae and the number of foci per larva, compared to control SKCM cells (Fig. 4b-d). We next examined whether Spint1a deficiency in the stroma, i.e. in a non-autonomous manner, also promoted SKCM aggressiveness. Spint1a wild type SKCMs were transplanted into the yolk sac of Spint1a-deficient and their wild type siblings larvae (Fig. 4e). Strikingly, it was found that Spint1a deficiency in the tumor microenvironment also promoted a significantly higher dissemination of SKCM compared to control tumor microenvironments (Fig. 4f-h).
To further confirm a role of Spint1a in both SKCM and tumor microenvironment cells, we next sorted tumor (eGFP+) and stromal (eGFP−) cells from both genotypes and then mixed in equal proportions (~ 90% of tumor and ~ 10% of stromal cells) in the 4 possible combinations (Fig. 5a), since it was found that all tumors had ~ 90% of tumor and ~ 10% of stromal cells (data not shown). Notably, both Spint1a-deficient tumor and stromal cells were able to increase SKCM cell invasion (Fig. 5b and c). Collectively, these results suggest that Spint1a deficiency enhances SKCM invasion by both cell autonomous and non-autonomous mechanisms.
Spint1a-deficient SKCM cells showed enhanced aggressiveness in adult zebrafish allotransplantation model
The results obtained in allotransplantation assay in larvae prompted us to analyze the role of Spint1a in SKCM aggressiveness and metastasis in adult casper zebrafish to directly visualize tumor cell proliferation and dissemination over time. spint1ahi2217Tg/hi2217Tg and spint1a+/hi2217Tg SKCMs were sampled, disaggregated and subcutaneously injected (300,000 cells) in the dorsal sinus of adult casper recipients previously irradiated with 30 Gy (Fig. 6a). Tumor engraftment was visible as early as 7 days post-transplantation in both genotypes. While 90% engraftment was obtained with wild type SKCM cells, Spint1a-deficient cells showed a significant enhancement of tumor engraftment rate, around 95% (Fig. 6b). In addition, adult zebrafish recipients transplanted with Spint1a-deficient SKCMs developed tumors with a significant higher growth rate than those injected with wild type SKCMs (Fig. 6c). Notably, Spint1a-deficient SKCM cells were able to invade the entire dorsal area, part of ventral cavity and the dorsal fin (Fig. 6c).
We next performed additional transplant assays following the same work-flow but injecting an increased number of cells (500,000 cells per recipient fish), that ensured a 100% of engraftment was for both genotypes (data not shown). From the first week of analysis, Spint1a-deficient SKCM tumor size was significantly larger than their control counterparts (Fig. 6d). In addition, the recipients injected with Spint1a-deficient SKCM cells developed larger tumors spanning the entire dorsal area and even exceed the notochord line and grew vertically, a clear aggressiveness signature of SKCM (Fig. 6d).
To further investigate the aggressiveness potential of Spint1a-deficient SKCMs, a serial dilution assay was performed following the work-flow previously described in Additional file 1: Figure S1. Cells from both Spint1a-deficient and wild type SKCMs were serially diluted and 3 different numbers of cells (30,000 cells, 100,000 cells and 300,000 cells) were transplanted into recipients as described above. Notably, while 30,000 and 100,000 Spint1a-deficient SKCM cells were able to engraft and the tumor grew over the time, wild type SKCMs hardly grew (Additional file 1: Figure S1A, S1B). However, injection of 300,000 Spint1a-deficient SKCM cells resulted in large tumor spanning the entire dorsal area and invading part of the ventral cavity (Additional file 1: Figure S1C), confirming previous results. Collectively, all these results confirm that Spint1a deficiency enhances SKCM aggressiveness.
Spint1a deficiency promotes SKCM dedifferentiation and inflammation
To understand the mechanisms involved in the Spint1a-mediated aggressiveness of SKCM, the expression of genes encoding important biomarkers was analyzed by RT-qPCR. The mRNA levels of sox10, tyr, dct and mitfa were lower in Spint1a-deficient SKCMs than in their wild type counterparts (Additional file 1: Figure S2). In addition, while the transcript levels of mmp9 and slug were similar in Spint1a-deficient and wild type SKCM, cdh1 levels significantly decreased in Spint1a-deficient compared to wild type SKCM (Additional file 1: Figure S2).
We next analyzed genes encoding key inflammatory molecules and immune cell markers, including the pro-inflammatory cytokine Il1b, the neutrophil markers Lyz and Mpx, the macrophage marker Mpeg1, and the ISGs B2m, Mxb and Pkz, in Spint1a-deficient and wild type SKCMs (Additional file 1: Figure S2). Curiously, it was found that while il1b, lyz and mpx mRNA levels were not affected by Spint1a deficiency, those of mpeg1 were elevated in Spint1a-deficient SKCMs. Furthermore, the ISGs b2m, mxb and pkz genes showed enhanced mRNA levels in Spint1a-deficient SKCMs. These results point out to altered immune surveillance and tumor cell dedifferentiation promoted by Spint1a-deficiency in SKCM.
The relationships between inflammation and cancer are ambiguous. Although it is estimated that underlying infections and inflammatory responses are linked to 15–20% of all deaths from cancer worldwide , immunosuppression is known to increase the risk for certain tumors . Furthermore, immunotherapy is considered the most promising cancer therapy for the next future . In this study, we have developed a preclinical zebrafish model to study the role of SPINT1-driven skin chronic inflammation in melanoma. We found that Spint1a deficiency is required at both cell autonomous and non-autonomous levels to enhance cell dissemination of SKCM by promoting tumor dedifferentiation and altered immune surveillance. These results may have important clinical impact, since genetic alterations of SPINT1 were found in 10% of SKCM patients and correlated with altered cell cycle, differentiation and innate and adaptive immune signaling pathways and, more importantly, with a poor prognosis. In addition, a positive correlation of both inflammation and macrophage marker with SPINT1 levels was found in tumor samples, fitting the increased number of TAM in SKCM with high SPINT1 levels. However, a negative correlation of SPINT1 and EMT marker levels was observed, suggesting that SPINT1 levels critically regulate tumor and immune cell crosstalk. Therefore, both high and low levels of SPINT1 result in an unbalanced crosstalk between tumor cells and their microenvironment promoting higher SKCM aggressiveness.
Another interesting observation is that SPINT1 transcript levels positively correlated with macrophage infiltration, but not neutrophil one, in SKCM tumor samples. Similarly, our gene expression studies also suggest that Spint1a regulates macrophage infiltration in the zebrafish SKCM model. Curiously, activated neutrophils in a condition of repeated wounding have been shown to interact with pre-neoplastic cells promoting their proliferation through the release of prostaglandin E2 and, more importantly, SKCM ulceration correlates with increased neutrophil infiltration and tumor cell proliferation, which are both associated with poor prognosis . Although we found a robust positive correlation between the transcript levels of SPINT1 and CXCR2, which encodes a major IL-8 receptor involved in SKCM neutrophil infiltration , CXCR2 has also been shown to promote tumor-induced angiogenesis and increased proliferation [37,38,39] and, therefore, the SPINT1/CXCR2 axis may regulate SKCM aggressiveness by neutrophil-independent pathways.
The zebrafish model developed in this study uncovered a role for Spint1a in facilitating oncogenic transformation which probably accelerates the SKCM onset. Curiously, Spint1a deficiency seems to have a prominent role in oncogenic transformation accelerating SKCM onset in vivo. However, the strong oncogenic activity of HRAS-G12V, which is even able to induce melanoma without the need of coadjuvating mutations in tumor suppressors, results in a similar tumor burden in wild type and Spint1a-deficient fish at later stages. It would be interesting, therefore, to analyze in future studies the impact of Spint1a in SKCM models driven by other major oncogenes, such as Braf, Kras, Nras and Nf1. Although SPINT1 is a serine protease inhibitor with several targets, including ST14 and HGFA, deregulation of the SPINT1/ST14 axis leads to spontaneous squamous cell carcinoma in mice  and keratinocyte hyperproliferation in zebrafish [19, 20] preceded by skin inflammation in both models. In addition, intestine-specific Spint1 deletion in mice induces the activation of the master inflammation transcription factor NF-κB and accelerated intestinal tumor formation . Strikingly, pharmacological inhibition of NF-κB activation reduces the formation of intestinal tumors in Spint1-deficient ApcMin/+ mice , unequivocally demonstrating that Spint1-driven inflammation promotes tumorigenesis.
The SKCM allotransplant assays in larvae revealed for the first time that Spint1a deficiency in both tumor and stromal cells increases SKCM invasiveness. In addition, Spint1a deficiency in both cell types does not show enhanced invasiveness compared to Spint1a deficiency in either cell type. This is an interesting observation, since SPINT1 is a membrane-bound protein that may, therefore, inhibit their targets in both tumor cell autonomous and non-autonomous manners. However, wild type Spint1a tumor microenvironment fails to compensate its loss in tumor cells, since Spint1a-deficient tumor cells show enhanced invasiveness in wild type larvae and adult recipients, and vice versa. Importantly, transplantation experiments of serial diluted SKCM cells revealed the crucial cell-autonomous role of Spint1a in inhibiting tumor aggressiveness. Similarly, loss of SPINT1 in human pancreatic cancer cells promotes ST14-dependent metastasis in nude mouse orthotopic xenograft models . We observed that genetic alterations in SPINT1 transcript levels in SKCM patient samples negatively correlated with EMT markers and that Spint1a-deficient zebrafish SKCM showed reduced cdh1 mRNA levels. EMT phenotype switching has been shown to be involved in acquisition of metastatic properties in the vertical growth phase of SKCM  and loss of E-cadherin, with gain of N-cadherin and osteonectin, was associated with SKCM metastasis . Importantly, the presence of aberrant E-cadherin expression in primary and metastatic SKCM is associated with a poor overall patient survival . Therefore, our results suggest that SPINT1 loss may facilitate metastatic invasion of human SKCM through EMT phenotype switching.
We have developed a preclinical model to study the role of altered expression of SPINT1 in early transformation, progression and metastatic invasion in SKCM. This model has revealed that Spint1a deficiency facilitates oncogenic transformation, regulates the tumor/immune microenvironment crosstalk and is associated to SKCM aggressiveness. In addition, Spint1a deficiency in either tumor or microenvironment compartment increases SKCM aggressiveness. The high prevalence of SPINT1 genetic alterations in SKCM patients and their association with a poor prognosis, point out to SPINT1 as a novel therapeutic target for the treatment of SKCM.
Availability of data and materials
Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Victoriano Mulero (email@example.com).
Analyzed by analysis of variance
Enhanced green fluorescent protein
Fetal bovine serum
Hepatocyte growth factor activator inhibitor 1
Hepatocyte growth factor activator
Skin cutaneous melanoma
Serine protease inhibitor, kunitz-type, 1
Suppression of tumorigenicity 14
The Cancer Genome Atlas
Wellbrock C, Arozarena I. The complexity of the ERK/MAP-kinase pathway and the treatment of melanoma skin cancer. Front Cell Dev Biol. 2016;4:33.
Schadendorf D, Hauschild A. Melanoma in 2013: melanoma--the run of success continues. Nat Rev Clin Oncol. 2014;11(2):75–6.
van Rooijen E, Fazio M, Zon LI. From fish bowl to bedside: the power of zebrafish to unravel melanoma pathogenesis and discover new therapeutics. Pigment Cell Melanoma Res. 2017;30(4):402–12.
NIH. Cancer stat facts: melanoma of the skin: National Cancer Institute; 2019. https://seer.cancer.gov/statfacts/html/melan.html
Schadendorf D, Fisher DE, Garbe C, Gershenwald JE, Grob JJ, Halpern A, et al. Melanoma. Nat Rev Dis Primers. 2015;1:15003.
Hoerter JD, Bradley P, Casillas A, Chambers D, Weiswasser B, Clements L, et al. Does melanoma begin in a melanocyte stem cell? J Skin Cancer. 2012;2012:571087.
Kaufman CK, Mosimann C, Fan ZP, Yang S, Thomas AJ, Ablain J, et al. A zebrafish melanoma model reveals emergence of neural crest identity during melanoma initiation. Science. 2016;351(6272):aad2197.
Santoriello C, Gennaro E, Anelli V, Distel M, Kelly A, Koster RW, et al. Kita driven expression of oncogenic HRAS leads to early onset and highly penetrant melanoma in zebrafish. PLoS One. 2010;5(12):e15170.
Nguyen N, Couts KL, Luo Y, Fujita M. Understanding melanoma stem cells. Melanoma Manag. 2015;2(2):179–88.
Shalapour S, Karin M. Immunity, inflammation, and cancer: an eternal fight between good and evil. J Clin Invest. 2015;125(9):3347–55.
Maru GB, Gandhi K, Ramchandani A, Kumar G. The role of inflammation in skin cancer. Adv Exp Med Biol. 2014;816:437–69.
Tang L, Wang K. Chronic inflammation in skin malignancies. J Mol Signal. 2016;11:2.
Lin CY, Anders J, Johnson M, Sang QA, Dickson RB. Molecular cloning of cDNA for matriptase, a matrix-degrading serine protease with trypsin-like activity. J Biol Chem. 1999;274(26):18231–6.
Benaud C, Dickson RB, Lin CY. Regulation of the activity of matriptase on epithelial cell surfaces by a blood-derived factor. Eur J Biochem. 2001;268(5):1439–47.
Tseng IC, Chou FP, Su SF, Oberst M, Madayiputhiya N, Lee MS, et al. Purification from human milk of matriptase complexes with secreted serpins: mechanism for inhibition of matriptase other than HAI-1. Am J Physiol Cell Physiol. 2008;295(2):C423–31.
Shimomura T, Denda K, Kitamura A, Kawaguchi T, Kito M, Kondo J, et al. Hepatocyte growth factor activator inhibitor, a novel Kunitz-type serine protease inhibitor. J Biol Chem. 1997;272(10):6370–6.
List K, Szabo R, Molinolo A, Sriuranpong V, Redeye V, Murdock T, et al. Deregulated matriptase causes ras-independent multistage carcinogenesis and promotes ras-mediated malignant transformation. Genes Dev. 2005;19(16):1934–50.
List K. Matriptase: a culprit in cancer? Future Oncol. 2009;5(1):97–104.
Mathias JR, Dodd ME, Walters KB, Rhodes J, Kanki JP, Look AT, et al. Live imaging of chronic inflammation caused by mutation of zebrafish Hai1. J Cell Sci. 2007;120(Pt 19):3372–83.
Carney TJ, von der Hardt S, Sonntag C, Amsterdam A, Topczewski J, Hopkins N, et al. Inactivation of serine protease Matriptase1a by its inhibitor Hai1 is required for epithelial integrity of the zebrafish epidermis. Development. 2007;134(19):3461–71.
Westerfield M. The zebrafish book. A guide for the laboratory use of zebrafish Danio* (Brachydanio) rerio. Eugene: University of Oregon Press; 2000.
White RM, Sessa A, Burke C, Bowman T, LeBlanc J, Ceol C, et al. Transparent adult zebrafish as a tool for in vivo transplantation analysis. Cell Stem Cell. 2008;2(2):183–9.
Dang M, Henderson RE, Garraway LA, Zon LI. Long-term drug administration in the adult zebrafish using oral gavage for cancer preclinical studies. Dis Model Mech. 2016;9(7):811–20.
Marques IJ, Weiss FU, Vlecken DH, Nitsche C, Bakkers J, Lagendijk AK, et al. Metastatic behaviour of primary human tumours in a zebrafish xenotransplantation model. BMC Cancer. 2009;9:128.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.
Coussens LM, Werb Z. Inflammation and cancer. Nature. 2002;420(6917):860–7.
Colegio OR, Chu NQ, Szabo AL, Chu T, Rhebergen AM, Jairam V, et al. Functional polarization of tumour-associated macrophages by tumour-derived lactic acid. Nature. 2014;513(7519):559–63.
Gordon SR, Maute RL, Dulken BW, Hutter G, George BM, McCracken MN, et al. PD-1 expression by tumour-associated macrophages inhibits phagocytosis and tumour immunity. Nature. 2017;545(7655):495–9.
Zaidi MR, Davis S, Noonan FP, Graff-Cherry C, Hawley TS, Walker RL, et al. Interferon-gamma links ultraviolet radiation to melanomagenesis in mice. Nature. 2011;469(7331):548–53.
Ronnstrand L, Phung B. Enhanced SOX10 and KIT expression in cutaneous melanoma. Med Oncol. 2013;30(3):648.
Ordonez NG. Value of melanocytic-associated immunohistochemical markers in the diagnosis of malignant melanoma: a review and update. Hum Pathol. 2014;45(2):191–205.
Wang H, Yang L, Wang D, Zhang Q, Zhang L. Pro-tumor activities of macrophages in the progression of melanoma. Hum Vaccin Immunother. 2017;13(7):1556–62.
Feng Y, Santoriello C, Mione M, Hurlstone A, Martin P. Live imaging of innate immune cell sensing of transformed cells in zebrafish larvae: parallels between tumor initiation and wound inflammation. PLoS Biol. 2010;8(12):e1000562.
Carreau N, Pavlick A. Revolutionizing treatment of advanced melanoma with immunotherapy. Surg Oncol. 2019. https://doi.org/10.1016/j.suronc.2019.01.002.
Antonio N, Bonnelykke-Behrndtz ML, Ward LC, Collin J, Christensen IJ, Steiniche T, et al. The wound inflammatory response exacerbates growth of pre-neoplastic cells and progression to cancer. EMBO J. 2015;34(17):2219–36.
Jablonska J, Wu CF, Andzinski L, Leschner S, Weiss S. CXCR2-mediated tumor-associated neutrophil recruitment is regulated by IFN-beta. Int J Cancer. 2014;134(6):1346–58.
Singh S, Sadanandam A, Nannuru KC, Varney ML, Mayer-Ezell R, Bond R, et al. Small-molecule antagonists for CXCR2 and CXCR1 inhibit human melanoma growth by decreasing tumor cell proliferation, survival, and angiogenesis. Clin Cancer Res. 2009;15(7):2380–6.
Gabellini C, Gomez-Abenza E, Ibanez-Molero S, Tupone MG, Perez-Oliva AB, de Oliveira S, et al. Interleukin 8 mediates bcl-xL-induced enhancement of human melanoma cell dissemination and angiogenesis in a zebrafish xenograft model. Int J Cancer. 2018;142(3):584–96.
Gabellini C, Trisciuoglio D, Desideri M, Candiloro A, Ragazzoni Y, Orlandi A, et al. Functional activity of CXCL8 receptors, CXCR1 and CXCR2, on human malignant melanoma progression. Eur J Cancer. 2009;45(14):2618–27.
Kawaguchi M, Yamamoto K, Kanemaru A, Tanaka H, Umezawa K, Fukushima T, et al. Inhibition of nuclear factor-kappaB signaling suppresses Spint1-deletion-induced tumor susceptibility in the ApcMin/+ model. Oncotarget. 2016;7(42):68614–22.
Ye J, Kawaguchi M, Haruyama Y, Kanemaru A, Fukushima T, Yamamoto K, et al. Loss of hepatocyte growth factor activator inhibitor type 1 participates in metastatic spreading of human pancreatic cancer cells in a mouse orthotopic transplantation model. Cancer Sci. 2014;105(1):44–51.
Bennett DC. How to make a melanoma: what do we know of the primary clonal events? Pigment Cell Melanoma Res. 2008;21(1):27–38.
Alonso SR, Tracey L, Ortiz P, Perez-Gomez B, Palacios J, Pollan M, et al. A high-throughput study in melanoma identifies epithelial-mesenchymal transition as a major determinant of metastasis. Cancer Res. 2007;67(7):3450–60.
Yan S, Holderness BM, Li Z, Seidel GD, Gui J, Fisher JL, et al. Epithelial-mesenchymal expression phenotype of primary melanoma and matched metastases and relationship with overall survival. Anticancer Res. 2016;36(12):6449–56.
We strongly acknowledge P. Martínez for his excellent technical assistance with zebrafish maintenance.
This work was supported by the Spanish Ministry of Science, Innovation and Universities (grants BIO2014–52655-R and BIO2017–84702-R to VM and PI13/0234 to MLC), all co-funded with Fondos Europeos de Desarrollo Regional/European Regional Development Funds), Fundación Séneca-Murcia (grant 19400/PI/14 to MLC), the University of Murcia (postdoctoral contracts to DGM), and the European Union Seventh Framework Programme-Marie Curie COFUND (FP7/2007–2013) under UMU Incoming Mobility Programme ACTion (U-IMPACT) Grant Agreement 267143. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The experiments complied with the Guidelines of the European Union Council (Directive 2010/63/EU) and the Spanish RD 53/2013. Experiments and procedures were performed as approved by the Consejería de Agua, Agricultura, Ganadería y Pesca de la CARM (authorization number # A13180602).
Consent for publication
L.I.Z. is a founder and stockholder of Fate Therapeutics, Inc., Scholar Rock and Camp4 Therapeutics.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Spint1a-deficient SKCM shows enhanced aggressiveness in adult zebrafish allotransplantation assays. Control and Spint1a deficient SKCMs were disaggregated and 30,000 (a), 100,000 (b) and 300,000 cells (c) were injected subcutaneously in pre-irradiated adult casper zebrafish. Fish were analyzed for average tumor size (pixels) from 1 to 4 weeks post-transplant. Representative images and quantification of the average tumor size are shown. Each dot corresponds to a recipient-transplanted fish and the mean ± SEM is also shown (n=2 SKCM tumors).*p<0.05, ***p<0.001 according to unpaired Student t test. Figure S2. Expression analysis of differentiation melanocyte, EMT, inflammation and immune cell markers in zebrafish SKCM. The mRNA levels of the genes encoding the differentiation melanocyte markers sox10, mitfa, tyr and dct, the EMT markers cdh1, slug and mmp9, the inflammation marker il1b, the neutrophil markers lyz and mpx, the macrophage marker mpeg1, and the ISGs b2m, mxb and pkz were analyzed by RT-qPCR in control and Spint1a-deficient SKCMs. *p < 0.05, **p<0.01 according to one-tailed Student t test. Table S1. Primers used in this study for RT-qPCR. The gene symbols followed the Zebrafish Nomenclature Guidelines (http://zfin.org/zf_info/nomen.html/). ENA, European Nucleotide Archive (http://www.ebi.ac.uk/ena/). (PDF 496 kb)
About this article
Cite this article
Gómez-Abenza, E., Ibáñez-Molero, S., García-Moreno, D. et al. Zebrafish modeling reveals that SPINT1 regulates the aggressiveness of skin cutaneous melanoma and its crosstalk with tumor immune microenvironment. J Exp Clin Cancer Res 38, 405 (2019). https://doi.org/10.1186/s13046-019-1389-3