Metabolomics unravels subtype-specific characteristics related to neoadjuvant therapy response in breast cancer patients.

INTRODUCTION: Breast cancer is the most diagnosed tumor and the leading cause of cancer death in women worldwide. Metabolomics allows the quantification of the entire set of metabolites in blood samples, making it possible to study differential metabolomics patterns related to neoadjuvant treatment in the breast cancer neoadjuvant setting. OBJECTIVES: Characterizing metabolic differences in breast cancer blood samples according to their response to neoadjuvant treatment. METHODS: One hundred and three plasma samples of breast cancer patients, before receiving neoadjuvant treatment, were analyzed through UPLC-MS/MS metabolomics. Then, metabolomics data were analyzed using probabilistic graphical models and biostatistics methods. RESULTS: Metabolomics data allowed the identification of differences between groups according to response to neoadjuvant treatment. These differences were specific to each breast cancer subtype. Patients with HER2+ tumors showed differences in metabolites related to amino acids and carbohydrates pathways between the two pathological response groups. However, patients with triple-negative tumors showed differences in metabolites related to the long-chain fatty acids pathway. Patients with Luminal B tumors showed differences in metabolites related to acylcarnitine pathways. CONCLUSIONS: It is possible to identify differential metabolomics patterns between complete and partial responses to neoadjuvant therapy, being this metabolomic profile specific for each breast cancer subtype.

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.

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.

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.

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.

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.

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.

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.

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.