Functional proteomics of colon cancer Consensus Molecular Subtypes.

BACKGROUND: The Colorectal Cancer Subtyping Consortium established four Consensus Molecular Subtypes (CMS) in colorectal cancer: CMS1 (microsatellite-instability [MSI], Immune), CMS2 (Canonical, epithelial), CMS3 (Metabolic), and CMS4 (Mesenchymal). However, only MSI tumour patients have seen a change in their disease management in clinical practice. This study aims to characterise the proteome of colon cancer CMS and broaden CMS's clinical utility. METHODS: One-hundred fifty-eight paraffin samples from stage II-III colon cancer patients treated with adjuvant chemotherapy were analysed through DIA-based mass-spectrometry proteomics. RESULTS: CMS1 exhibited overexpression of immune-related proteins, specifically related to neutrophils, phagocytosis, antimicrobial response, and a glycolytic profile. These findings suggested potential therapeutic strategies involving immunotherapy and glycolytic inhibitors. CMS3 showed overexpression of metabolic proteins. CMS2 displayed a heterogeneous protein profile. Notably, two proteomics subtypes within CMS2, with different protein characteristics and prognoses, were identified. CMS4 emerged as the most distinct group, featuring overexpression of proteins related to angiogenesis, extracellular matrix, focal adhesion, and complement activation. CMS4 showed a high metastatic profile and suggested possible chemoresistance that may explain its worse prognosis. CONCLUSIONS: DIA proteomics revealed new features for each colon cancer CMS subtype. These findings provide valuable insights into potential therapeutic targets for colorectal cancer subtypes in the future.

Utility of CYP2D6 copy number variants as prognostic biomarker in localized anal squamous cell carcinoma.

BACKGROUND: Anal squamous cell carcinoma (ASCC) is an infrequent tumor whose treatment has not changed since the 1970s. The aim of this study is the identification of biomarkers allowing personalized treatments and improvement of therapeutic outcomes. METHODS: Forty-six paraffin tumor samples from ASCC patients were analyzed by whole-exome sequencing. Copy number variants (CNVs) were identified and their relation to disease-free survival (DFS) was studied and validated in an independent retrospective cohort of 101 ASCC patients from the Multidisciplinary Spanish Digestive Cancer Group (GEMCAD). GEMCAD cohort proteomics allowed assessing the biological features of these tumors. RESULTS: On the discovery cohort, the median age was 61 years old, 50% were males, stages I/II/III: 3 (7%)/16 (35%)/27 (58%), respectively, median DFS was 33 months, and overall survival was 45 months. Twenty-nine genes whose duplication was related to DFS were identified. The most representative was duplications of the CYP2D locus, including CYP2D6, CYP2D7P, and CYP2D8P genes. Patients with CYP2D6 CNV had worse DFS at 5 years than those with two CYP2D6 copies (21% vs. 84%; p < .0002, hazard ratio [HR], 5.8; 95% confidence interval [CI], 2.7-24.9). In the GEMCAD validation cohort, patients with CYP2D6 CNV also had worse DFS at 5 years (56% vs. 87%; p = .02, HR = 3.6; 95% CI, 1.1-5.7). Mitochondria and mitochondrial cell-cycle proteins were overexpressed in patients with CYP2D6 CNV. CONCLUSIONS: Tumor CYP2D6 CNV identified patients with a significantly worse DFS at 5 years among localized ASCC patients treated with 5-fluorouracil, mitomycin C, and radiotherapy. Proteomics pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets for these high-risk patients. PLAIN LANGUAGE SUMMARY: Anal squamous cell carcinoma is an infrequent tumor whose treatment has not been changed since the 1970s. However, disease-free survival in late staged tumors is between 40% and 70%. The presence of an alteration in the number of copies of CYP2D6 gene is a biomarker of worse disease-free survival. The analysis of the proteins in these high-risk patients pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets. Therefore, the determination of the number of copies of CYP2D6 allows the identification of anal squamous carcinoma patients with a high-risk of relapse that could be redirected to a clinical trial. Additionally, this study may be useful to suggest new treatment strategies to increase current therapy efficacy.

