Proteomics analysis of melanoma metastases: association between S100A13 expression and chemotherapy resistance.

BACKGROUND: Disseminated cutaneous malignant melanoma (CMM) is commonly unresponsive to standard chemotherapies, and there are as yet no predictive markers of therapy response.
METHODS: In the present study we collected fresh-frozen pretreatment lymph-node metastasis samples (n=14) from melanoma patients with differential response to dacarbazine (DTIC) or temozolomide (TMZ) chemotherapy, to identify proteins with an impact on treatment response. We performed quantitative protein profiling using tandem mass spectrometry and compared the proteome differences between responders (R) and non-responders (NR), matched for age, gender and histopathological type of CMM.
RESULTS: Biological pathway analyses showed several signalling pathways differing between R vs NR, including Rho signalling. Gene expression profiling data was available for a subset of the samples, and the results were compared with the proteomics data. Four proteins with differential expression between R and NR were selected for technical validation by immunoblotting (ISYNA1, F13A1, CSTB and S100A13), and CSTB and S100A13 were further validated on a larger sample set by immunohistochemistry (n=48). The calcium binding protein S100A13 was found to be significantly overexpressed in NR compared with R in all analyses performed.
CONCLUSIONS: Our results suggest that S100A13 is involved in CMM resistance to DTIC/TMZ.

Cutaneous malignant melanoma (CMM) is one of the most aggressive skin cancer types, with a rapid increase in prevalence, especially among western countries (Erdmann ). So far, curative treatment is achieved only by surgical resection of primary tumours at an early stage (Balch ; Balch ). Chemotherapy with the alkylating agents dacarbazine (DTIC) and temozolomide (TMZ) has not shown any effect on median overall survival of patients (Hill ; Gogas ). However, ∼13% of the patients have an initial partial or complete response to these drugs and rare long-time responses are seen (Middleton ; Kim ; Patel ). Several chemoresistance mechanisms have been proposed, but the clinically relevant mechanisms are to a large extent still unknown.

DNA repair proteins, such as O-6-methylguanine-DNA methyltransferase (MGMT), which counteract the alkylating effects of DTIC/TMZ, have been suggested to confer resistance to these drugs (for a review, refer to the study by Zhang ). The relationship between MGMT and response to therapy has mainly been reported for TMZ treatment in glioma, but association between MGMT expression and resistance to DTIC/TMZ has been found also in melanoma (Ma ; Busch ). In addition, studies using melanoma cell lines have suggested a role for melanosomes and related protein-trafficking pathways in melanoma chemoresistance (Chen ; Xie ), and we have shown higher expression of the key melanogenesis regulator microphthalmia-associated transcription factor (MITF) and the melanosomal G-protein-coupled receptor 143 (GPR143) in tumours from non-responders to DTIC/TMZ therapy, compared with responders (Hertzman Johansson ).

Recently, the treatment options for melanoma have expanded rapidly, and several novel immunotherapies and small-molecule drugs targeting the MAPK pathway have been approved or are in clinical trials (Hodi ; Chapman ; Hauschild ; Flaherty , 2012b; Menzies and Long, 2013). Unfortunately, primary resistance is widely encountered for immunotherapies, and despite many impressive initial responses to MAPK signalling inhibitors, acquired drug resistance eventually develops. Thus, chemotherapy could have a role in combination therapies by inhibiting the growth of tumour cells that are resistant to specific inhibitors. In addition, standard chemotherapy is still relevant for patients not responding to immunotherapy, and for those whose tumours lack activated druggable targets.

To identify proteins associated with chemotherapy resistance, we carried out protein expression profiling of tumour biopsies from patients with differential response to DTIC/TMZ treatment. Gene expression profiling data was available for a majority of the tumours (Hertzman Johansson ), and several candidate genes and biological pathways were identified by both methods. Four of the top candidates among the identified proteins were validated in the same tumour samples by immunoblotting. In addition, two of the candidate proteins with higher expression in non-responders, the calcium binding protein S100A13 and the cathepsin inhibitor cystatin B (CSTB), were further studied by immunohistochemistry in an extended set of pre-treatment tumour samples. The significantly higher expression of S100A13 in non-responders was confirmed also in this sample set.

