The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD011175.

Introduction {#sec001}
============

Malignant melanoma (MM) is a type of cancer that has a high rate of incidence in many countries \[[@pone.0208562.ref001]\]. Along with an unfavorable prognosis, MM lacks any well-established biomarkers to aid in early detection, progression and treatment \[[@pone.0208562.ref001],[@pone.0208562.ref002]\]. In this context, the use of -*omic* approaches is essential to further understand the disease. In particular, proteomics has led to the identification of disease-related mechanisms by exploring the proteomes from melanoma-derived cell lines, blood and patient tumors \[[@pone.0208562.ref003]\]. To deeply cover several aspects of MM including the discovery of biomarkers of clinical relevance, current clinical proteomic studies primarily consider large MM patient cohorts where samples have been correctly collected, stored and documented \[[@pone.0208562.ref004]\].

The results acquired from cell line-based studies do not necessarily reflect the complexity of the disease. Exploring MM in patient tumors is crucial because the heterogeneity of the tissue and the microenvironment of the surrounding cells is taken into consideration \[[@pone.0208562.ref005]--[@pone.0208562.ref007]\]. In addition, the analysis of individual patients could be integrated with clinical data to acquire more reliable information on biomarker candidates \[[@pone.0208562.ref004]\].

Phosphorylation is highly-relevant in many biological processes (signaling pathways, metabolism, cell growth, *etc*.) \[[@pone.0208562.ref008]\] and in particular, is central to the role of the RAS/RAF/MAPK pathway in MM proliferation and survival \[[@pone.0208562.ref003]\]. Thus, proteomics of MM has focused on the analysis of the phosphoproteome. In most cases, these studies have been performed using cell lines because of ease of handling, replication power and large quantities of available protein to perform proteomic experiments \[[@pone.0208562.ref009]--[@pone.0208562.ref011]\]. The main phosphoproteomics strategies have taken advantage of these characteristics and have been shown to be reliable in detecting and identifying thousands of phosphoproteins in a single experiment \[[@pone.0208562.ref012]\].

Amongst the most-common methodologies used to enrich phosphopeptides from samples are IMAC \[[@pone.0208562.ref013]\] and titanium dioxide (TiO~2~) \[[@pone.0208562.ref014]\], and combinations thereof with fractionation techniques such as strong-cation exchange (SCX) \[[@pone.0208562.ref015]\], hydrophilic interaction liquid chromatography (HILIC) \[[@pone.0208562.ref016]\] or basic reversed-phase chromatography \[[@pone.0208562.ref017]\]. Although these approaches have been successful, some limitations should be considered. These include: reproducibility, the large quantities of starting material required, very labor-intensive manual steps in such a complex experiment; and the lack of automation.

Such limitations are particularly challenging when investigating MM-derived patient tissue. As such, the enrichment protocols that are often utilized are susceptible to sample variability and material loss. In this context, the Agilent AssayMAP Bravo platform has been reported to have the ability to conduct protein digestion, peptide desalting and phosphopeptide enrichment as an automated workflow \[[@pone.0208562.ref018]--[@pone.0208562.ref021]\]. This automation enables the simultaneous processing of multiple samples and thereby increases both efficiency and reproducibility. In particular, the Fe(III)-IMAC-based phosphorylation enrichment protocol has shown outstanding performance in detecting thousands of phosphopeptides from only a few micrograms of protein sample with high selectivity \[[@pone.0208562.ref019],[@pone.0208562.ref021]\]. This technology thus represents an alternative means to study MM tissues.

Our group recently reported a study on protein extraction, digestion, peptide enrichment and fractionation of MM tumor samples using the Agilent AssayMap Bravo platform (submitted manuscript). In the present work, the overall performance of this instrument for the phosphoproteomic analysis of MM lymph node tumors deposited in a Biobank was evaluated. Automated phosphopeptide enrichment on the Bravo AssayMap platform revealed a very high selectivity (85--98%) for all the input material that was assessed. Additionally, reproducible and consistent results were obtained as shown by the high correlation between experimental replicates. More importantly, the high sensitivity achieved with this technology led to the generation of results that were comparable to studies where at least 12-fold larger quantities of starting material from melanoma cultured cells were used. Finally, specific biological pathways that have not been previously described in melanoma cultured cells were identified in the MM-derived tissue.

Materials and methods {#sec002}
=====================

This study was approved by the Regional Ethical Committee at Lund University, Southern Sweden, approval numbers: DNR 191/2007, 101/2013 and 2015/266, 2015/618. All patients within the study provided written informed consent. The MM tumor tissues used had been deposited in the Biobank located at the Department of Biomedical Engineering (Lund University, Sweden). Samples corresponded to lymph node metastases from melanoma patients undergoing surgery at Lund University Hospital, Sweden.

Protein extraction, digestion and automated C18 desalting workflow {#sec003}
------------------------------------------------------------------

Protein extraction was performed on sectioned, fresh-frozen tissues from human MM lymph node metastases tissues (10 μm) using the Bioruptor plus, model UCD-300 (Dieagenode). A total of 20 tissue sections were lysed in 100 μL lysis buffer containing 4 M urea and 100 mM ammonium bicarbonate. After briefly vortexing, samples were sonicated in the Bioruptor for 40 cycles. Each cycle consisted of 15 s at high power and 15 s without sonication at 4°C. The samples were then centrifuged at 10,000 ×*g* for 10 min at 4°C. According to the instructions provided by the manufacturer, the protein content in the supernatant was determined using the colorimetric micro BCA Protein Assay kit (Thermo Fisher Scientific, Rockford, IL).

Proteins were reduced with 10 mM DTT for 1 h at room temperature (RT) and sequentially alkylated with 20 mM iodoacetamide for 30 min in the dark at RT. To decrease the urea concentration, the samples were then diluted approximately seven times with 100 mM ammonium bicarbonate. Digestion was performed in two steps. Proteins were firstly incubated with Lys-C at a 1:50 (w/w) ratio (enzyme:protein) for 5 h and then trypsin was then added at a 1:50 (w/w) ratio and the mixture incubated overnight at RT. The reaction was quenched by adding 20% TFA to a final concentration of \~1%. Peptides were desalted on the AssayMAP Bravo (Agilent Technologies) platform using the peptide cleanup v2.0 protocol. Briefly, C18 cartridges (Agilent, 5 μL bed volume) were primed with 100 μL 90% acetonitrile (ACN) and equilibrated with 70 μL 0.1% TFA at a flow rate of 10 μL/min. The samples were loaded at 5 μL/min, followed by an internal cartridge wash with 0.1% TFA at a flow rate of 10 μL/min. Peptides were eluted with 30 μL 80% ACN, 0.1% TFA and diluted to 200 μL with the same buffer for subsequent phosphopeptide enrichment.

Automated Fe(III)-IMAC-based workflow {#sec004}
-------------------------------------

Following the Phospho Enrichment v2.0 protocol, phosphorylated peptides were enriched using 5 μL Fe(III)-NTA cartridges on the AssayMAP Bravo platform. The cartridges were primed with 100 μL 50% ACN, 0.1% TFA at a flow rate of 300 μL/min and equilibrated with 50 μL loading buffer (80% ACN, 0.1% TFA) at 10 μL/min. Samples eluted from the C18 material were loaded onto the cartridge at 3.5 μL/min. The columns were washed with 50 μL loading buffer and the phosphorylated peptides were eluted with 25 μL 5% ammonia directly into 10 μL 50% formic acid. Samples were lyophilized in a vacuum concentrator and stored at -80°C until analysis by LC-MS/MS.

Automated off line high pH fractionation {#sec005}
----------------------------------------

Reversed-phase high pH (HpH) fractionation of the MM digests were performed on the Bravo AssayMAP system using the Fractionation v1.0 protocol. The pH of the digested samples were adjusted to 10 with 200 mM ammonium formate. RPS cartridges (Agilent, 5 μL bed volume) were primed with 100 μL 90% ACN and equilibrated with 70 μL 20 mM ammonium formate (pH 10) at a flow rate of 10 μL/min. Samples were loaded at 5 μL/min followed by an internal cartridge wash with 20 mM ammonium formate at 10 μL/min. Elution was performed in six steps with increasing concentrations of ACN in ammonium formate: 4% (F1), 10% (F2), 15% (F3), 20% (F4), 25% (F5) and 80% (F6); 30 μL was eluted into 50 μL of pre-existing buffer containing 50% ACN, 0.1% TFA. The flow through (FT) was collected, processed on C18 cartridges and analyzed as an additional fraction.

LC-MS/MS analysis {#sec006}
-----------------

LC-MS/MS analysis was performed on a Thermo Easy-nLC 1000 HPLC coupled on-line to a Q Exactive mass spectrometer (Thermo Scientific, San Jose, CA) or on an Ultimate 3000 HPLC coupled to a Q Exactive HF-X mass spectrometer (Thermo Scientific, San Jose, CA). Peptides were loaded onto a trap column (Acclaim1 PepMap 100 pre-column, 75 μm, 2 cm, C18, 3 mm, 100 Å, Thermo Scientific, San José, CA) and then separated on an analytical column (EASY-Spray column, 25 cm, 75 μm i.d., PepMap RSLC C18, 2 mm, 100Å, Thermo Scientific, San José, CA) using a 120 min ACN gradient in 0.1% formic acid at a flow rate of 300 nL/min and a column temperature of 45°C. Settings on the Q Exactive mass spectrometer were as follows: data-dependent analysis (DDA) selection of the 10 most intense ions for fragmentation, full MS scans at *m/z* 375--1,750 with a resolution of 60,000 at *m/z* 200, target AGC value of 3×10^6^ and an injection time (IT) of 60 ms, fragmentation in the HCD collision cell with an NCE of 25 and MS/MS spectra acquisition in the Orbitrap analyzer at a resolution of 35,000 (at *m/z* 200) with a maximum IT of 120 ms and dynamic exclusion of 30 s. Settings on the Q Exactive HF-X mass spectrometer were as follows: DDA selection of the 15 most intense ions for fragmentation, full MS scans at *m/z* 375--1,750 with a resolution of 120,000 at *m/z* 200, a target AGC value of 3×10^6^ and IT of 100 ms, fragmentation in HCD collision cell with an NCE of 25 and MS/MS spectra acquisition in the Orbitrap analyzer at a resolution of 60,000 (at *m/z* 200) with a maximum IT of 120 ms and dynamic exclusion of 30 s.

Data analysis {#sec007}
-------------

Raw files were processed with Proteome Discoverer 2.1 (Thermo Scientific). The search was performed against the *Homo sapiens* UniProt revised database with the Sequest HT search engine. Cysteine carbamidomethylation was set as fixed modification and methionine oxidation, protein N-terminal acetylation and phosphorylation of serine, threonine and tyrosine were set as variable modifications; peptide mass tolerance for the precursor ions and MS/MS spectra were 10 ppm and 0.02 Da, respectively. A maximum of two missed cleavage sites were accepted. The ptmRS algorithm was used to score phosphorylation sites with a site probability threshold \>75. Filtering was performed at a 1% false discovery rate (FDR) for both peptides and proteins. For statistical analyses, label-free quantification (LFQ) data was log~2~-transformed and normalized to column median as standard methods. Analyses such as Pearson correlation, ANOVA and clustering were performed using the Perseus software.

Results and discussion {#sec008}
======================

Experiments were conducted on proteins obtained from the lysis of fresh, frozen MM tissues (20 sections per tissue). Three experiments were performed: (A) sensitivity assessment of automated enrichment from pooled MM lysates; (B) analysis of the MM phosphoproteome from individual tumor samples; and (C) deep coverage of the MM phosphoproteome using HpH fractionation ([Fig 1](#pone.0208562.g001){ref-type="fig"}).

![Overview of the workflows performed for the phosphoproteomic analyses.\
(A) Phosphoproteome analysis with varying quantities of MM peptides extracted from pooled tissues; each enrichment was acquired in technical triplicates. (B) Phosphopeptide enrichment of digested MM tissues performed on three different tissues in biological duplicates. (C) Phosphopeptide enrichment combined with HpH fractionation for a single MM sample (one biological replicate).](pone.0208562.g001){#pone.0208562.g001}

Sensitivity analysis of automated phosphopeptide enrichment {#sec009}
-----------------------------------------------------------

To evaluate the sensitivity of the protocol, the phosphopeptide enrichment was assessed on increasing quantities of MM protein digests (12.5--200 μg) taken from a pool of eight MM tumor lysates (characteristics of the tissue samples used in the pool are described in [S1 Table](#pone.0208562.s005){ref-type="supplementary-material"}). A pooled sample was generated to avoid a bias towards the characteristics of a single tissue characteristic and to obtain data that could be extrapolated to different MM tissues. As shown in [Fig 2](#pone.0208562.g002){ref-type="fig"}, the method was sufficiently sensitive to identify more than 1,000 phosphopeptides from 12.5 μg of input material. As expected, the number of identifications increased when more input material was used. The samples displayed linearity up to 100 μg of loaded peptides. When higher quantities were injected onto the cartridges, however, linearity was compromised suggesting that the overall capacity was affected ([S1 Fig](#pone.0208562.s001){ref-type="supplementary-material"}). The same trend has been previously observed with cell lines, thus reflecting the credibility of our results and the robustness of the methodology \[[@pone.0208562.ref019]\]. For all peptide amounts analyzed, the enrichment selectivity was \>90%, however, this value may be compromised when lower quantities of material are analyzed. An assumption such as this can be explained partly by the fact that a poor MS signal-to-noise ratio markedly affects phosphopeptide identification \[[@pone.0208562.ref019]\]. The enrichment from 100 μg of the MM peptides resulted in the best data. Over 5,226 unique phosphopeptides and 5,524 phosphosites class I were identified from a single 120 min LC-MS/MS analysis. Excluding slight differences in phosphopeptide ratios, for all peptide quantities loaded onto the Fe(III)-NTA cartridges the enrichment of phosphorylated serine (pS), threonine (pT) and tyrosine (pY) followed the expected pattern whereby pS\>pT\>\>pY. The same observations were made with respect to the number of phosphopeptides that were identified containing multiple phosphorylation sites ([S2 Table](#pone.0208562.s006){ref-type="supplementary-material"}).

![Sensitivity curve of the automated phosphoproteome analysis.\
Phosphopeptide enrichment was performed with a pooled MM digest using 12.5 to 200 ug of input material. Values expressed at the top of the graph represent the ratio pS:pT:pY obtained in each enrichment; left axis represents the number of phosphopeptides and phosphoproteins in each experiment (from three technical replicates); right axis expresses the enrichment selectivity of each experiment. LC-MS/MS analysis was performed on a Q Exactive mass spectrometer.](pone.0208562.g002){#pone.0208562.g002}

In summary, the factors discussed above suggest that the optimal performance of the cartridges is in the range of 12.5--100 μg of input material without any compromise in the overall enrichment efficiency. Additionally, quantities less than 100 μg still produced high-quality results and are ideal for cases where very limited quantity of starting material are available.

Analysis of the MM phosphoproteome from individual tumor samples {#sec010}
----------------------------------------------------------------

On the basis of the results obtained from the sensitivity analysis, phosphopeptide enrichment from three different MM tissues was performed (two independent experiments in each tissue). The description of the MM tumors used here are provided in [S3 Table](#pone.0208562.s007){ref-type="supplementary-material"}. From 60 μg of input material, the data was outstanding with more than 5,600 phosphopeptides and 9,000 phosphosites identified per sample. The three enrichments reflected similar trends in selectivity, pS:pT:pY ratio, and proportion of phosphosite class I and phosphosite multiplicity ([Table 1](#pone.0208562.t001){ref-type="table"}). The reproducibility of the enrichments was assessed by correlation analyses of the quantitated phosphopeptide intensities across independent replicates. Good correlations (Pearson correlation, r\>0.88) and a uniform distribution was observed for each comparison, indicating a high reproducibility between experiments ([Fig 3A](#pone.0208562.g003){ref-type="fig"} and [S3 Fig](#pone.0208562.s003){ref-type="supplementary-material"}). These data reflect the consistency of the enrichment across the samples. In this sense, the physical characteristics of the MM tissues appeared to not have any significant influence on the entire sample preparation process.

![Correlation and overlapping between the three MM samples analyzed.\
(A) Pearson correlation between independently-enriched samples for the three different MM tissues. (B) Overlap between the phosphoproteomes of the three MM samples.](pone.0208562.g003){#pone.0208562.g003}

10.1371/journal.pone.0208562.t001

###### Phosphopeptide enrichment for each of the analyzed MM samples.

Data represented here was obtained from merging two independent biological replicates.

![](pone.0208562.t001){#pone.0208562.t001g}

  --------------------------------------------------------------------------
  Description                Sample                           
  -------------------------- -------------- ----------------- --------------
  Total peptides             5,816          12,512            8,614

  Total phosphopeptides      5,638          12,201            8,242

  \% Phosphopeptides         96.9           97.5              95.6

  Total proteins (groups)    2,268          3,599             2,831

  Phosphosite multiplicity   Mono: 5,076\   Mono: 10,370\     Mono: 7,100\
                             Di: 521\       Di: 1,561\        Di: 1,016\
                             Tri: 39\       Tri: 245\         Tri: 116\
                             Tetra: 2       Tetra: 25         Tetra: 10

  Ratio phosphosites\        38: 7: 1       63: 9: 1          58: 9 1
  (pS: pT: pY)                                                

  Total phosphosites         9,082          17,669            12,164

  \% Phosphosites Class I    71             71                75

  Starting material (μg)     60                               

  Instrument                 Q Exactive     Q Exactive HF-X   
  --------------------------------------------------------------------------

In general, the overlap of the identified phosphopeptides between replicates and between patients was higher than 73% and 45%, respectively ([S2 Fig](#pone.0208562.s002){ref-type="supplementary-material"} and [Fig 3B](#pone.0208562.g003){ref-type="fig"}). In principle, the variation amongst the samples could be explained by the differences within the tissues and the performance of the mass spectrometer that was used. Despite these variations, the high overlap between the three samples enabled both pathway enrichment analysis and a quantitative comparison of the phosphoproteomes to detect protein clusters with biological significance. As summarized in [Fig 4A](#pone.0208562.g004){ref-type="fig"}, the main pathways that were commonly-enriched for the identified phosphoproteins in the three samples were related to aspects of the melanoma phenotype. Namely, mitogen-activated protein kinases (MAPK) signaling pathways, control of gene expression by vitamin D receptor, integrin signaling, melanocyte development and pigmentation *etc*.

![Enriched pathways and clustering of common proteins in MM samples.\
(A) Enriched pathways related to common proteins identified in MM1, MM2 and MM3. (B) Heat map clustering of proteins with high significance amongst the MM samples obtained from ANOVA analysis (permutation-based FDR \>0.01).](pone.0208562.g004){#pone.0208562.g004}

MAPK pathways are essential in the development of MM because of the involvement in cell proliferation and evasion of apoptosis \[[@pone.0208562.ref003]\]. Proteins that belong to the classical MAPK pathway, such as v-Raf murine sarcoma viral oncogene product (BRAF) and the downstream partners extracellular signal-regulated kinases 1 and 2 (ERK1 and ERK2) were all observed in our data set. These proteins are extremely relevant because of the role played in the activation of several transcription factors that leads to cell proliferation and cell development \[[@pone.0208562.ref022]\]. Considering that BRAF is mutated in 50--60% of melanomas \[[@pone.0208562.ref023]\], it is not unreasonable that major drug development has focused on targeting these protein variants. Hence, vemurafenib, trametinib \[[@pone.0208562.ref024]\], dabrafenib \[[@pone.0208562.ref025]\], and ipilimumab \[[@pone.0208562.ref026]\] have been approved or considered as potential cancer therapies against MM. Moreover, MAPK variants were also observed in our samples. In the p38 MAPK cascade, the phosphorylation of AKT, PTP, MEKK1, MKK6, 7, 3 and 4 and p38 are relevant in regulating inflammation, apoptosis and cell cycle steps \[[@pone.0208562.ref022]\]. In particular, it is known that pAKT acts as an oncogenic signal that promotes tumor progression \[[@pone.0208562.ref027]\]. In addition, this cascade allows the activation of p53, which is one of the best known cancer suppressors \[[@pone.0208562.ref028]\]. This protein remains controversial in MM as it is widely present in MM patients but only with very rare mutations. In this sense, the protein remain silent in 90% of MM patients but reactivation combined with BRAF^V600E^ inhibition leads to apoptosis and suppression of melanoma growth \[[@pone.0208562.ref028]\].

Other enriched MAPK-related pathways were related to proliferation, *e*.*g*., ERK5 and some scaffold proteins (HGK, HPK1 and MEKK1) that are responsible for activating the Wnt pathway \[[@pone.0208562.ref003]\]. Synergy between MAPK and Wnt leads to the formation and development of melanoma \[[@pone.0208562.ref003]\]. In addition, non-canonical Wnt proteins induce the release of exosomes from tumor cells containing pro-angiogenic and immunosuppressive factors \[[@pone.0208562.ref029]\].

Together with MAPK, another pathway that acts synergistically in the development of MM is the melanocyte development and pigmentation pathway. Key proteins such as proto-oncogene receptor tyrosine kinase KIT were observed. KIT is involved in migration and survival of tumor cells, but in higher metastatic states expression is lost \[[@pone.0208562.ref030]\]. Other phosphoproteins that belong to this cascade were the melanogenesis associated transcriptor factor (MITF), apoptosis regulator BCL2 and cAMP responsive element binding protein 1 (CREB1) \[[@pone.0208562.ref003]\].

The integrin β pathway was also identified in our data set. The expression of integrin β1 facilitates metastases (migration and proliferation) in lymph nodes \[[@pone.0208562.ref031]\]. Several downstream proteins activated by phosphorylated integrin β1 such as tensin, zyxin, talin, vinculin, SRC proto-oncogene and paxillin (all involved in F-actin dynamics during migration and cell division) were also identified \[[@pone.0208562.ref032]\]. Activation of caveolin by interaction with collagen also activates the classical MAPK pathway by the phosphorylation of RAS/RAF. In this case, SRC acts as a redundant protein recruiting actin stress fibers and phosphorylating ERK1/ERK2 \[[@pone.0208562.ref022],[@pone.0208562.ref033],[@pone.0208562.ref034]\].

ANOVA analysis of the label-free quantitated phosphoproteins that were common to the patient samples reflected interesting results ([Fig 4B](#pone.0208562.g004){ref-type="fig"}). The main pathways that were previously mentioned appeared as key clusters with differential expression in the three samples, thus demonstrating the relevance of these proteins in MM phenotypes. The differences in expression could potentially be explained by several factors related to the characteristic of each patient (sex, age, percentage of tumor in the tissues, treatment, survival, *etc*.). Unfortunately, the number of samples used here is not sufficient to perform a truly robust statistical analysis to correlate proteins or clusters with clinical insights. Methodologically speaking however, the results are valuable because identification/quantitation of significant differences in phosphoproteomic profile of MM samples regardless of the clinical or physical characteristics of the samples was feasible.

Deep coverage of the MM phosphoproteome via basic reversed-phase fractionation {#sec011}
------------------------------------------------------------------------------

To increase the depth of the analysis with respect to the number of phosphopeptide identifications, high pH reversed-phase peptide pre-fractionation of the MM2 sample was also performed on the AssayMAP Bravo platform. Application of an optimized step elution method (see [Materials and methods](#sec002){ref-type="sec"}) resulted in the collection of seven fractions that included the flow-through. Using the Fe(III)-NTA cartridges, phosphopeptides were enriched from each fraction. [Fig 5](#pone.0208562.g005){ref-type="fig"} and [S4 Table](#pone.0208562.s008){ref-type="supplementary-material"} show the number of phosphopeptides that were identified from each individual fraction, the overall analysis and in the sample without fractionation. From 60 ug of starting material, 13,036 phosphopeptides and 4,359 phosphoproteins were identified. Compared to the unfractionated sample (8,723 ± 600 phosphopeptides per replicate), this represents an increase of 49% in the identification of phosphopeptides and phosphoproteins.

![Profile of the phosphopeptide enrichment from the six fractions and the FT obtained from HpH fractionation on the Bravo AssayMap platform for sample MM2.\
LC-MS/MS analysis was performed on a Q Exactive HF-X mass spectrometer. Results are compared with the overall analysis and the sample without fractionation. Data corresponds to a single replicate.](pone.0208562.g005){#pone.0208562.g005}

This data set was compared with other available high-throughput phosphoproteomic analyses involving MM cell lines \[[@pone.0208562.ref009]--[@pone.0208562.ref011]\]. The numbers obtained here are relatively close to those obtained in previous studies ([Fig 6A](#pone.0208562.g006){ref-type="fig"} and [S5 Table](#pone.0208562.s009){ref-type="supplementary-material"}), with approximately 90% of the proteins identified in this experiment observed in prior studies utilizing various approaches and methods. This result is beyond expectations because 12.5- to 100-fold less starting material was used. Similarly, the fractionation method used here was not as extensive as has been shown in earlier studies. Interestingly, 672 proteins were exclusively detected in our data set ([Fig 6A and 6B](#pone.0208562.g006){ref-type="fig"}). These proteins are involved in pathways and biological processes related to immunological responses from T cells, B cells and NK cells plus interleukin (IL) pathways ([Fig 6C](#pone.0208562.g006){ref-type="fig"}). As the phosphoproteomic analysis in our study was performed on a tissue sample, it was possible to enrich for pathways closely-related to lymph node functions. By using tissues rather than cultured cells, processes that cannot be fully explored with cells were immediately apparent. This is simply because cell culture do not truly reflect the *in vivo* environment of the disease.

![Comparison with previously reported phosphoproteome data sets.\
(A) Venn diagram comparison of three different reported phosphoproteome data sets from melanoma-derived cell lines with data obtained from HpH-fractionated sample MM2 (single biological replicate). Each comparison is accompanied by the quantity of starting material and the reported fractionation method. (B) Overall comparison between the phosphoproteomic data sets and HpH-fractionated MM2. (C) Enriched pathways related to proteins exclusively identified in our study.](pone.0208562.g006){#pone.0208562.g006}

As lymph nodes are the major sites of B and T lymphocytes and represent the most frequent early metastatic sites in MM, it is not unreasonable to assume that the antitumor immune response is relevant for the reduction of metastases \[[@pone.0208562.ref035]\]. Several molecules involved in the progression of antitumor response were identified, *i*.*e*., CD8, CD2, CD11a, CD18, CD45, protein tyrosine phosphatase receptor type C (in T cytotoxic and T helper cells), JAK and STAT proteins, CD3, CXCR3, LCK proto-oncogene, MEK and IL-2R (IL-2 response and NO-dependent IL-12 activation of NK cells). These immune cells play an important role in combatting MM as these assist in establishing a first line anti-tumor immune response in lymph nodes \[[@pone.0208562.ref036]\]. Immunological treatment of MM that targets some of these proteins is under investigation. Interleukins such as IL-2 or IL-15 are candidates for therapeutic applications based on cytokine action, but severe side effects are associated with systemic cytokine treatment \[[@pone.0208562.ref037]\]. This 'exclusive' data set provides complementary information regarding the biology of MM and enables the exploration of alternative/potential targets in immunological treatment of MM.

Despite a lower number of identified phosphopeptides in the unfractionated sample, sufficient information was still obtained to provide satisfactory biological relevance and context. The comparison with other studies revealed that a significant number of proteins (468 proteins, [S4 Fig](#pone.0208562.s004){ref-type="supplementary-material"}) were still exclusively identified by this approach. In this context, processing samples without fractionation is a reasonable option to study the phosphoproteome of MM. In particular, such an option would be more practical when many samples, *e*.*g*., from a large patient cohort, would need to be simultaneously processed and analyzed.

Conclusion {#sec012}
==========

Clinical proteomic studies require large sample cohorts to discover disease-related mechanisms and biomarkers of clinical significance. Thus, it is essential that the field progresses towards a high degree of sensitivity plus robustness in sample preparation. Realistically, this can only be achieved with automated technologies. Hence, it is our belief that instruments such as the Bravo AssayMap platform could play a key role in the analysis of large cohorts of challenging biological samples such as MM. An automated phosphorylation enrichment protocol was assessed. Here, Fe(III)-NTA IMAC cartridges were used in the Bravo AssayMap system and combined with high-performance MS instruments. This method proved to be efficient, sensitive and reproducible in tissue samples obtained from patients with MM. The protocol was robust and enabled the identification of thousands of phosphorylated peptides from a limited quantity of starting material. The MM tissue-based phosphoproteome covered essential features of MM biology and was comparable to reports based on MM-derived cell lines. More importantly, additional MM-related pathways were revealed that can only be detected in human tissue. Overall, our study represents a significant step towards a deeper understanding of malignant melanoma and the opportunity to integrate clinical and phosphoproteomic data in biomarker discovery.

Supporting information {#sec013}
======================

###### Linearity of the Fe(III)-NTA IMAC-enriched and quantitated phosphopeptides on the Bravo AsssayMap platform.

Only phosphopeptides identified in all replicates were considered for this analysis. All values were normalized with respect to the initial 200 μg protein quantity (intensity = 1). The enrichment displayed quantitative linear performance up to 100 μg of input material.

(TIF)

###### 

Click here for additional data file.

###### Venn diagram comparison of identified phosphoproteins in each replicate for the three patient tissues that were analyzed.

(TIF)

###### 

Click here for additional data file.

###### Normalized protein intensities for each biological replicate of the three MM tumors (only including common proteins identified in all replicates).

(TIF)

###### 

Click here for additional data file.

###### Venn diagram comparison of melanoma-derived cell line data obtained from major published phosphoproteomic data sets with the data obtained from sample MM2 without fractionation.

(TIF)

###### 

Click here for additional data file.

###### Overall description of MM tissues used in the pool to perform sensitivity analysis of phosphopeptide enrichment.

Some patient clinical information is also provided.

(DOCX)

###### 

Click here for additional data file.

###### Detailed results of the sensitivity analysis of Fe(III)-NTA IMAC phosphoenrichment performed on the Bravo AssayMap platform.

(DOCX)

###### 

Click here for additional data file.

###### Overall description of MM tissues used in the phosphopeptide enrichment of individual samples.

Some patient clinical information is also provided.

(DOCX)

###### 

Click here for additional data file.

###### Detailed description of the phosphoproteome obtained from Fe(III)-NTA IMAC phosphoenrichment on the Bravo AssayMap platform combined with HpH fractionation.

(DOCX)

###### 

Click here for additional data file.

###### Overall comparison of other available high-throughput phosphoproteomic analyses involving MM cell lines with this study.

(DOCX)

###### 

Click here for additional data file.

The authors acknowledge Agilent Technologies for the supply of the instrument and The Royal Physiographic Society of Lund (2017) for the support with the experiments.

[^1]: **Competing Interests:**The authors (including those from AstraZeneca) have declared that no competing interests exist. AstraZeneca affiliation does not alter our adherence to PLOS ONE policies on sharing data and materials.