Sorting Transcriptomics Immune Information from Tumor Molecular Features Allows Prediction of Response to Anti-PD1 Therapy in Patients with Advanced Melanoma.

Immunotherapy based on anti-PD1 antibodies has improved the outcome of advanced melanoma. However, prediction of response to immunotherapy remains an unmet need in the field. Tumor PD-L1 expression, mutational burden, gene profiles and microbiome profiles have been proposed as potential markers but are not used in clinical practice. Probabilistic graphical models and classificatory algorithms were used to classify melanoma tumor samples from a TCGA cohort. A cohort of patients with advanced melanoma treated with PD-1 inhibitors was also analyzed. We established that gene expression data can be grouped in two different layers of information: immune and molecular. In the TCGA, the molecular classification provided information on processes such as epidermis development and keratinization, melanogenesis, and extracellular space and membrane. The immune layer classification was able to distinguish between responders and non-responders to immunotherapy in an independent series of patients with advanced melanoma treated with PD-1 inhibitors. We established that the immune information is independent than molecular features of the tumors in melanoma TCGA cohort, and an immune classification of these tumors was established. This immune classification was capable to determine what patients are going to respond to immunotherapy in a new cohort of patients with advanced melanoma treated with PD-1 inhibitors Therefore, this immune signature could be useful to the clinicians to identify those patients who will respond to immunotherapy.

Molecular characterization of triple negative breast cancer formaldehyde-fixed paraffin-embedded samples by data-independent acquisition proteomics.

Triple negative breast cancer accounts for 15%-20% of all breast carcinomas and is clinically characterized by an aggressive phenotype and poor prognosis. Triple negative tumors do not benefit from targeted therapies, so further characterization is needed to define subgroups with potential therapeutic value. In this work, the proteomes of 125 formalin-fixed paraffin-embedded samples from patients diagnosed with non-metastatic triple negative breast cancer were analyzed using data-independent acquisition + in a LTQ-Orbitrap Fusion Lumos mass spectrometer coupled to an EASY-nLC 1000. 1206 proteins were identified in at least 66% of the samples. Hierarchical clustering, probabilistic graphical models and Significance Analysis of Microarrays were combined to characterize proteomics-based molecular groups. Two molecular groups were defined with differences in biological processes such as glycolysis, translation and immune response. These two molecular groups showed also several differentially expressed proteins. This clinically homogenous dataset may serve to design new therapeutic strategies in the future.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification.

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

Computational models applied to metabolomics data hints at the relevance of glutamine metabolism in breast cancer.

BACKGROUND: Metabolomics has a great potential in the development of new biomarkers in cancer and it has experiment recent technical advances. METHODS: In this study, metabolomics and gene expression data from 67 localized (stage I to IIIB) breast cancer tumor samples were analyzed, using (1) probabilistic graphical models to define associations using quantitative data without other a priori information; and (2) Flux Balance Analysis and flux activities to characterize differences in metabolic pathways. RESULTS: On the one hand, both analyses highlighted the importance of glutamine in breast cancer. Moreover, cell experiments showed that treating breast cancer cells with drugs targeting glutamine metabolism significantly affects cell viability. On the other hand, these computational methods suggested some hypotheses and have demonstrated their utility in the analysis of metabolomics data and in associating metabolomics with patient's clinical outcome. CONCLUSIONS: Computational analyses applied to metabolomics data suggested that glutamine metabolism is a relevant process in breast cancer. Cell experiments confirmed this hypothesis. In addition, these computational analyses allow associating metabolomics data with patient prognosis.

Computational metabolism modeling predicts risk of distant relapse-free survival in breast cancer patients.

Aim: Differences in metabolism among breast cancer subtypes suggest that metabolism plays an important role in this disease. Flux balance analysis is used to explore these differences as well as drug response. Materials & methods: Proteomics data from breast tumors were obtained by mass-spectrometry. Flux balance analysis was performed to study metabolic networks. Flux activities from metabolic pathways were calculated and used to build prognostic models. Results: Flux activities of vitamin A, tetrahydrobiopterin and ß-alanine metabolism pathways split our population into low- and high-risk patients. Additionally, flux activities of glycolysis and glutamate metabolism split triple negative tumors into low- and high-risk groups. Conclusion: Flux activities summarize flux balance analysis data and can be associated with prognosis in cancer.