Materials and methods

Tumour samples

Fresh-frozen pre-treatment lymph-node biopsies were excised from patients with metastatic cutaneous melanoma before chemotherapy with DTIC/TMZ. All tumour pieces included had a tumour cell burden of at least 50% (the majority of tumours >70%). Although there were individual differences between the tumour pieces with regard to cell composition, vascularisation and infiltration of leukocytes, there were no major differences between the groups of responders and non-responders (see Supplementary Figure 1 for histological images).

Five patients were responders (defined as a reduction of tumour size of at least 50%), and five non-responders (with stable or progressive disease) were matched to the responders with regard to gender, histopathology (superficial spreading or nodular melanoma) and, when possible, age. Four of the responders had previously been included in a gene expression profiling study (Hertzman Johansson ) and four non-responders from this study material were also included in the present analysis, in addition to the five matched non-responders. For details on the clinical material included in the proteomics analysis, see Supplementary Table 1 and Supplementary Figure 1.

For the immunohistochemistry, an extended material of formalin-fixed paraffin-embedded pre-treatment biopsies from 16 DTIC/TMZ responders and 34 non-responders were analysed.

Patients have given informed consent and the study was approved by the Regional Ethics Committee of Stockholm.

Protein sample preparation

Tissue disruption was performed using a Mixer mill MM200 (Retsch, Haan, Germany) with a 1-cm ø Teflon-coated tungsten ball. Teflon vials, caps and the ball were first pre-cooled in −80 °C and then in liquid nitrogen. Frozen tumour samples (∼3 × 3 × 3 mm) were cut into smaller pieces (1 × 1 × 1 mm) on a block of dry ice. Samples were then transferred to the pre-cooled vials and cooled again for 2 min in liquid nitrogen. The samples were then homogenised in the Mixer mill for 2 min at maximum speed after which they were put back into the liquid nitrogen. This homogenisation/cooling procedure was repeated three times until a frozen tumour powder was achieved. The powder was dissolved in 1 ml 20 mM HEPES (pH 7.6), transferred to Eppendorf tubes and kept in −80 °C.

The samples were thawed and DTT and SDS were added to an end concentration of 1 mM and 3.75%, respectively. Samples were then lysed by heating (at 90 °C for 5 min on a thermomixer) and sonication (at room temperature for 5 min). The samples were then centrifuged at 10 000 g for 10 min and the supernatant was transferred to a new tube. Protein concentrations of the lysates were measured using the DC protein assay kit 2 (BioRad, Hercules, CA, USA). The protein yield was at least 1.6 mg per tumour piece.

For acetone precipitation, 120 mg of protein from each sample was mixed with four volumes of ice-cold acetone. The samples were incubated in +4 °C for 2 hours until a flocculent formed. The samples were then centrifuged for 10 min at 10 000 g. The supernatant was discarded and the protein was allowed to air dry. This step was performed for the MS/MS analysis, but not for immunoblotting samples.

Digestion, iTRAQ labelling and clean-up of protein samples

Samples were digested and labelled according to a standard iTRAQ eight-plex protocol (Thermo Fisher Scientific, Waltham, MA, USA). The trypsin-to-protein ratio was 1 : 20. One hundred micrograms of each sample were labelled and two pools were made each with seven samples and one internal standard for normalisation.

