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.

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.

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.

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.

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.

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.

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.

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 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.

Biological molecular layer classification of muscle-invasive bladder cancer opens new treatment opportunities.

BACKGROUND: Muscle-invasive bladder tumors are associated with a high risk of relapse and metastasis even after neoadjuvant chemotherapy and radical cystectomy. Therefore, further therapeutic options are needed and molecular characterization of the disease may help to identify new targets. The aim of this study was to characterize muscle-invasive bladder tumors at the molecular level using computational analyses. METHODS: The TCGA cohort of muscle-invasive bladder cancer patients was used to describe these tumors. Probabilistic graphical models, layer analyses based on sparse k-means coupled with Consensus Cluster, and Flux Balance Analysis were applied to characterize muscle-invasive bladder tumors at a functional level. RESULTS: Luminal and Basal groups were identified, and an immune molecular layer with independent value was also described. Luminal tumors showed decreased activity in the nodes of epidermis development and extracellular matrix, and increased activity in the node of steroid metabolism leading to a higher expression of the androgen receptor. This fact points to the androgen receptor as a therapeutic target in this group. Basal tumors were highly proliferative according to Flux Balance Analysis, which makes these tumors good candidates for neoadjuvant chemotherapy. The Immune-high group showed a higher degree of expression of immune biomarkers, suggesting that this group may benefit from immune therapy. CONCLUSIONS: Our approach, based on layer analyses, established a Luminal group candidate for therapy with androgen receptor inhibitors, a proliferative Basal group which seems to be a good candidate for chemotherapy, and an immune-high group candidate for immunotherapy.

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.

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.

Probabilistic graphical models relate immune status with response to neoadjuvant chemotherapy in breast cancer.

Breast cancer is the most frequent tumor in women and its incidence is increasing. Neoadjuvant chemotherapy has become standard of care as a complement to surgery in locally advanced or poor-prognosis early stage disease. The achievement of a complete response to neoadjuvant chemotherapy correlates with prognosis but it is not possible to predict who will obtain an excellent response. The molecular analysis of the tumor offers a unique opportunity to unveil predictive factors. In this work, gene expression profiling in 279 tumor samples from patients receiving neoadjuvant chemotherapy was performed and probabilistic graphical models were used. This approach enables addressing biological and clinical questions from a Systems Biology perspective, allowing to deal with large gene expression data and their interactions. Tumors presenting complete response to neoadjuvant chemotherapy had a higher activity of immune related functions compared to resistant tumors. Similarly, samples from complete responders presented higher expression ​​of lymphocyte cell lineage markers, immune-activating and immune-suppressive markers, which may correlate with tumor infiltration by lymphocytes (TILs). These results suggest that the patient's immune system plays a key role in tumor response to neoadjuvant treatment. However, future studies with larger cohorts are necessary to validate these hypotheses.

Molecular characterization of breast cancer cell response to metabolic drugs.

Metabolic reprogramming is a hallmark of cancer. It has been described that breast cancer subtypes present metabolism differences and this fact enables the possibility of using metabolic inhibitors as targeted drugs in specific scenarios. In this study, breast cancer cell lines were treated with metformin and rapamycin, showing a heterogeneous response to treatment and leading to cell cycle disruption. The genetic causes and molecular effects of this differential response were characterized by means of SNP genotyping and mass spectrometry-based proteomics. Protein expression was analyzed using probabilistic graphical models, showing that treatments elicit various responses in some biological processes such as transcription. Moreover, flux balance analysis using protein expression values showed that predicted growth rates were comparable with cell viability measurements and suggesting an increase in reactive oxygen species response enzymes due to metformin treatment. In addition, a method to assess flux differences in whole pathways was proposed. Our results show that these diverse approaches provide complementary information and allow us to suggest hypotheses about the response to drugs that target metabolism and their mechanisms of action.

Urothelial cancer proteomics provides both prognostic and functional information.

Traditionally, bladder cancer has been classified based on histology features. Recently, some works have proposed a molecular classification of invasive bladder tumors. To determine whether proteomics can define molecular subtypes of  muscle invasive urothelial cancer (MIUC) and allow evaluating the status of biological processes and its clinical value. 58 MIUC patients who underwent curative surgical resection at our institution between 2006 and 2012 were included. Proteome was evaluated by high-throughput proteomics in routinely archive FFPE tumor tissue. New molecular subgroups were defined. Functional structure and individual proteins prognostic value were evaluated and correlated with clinicopathologic parameters. 1,453 proteins were quantified, leading to two MIUC molecular subgroups. A protein-based functional structure was defined, including several nodes with specific biological activity. The functional structure showed differences between subtypes in metabolism, focal adhesion, RNA and splicing nodes. Focal adhesion node has prognostic value in the whole population. A 6-protein prognostic signature, associated with higher risk of relapse (5 year DFS 70% versus 20%) was defined. Additionally, we identified two MIUC subtypes groups. Prognostic information provided by pathologic characteristics is not enough to understand MIUC behavior. Proteomics analysis may enhance our understanding of prognostic and classification. These findings can lead to improving diagnosis and treatment selection in these patients.

Functional proteomics outlines the complexity of breast cancer molecular subtypes.

Breast cancer is a heterogeneous disease comprising a variety of entities with various genetic backgrounds. Estrogen receptor-positive, human epidermal growth factor receptor 2-negative tumors typically have a favorable outcome; however, some patients eventually relapse, which suggests some heterogeneity within this category. In the present study, we used proteomics and miRNA profiling techniques to characterize a set of 102 either estrogen receptor-positive (ER+)/progesterone receptor-positive (PR+) or triple-negative formalin-fixed, paraffin-embedded breast tumors. Protein expression-based probabilistic graphical models and flux balance analyses revealed that some ER+/PR+ samples had a protein expression profile similar to that of triple-negative samples and had a clinical outcome similar to those with triple-negative disease. This probabilistic graphical model-based classification had prognostic value in patients with luminal A breast cancer. This prognostic information was independent of that provided by standard genomic tests for breast cancer, such as MammaPrint, OncoType Dx and the 8-gene Score.

Prediction of adjuvant chemotherapy response in triple negative breast cancer with discovery and targeted proteomics.

BACKGROUND: Triple-negative breast cancer (TNBC) accounts for 15-20% of all breast cancers and usually requires the administration of adjuvant chemotherapy after surgery but even with this treatment many patients still suffer from a relapse. The main objective of this study was to identify proteomics-based biomarkers that predict the response to standard adjuvant chemotherapy, so that patients at are not going to benefit from it can be offered therapeutic alternatives. METHODS: We analyzed the proteome of a retrospective series of formalin-fixed, paraffin-embedded TNBC tissue applying high-throughput label-free quantitative proteomics. We identified several protein signatures with predictive value, which were validated with quantitative targeted proteomics in an independent cohort of patients and further evaluated in publicly available transcriptomics data. RESULTS: Using univariate Cox analysis, a panel of 18 proteins was significantly associated with distant metastasis-free survival of patients (p<0.01). A reduced 5-protein profile with prognostic value was identified and its prediction performance was assessed in an independent targeted proteomics experiment and a publicly available transcriptomics dataset. Predictor P5 including peptides from proteins RAC2, RAB6A, BIEA and IPYR was the best performance protein combination in predicting relapse after adjuvant chemotherapy in TNBC patients. CONCLUSIONS: This study identified a protein combination signature that complements histopathological prognostic factors in TNBC treated with adjuvant chemotherapy. The protein signature can be used in paraffin-embedded samples, and after a prospective validation in independent series, it could be used as predictive clinical test in order to recommend participation in clinical trials or a more exhaustive follow-up.