A Novel Molecular Analysis Approach in Colorectal Cancer Suggests New Treatment Opportunities.

Colorectal cancer (CRC) is a molecular and clinically heterogeneous disease. In 2015, the Colorectal Cancer Subtyping Consortium classified CRC into four consensus molecular subtypes (CMS), but these CMS have had little impact on clinical practice. The purpose of this study is to deepen the molecular characterization of CRC. A novel approach, based on probabilistic graphical models (PGM) and sparse k-means-consensus cluster layer analyses, was applied in order to functionally characterize CRC tumors. First, PGM was used to functionally characterize CRC, and then sparse k-means-consensus cluster was used to explore layers of biological information and establish classifications. To this aim, gene expression and clinical data of 805 CRC samples from three databases were analyzed. Three different layers based on biological features were identified: adhesion, immune, and molecular. The adhesion layer divided patients into high and low adhesion groups, with prognostic value. The immune layer divided patients into immune-high and immune-low groups, according to the expression of immune-related genes. The molecular layer established four molecular groups related to stem cells, metabolism, the Wnt signaling pathway, and extracellular functions. Immune-high patients, with higher expression of immune-related genes and genes involved in the viral mimicry response, may benefit from immunotherapy and viral mimicry-related therapies. Additionally, several possible therapeutic targets have been identified in each molecular group. Therefore, this improved CRC classification could be useful in searching for new therapeutic targets and specific therapeutic strategies in CRC disease.

Proteomics Characterization of Clear Cell Renal Cell Carcinoma.

PURPOSE: To explore the tumor proteome of patients diagnosed with localized clear cell renal cancer (ccRCC) and treated with surgery. MATERIAL AND METHODS: A total of 165 FFPE tumor samples from patients diagnosed with ccRCC were analyzed using DIA-proteomics. Proteomics ccRCC subtypes were defined using a consensus cluster algorithm (CCA) and characterized by a functional approach using probabilistic graphical models and survival analyses. RESULTS: We identified and quantified 3091 proteins, including 2026 high-confidence proteins. Two proteomics subtypes of ccRCC (CC1 and CC2) were identified by CC using the high-confidence proteins only. Characterization of molecular differences between CC1 and CC2 was performed in two steps. First, we defined 514 proteins showing differential expression between the two subtypes using a significance analysis of microarrays analysis. Proteins overexpressed in CC1 were mainly related to translation and ribosome, while proteins overexpressed in CC2 were mainly related to focal adhesion and membrane. Second, a functional analysis using probabilistic graphical models was performed. CC1 subtype is characterized by an increased expression of proteins related to glycolysis, mitochondria, translation, adhesion proteins related to cytoskeleton and actin, nucleosome, and spliceosome, while CC2 subtype showed higher expression of proteins involved in focal adhesion, extracellular matrix, and collagen organization. CONCLUSIONS: ccRCC tumors can be classified in two different proteomics subtypes. CC1 and CC2 present specific proteomics profiles, reflecting alterations of different molecular pathways in each subtype. The knowledge generated in this type of studies could help in the development of new drugs targeting subtype-specific deregulated pathways.

Multi-omics Characterization of Response to PD-1 Inhibitors in Advanced Melanoma.

Immunotherapy improves the survival of patients with advanced melanoma, 40% of whom become long-term responders. However, not all patients respond to immunotherapy. Further knowledge of the processes involved in the response and resistance to immunotherapy is still needed. In this study, clinical paraffin samples from fifty-two advanced melanoma patients treated with anti-PD-1 inhibitors were assessed via high-throughput proteomics and RNA-seq. The obtained proteomics and transcriptomics data were analyzed using multi-omics network analyses based on probabilistic graphical models to identify those biological processes involved in the response to immunotherapy. Additionally, proteins related to overall survival were studied. The activity of the node formed by the proteins involved in protein processing in the endoplasmic reticulum and antigen presentation machinery was higher in responders compared to non-responders; the activity of the immune and inflammatory response node was also higher in those with complete or partial responses. A predictor for overall survival based on two proteins (AMBP and PDSM5) was defined. In summary, the response to anti-PD-1 therapy in advanced melanoma is related to protein processing in the endoplasmic reticulum, and also to genes involved in the immune and inflammatory responses. Finally, a two-protein predictor can define survival in advanced disease. The molecular characterization of the mechanisms involved in the response and resistance to immunotherapy in melanoma leads the way to establishing therapeutic alternatives for patients who will not respond to this treatment.