Pooled iTRAQ-labelled digests were applied to 1 ml Strata X-C 33 μm polymeric strong cation exchange (SCX) microcolumns (Phenomenex, Torrance, CA, USA). The microcolumns were initially washed with 1 ml 100% methanol followed by 1 ml MilliQ grade water (EMD Millipore Corporation, Billerica, MA, USA). The sample was adjusted to 500 μl 0.1% formic acid and then applied to the columns. After washing with 1 ml 30% methanol and 0.1% formic acid the samples were eluted with 30% methanol and 5% ammonium hydroxide. Samples were then dried in a SpeedVac system (EMD Millipore Corporation).

Narrow-range IEF

Briefly, samples were added to pH 3.3 loading strips, kindly provided by GE healthcare (Uppsala, Sweden), and allowed to re-swell overnight. Loading strips were then put on re-swelled 24-cm narrow-range isoelectric focusing strips (pH 3.4–4.7) also provided by GE Healthcare. Samples were focused until 100 kVhs was reached. The peptides were passively eluted in 72 fractions of 150 μl milliQ water using an in-house robot. Eluted peptide samples were then dried in a SpeedVac system.

Mass spectrometry

The 72 fractions were re-suspended in 10 μl of 3% acetonitrile in 0.1% formic acid. From each IPG fraction 3 μl was injected into online HPLC/MS performed on a hybrid LTQ-Orbitrap Velos mass spectrometer (Thermo Fischer Scientific, San Jose, CA, USA). An Agilent HPLC 1200 system (Agilent Technologies, Santa Clara, CA, USA) was used to provide the gradient for online reversed-phase nano-LC at a flow of 0.4 μl min−1. Solvent A was 97% water, 3% ACN and 0.1% formic acid, and solvent B was 5% water, 95% ACN and 0.1% formic acid. The curved gradient went from 2% B up to 40% B in 45 min, followed by a steep increase to 100% B in 5 min. The sample was injected into a C18 guard-desalting column (Agilent Technologies) prior to a 15 cm long C18 picofrit column (100 μm internal diameter, 5 μm bead size, Nikkyo Technos Co., Tokyo, Japan) installed on to the nano electrospray ionisation (NSI) source of the Orbitrap Velos instrument (Thermo Fisher, Waltham, MA, USA). Acquisition proceeded in ∼3.5-s scan cycles, starting by a single full-scan MS at 30 000 resolution (profile mode), followed by two stages of data-dependent tandem MS (centroid mode): the top five ions from the full scan MS were selected first for collision-induced dissociation (CID, at 35% energy) with MS/MS detection in the ion trap, and finally for higher-energy collision dissociation (HCD, at 45% energy) with MS/MS detection in the orbitrap. Precursors were isolated with a 2 m/z width and dynamic exclusion was used with 60 s duration.

Database search

Orbitrap data was searched by Mascot 2.2 (Matrix Science Limited, London, UK) under the software platform Proteome Discoverer 1.1 (Thermo) against the human Uniref100 protein sequence database. A precursor mass tolerance of 10 p.p.m., and product mass tolerances of 0.02 Da for HCD-FTMS and 0.8 Da for CID-ITMS were used. Further settings used were the following: trypsin with 1 missed cleavage; carbamidomethylation on cysteine and iTRAQ-8plex on lysine and N-terminal as fixed modifications; and oxidation of methionine as variable modification. Quantification of iTRAQ-8plex reporter ions was performed by using Proteome Discoverer on HCD-FTMS tandem mass spectra using an integration window tolerance of 20 p.p.m. In total, 5873 proteins were identified with ⩾1 peptide 95% confidence. Out of these 3029 proteins were significantly (95% confidence) detected in both of the sample pools. The downstream data analysis was performed on these 3029 overlapping proteins (for details, see Supplementary Table 2).

Immunoblotting