MiRNAs in renal cell carcinoma.

MicroRNAs (miRNAs) are small RNA sequences that act as post-transcriptional regulatory genes to control many cellular processes through pairing bases with a complementary messenger RNA (mRNA). A single miRNA molecule can regulate more than 200 different transcripts and the same mRNA can be regulated by multiple miRNAs. In this review, we highlight the importance of miRNAs and collect the existing evidence on their relationship with kidney cancer.

Synergistic effect of antimetabolic and chemotherapy drugs in triple-negative breast cancer.

The triple-negative breast cancer (TNBC) subtype comprises approximately 15% of all breast cancers and is associated with poor long-term outcomes. Classical chemotherapy remains the standard of treatment, with toxicity and resistance being major limitations. TNBC is a high metabolic group, and antimetabolic drugs are effective in inhibiting TNBC cell growth. We analyzed the combined effect of chemotherapy and antimetabolic drug combinations in MDA-MB-231, MDA-MB-468 and HCC1143 human TNBC cell lines. Cells were treated with each drug or with drug combinations at a range of concentrations to establish the half-maximal inhibitory concentrations (IC50). The dose-effects of each drug or drug combination were calculated, and the synergistic or antagonistic effects of drug combinations were defined. Chemotherapy and antimetabolic drugs exhibited growth inhibitory effects on TNBC cell lines. Antimetabolic drugs targeting the glycolysis pathway had a synergistic effect with chemotherapy drugs, and antiglycolysis drug combinations also had a synergistic effect. The use of these drug combinations could lead to new therapeutic strategies that reduce chemotherapy drug doses, decreasing their toxic effect, or that maintain the doses but enhance their efficacy by their synergistic effect with other drugs.

Identification of Carcinogenesis and Tumor Progression Processes in Pancreatic Ductal Adenocarcinoma Using High-Throughput Proteomics.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease with an overall 5-year survival rate of just 5%. A better understanding of the carcinogenesis processes and the mechanisms of the progression of PDAC is mandatory. Fifty-two PDAC patients treated with surgery and adjuvant therapy, with available primary tumors, normal tissue, preneoplastic lesions (PanIN), and/or lymph node metastases, were selected for the study. Proteins were extracted from small punches and analyzed by LC-MS/MS using data-independent acquisition. Proteomics data were analyzed using probabilistic graphical models, allowing functional characterization. Comparisons between groups were made using linear mixed models. Three proteomic tumor subtypes were defined. T1 (32% of patients) was related to adhesion, T2 (34%) had metabolic features, and T3 (34%) presented high splicing and nucleoplasm activity. These proteomics subtypes were validated in the PDAC TCGA cohort. Relevant biological processes related to carcinogenesis and tumor progression were studied in each subtype. Carcinogenesis in the T1 subtype seems to be related to an increase of adhesion and complement activation node activity, whereas tumor progression seems to be related to nucleoplasm and translation nodes. Regarding the T2 subtype, it seems that metabolism and, especially, mitochondria act as the motor of cancer development. T3 analyses point out that nucleoplasm, mitochondria and metabolism, and extracellular matrix nodes could be involved in T3 tumor carcinogenesis. The identified processes were different among proteomics subtypes, suggesting that the molecular motor of the disease is different in each subtype. These differences can have implications for the development of future tailored therapeutic approaches for each PDAC proteomics subtype.

Description of the genetic variants identified in a cohort of patients diagnosed with localized anal squamous cell carcinoma and treated with panitumumab.

Squamous cell carcinoma is the most frequent histologic type of anal carcinoma. The standard of care since the 1970s has been a combination of 5-fluorouracil, mitomycin C, and radiotherapy. This treatment is very effective in T1/T2 tumors (achieving complete regression in 80-90% of tumors). However, in T3/T4 tumors, the 3-year relapse free survival rate is only 50%. The VITAL trial aimed to assess the efficacy and safety of panitumumab in combination with this standard treatment. In this study, 27 paraffin-embedded samples from the VITAL trial and 18 samples from patients from daily clinical practice were analyzed by whole-exome sequencing and the influence of the presence of genetic variants in the response to panitumumab was studied. Having a moderate- or high-impact genetic variant in PIK3CA seemed to be related to the response to panitumumab. Furthermore, copy number variants in FGFR3, GRB2 and JAK1 were also related to the response to panitumumab. These genetic alterations have also been studied in the cohort of patients from daily clinical practice (not treated with panitumumab) and they did not have a predictive value. Therefore, in this study, a collection of genetic alterations related to the response with panitumumab was described. These results could be useful for patient stratification in new anti-EGFR clinical trials.

Bayesian networks established functional differences between breast cancer subtypes.

Breast cancer is a heterogeneous disease. In clinical practice, tumors are classified as hormonal receptor positive, Her2 positive and triple negative tumors. In previous works, our group defined a new hormonal receptor positive subgroup, the TN-like subtype, which had a prognosis and a molecular profile more similar to triple negative tumors. In this study, proteomics and Bayesian networks were used to characterize protein relationships in 96 breast tumor samples. Components obtained by these methods had a clear functional structure. The analysis of these components suggested differences in processes such as mitochondrial function or extracellular matrix between breast cancer subtypes, including our new defined subtype TN-like. In addition, one of the components, mainly related with extracellular matrix processes, had prognostic value in this cohort. Functional approaches allow to build hypotheses about regulatory mechanisms and to establish new relationships among proteins in the breast cancer context.

miRNA profiling in renal carcinoma suggest the existence of a group of pro-angionenic tumors in localized clear cell renal carcinoma.

Renal cell carcinoma comprises a variety of entities, the most common being the clear-cell, papillary and chromophobe subtypes. These subtypes are related to different clinical evolution; however, most therapies have been developed for clear-cell carcinoma and there is not a specific treatment based on different subtypes. In this study, one hundred and sixty-four paraffin samples from primary nephrectomies for localized tumors were analyzed. MiRNAs were isolated and measured by microRNA arrays. Significance Analysis of Microarrays and Consensus Cluster algorithm were used to characterize different renal subtypes. The analyses showed that chromophobe renal tumors are a homogeneous group characterized by an overexpression of miR 1229, miR 10a, miR 182, miR 1208, miR 222, miR 221, miR 891b, miR 629-5p and miR 221-5p. On the other hand, clear cell renal carcinomas presented two different groups inside this histological subtype, with differences in miRNAs that regulate focal adhesion, transcription, apoptosis and angiogenesis processes. Specifically, one of the defined groups had an overexpression of proangiogenic microRNAs miR185, miR126 and miR130a. In conclusion, differences in miRNA expression profiles between histological renal subtypes were established. In addition, clear cell renal carcinomas had different expression of proangiogenic miRNAs. With the emergence of antiangiogenic drugs, these differences could be used as therapeutic targets in the future or as a selection method for tailoring personalized treatments.

Comprehensive Characterization of the Mutational Landscape in Localized Anal Squamous Cell Carcinoma.

Anal squamous cell carcinoma (ASCC) is a rare neoplasm. Chemoradiotherapy is the standard of care, with no therapeutic advances achieved over the past three decades. Thus, a deeper molecular characterization of this disease is still necessary. We analyzed 46 paraffin-embedded tumor samples from patients diagnosed with primary ASCC by exome sequencing. A bioinformatics approach focused in the identification of high-impact genetic variants, which may act as drivers of oncogenesis, was performed. The relation between genetics variants and prognosis was also studied. The list of high-impact genetic variants was unique for each patient. However, the pathways in which these genes are involved are well-known hallmarks of cancer, such as angiogenesis or immune pathways. Additionally, we determined that genetic variants in BRCA2, ZNF750, FAM208B, ZNF599, and ZC3H13 genes are related with poor disease-free survival in ASCC. This may help to stratify the patient's prognosis and open new avenues for potential therapeutic intervention. In conclusion, sequencing of ASCC clinical samples appears an encouraging tool for the molecular portrait of this disease.