For technical validation of selected biomarker candidates, the same tumour protein extracts as in the protein profiling were analysed by immunoblotting, using NuPAGE Novex Bis-Tris Gel (Life Technologies, Carlsbad, CA, USA) and PVDF membranes (Thermo Scientific, Rockford, IL, USA), according to the manufacturer's standard protocol and the following primary antibodies from Sigma Aldrich Sweden AB (Stockholm, Sweden): ISYNA1 produced in rabbit (1 : 500; catalogue no. AV53716); CSTB produced in mouse (1 : 1000. catalogue no. C5243); F13A1 antibody produced in rabbit (1 : 250; catalogue no. HPA001804); S100A13 produced in rabbit (1 : 250. catalogue no. HPA019592); and β-actin produced in mouse (1 : 5000; catalogue no. A544). The primary antibody for β-tubulin (sc-9104 from Santa Cruz Biotechnology, Inc., Santa Cruz, CA,, USA) was diluted 1 : 200. The HRP-conjugated secondary antibodies used were bovine anti-rabbit antibody (1 : 1000; Santa Cruz Biotechnology catalogue no. SC-2379) and goat anti-mouse antibody (1 : 2000; Life Technologies catalogue no. 626520). The Pierce ECL Western Blotting Substrate (Thermo Scientific) was used for detection. Quantification of the protein expression was performed using Image J (http://rsbweb.nih.gov/ij/index.html)

Immunohistochemistry

Immunohistochemistry was performed on formalin-fixed, paraffin-embedded 3–4 μm thick tumour sample sections. Antigen retrieval was induced by heating the sections in humid decloaking heat chambers (Biocare, Concord, CA, USA) in citrate buffer (pH 6.0) according to the manufacturer's instruction followed by 10 min incubation with 3% hydrogen peroxide at room temperature and washing with 1 × TBS buffer. To avoid unspecific binding, the sections were incubated with 2.5% horse serum for 20 min at room temperature followed by overnight incubation at 4 °C with primary antibodies against CSTB (1 : 1000 Sigma Aldrich Sweden AB catalogue no. C5243) or S100A13 (1 : 750 Sigma Aldrich Sweden AB catalogue no. HPA019592) in TBS buffer containing 1.5% horse serum. Negative controls were incubated with the same TBS buffer without any primary antibody.

Incubation with the secondary antibody plus streptavidin/peroxidase was performed according to the manufacturer's instructions (Vectastain Universal Quick Kit, Vector Laboratories Inc., Burlingame, CA, USA) and were developed with 3,3′-diaminobenzidine (DAB Kit, Vector Laboratories Inc.). As the final step, slides are rinsed in water and counterstained with Mayer‘s haematoxylin, rinsed with water, dehydrated and sealed with glass lamella and mounting solution Mountex (Histolab, Spånga, Sweden).

Independent assessment of all slides was performed by at least three observers (AA, CHJ, MFS or SEB). The intensity was scored 0–3 (0=no staining, 1=weak, 2=moderate, 3=strong staining) and samples where >30% of the tumour cells had an intensity score of at least 2 were considered to have positive staining.

Immunohistochemistry for ISYNA1 and F13A1 was not completed due to unspecific staining for the antibodies available.

Statistical analysis

Quality assessment of the LC/MS/MS data was done using principal component analysis (PCA) to check for outliers (SIMCA statistical analysis software, Umetrics, Umeå, Sweden). Multivariate orthogonal-partial least squares (O-PLS) analysis was performed, also using SIMCA. For optimisation of O-PLS models we used VIP (Variable Importance in Projection) value to judge protein importance in the model. The O-PLS models were validated by sevenfold cross validation. Proteins with significant VIP throughout the cross validation of the model were selected for the optimised model. CV-ANOVA was used to judge the model validity (Eriksson ). Biological pathways analyses were performed using the Ingenuity IPA software (Ingenuity Systems, Qiagen, Redwood City, CA, USA).

For the immunohistochemistry results, Fisher's exact test was used to analyse the difference in S100A13 and CSTB protein expression between responders and non-responders.

Results