Proteomics characterisation of central nervous system metastasis biomarkers in triple negative breast cancer.

BACKGROUND: Breast cancer (BC) is the most frequent tumour in women. Triple negative tumours (TNBC)-which are associated with minor survival rates-lack markers predictive of response to anticancer drugs. Triple negative tumours frequently metastasise to the central nervous system (CNS). OBJECTIVE: The main objective of this study was to study differences in tumour protein expression between patients with CNS metastases and those without this kind of spread, and propose new biomarkers. METHODS: A retrospective study was performed. Targeted proteomics and statistical analyses were used to identify possible biomarkers. RESULTS: Proteins were quantified by a targeted proteomics approach and protein expression data were successfully obtained from 51 triple negative formalin-fixed paraffin-embedded samples. ISG15, THBS1 and AP1M1 were identified as possible biomarkers related with CNS metastasis development. CONCLUSIONS: Three possible biomarkers associated with CNS metastases in TNBC tumours were identified: ISG15, THBS1 and AP1M1. They may become markers predicting the appearance of CNS infiltration in triple negative BC.

A novel approach to triple-negative breast cancer molecular classification reveals a luminal immune-positive subgroup with good prognoses.

Triple-negative breast cancer is a heterogeneous disease characterized by a lack of hormonal receptors and HER2 overexpression. It is the only breast cancer subgroup that does not benefit from targeted therapies, and its prognosis is poor. Several studies have developed specific molecular classifications for triple-negative breast cancer. However, these molecular subtypes have had little impact in the clinical setting. Gene expression data and clinical information from 494 triple-negative breast tumors were obtained from public databases. First, a probabilistic graphical model approach to associate gene expression profiles was performed. Then, sparse k-means was used to establish a new molecular classification. Results were then verified in a second database including 153 triple-negative breast tumors treated with neoadjuvant chemotherapy. Clinical and gene expression data from 494 triple-negative breast tumors were analyzed. Tumors in the dataset were divided into four subgroups (luminal-androgen receptor expressing, basal, claudin-low and claudin-high), using the cancer stem cell hypothesis as reference. These four subgroups were defined and characterized through hierarchical clustering and probabilistic graphical models and compared with previously defined classifications. In addition, two subgroups related to immune activity were defined. This immune activity showed prognostic value in the whole cohort and in the luminal subgroup. The claudin-high subgroup showed poor response to neoadjuvant chemotherapy. Through a novel analytical approach we proved that there are at least two independent sources of biological information: cellular and immune. Thus, we developed two different and overlapping triple-negative breast cancer classifications and showed that the luminal immune-positive subgroup had better prognoses than the luminal immune-negative. Finally, this work paves the way for using the defined classifications as predictive features in the neoadjuvant scenario.

Melanoma proteomics suggests functional differences related to mutational status.

Melanoma is the most lethal cutaneous cancer. New drugs have recently appeared; however, not all patients obtain a benefit of these new drugs. For this reason, it is still necessary to characterize melanoma at molecular level. The aim of this study was to explore the molecular differences between melanoma tumor subtypes, based on BRAF and NRAS mutational status. Fourteen formalin-fixed, paraffin-embedded melanoma samples were analyzed using a high-throughput proteomics approach, combined with probabilistic graphical models and Flux Balance Analysis, to characterize these differences. Proteomics analyses showed differences in expression of proteins related with fatty acid metabolism, melanogenesis and extracellular space between BRAF mutated and BRAF non-mutated melanoma tumors. Additionally, probabilistic graphical models showed differences between melanoma subgroups at biological processes such as melanogenesis or metabolism. On the other hand, Flux Balance Analysis predicts a higher tumor growth rate in BRAF mutated melanoma samples. In conclusion, differential biological processes between melanomas showing a specific mutational status can be detected using combined proteomics and computational approaches.