Fourteen fresh-frozen melanoma lymph-node metastasis samples (see Supplementary Table 1), from five chemotherapy responders and nine non-responders, were profiled by liquid chromatography tandem mass spectrometry (LC/MS-MS). In total, 3029 proteins were detected in all samples, and the downstream data analyses were performed on this data set. Detailed protein data are provided in Supplementary Table 2. To check for outliers, unsupervised assessment was performed using PCA and no significant outlier was identified. Multivariate O-PLS was performed to create models that could separate the responders and non-responders.

A significant O-PLS model consisting of ninety-four proteins was created when comparing pretreatment tumours from five responders and five non-responders that were matched for age, gender and histopathological type (CV-ANOVA, P=0.0007, Q2=0.874, R2Y=0.906 and R2X=0.619, see Supplementary Table 3). The enriched top molecular and cellular functions and canonical pathways of this selected group of 94 proteins was defined by ingenuity pathway analysis (see Table 1), and include several key cell signalling pathways such as Rho and Rac signalling.

Table 1

Top cellular and molecular functions and canonical pathways identified when comparing the tumour proteome profile of pre-treatment biopsies from five DTIC/TMZ responders and five matched non-responders

Top cellular and molecular functions	P-value
RNA post-transcriptional modification	1.10E−05–3.29E−02
Cell morphology	2.22E−04–3.75E−02
Cell-to-cell signalling and interaction	2.22E−04–3.75E−02
Carbohydrate metabolism	3.32E−04–3.75E−02
Cellular movement	3.47E−04–3.75E−02
Top canonical pathways
Signalling by Rho Family GTPases	0.000108
Actin nucleation by ARP–WASP complexes	0.000145
Rac signalling	0.00147
fMLP signalling in neutrophils	0.00175
Clathrin-mediated endocytosis signalling	0.00189

For four of the responders in the present study, gene expression profiling data from the same tumours were available (Hertzman Johansson ). The proteomic results from these samples and four other (unmatched) non-responders, which were also included in the gene expression profiling study, were analysed separately with the purpose of comparing the differences on protein level to the transcriptome data (see Supplementary Table 1 for details regarding the samples included). An O-PLS model consisting of 146 proteins was created (CV-ANOVA P=0.006, Q2=0.998, R2Y=1 and R2X=0,783, see Supplementary Table 3). Compared with the first O-PLS model, with the matched non-responders, 11 proteins overlapped and 76 canonical pathways were identified in both models (see Supplementary Table 3).

When comparing the proteomics results to the transcriptome data for the samples where this was available, 13 candidates were identified by both methods, 6 with concordant expression levels (see Table 2 for details). We could also detect 16 overlapping biological pathways in the gene expression profiling and the corresponding proteomics study (see Supplementary Table 4). Four of these pathways were immune system-related, and the remaining 12 were signalling cascades, including actin cytoskeleton signalling, RhoA signalling, PPAR signalling and IL-8 signalling.

Table 2

Candidates identified as differentially expressed in DTIC/TMZ responders and non-responders by both proteomics and gene expression profiling

Protein accession	Symbol	Protein name	Proteomics GE profiled samples	Gene expression profiling
Q9Y2Q5	ROBLD3	Roadblock domain containing 3	−1.2	0.8
P59998	ARPC4	Actin-related protein 2/3 complex, subunit 4	−0.8	0.4
P07203	GPX1	Glutathione peroxidase	−0.6	0.9
Q2NLC8	GSTK1	Glutathione S−transferase kappa 1	−0.6	0.7
P04080	CSTB	Cystatin B (stefin B)	−0.5	−1
Q6FG43	FLOT2	Flotilin 2	−0.5	0.7
B4DN37	DEK	DEK oncogene	−0.4	0.5
P62136	PPP1CA	Protein phosphatase 1, catalytic subunit, alpha isoenzyme	−0.4	1.4
B4DJ89	POLR2E	Polymerase (RNA) II (DNA-directed) polypeptide E	0.2	0.9
A2RRE5	GRLF1	Glucocorticoid receptor DNA binding factor 1	0.4	0.7
Q96G21	IMP4	IMP4, U3 small nucleolar ribonucleoprotein	0.5	0.6
B4E1U1	PITPNA	Phosphatidylinositol transfer protein, alpha	0.6	0.9
B7Z3K3	ISYNA	Inositol-3-phosphate synthase 1	1.4	1.2

The results are presented as log2 values of the fold change protein or mRNA expression responders/non-responders.

The proteins in the two O-PLS models were also compared with known targets of the transcription factor MITF (Hoek ; Strub ), a potential chemoresistance candidate (Hertzman Johansson ). About 40% of the proteins identified in the present study are previously reported to be regulated by MITF (see Supplementary Table 3 for details).

Four of the identified chemotherapy response candidates, CSTB, F13A1, ISYNA1 and S100A13, were selected for technical validation by immunoblotting, using the same protein extract as in the MS analyses. The selection was based on the strength of the proteomics and gene expression profiling results, as well as the biological function of the protein. Among the proteins that were investigated further, F13A1 and S100A13 were among the top 10 most upregulated in non-responders in the proteomics analysis of the matched samples (see Supplementary Table 3). F13A1 activation is the last event in the blood coagulation cascade and high abundance of the protein is expected in an advanced tumour with the hallmark of blunt and leaky vessels. Presence of high levels of this coagulation factor in the tumours can also lead to more blood clotting events in the tumour vessels and consequently, less efficiency in drug penetration in inner layers of the metastasis. S100A13 is suggested to have a role in melanoma progression and angiogenesis (Massi ), and calcium-binding proteins, including S100A13, have been reported to be elevated in cisplatin-resistant melanoma cell lines in a recent proteomics study (Paulitschke ).

The two other selected candidates, CSTB and ISYNA1, had been identified by both proteome and gene expression profiling and showed similar results in both analyses (see Table 2). The three selected candidates with higher expression in resistant tumours, CSTB, F13A1 and S100A13, have all been reported to be targets of MITF (Hoek ; Strub ).

All four proteins analysed showed similar results in the proteomics and immunoblotting analyses (see Figure 1), thereby confirming the quantitative proteomics analysis.

Figure 1

Immunoblot results. (A) Comparison of proteomics and immunobloting results. Data are shown as the log2 of the mean of the fold change in responders over the value in non-responders. (B–D) Immunoblot images for four candidate proteins.

To evaluate the clinical relevance of the identified chemoresistance candidates we performed immunohistochemistry on an extended set of tumours from 16 responders and 32–34 non-responders. The antibodies for ISYNA1 and F13A1 were tested, but did not work well for immunohistochemistry and were not included in this analysis. Representative images of CSTB and S100A13 staining are shown in Figure 2.

Figure 2

Immunohistochemistry. Representative images of metastases from DTIC/TMZ responders and non-responders, stained with antibodies for S100A13, CSTB or negative control without any primary antibody.

Analysis of the cytoplasmic expression showed that, similar to the proteomics results, significant S100A13 staining was more often found in non-responders compared with responders (P=0.013), while CSTB did not deviate between the groups (see Table 3). No difference in nuclear expression was observed for either protein (data not shown).

Table 3

Cytoplasmic expression of S100A13 and CSTB in pre-treatment tumour biopsies from DTIC/TMZ responders (R) and non-responders (NR)

	Positive cytoplasmic staining (%)	Fisher's exact test (P)
S100A13
R	4/16 (25)	0.013
NR	21/32 (66)
CSTB
R	5/16 (31)	0.49
NR	7/34 (21)

Discussion

The aim of the present study was to unravel a protein signature predictive of the efficacy of DTIC/TMZ therapy in pre-treatment melanoma tumour biopsies, using mass spectrometry-based proteomics. We compared the proteome profile of tumour samples from chemotherapy non-responders with that of responders. For four of the identified proteins (the cathepsin inhibitor cystatin B (CSTB), myo-inositol 1-phosphate synthase (ISYNA1, that catalyses the de novo synthesis of inositol 1-phosphate from glucose 6-phosphate), the coagulation factor F13A1 and the calcium binding protein S100A13) the expression differences between tumours from responders and non-responders were validated by immunoblotting, indicating the robustness of the MS analysis.

To our knowledge, this is the first LC/MS-MS study of drug resistance in human melanoma tumours. Resistance to TMZ and other alkylating agents has been studied by other proteomics methods mainly in glioma samples and cell lines (refer to the study by Suk (2012) for a review). Resistance to TMZ or other chemotherapeutic agents (vindesine, cisplatin, fotemustine or etoposide) has also been investigated by gene expression microarray analyses (Augustine ) and proteomics (Sinha ; Paulitschke ) in melanoma cell lines. Interestingly, members of the top canonical pathway identified in our study, ‘Signalling by Rho Family GTPases', are among the strongest chemoresistance candidates also in melanoma cell lines treated with other anticancer agents (Augustine ). Additional functional categories from other studies that overlap with our results include potential druggable candidate pathways, namely heat shock proteins and proteasomal proteins (Sinha ; Suk, 2012).

The correlation between S100A13 expression and DTIC/TMZ resistance detected in the global proteomics analysis was confirmed by immunohistochemistry in an extended set of pre-treatment melanoma tumours. The calcium binding S100 family proteins are involved in several intracellular processes such as phosphorylation regulation of proteins, cell differentiation and cell growth and cytoskeleton dynamics (Donato ). They are reported to be over-expressed in several cancer types, including melanoma tumour cells (Hardesty ) and in gliomas (Schwartz ). S100A13 has previously been suggested to be an angiogenic marker for melanoma and astrocytic gliomas (Landriscina ; Massi ) and also to be involved in the invasiveness of lung cancer cell lines (Pierce ).

There could be several potential mechanisms by which S100A13 influences chemoresistance: in tumours, S100A13 could be a marker for vessel density and on a cellular level, S100A13 regulates secretion of FGF1 (Cao ) and IL1A (Mohan and Yu, 2011). In addition to its proposed role in angiogenesis and invasion, S100A13 is also involved in several key signalling pathways with relevance to cell cycle progression and differentiation, including cytokine and NFKB signalling, supporting that S100A13 may be related to increased aggressiveness of melanoma tumours (Hsieh ; Massi ). Our results are also in line with the reported positive correlation between S100A13 mRNA levels, risk of relapse and status of melanoma patients at follow-up (Massi ). The elevated expression of S100A proteins, including S100A13, in melanoma cell lines resistant to cisplatin (Paulitschke ) indicate that S100A13 could have a more general role in melanoma chemoresistance, not specific for DTIC/TMZ.

The second candidate protein that was investigated in the larger material was CSTB. This endogenous inhibitor of cathepsins l, h, b and papain is thought to be involved in the protection against lysosomal protease leakage, and has been suggested to be a marker for more advanced disease and bad prognosis in various forms of cancer (Kos ; Werle ; Ghashenko ). However, our immunohistochemistry results did not support the hypothesis of CSTB playing a significant role in DTIC/TMZ resistance.

We have previously shown that the main melanogenesis regulator, the transcription factor MITF, is associated with resistance to DTIC/TMZ (Hertzman Johansson ) and the relatively high proportion of MITF-regulated proteins among the candidates identified in this study supports this finding.

In summary, our results showed significantly higher cytoplasmic levels of the S100A13 protein in pre-treatment tumour biopsies from DTIC/TMZ non-responders compared with responders. Thus, low levels or no expression of S100A13 may be one of the key predictive markers to identify melanoma patients responding to DTIC/TMZ therapy. Further in-depth studies on how S100A13 regulates response to DTIC/TMZ are needed and the finding should be validated in larger independent sets of tumour samples from melanoma patients.