Next Article in Journal
NSCLC Biomarkers to Predict Response to Immunotherapy with Checkpoint Inhibitors (ICI): From the Cells to In Vivo Images
Next Article in Special Issue
Metabolic Flexibility Is a Determinant of Breast Cancer Heterogeneity and Progression
Previous Article in Journal
Facile Generation of Potent Bispecific Fab via Sortase A and Click Chemistry for Cancer Immunotherapy
Previous Article in Special Issue
Targeting Amino Acid Metabolic Reprogramming via L-Type Amino Acid Transporter 1 (LAT1) for Endocrine-Resistant Breast Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Multi-Omic Approaches to Breast Cancer Metabolic Phenotyping: Applications in Diagnosis, Prognosis, and the Development of Novel Treatments

by
Nuria Gómez-Cebrián
1,
Inés Domingo-Ortí
1,2,
José Luis Poveda
3,
María J. Vicent
2,
Leonor Puchades-Carrasco
1,* and
Antonio Pineda-Lucena
4,*
1
Drug Discovery Unit, Instituto de Investigación Sanitaria La Fe, 46026 Valencia, Spain
2
Polymer Therapeutic Laboratory, Centro de Investigación Príncipe Felipe, 46012 Valencia, Spain
3
Pharmacy Department, Hospital Universitario y Politécnico La Fe, 46026 Valencia, Spain
4
Molecular Therapeutics Program, Centro de Investigación Médica Aplicada, 31008 Navarra, Spain
*
Authors to whom correspondence should be addressed.
Cancers 2021, 13(18), 4544; https://doi.org/10.3390/cancers13184544
Submission received: 19 July 2021 / Revised: 1 September 2021 / Accepted: 8 September 2021 / Published: 10 September 2021
(This article belongs to the Special Issue Understanding and Modelling Metabolic Reprogramming in Breast Cancer)

Abstract

:

Simple Summary

Breast cancer (BC) is a heterogeneous tumor type and has become the leading cause of cancer worldwide, with 685,000 deaths forecast in 2020. The clinical management of BC patients remains challenging, and there exists an urgent need for improved diagnostic, prognostic, and therapeutic strategies. Multi-omics platforms represent a promising tool for discovering novel biomarkers and identifying new therapeutic targets. In addition, the ongoing development of multi-omics approaches may foster the identification of more robust and accurate algorithms for data analysis. This review aims to summarize the results of recent multi-omics-based studies focused on the characterization of the metabolic phenotype of BC.

Abstract

Breast cancer (BC) is characterized by high disease heterogeneity and represents the most frequently diagnosed cancer among women worldwide. Complex and subtype-specific gene expression alterations participate in disease development and progression, with BC cells known to rewire their cellular metabolism to survive, proliferate, and invade. Hence, as an emerging cancer hallmark, metabolic reprogramming holds great promise for cancer diagnosis, prognosis, and treatment. Multi-omics approaches (the combined analysis of various types of omics data) offer opportunities to advance our understanding of the molecular changes underlying metabolic rewiring in complex diseases such as BC. Recent studies focusing on the combined analysis of genomics, epigenomics, transcriptomics, proteomics, and/or metabolomics in different BC subtypes have provided novel insights into the specificities of metabolic rewiring and the vulnerabilities that may guide therapeutic development and improve patient outcomes. This review summarizes the findings of multi-omics studies focused on the characterization of the specific metabolic phenotypes of BC and discusses how they may improve clinical BC diagnosis, subtyping, and treatment.

1. Introduction

Breast cancer (BC) is the most frequently diagnosed tumor and the leading cause of cancer deaths in women worldwide [1]. The year 2020 saw an estimated 2.3 million new cases of BC (11.7% of all cancer cases), with 685,000 deaths worldwide [1]; however, advances in population screening and early treatment (among other factors) have supported a steady decrease in BC mortality [2,3]. Unfortunately, figures from the American Cancer Society place the five-year survival rate after diagnosis of metastatic BC at 27%, a low value considering the 99% five-year survival rate for localized disease [4]. Therefore, early BC detection represents a crucial step in reducing disease mortality [5]. BC screening currently relies on mammography, a non-invasive strategy primarily performed in women between 50 to 69 years of age that has prompted a reduction in BC mortality [6]. Nevertheless, this approach suffers from several limitations, including false-positive reporting and overdiagnosis [6,7,8,9,10]. Ultrasound, magnetic resonance imaging, and computed tomography can overcome such problems thanks to their high sensitivity; however, the elevated costs associated with these tools make this approach less accessible. Thus, we still lack alternative methods for the accurate, non-invasive, and low-cost diagnosis of early-stage BC.
BC is a highly heterogeneous disease from a molecular perspective and is primarily characterized by the overexpression of the HER2 growth factor, estrogen receptor (ER), and progesterone receptor (PR) and mutations in the BRCA1/2 genes, with the latter associated with a higher risk of developing BC [11]. Former classifications of BC tumors employed tumor size, histological grade, immunohistochemistry of ER/PR status, and the amplification of HER2. The addition of gene expression profiling to these molecular features has resulted in the classification currently used by the European Society for Medical Oncology as a clinical guideline for BC diagnosis, follow-up, and treatment [12,13]. This system classifies BC tumors into four major intrinsic molecular subgroups: luminal A (ER+ and/or PR+, HER2-, low Ki67), luminal B (ER+ and/or PR+, HER2+ or HER2- with high Ki67), basal-like (ER/PR-, HER2-), and HER2-enriched (ER/PR-, HER2+). Luminal A tumors (the low-grade group) are the most common BC subtype, comprising over 60–70% of all cases. Meanwhile, basal-like tumors, with an ~80% overlap with highly proliferative triple-negative breast cancer (TNBC) [14], exhibit aggressive behavior and suffer from poor prognosis [15]. Although ER+ tumors present lower recurrence rates within the first five years, over 50% of tumor recurrences occur after this time and cause most BC-related deaths [16,17]. Each BC subtype has a characteristic biological profile, prognosis, and treatment strategy [18,19,20,21,22], and several scoring systems aid prognosis and treatment decision-making processes. Unfortunately, systems based on different molecular features related to tumor biology, including histological type, grade, lymphovascular invasion, and marker status, do not accurately reflect BC subtype heterogeneity or specific patient subtypes [23]. Thus, enormous efforts have been devoted to classifying heterogeneous BC subtypes into molecular subtypes that guide treatment decisions [24,25,26,27,28].
Current BC treatment strategies are applied according to BC subtype and prognosis. BC subtypes expressing hormone receptors (Luminal A, Luminal B, and HER2+) associate with good/intermediate prognosis, and their clinical management includes endocrine therapy alone or in combination with chemotherapy if there exists the risk of recurrence. Non-metastatic BC treatment also includes surgical resection and postoperative radiation therapy. Nevertheless, we lack effective therapeutic options for TNBC, the most aggressive BC subtype [29,30,31], and we urgently require the development of novel treatment strategies specifically targeting TNBC tumors.
Importantly, a range of studies employing multi-omics-based approaches, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics, have noted the existence of characteristic metabolic profiles for BC subtypes. Metabolomics-based studies performed in BC preclinical models [32,33] and patient tissues [34,35] have identified alterations to metabolites associated with glutamine, glutathione, and choline (Cho) metabolism. Similar studies have been carried out in biofluids such as urine [36,37], blood [38,39], and saliva [40], while other studies focusing on the metabolic characterization of normal breast epithelial and BC cell lines have revealed specific metabolomic differences between TNBC, Luminal B, and HER2+ subtypes [41].
Alterations to metabolic pathways associated with BC tumors and disease progression have been broadly explored at the genomic level [42,43,44]. Transcriptomic-based approaches have characterized specific metabolic changes in different BC subtypes [41] and evaluated the impact of various treatments on metabolic processes [45]. Furthermore, proteomics studies have revealed alterations in metabolism-associated protein expression in BC patients and correlations to overall and recurrence-free survival [46].
The combination of omic approaches has recently emerged as a promising strategy for generating information on the interrelationship between genetic aberrations, epigenetic alterations, changes in gene transcription and signaling pathways, and metabolic alterations that may contribute to the progression of a specific condition (Figure 1).
Interestingly, a range of studies has applied multi-omics approaches to the study of BC metabolism and has provided a more accurate understanding of the disease and its sub-classification [47,48,49]. These studies enabled the identification of characteristic disease biomarkers [50] or metabolic pathways related to specific subgroups of BC patients [51]. This review summarizes the results obtained in recently published multi-omics-based studies focused on analyzing those metabolic changes related to BC progression.

2. BC Diagnosis

Potential metabolic biomarkers for the detection and stratification of BC patients have been described in genomic [42], epigenomic [52,53], transcriptomic [54], proteomic [55], and metabolomic [34,35,36,37,38,39,40,41,56,57,58] studies. However, multi-omics-based analyses have been carried out in a reduced number of studies, which probably reflects limited sample availability and/or technical difficulties associated with generating complete multi-omics datasets. Thirteen studies published between 2010 and 2021 have integrated two different omics datasets to identify disease-associated metabolic alterations that could improve the diagnosis and subclassification of BC patients [48,57,59,60,61,62,63,64,65,66,67,68,69]. Tissue samples represented the most frequently used sample type, with only three studies collecting biofluids (primarily blood). Most cases evaluated a combination of transcriptomics and metabolomics, while one study included lipidomics, a recently emerged discipline in the omics field [70].

2.1. BC Metabolic Markers

The identification of novel biomarkers useful for early detection remains a critical clinical challenge in BC. A range of studies has compared the metabolic profiles of BC patients and healthy individuals using multi-omics approaches to identify potential metabolic markers with clinical utility in early BC diagnosis (Table 1). Integration of transcriptomics and metabolomics datasets represent the primary strategy followed in these studies, with tissue samples the preferable choice for analysis. Overall, the characterization of metabolic differences between BC and healthy tissues has revealed alterations in metabolites and/or enzymes involved in glycolysis, amino acid, lipid, and nucleotide metabolism.
A multi-omics study conducted by Iqbal et al. described the role of chromobox (CBX) proteins in BC metabolic reprogramming [59]. The authors identified a significant correlation between glycolytic activity and changes in CBX2 and CBX7 isoforms levels. Glycolytic metabolites, including glucose-6-phosphate, 3-phosphoglycerate, and fructose-6-phosphate, exhibited upregulated expression in high CBX2 and low CBX7 tumor samples. A study focused on the characterization of specific dysregulations in lipid metabolism conducted by Hilvo et al. reported increased levels of phospholipids, sphingomyelins, and ceramides, and the overexpression of four enzymes involved in de novo lipid metabolism (acetyl-CoA carboxylase alpha (ACACA), fatty acid synthase (FASN), insulin-induced gene 1 (INSIG1), and sterol regulatory element-binding transcription factor 1 (SREBP1)) in BC compared to healthy tissues [60]. Enhanced phospholipid metabolism in BC tissues has been previously observed, reflected by high phosphatidylcholine (PtdCho), phosphatidylethanolamine, or phosphocholine (PCho) levels [71,72,73]. In addition, reports suggest that certain of the above-mentioned lipid metabolism enzymes play a significant role in BC progression. For instance, acetyl-CoA carboxylase alpha (ACACA) overexpression associates with increased cancer risk [74] and can be detected in early-stage BC [75], while fatty acid synthase (FASN) promoted BC metastasis [76,77] and was proposed as a potential therapeutic target for BC treatment [78].
Additional multi-omics studies have focused on analyzing metabolic alterations in plasma samples from BC patients and correlating them with transcriptomic changes in BC tumors. Luo and coworkers identified significant changes in purine metabolism when combining data from transcriptomics and metabolomics analyses of BC [61]. Enrichment-pathway analysis demonstrated the lower expression of genes related to purine nucleotides degradation. Notably, the authors also reported the downregulation of purine nucleoside phosphorylase (PNP) and hypoxanthine phosphoribosyltransferase (HPRT1) expression, two enzymes involved in the purine salvage pathway, as alterations in this pathway. Furthermore, the authors noted lower guanine and hypoxanthine levels in BC tissues. Although studies have suggested that de novo nucleotide synthesis can promote BC metastasis [79], recent studies have reported the downregulation of genes involved in salvage synthesis concurrent with the upregulated expression of nucleotide de novo synthesis genes in metastatic BC cells [79]. Related studies suggested the overexpression of ecto-5’-nucleotidase (CD73, which catalyzes the dephosphorylation of adenosine monophosphates into adenosine) in BC cells and an associated modulation of tumor growth and metastasis [80,81]; meanwhile, the observed increased activity of adenosine deaminase (ADA) in BC tumors [82,83,84] may represent a diagnostic tool for BC [83,84]. Additionally, enzymes involved in purine metabolism have also been proposed as potential drug targets for BC, and these include the ribonucleotide reductase regulatory subunit M2 (RRM2), natriuretic peptide receptor 1 (NPR1,) and phosphodiesterases (PDE) enzymes [54,85].
Multi-omics approaches have also been explored in the development of novel biomarker combinations that discriminate healthy tissue from BC at different stages [62]. In a multi-omics study aiming to discover crucial metabolic pathway signatures for more accurate BC diagnosis, Huang et al. reported an eight pathway-based classification model that provided highly accurate predictions in all-stage and early-stage BC patients [57]. Taurine and hypotaurine metabolism and the alanine, aspartate, and glutamate pathway were identified as critical to the early diagnosis of BC, with consistent changes occurring in plasma and serum metabolomics analyses and the transcriptomics-based analyses of BC tumors and adjacent normal tissue. The study also reported significant correlations between alterations to metabolite concentrations and gene expression levels when combining both omics approaches. For instance, increased asparagine, cysteine, and hypotaurine concentrations in early-stage BC correlated with the overexpression of the relevant synthetic enzymes, while lower levels of oxoglutarate associated with a significant decrease in the expression of the enzymes involved in its generation. Recent metabolomics studies have also reported altered taurine levels in plasma and serum samples of BC patients [86] and altered asparagine and cysteine levels in blood [87,88] and urine [36] of BC patients.
Dowling et al. evaluated protein and metabolite levels in BC serum samples using high-throughput antibody-based screening and mass spectrometry platforms, respectively [62]. A combined analysis demonstrated that the addition of glutamate, 12-hydroxyeicosatetraenoic acid, β-hydroxybutyrate, coagulation factor V, and matrix metalloproteinase-1 to cancer antigen 15-3 (CA15-3) values, a widely used serum marker in BC, permitted accurate discrimination when comparing non-malignant breast disease and BC at different stages. This combination resulted in noticeably improved prediction values compared with a classification based only on CA15-3 alone, confirming the potential of multi-omics approaches to provide clinically relevant information.

2.2. BC Metabolic Subtyping

Multi-omics approaches have also shown great potential for redefining cancer subtypes [89,90]. Recent studies have focused on exploiting the information derived from multi-omics strategies to improve tumor subtyping (Table 2). Most of these studies, which relied on tissue sample analysis and the integration of transcriptomics and metabolomics data, focused on the characterization of ER- vs. ER+ BC tumors and basal- vs. luminal-like BC tumors, while a few cases explored unsupervised BC subtyping. Overall, metabolic alterations in glycolysis, the tricarboxylic acid (TCA) cycle, the pentose phosphate pathway (PPP), and lipid, nucleotide, and amino acid metabolism distinguished different BC tumor subtypes.
The results of different multi-omics studies reveal common alterations to lipid, nucleotide, and amino acid metabolism when comparing basal- and luminal-like BC subtypes, and Moestue and coworkers focused on differences in Cho metabolism profiles between these two BC subtypes [63]. In this study, basal-like BC models displayed higher concentrations of glycerophosphocholine (GPC) than PCho and higher glycine concentrations, which could be explained by lower choline kinase (CHKA, CHKB) expression and higher PtdCho degradation mediated by the overexpression of phospholipase A2 group 4A (PLA2G4A) and phospholipase B1 (PLB1). This group of tumors also exhibited a characteristic metabolic shift from PtdCho synthesis to glycine formation mediated by choline dehydrogenase (CHDH) and sarcosine dehydrogenase (SARDH) expression. A related metabolomics study established increased GPC and decreased PCho concentration levels in TNBC samples compared with luminal A subtypes [35]. In addition, BC patients exhibiting better prognoses exhibited lower concentrations of glycine than those with worse prognoses [91]. In another study, Grinde and colleagues discovered that the regulation of Cho metabolism varied between different BC molecular subgroups [64]. The authors analyzed the metabolomic and transcriptomic profiles of thirty-four patient-derived xenograft BC models to explore the differences between basal- and luminal-like subgroups and revealed a higher PCho/GPC ratio in luminal B tumors. The Cho metabolic profiles and the expression of genes involved in Cho metabolism, including CHKA and glycerophosphodiester phosphodiesterase domain containing 5 (GDPD5), displayed differences in luminal B and basal-like xenografts. Overall, the results of these and other studies highlight the heterogeneity in Cho metabolism in different BC molecular subtypes [92,93].
The luminal- and basal-like BC subgroups also display additional characteristic metabolic differences. Putluri et al. used an integrated analysis of the metabolome and transcriptome of these subgroups to identify eleven differentially expressed metabolic pathways between the luminal- and basal-like BC subgroups, including nucleotide (purine and pyrimidine) synthesis and different amino acid metabolism pathways, such as lysine degradation or branched-chain amino acid (BCAA) and glutamate metabolism [65]. The increased glutamate and nucleotide metabolism observed in basal-like tumors agree with previous studies describing a higher susceptibility to glutaminolysis-targeting therapies in TNBC [94,95] and an enhanced de novo nucleotide synthesis in metastatic BC [79]. Enhanced pyrimidine metabolism in TNBC has been reported [35], and inhibition of this metabolic route represents a potential therapeutic strategy to promote the sensitivity of TNBC cells to chemotherapy by exacerbating DNA damage [96].
Metabolic differences between basal and luminal-like BC phenotypes were recently explored by Mahendralingam and colleagues in a study based on a combined proteomic and transcriptomic analysis of normal mammary epithelial cells (MECs) [66]. The authors identified distinct metabolic phenotypes corresponding to basal, luminal progenitor, and mature luminal cells. Specifically, basal-like cells exhibited a glycolytic phenotype characterized by the abundance of glycolytic enzymes (phosphofructokinase (muscle) (PFKM), aldolase, fructose-bisphosphate C (ALDOC), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and pyruvate kinase M1/2 (PKM)), while luminal progenitors exhibited enhanced oxidative phosphorylation (OXPHOS). Bioinformatics analyses indicated that BC subtypes retain metabolic characteristics of their cell of origin. These findings agree with recent multi-omics studies that observed increased OXPHOS activity in TNBC, whose cells of origin are luminal progenitors [69].
Multiple studies have revealed relevant metabolic differences between the ER- and ER+ BC subtypes using multi-omics approaches. Using a combination of genomics, transcriptomics, and metabolomics, Tang and coworkers reported correlations between gene mutations/gene expression levels and tumor metabolism [67]. Cancers with TP53 alterations exhibited decreased levels of lipid glycerophosphocholines, tumors with ERBB2 amplification associated with changes in docosapentaenoate, fucose, and 1-oleoylglycerophosphoethanolamine levels, and PIK3CA mutations correlated with altered malonylcarnitine levels. Among other significant associations, indoleamine 2,3-dioxygenase 1 (IDO1) overexpression correlated with increased kynurenine levels. IDO1 functions by generating kynurenine from tryptophan, and a higher kynurenine/tryptophan ratio was observed during the progression of tumors, including BC [97,98]. Overall, Tang et al. established an increase in eight metabolites in ER+ tumors (including carnitine derivates and short- and medium-chain fatty acids) and an increase in long-chain fatty acid and monoacylglycerol levels in ER- tumors, indicating that differences in lipolysis and fatty acid oxidation (FAO) correlate with hormone receptor status [67]. Furthermore, higher BRCA1 expression levels associated with changes in a set of metabolites reflecting FAO activation in these tumors. FAO-related genes appear to promote cell proliferation and survival [99,100], and FAO inhibition has been proposed as a potential therapeutic strategy for BC treatment [101,102,103].
In a more recent study, Barupal et al. used reactome pathway mapping analyses to associate alterations to glycolysis, PPP, TCA cycle, nucleotide salvage, glutathione conjugation, steroid metabolism, fatty acyl-CoA biosynthesis, serine biosynthesis, and metabolism of aromatic amino acids with the ER- phenotype [68]. Increased PPP activity represented the most significant change supported by the multi-omics analysis. The PPP intermediates ribose-5-phosphate and ribulose-5-phosphate increased in ER- tumors and associated with increased protein and mRNA expression for the critical oxidation enzymes glucose-6-phosphate dehydrogenase (G6PD) and 6-phosphogluconate dehydrogenase (PGD). Furthermore, ER- tumors also displayed increased levels of transketolase (TKT), a key enzyme in the non-oxidative branch of the PPP, phosphoglucomutase 1 (PGM1), ribose-5-phosphate isomerase (RPIA), and deoxyribose-phosphate aldolase (DERA). The PPP delivers the nicotinamide adenine nucleotide phosphate (NADPH) and pentose phosphates required for nucleotide and fatty acid synthesis during cell division and tumor proliferation [104,105]. Consistent with PPP activation, the study discovered elevated levels of five purine metabolites (adenine, guanosine, guanine, xanthine, and hypoxanthine). In addition, β-alanine (an intermediate of the pyrimidine salvage pathway) and uracil, pseudo-uridine, uridine monophosphate (UMP), and cytidine monophosphate (CMP) levels also displayed increased expression in ER- tumors [68]. Interestingly, Budczies et al. previously described β-alanine as a significant marker differentiating ER- and ER+ BC [106]. Additional studies reported significant increases in lipid content, including phospholipids, in ER- compared with ER+ tumors [60,107]. Based on the analysis of lipidomics and gene expression data from 257 human BC tissues, Hilvo et al. associated alterations in phospholipid metabolism with ER status and tumor grade [60]. In this study, both ER- and high-grade tumors exhibited increased levels of phospholipids that incorporated de novo synthesized palmitate and myristic acid, while only tumor grade (and not ER- status) associated with changes in lipid-related gene expression.
Related studies employed integrated approaches to metabolically identify potential BC subgroups following an unsupervised approach. Haukaas et al. applied hierarchical clustering on 228 primary BC samples and identified three different metabolic subgroups (MC1, MC2, and MC3) [48]. Integrated pathway analysis of metabolite and gene expression data uncovered differences in glycolysis/gluconeogenesis and glycerophospholipid metabolism between the clusters. MC1 exhibited a lipidomic phenotype with higher levels of GPC and PCho; MC2 displayed a low glycolytic profile exhibiting increased glucose concentrations; and MC3 exhibited elevated lactate and alanine levels. Compared with the other two subgroups, MC1 possessed increased levels of GPC and PCho due to the downregulation of genes involved in PtdCho degradation and higher CHKA levels. In addition, the MC1 cluster exhibited lower acetate and glutamine levels due to the lower expression levels of aldehyde dehydrogenase (ALDH) and glutaminase (GLS) enzymes, respectively. However, the study failed to encounter differences in mRNA levels between MC2 and MC3 subgroups.
More recently, Gong et al. applied a similar approach to identify metabolic pathway-based subtypes (MPSs) in TNBC and identified three different metabolic phenotypes (lipogenic, glycolytic, and mixed subtypes) with distinct molecular features and sensitivities to various metabolic inhibitors [69]. In particular, the lipogenic phenotype (MPS1) associated with the high expression of genes involved in cholesterol and de novo lipid metabolism and the accumulation of several lipids, while the glycolytic group (MPS2) exhibited the upregulation of genes involved in glycolysis and increased concentrations of metabolic intermediates of glycolysis and nucleotide metabolism. The authors also evaluated the mRNA expression of metabolic enzymes and corresponding metabolite abundance to assess their impact on the MPS classification. Tumors corresponding to the MPS1 subtype displayed increased expression of acetyl-CoA carboxylase alpha (ACACA), HMG-CoA reductase (HMGCR), FASN, and stearoyl-CoA desaturase (SCD), and a significant accumulation of various lipids, such as myristic, palmitoleic, oleic, and arachidonic acid. In contrast, the MSP2 subtype displayed the elevated expression of metabolic enzymes associated with glycolysis and nucleotide metabolism, with the most significant changes associated with phosphofructokinase (PFKP), enolase 2 (ENO2), thymidylate synthetase (TYMS), CTP synthase 1 (CTPS1), but also with glucose transporter (GLUD), solute carrier family 2 member 1 (SLC2A1), and lactate transporter SLC16A1. The MSP2 subtype also displayed lower levels of glucose and the accumulation of intermediates in glycolysis and nucleotide metabolism, including glucose-1-phosphate, dihydroxyacetone phosphate, lactate, and adenosine-3′-5′-diphosphate. MPS3, the mixed subtype, showed only partial pathway dysregulation.

3. Multi-Omics Studies of BC Prognosis

Metabolic deregulation can impact various molecular processes (e.g., cell proliferation, apoptosis, migration, and invasion) that contribute to tumor progression [44,108,109,110] and influence cancer patient survival (Figure 2) [111].
Several studies have reported associations between metabolic alterations and BC patient survival based on single omic analysis, including genomic [112], transcriptomic [113], proteomic [46], and metabolomic [56,114,115] studies in tissue [46,112,113] and serum [56,114,115] samples. However, additional studies based on the integration of data from multi-omics analyses have provided more accurate information regarding the molecules involved in metabolic rewiring associated with BC progression. Of the seven studies following a multi-omics approach to identify metabolic alterations associated with BC prognosis (Table 3), most relied on the analysis of tissue sample analysis and the integration of transcriptomic and metabolomic datasets. Overall, BC patient survival associated with the altered expression of enzymes involved in nucleotide, lipid, and amino acid metabolism.
Previous studies have demonstrated that nucleotide biosynthesis plays a vital role in BC [79,80,83,84,117], and could represent a promising therapeutic strategy [54,85,118]. Notably, two of the multi-omics studies included in Table 3 established an inverse correlation between the expression levels of genes involved in de novo purine and pyrimidine syntheses and BC patient survival [61,65]. Putluri et al. performed an in silico analysis to evaluate the association between omics-based enrichment and patient survival, using ten independent gene expression data sets to select clinically relevant prognostic biomarkers [65]. Kaplan–Meyer curves revealed an association between increased expression of pyrimidine metabolism-related genes and shorter metastasis-free survival across all BC and within the subset of ER + tumors. RRM2, a critical gene in pyrimidine metabolism displaying elevated expression in aggressive BC [119], has prognostic relevance in BC [120,121] when combined with proliferation markers [122]. RRM2 expression distinguished good vs. poor survival within the entire BC patient group in this multi-omics study, including a significant proportion of luminal A subtype typically considered to have better survival outcomes [80]. Luo and colleagues integrated metabolomic and transcriptomic analysis to confirm the association between alterations in nucleotide metabolism and BC patient survival in a TCGA cohort of patients [61]. The authors observed a significant correlation between poor survival of BC patients and changes to the expression of RRM2 and adenosine monophosphate deaminase 1 (AMPD1), a key enzyme in de novo purine synthesis. These enzymes have been postulated as promising therapeutic targets in different tumor types, including BC [65,85,123,124,125,126].
Iqbal et al. established antagonistic roles of CBX2 and CBX7 in metabolic reprogramming of BC and an association with BC patient survival [59]. The authors described a significant correlation between higher CBX2 and lower CBX7 mRNA levels and worse BC prognosis, which agrees with previous findings that correlated the CBX2 or CBX7 expression with overall patient survival [127,128,129,130].
As for the alterations in lipid metabolism, Camarda and colleagues followed a targeted metabolomics approach and reported the dramatic upregulation of FAO intermediates in a MYC-driven model of TNBC [102]. To characterize a potential association between FAO gene expression and prognosis in TNBC, the authors performed a univariate analysis of 336 fatty acid metabolism genes on a patient cohort with long-term distant recurrence-free survival data. The analysis revealed that decreased ACACB (acetyl-CoA carboxylase 2, ACC2) expression levels associated with worse prognoses in all BC and TNBC patients. Subsequent studies also described significant associations between increased levels of ACC2 and better BC prognosis [131,132,133,134]. Kang et al. conducted a multi-layered lipidomics and transcriptomics analysis to describe the rewiring of the BC lipidome during malignant transformation [116]. Analyses in a spheroid-induced epithelial-mesenchymal transition (EMT) model demonstrated a dramatic reduction in the ratio of C22:6n3 (docosahexaenoic acid, DHA) to C22:5n3 in spheroid cells, similarly to the down-regulation of ELOVL2, a process associated with the induction of metastatic characteristics in BC cells. The authors examined the relationship between ELOVL2 expression and metastatic relapse-free in a BC cohort with a follow-up of ten years, resulting in the discovery of an association between lower ELOVL2 expression levels and shorter metastasis-free survival and higher tumor grade [116]. A recent study investigating the molecular mechanisms of tamoxifen resistance in BC confirmed these findings and described lower ELOVL2 expression in tamoxifen-resistant models and ELOVL2 downregulation in patients with tamoxifen resistance [135].
Additional studies have revealed a correlation between alterations in the levels of genes and metabolites involved in amino acid metabolism and BC prognosis. Terunuma and coworkers identified a subset of BC tumors accumulating high levels of 2-hydroxyglutarate (2HG). Further analyses revealed the presence of a subgroup of BC patients with significantly decreased survival characterized by an exceptionally high accumulation of 2HG, reduced DNA methylation at the isocitrate dehydrogenase (IDH2) locus, increased IDH2 expression, and increased levels of S-adenosyl- methionine (SAM) and S-adenosylhomocysteine (SAH) [107]. Previous studies reported an accumulation of the oncometabolite 2HG in different tumor types, including BC [67], glioma [136], and leukemia [137]. In another study, concerning amino acid metabolism, Budczies et al. described an association between alterations in the metabolism of β-alanine and shorter recurrence-free survival of BC patients [106]. Specifically, the authors demonstrated that lower expression levels of 4-aminobutyrate aminotransferase (ABAT), which negatively correlated with the concentration of β-alanine, indicated worse prognoses in BC patients. A similar study reported decreased ABAT expression in more aggressive BC subtypes, which correlated with an increased risk of metastasis and shorter overall, relapse-free, and distant metastasis-free survival [138]. Finally, Jansen et al. correlated ABAT downregulation with poor progression-free and metastasis-free survival in tamoxifen-treated patients [139].

4. Multi-Omics Studies and Novel BC Treatment Strategies

Omics-based technologies have also been used to identify novel therapeutic targets and monitor biological alterations related to BC metabolism following treatment [32,45,140,141,142,143,144,145,146]. The majority of studies relied on the application of metabolomics-based approaches in BC tissue [141,143] and serum [140,142,144] samples, although groups have evaluated transcriptomic [33,45] and proteomic [146] profiles in tissue samples. Various studies have described how the combination of omics approaches could characterize specific targets and foster the development of novel therapeutic strategies for specific subgroups of BC patients [147,148,149]. In particular, multi-omics studies have focused on identifying and validating metabolic enzymes as promising therapeutic strategies for the treatment of different BC tumors (Table 4). Figure 3 illustrates those metabolic-related genes proposed as potential therapeutic targets for treating BC patients in these studies. Overall, these findings suggest the therapeutic potential of inhibiting specific metabolic enzymes associated with glycolysis or involved in nucleotide, amino acid, and lipid metabolism in BC patients.
Iqbal and coworkers demonstrated that silencing CBX2 and CBX7 exhibited inverse effects on glycolysis, ATP production, viability, and proliferation [59]. CBX7 overexpression provided comparable results to CBX2 knockdown, which included decreased biomass production and reduced cell viability and proliferation. These in vitro results agreed with the findings of the transcripto-metabolomic analyses performed on BC patients and validated the roles of CBX2 and CBX7 in metabolic reprogramming of BC, highlighting the potential of these targets for the development of therapeutic strategies in BC. Of note, additional studies have provided similar results for CBX2 and CBX7 in BC [130,150,151], pancreatic adenocarcinoma [128] and metastatic prostate cancer [152]. Gong and coworkers evaluated the sensitivity of different BC metabolic phenotypes to metabolic inhibitors targeting glycolysis or de novo fatty acid synthesis [69]. The glycolytic BC phenotype displayed greater sensitivity to glycolytic inhibitors (oxamate, lactate dehydrogenase (LDH) inhibitors, and 2-deoxy-D-glucose), while inhibitors of lipid synthesis (cerulin and FASN inhibitor) exhibited higher efficacy against the lipogenic phenotype. Significantly, in vivo LDH inhibition enhanced tumor response to anti-PD-1 immunotherapy in the BC glycolytic phenotype. Previous studies have shown that LDH inhibition can suppress glycolysis [153,154] and cell proliferation [155] in BC cell lines. Differences in sensitivity to pharmacological inhibitors targeting glycolysis or electron transport chain (ECT) subunits were evaluated in different metabolic MECs phenotypes by Mahendralingam and colleagues [66]. In this study, the glycolytic phenotype displayed greater sensitivity to inhibitors targeting glucose transporter 1 (GLUT1), hexokinase (HK), LDH, and pyruvate dehydrogenase kinase (PDC), which agrees with results reported by Gong et al. [69]. Furthermore, studies have demonstrated that HK inhibition prevents BC growth [156,157].
The metabolic enzyme RRM2 has also been proposed as a potential therapeutic target for the treatment of BC. Putluri et al. observed that inhibiting RRM2 in BC cells significantly decreased proliferation and the expression of cell cycle genes and sensitized cells to tamoxifen treatment [65]. In agreement with the potential relevance of RRM2, additional studies have reported a reduction in proliferation [85,126] and tamoxifen resistance [158] in BC cell lines following RRM2 inhibition. Furthermore, associations between RRM2 overexpression and deterioration in BC survival have been widely reported, strongly suggesting a role as a targeted therapy for BC [61,120,121].
A multi-omics-based study by Terunuma et al. identified alterations associated with glutamine metabolism in a subset of BC tumors [107]. The authors described a subset of BC tumors with high 2HG levels and a distinct DNA methylation pattern associated with worse prognoses. Overall, studies have underscored the critical role of epigenetic-metabolomic interplay in promoting tumorigenesis [159]; in particular, high 2-HG levels induce epigenetic reprogramming associated with progression in different tumors [160,161]. Interestingly, the silencing of IDH2 and alcohol dehydrogenase iron containing 1 (ADHFE1), two enzymes implicated in the mitochondria-associated α-ketoglutarate–dependent production of 2HG [137,162], prompted a marked reduction of endogenous 2HG in BC cells. Furthermore, ADHFE1 loss resulted in a moderate but significant inhibition of cell cycle kinetics and reduced migration and invasion, suggesting an oncogenic role for ADHFE1 in BC. In agreement, several studies have associated high ADHFE1 expression levels with increased synthesis of 2HG and worse patient prognosis in BC [163,164].
In the context of a multi-omics study focused on lipid metabolism, Hilvo et al. conducted gene silencing experiments on seven enzymes involved in phospholipid remodeling and de novo lipid synthesis [60]. The results established that the individual inhibition of multiple lipid metabolism-regulating genes reduced the growth and viability of BC cell lines, which agrees with studies reporting reduced cell migration, invasion, and tumor proliferation in BC [165,166,167] and other tumor types [168,169,170] following the inhibition of specific lipid metabolism-related enzymes. Interestingly, a more recent study described FASN, another enzyme involved in de novo lipid metabolism, as a promising therapeutic target for BC treatment [78].
The results of a multi-omics-based study conducted by Kang et al. revealed lipid composition alterations during the EMT in BC [116]. The inhibition of ELOVL2 increased malignant potential, higher migration rate, and elevated colony formation. Mechanistically, downregulation of ELOVL2 increased sterol regulatory element-binding transcription factor 1 (SREBP1) expression in BC cells and activated lipogenesis, a process associated with the promotion of malignant BC phenotypes. SREBP1 is a crucial regulator of fatty acid metabolism and plays a pivotal role in the transcriptional regulation of different lipogenic genes mediating lipid synthesis [171,172]. SREBP1 overexpression has been observed in different tumor types, including BC [60,173,174], and supports the malignant BC phenotype [167]. Finally, Camarda et al. suggested the inhibition of FAO as a novel therapeutic target in a MYC-driven model of TNBC [102]. Based on the characterization of the effects of small-molecule inhibition and knock-down of carnitine palmitoyltransferase 1A (CPT1) and carnitine palmitoyltransferase 2 (CPT2) in BC cell lines, the authors demonstrated that FAO plays an essential role in this BC model. Furthermore, in vivo experiments demonstrated that treatment with etomoxir, a CPT1 inhibitor, significantly attenuated tumor growth in various BC models. Additional studies have also demonstrated that individual knockdown of both CPT enzymes reduces FAO metabolism and cell proliferation in different tumors [175,176,177]. In addition, more recent studies suggested CPT1 as a potential BC tumor target [103,178].

5. Future Perspectives and Conclusions

The application of single-omic approaches has significantly contributed to the characterization of BC tumors’ metabolic profile and the identification of BC-specific metabolic biomarkers. Individual omics have utility in clinical practice and can aid early diagnosis and identify indicators of disease progression. Among these experimental approaches, metabolomics represents a non-invasive strategy that could facilitate the development of non-invasive biomarkers with enormous potential for high-risk population screening, patient stratification, and treatment follow-up. However, the metabolome coverage remains limited compared to the genomic or transcriptomic layer, which limits the interpretation of the final results [179]. Moreover, compared with other omics platforms, metabolomics datasets suffer from a lower level of standardization, and we lack extensive data repositories [180,181,182]. Combining different omics datasets and taking advantage of the associated enhancement of statistical power may provide a more powerful strategy to characterize robust and consistent metabolic changes underlying disease states and identify novel therapeutic strategies [183]. Overall, studies describing metabolic alterations in glycolysis, fatty acid, nucleotide, lipid, and amino acid metabolism, could be employed to develop BC biomarkers and identify targets of great interest in clinical practice. This review focused on those multi-omic studies leading to the identification of metabolic markers with clinical utility for the early diagnosis, prognosis, and subtyping of BC tumors and the development of novel therapeutic strategies for BC based on metabolic-associated alterations. Multi-omics approaches have also been used to analyze the potential relationship between particular clinical variables (e.g., obese/non-obese [51]) and higher risk of BC. These studies have paved the way for the application of multi-omic approaches to the exploration of how patient-specific data affect the risk of BC (i.e., diabetes, metabolic syndrome, insulin resistance). Nevertheless, while certain studies included a vast number of samples [59,60,68,69,106], only a few include external independent cohorts of patients to assess the clinical significance of the results. Thus, future studies should analyze larger integrated sample cohorts that could generate sufficiently well-powered datasets and produce accurate and robust findings to validate the clinical utility of these findings. In these studies, coordinated sample processing poses technical challenges, given that each omics platform has specific requirements for sample treatment, and accessibility of available patient material is usually limited [184]. The lack of gold standard unified sample processing workflows and data analysis protocols for normalization, transformation, and scaling to ensure robustness and reproducibility of the results represent challenges to the integration of different datasets [185]. Moreover, in most reviewed studies, analyses of the different omics datasets were conducted individually before the combination of results, reflecting the lack of implementation of specific integration tools directed to the analysis of metabolic alterations. Therefore, developing new computational approaches that could facilitate the integrated analysis of different multi-omics datasets focusing on cancer metabolism would be incredibly beneficial [186]. In this scenario, the application of machine learning approaches to analyze multi-omics datasets could provide novel insights into the characterization of specific BC metabolic phenotypes [187].

Author Contributions

Conceptualization, N.G.-C., L.P.-C. and A.P.-L.; writing—original draft preparation, N.G.-C. and I.D.-O. with input from L.P.-C.; writing—review and editing, J.L.P., M.J.V., L.P.-C. and A.P.-L. All authors have read and agreed to the published version of the manuscript.

Funding

I.D.-O. is supported by an FPU grant (FPU19/03729) from the Spanish Ministry of Science and Innovation (PID2019-108806RB-I00). N.G.-C. and L.P.-C. are supported by the European Regional Development Fund (FEDER). N.G.-C., L.P.-C., and A.P.-L. are supported by the Spanish Ministry of Economy and Competitiveness (PID2020-115875RB-I00).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Stuart P. Atkinson for English editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Berry, D.A.; Cronin, K.A.; Plevritis, S.K.; Fryback, D.G.; Clarke, L.; Zelen, M.; Mandelblatt, J.S.; Yakovlev, A.Y.; Habbema, J.D.F.; Feuer, E.J. Effect of screening and adjuvant therapy on mortality from breast cancer. N. Engl. J. Med. 2005, 353, 1784–1792. [Google Scholar] [CrossRef]
  3. Plevritis, S.K.; Munoz, D.; Kurian, A.W.; Stout, N.K.; Alagoz, O.; Near, A.M.; Lee, S.J.; Van den Broek, J.J.; Huang, X.; Schechter, C.B.; et al. Association of screening and treatment with breast cancer mortality by molecular subtype in US women, 2000–2012. JAMA 2018, 319, 154–164. [Google Scholar] [CrossRef]
  4. Survival Rates for Breast Cancer. Available online: https://www.cancer.org/cancer/breast-cancer/understanding-a-breast-cancer-diagnosis/breast-cancer-survival-rates.html (accessed on 23 June 2021).
  5. Li, J.; Guan, X.; Fan, Z.; Ching, L.-M.; Li, Y.; Wang, X.; Cao, W.-M.; Liu, D.-X. Non-invasive biomarkers for early detection of breast cancer. Cancers 2020, 12, 2767. [Google Scholar] [CrossRef]
  6. Independent UK panel on breast cancer screening the benefits and harms of breast cancer screening: An independent review. Lancet 2012, 380, 1778–1786. [CrossRef]
  7. Lauby-Secretan, B.; Scoccianti, C.; Loomis, D.; Benbrahim-Tallaa, L.; Bouvard, V.; Bianchini, F.; Straif, K. International agency for research on cancer handbook working group breast-cancer screening—Viewpoint of the IARC working group. N. Engl. J. Med. 2015, 372, 2353–2358. [Google Scholar] [CrossRef] [Green Version]
  8. Houssami, N. Overdiagnosis of breast cancer in population screening: Does it make breast screening worthless? Cancer Biol. Med. 2017, 14, 1–8. [Google Scholar] [CrossRef] [Green Version]
  9. Nelson, H.D.; Pappas, M.; Cantor, A.; Griffin, J.; Daeges, M.; Humphrey, L. Harms of breast cancer screening: Systematic review to update the 2009 U.S. preventive services task force recommendation. Ann. Intern. Med. 2016, 164, 256. [Google Scholar] [CrossRef] [Green Version]
  10. Carter, J.L.; Coletti, R.J.; Harris, R.P. Quantifying and monitoring overdiagnosis in cancer screening: A systematic review of methods. BMJ 2015, 350, g7773. [Google Scholar] [CrossRef] [Green Version]
  11. Brewer, H.R.; Jones, M.E.; Schoemaker, M.J.; Ashworth, A.; Swerdlow, A.J. Family history and risk of breast cancer: An analysis accounting for family structure. Breast Cancer Res. Treat. 2017, 165, 193–200. [Google Scholar] [CrossRef] [Green Version]
  12. Cardoso, F.; Kyriakides, S.; Ohno, S.; Penault-Llorca, F.; Poortmans, P.; Rubio, I.T.; Zackrisson, S.; Senkus, E. ESMO Guidelines Committee. Early breast cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 2019, 30, 1194–1220. [Google Scholar] [CrossRef] [Green Version]
  13. Cardoso, F.; Senkus, E.; Costa, A.; Papadopoulos, E.; Aapro, M.; André, F.; Harbeck, N.; Aguilar Lopez, B.; Barrios, C.H.; Bergh, J.; et al. 4th ESO-ESMO international consensus guidelines for advanced breast cancer (ABC 4). Ann. Oncol. 2018, 29, 1634–1657. [Google Scholar] [CrossRef] [PubMed]
  14. Network, T.C.G.A. Cancer genome atlas network comprehensive molecular portraits of human breast tumours. Nature 2012, 490, 61–70. [Google Scholar] [CrossRef] [Green Version]
  15. Milioli, H.H.; Tishchenko, I.; Riveros, C.; Berretta, R.; Moscato, P. Basal-like breast cancer: Molecular profiles, clinical Features and survival outcomes. BMC Med. Genom. 2017, 10, 19. [Google Scholar] [CrossRef] [Green Version]
  16. Anderson, W.F.; Chen, B.E.; Jatoi, I.; Rosenberg, P.S. Effects of estrogen receptor expression and histopathology on annual hazard rates of death from breast cancer. Breast Cancer Res. Treat. 2006, 100, 121–126. [Google Scholar] [CrossRef] [PubMed]
  17. Goncalves, R.; Ma, C.; Luo, J.; Suman, V.; Ellis, M.J. Use of neoadjuvant data to design adjuvant endocrine therapy trials for breast cancer. Nat. Rev. Clin. Oncol. 2012, 9, 223–229. [Google Scholar] [CrossRef] [PubMed]
  18. Rouzier, R.; Perou, C.M.; Symmans, W.F.; Ibrahim, N.; Cristofanilli, M.; Anderson, K.; Hess, K.R.; Stec, J.; Ayers, M.; Wagner, P.; et al. Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Clin. Cancer Res. 2005, 11, 5678–5685. [Google Scholar] [CrossRef] [Green Version]
  19. Sotiriou, C.; Neo, S.-Y.; McShane, L.M.; Korn, E.L.; Long, P.M.; Jazaeri, A.; Martiat, P.; Fox, S.B.; Harris, A.L.; Liu, E.T. Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc. Natl. Acad. Sci. USA 2003, 100, 10393–10398. [Google Scholar] [CrossRef] [Green Version]
  20. Van’t Veer, L.J.; Dai, H.; Van de Vijver, M.J.; He, Y.D.; Hart, A.A.M.; Mao, M.; Peterse, H.L.; Van der Kooy, K.; Marton, M.J.; Witteveen, A.T.; et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature 2002, 415, 530–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Carey, L.A.; Dees, E.C.; Sawyer, L.; Gatti, L.; Moore, D.T.; Collichio, F.; Ollila, D.W.; Sartor, C.I.; Graham, M.L.; Perou, C.M. The triple negative paradox: Primary tumor chemosensitivity of breast cancer subtypes. Clin. Cancer Res. 2007, 13, 2329–2334. [Google Scholar] [CrossRef] [Green Version]
  22. Huang, E.; Cheng, S.H.; Dressman, H.; Pittman, J.; Tsou, M.H.; Horng, C.F.; Bild, A.; Iversen, E.S.; Liao, M.; Chen, C.M.; et al. Gene expression predictors of breast cancer outcomes. Lancet 2003, 361, 1590–1596. [Google Scholar] [CrossRef]
  23. Phung, M.T.; Tin Tin, S.; Elwood, J.M. Prognostic models for breast cancer: A systematic review. BMC Cancer 2019, 19, 230. [Google Scholar] [CrossRef] [Green Version]
  24. Lehmann, B.D.; Bauer, J.A.; Chen, X.; Sanders, M.E.; Chakravarthy, A.B.; Shyr, Y.; Pietenpol, J.A. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J. Clin. Investig. 2011, 121, 2750–2767. [Google Scholar] [CrossRef] [Green Version]
  25. Burstein, M.D.; Tsimelzon, A.; Poage, G.M.; Covington, K.R.; Contreras, A.; Fuqua, S.A.W.; Savage, M.I.; Osborne, C.K.; Hilsenbeck, S.G.; Chang, J.C.; et al. Comprehensive genomic analysis identifies novel subtypes and targets of triple-negative breast cancer. Clin. Cancer Res. 2015, 21, 1688–1698. [Google Scholar] [CrossRef] [Green Version]
  26. Bareche, Y.; Venet, D.; Ignatiadis, M.; Aftimos, P.; Piccart, M.; Rothe, F.; Sotiriou, C. Unravelling triple-negative breast cancer molecular heterogeneity using an integrative multiomic analysis. Ann. Oncol. 2018, 29, 895–902. [Google Scholar] [CrossRef] [PubMed]
  27. Jiang, Y.-Z.; Ma, D.; Suo, C.; Shi, J.; Xue, M.; Hu, X.; Xiao, Y.; Yu, K.-D.; Liu, Y.-R.; Yu, Y.; et al. Genomic and transcriptomic landscape of triple-negative breast cancers: Subtypes and treatment strategies. Cancer Cell 2019, 35, 428–440.e5. [Google Scholar] [CrossRef] [Green Version]
  28. Garrido-Castro, A.C.; Lin, N.U.; Polyak, K. Insights into molecular classifications of triple-negative breast cancer: Improving patient selection for treatment. Cancer Discov. 2019, 9, 176–198. [Google Scholar] [CrossRef] [Green Version]
  29. Dignam, J.J.; Dukic, V.; Anderson, S.J.; Mamounas, E.P.; Wickerham, D.L.; Wolmark, N. Hazard of Recurrence and Adjuvant Treatment Effects over Time in Lymph Node-Negative Breast Cancer. Breast Cancer Res. Treat. 2009, 116, 595–602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Denkert, C.; Liedtke, C.; Tutt, A.; Von Minckwitz, G. Molecular alterations in triple-negative breast cancer-the road to new treatment strategies. Lancet 2017, 389, 2430–2442. [Google Scholar] [CrossRef] [Green Version]
  31. Harbeck, N.; Penault-Llorca, F.; Cortes, J.; Gnant, M.; Houssami, N.; Poortmans, P.; Ruddy, K.; Tsang, J.; Cardoso, F. Breast cancer. Nat. Rev. Dis. Primers 2019, 5, 66. [Google Scholar] [CrossRef]
  32. Armiñán, A.; Palomino-Schätzlein, M.; Deladriere, C.; Arroyo-Crespo, J.J.; Vicente-Ruiz, S.; Vicent, M.J.; Pineda-Lucena, A. Metabolomics facilitates the discrimination of the specific anti-cancer effects of free- and polymer-conjugated doxorubicin in breast cancer models. Biomaterials 2018, 162, 144–153. [Google Scholar] [CrossRef]
  33. Arroyo-Crespo, J.J.; Armiñán, A.; Charbonnier, D.; Deladriere, C.; Palomino-Schätzlein, M.; Lamas-Domingo, R.; Forteza, J.; Pineda-Lucena, A.; Vicent, M.J. Characterization of triple-negative breast cancer preclinical models provides functional evidence of metastatic progression. Int. J. Cancer 2019, 145, 2267–2281. [Google Scholar] [CrossRef] [Green Version]
  34. Cao, M.D.; Lamichhane, S.; Lundgren, S.; Bofin, A.; Fjøsne, H.; Giskeødegård, G.F.; Bathen, T.F. Metabolic characterization of triple negative breast cancer. BMC Cancer 2014, 14, 941. [Google Scholar] [CrossRef] [Green Version]
  35. Tayyari, F.; Gowda, G.A.N.; Olopade, O.F.; Berg, R.; Yang, H.H.; Lee, M.P.; Ngwa, W.F.; Mittal, S.K.; Raftery, D.; Mohammed, S.I. Metabolic profiles of triple-negative and luminal A breast cancer subtypes in African-American identify key metabolic differences. Oncotarget 2018, 9, 11677–11690. [Google Scholar] [CrossRef] [Green Version]
  36. Slupsky, C.M.; Steed, H.; Wells, T.H.; Dabbs, K.; Schepansky, A.; Capstick, V.; Faught, W.; Sawyer, M.B. Urine metabolite analysis offers potential early diagnosis of ovarian and breast cancers. Clin. Cancer Res. 2010, 16, 5835–5841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cala, M.; Aldana, J.; Sánchez, J.; Guio, J.; Meesters, R.J.W. Urinary metabolite and lipid alterations in colombian hispanic women with breast cancer: A pilot study. J. Pharm. Biomed. Anal. 2018, 152, 234–241. [Google Scholar] [CrossRef] [PubMed]
  38. Poschke, I.; Mao, Y.; Kiessling, R.; de Boniface, J. Tumor-dependent increase of serum amino acid levels in breast cancer patients has diagnostic potential and correlates with molecular tumor subtypes. J. Transl. Med. 2013, 11, 290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Fan, Y.; Zhou, X.; Xia, T.-S.; Chen, Z.; Li, J.; Liu, Q.; Alolga, R.N.; Chen, Y.; Lai, M.-D.; Li, P.; et al. Human plasma metabolomics for identifying differential metabolites and predicting molecular subtypes of breast cancer. Oncotarget 2016, 7, 9925–9938. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Sugimoto, M.; Wong, D.T.; Hirayama, A.; Soga, T.; Tomita, M. Capillary electrophoresis mass spectrometry-based saliva metabolomics identified oral, breast and pancreatic cancer-specific profiles. Metabolomics 2010, 6, 78–95. [Google Scholar] [CrossRef] [Green Version]
  41. Willmann, L.; Schlimpert, M.; Halbach, S.; Erbes, T.; Stickeler, E.; Kammerer, B. Metabolic profiling of breast cancer: Differences in central metabolism between subtypes of breast cancer cell lines. J. Chromatogr. B 2015, 1000, 95–104. [Google Scholar] [CrossRef]
  42. Gaude, E.; Frezza, C. Tissue-specific and convergent metabolic transformation of cancer correlates with metastatic potential and patient survival. Nat. Commun. 2016, 7, 13041. [Google Scholar] [CrossRef] [Green Version]
  43. Rosario, S.R.; Long, M.D.; Affronti, H.C.; Rowsam, A.M.; Eng, K.H.; Smiraglia, D.J. Pan-cancer analysis of transcriptional metabolic dysregulation using the cancer genome atlas. Nat. Commun. 2018, 9, 5330. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, L.; Zhang, S.; Wang, X. The Metabolic mechanisms of breast cancer metastasis. Front. Oncol. 2021, 10, 602416. [Google Scholar] [CrossRef]
  45. Lord, S.R.; Collins, J.M.; Cheng, W.-C.; Haider, S.; Wigfield, S.; Gaude, E.; Fielding, B.A.; Pinnick, K.E.; Harjes, U.; Segaran, A.; et al. Transcriptomic analysis of human primary breast cancer identifies fatty acid oxidation as a target for metformin. Br. J. Cancer 2020, 122, 258–265. [Google Scholar] [CrossRef] [Green Version]
  46. Bernhardt, S.; Bayerlová, M.; Vetter, M.; Wachter, A.; Mitra, D.; Hanf, V.; Lantzsch, T.; Uleer, C.; Peschel, S.; John, J.; et al. Proteomic profiling of breast cancer metabolism identifies SHMT2 and ASCT2 as prognostic factors. Breast Cancer Res. 2017, 19, 112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Borgan, E.; Sitter, B.; Lingjærde, O.C.; Johnsen, H.; Lundgren, S.; Bathen, T.F.; Sørlie, T.; Børresen-Dale, A.-L.; Gribbestad, I.S. Merging transcriptomics and metabolomics—Advances in breast cancer profiling. BMC Cancer 2010, 10, 628. [Google Scholar] [CrossRef] [Green Version]
  48. The Oslo Breast Cancer Consortium (OSBREAC); Haukaas, T.H.; Euceda, L.R.; Giskeødegård, G.F.; Lamichhane, S.; Krohn, M.; Jernström, S.; Aure, M.R.; Lingjærde, O.C.; Schlichting, E.; et al. Metabolic clusters of breast cancer in relation to gene- and protein expression subtypes. Cancer Metab. 2016, 4, 12. [Google Scholar] [CrossRef] [Green Version]
  49. OSBREAC; Aure, M.R.; Vitelli, V.; Jernström, S.; Kumar, S.; Krohn, M.; Due, E.U.; Haukaas, T.H.; Leivonen, S.-K.; Vollan, H.K.M.; et al. Integrative clustering reveals a novel split in the luminal A subtype of breast cancer with impact on outcome. Breast Cancer Res. 2017, 19, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Nam, H.; Chung, B.C.; Kim, Y.; Lee, K.; Lee, D. Combining tissue transcriptomics and urine metabolomics for breast cancer biomarker identification. Bioinformatics 2009, 25, 3151–3157. [Google Scholar] [CrossRef] [Green Version]
  51. Hassan, M.A.; Al-Sakkaf, K.; Shait Mohammed, M.R.; Dallol, A.; Al-Maghrabi, J.; Aldahlawi, A.; Ashoor, S.; Maamra, M.; Ragoussis, J.; Wu, W.; et al. Integration of transcriptome and metabolome provides unique insights to pathways associated with obese breast cancer patients. Front. Oncol. 2020, 10, 804. [Google Scholar] [CrossRef]
  52. Mastoraki, S.; Strati, A.; Tzanikou, E.; Chimonidou, M.; Politaki, E.; Voutsina, A.; Psyrri, A.; Georgoulias, V.; Lianidou, E. ESR1 Methylation: A liquid biopsy-based epigenetic assay for the follow-up of patients with metastatic breast cancer receiving endocrine treatment. Clin. Cancer Res. 2018, 24, 1500–1510. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Schröder, R.; Illert, A.-L.; Erbes, T.; Flotho, C.; Lübbert, M.; Duque-Afonso, J. The epigenetics of breast cancer—Opportunities for diagnostics, risk stratification and therapy. Epigenetics 2021, 1–13. [Google Scholar] [CrossRef] [PubMed]
  54. Serrano-Carbajal, E.A.; Espinal-Enríquez, J.; Hernández-Lemus, E. Targeting metabolic deregulation landscapes in breast cancer subtypes. Front. Oncol. 2020, 10, 97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Bernhardt, S.; Tönsing, C.; Mitra, D.; Erdem, N.; Müller-Decker, K.; Korf, U.; Kreutz, C.; Timmer, J.; Wiemann, S. Functional proteomics of breast cancer metabolism identifies GLUL as responder during hypoxic adaptation. J. Proteome Res. 2019, 18, 1352–1362. [Google Scholar] [CrossRef] [PubMed]
  56. Asiago, V.M.; Alvarado, L.Z.; Shanaiah, N.; Gowda, G.A.N.; Owusu-Sarfo, K.; Ballas, R.A.; Raftery, D. Early detection of recurrent breast cancer using metabolite profiling. Cancer Res. 2010, 70, 8309–8318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Huang, S.; Chong, N.; Lewis, N.E.; Jia, W.; Xie, G.; Garmire, L.X. Novel personalized pathway-based metabolomics models reveal key metabolic pathways for breast cancer diagnosis. Genome Med. 2016, 8, 34. [Google Scholar] [CrossRef] [Green Version]
  58. Wei, Y.; Jasbi, P.; Shi, X.; Turner, C.; Hrovat, J.; Liu, L.; Rabena, Y.; Porter, P.; Gu, H. Early breast cancer detection using untargeted and targeted metabolomics. J. Proteome Res. 2021, 20, 3124–3133. [Google Scholar] [CrossRef]
  59. Iqbal, M.A.; Siddiqui, S.; Ur Rehman, A.; Siddiqui, F.A.; Singh, P.; Kumar, B.; Saluja, D. Multiomics integrative analysis reveals antagonistic roles of CBX2 and CBX7 in metabolic reprogramming of breast cancer. Mol. Oncol. 2021, 15, 1450–1465. [Google Scholar] [CrossRef]
  60. Hilvo, M.; Denkert, C.; Lehtinen, L.; Müller, B.; Brockmöller, S.; Seppänen-Laakso, T.; Budczies, J.; Bucher, E.; Yetukuri, L.; Castillo, S.; et al. Novel theranostic opportunities offered by characterization of altered membrane lipid metabolism in breast cancer progression. Cancer Res. 2011, 71, 3236–3245. [Google Scholar] [CrossRef] [Green Version]
  61. Luo, X.; Yu, H.; Song, Y.; Sun, T. Integration of metabolomic and transcriptomic data reveals metabolic pathway alteration in breast cancer and impact of related signature on survival. J. Cell. Physiol. 2019, 234, 13021–13031. [Google Scholar] [CrossRef]
  62. Dowling, P.; Henry, M.; Meleady, P.; Clarke, C.; Gately, K.; O’Byrne, K.; Connolly, E.; Lynch, V.; Ballot, J.; Gullo, G.; et al. Metabolomic and proteomic analysis of breast cancer patient samples suggests that glutamate and 12-HETE in combination with CA15-3 may be useful biomarkers reflecting tumour burden. Metabolomics 2015, 11, 620–635. [Google Scholar] [CrossRef]
  63. Moestue, S.A.; Borgan, E.; Huuse, E.M.; Lindholm, E.M.; Sitter, B.; Børresen-Dale, A.-L.; Engebraaten, O.; Mælandsmo, G.M.; Gribbestad, I.S. Distinct choline metabolic profiles are associated with differences in gene expression for basal-like and luminal-like breast cancer xenograft models. BMC Cancer 2010, 10, 433. [Google Scholar] [CrossRef] [Green Version]
  64. Grinde, M.T.; Skrbo, N.; Moestue, S.A.; Rødland, E.A.; Borgan, E.; Kristian, A.; Sitter, B.; Bathen, T.F.; Børresen-Dale, A.-L.; Mælandsmo, G.M.; et al. Interplay of choline metabolites and genes in patient-derived breast cancer xenografts. Breast Cancer Res 2014, 16, R5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Putluri, N.; Maity, S.; Kommagani, R.; Creighton, C.J.; Putluri, V.; Chen, F.; Nanda, S.; Bhowmik, S.K.; Terunuma, A.; Dorsey, T.; et al. Pathway-centric integrative analysis identifies RRM2 as a prognostic marker in breast cancer associated with poor survival and tamoxifen resistance. Neoplasia 2014, 16, 390–402. [Google Scholar] [CrossRef] [Green Version]
  66. Mahendralingam, M.J.; Kim, H.; McCloskey, C.W.; Aliar, K.; Casey, A.E.; Tharmapalan, P.; Pellacani, D.; Ignatchenko, V.; Garcia-Valero, M.; Palomero, L.; et al. Mammary epithelial cells have lineage-rooted metabolic identities. Nat. Metab. 2021, 3, 665–681. [Google Scholar] [CrossRef]
  67. Tang, X.; Lin, C.-C.; Spasojevic, I.; Iversen, E.S.; Chi, J.-T.; Marks, J.R. A joint analysis of metabolomics and genetics of breast cancer. Breast Cancer Res. 2014, 16, 415. [Google Scholar] [CrossRef] [Green Version]
  68. Barupal, D.K.; Gao, B.; Budczies, J.; Phinney, B.S.; Perroud, B.; Denkert, C.; Fiehn, O. Prioritization of metabolic genes as novel therapeutic targets in estrogen-receptor negative breast tumors using multi-omics data and text mining. Oncotarget 2019, 10, 3894–3909. [Google Scholar] [CrossRef] [Green Version]
  69. Gong, Y.; Ji, P.; Yang, Y.-S.; Xie, S.; Yu, T.-J.; Xiao, Y.; Jin, M.-L.; Ma, D.; Guo, L.-W.; Pei, Y.-C.; et al. Metabolic-pathway-based subtyping of triple-negative breast cancer reveals potential therapeutic targets. Cell Metab. 2021, 33, 51–64.e9. [Google Scholar] [CrossRef]
  70. Yang, K.; Han, X. Lipidomics: Techniques, applications, and outcomes related to biomedical sciences. Trends Biochem. Sci. 2016, 41, 954–969. [Google Scholar] [CrossRef] [Green Version]
  71. Podo, F. Tumour phospholipid metabolism. NMR Biomed. 1999, 12, 413–439. [Google Scholar] [CrossRef]
  72. Menendez, J.A.; Lupu, R.; Colomer, R. Targeting fatty acid synthase: Potential for therapeutic intervention in her-2/neu-overexpressing breast cancer. Drug News Perspect. 2005, 18, 375–385. [Google Scholar] [CrossRef] [PubMed]
  73. Cairns, R.A.; Harris, I.S.; Mak, T.W. Regulation of cancer cell metabolism. Nat. Rev. Cancer 2011, 11, 85–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bhattacharjee, K.; Nath, M.; Choudhury, Y. Fatty acid synthesis and cancer: Aberrant expression of the ACACA and ACACB genes increases the risk for cancer. Meta Gene 2020, 26, 100798. [Google Scholar] [CrossRef]
  75. Milgraum, L.Z.; Witters, L.A.; Pasternack, G.R.; Kuhajda, F.P. Enzymes of the fatty acid synthesis pathway are highly expressed in in situ breast carcinoma. Clin. Cancer Res. 1997, 3, 2115–2120. [Google Scholar]
  76. Li, J.; Dong, L.; Wei, D.; Wang, X.; Zhang, S.; Li, H. Fatty acid synthase mediates the epithelial-mesenchymal transition of breast cancer cells. Int. J. Biol. Sci. 2014, 10, 171–180. [Google Scholar] [CrossRef] [Green Version]
  77. Xu, S.; Chen, T.; Dong, L.; Li, T.; Xue, H.; Gao, B.; Ding, X.; Wang, H.; Li, H. Fatty acid synthase promotes breast cancer metastasis by mediating changes in fatty acid metabolism. Oncol. Lett. 2020, 21, 27. [Google Scholar] [CrossRef]
  78. Menendez, J.A.; Lupu, R. Fatty acid synthase (FASN) as a therapeutic target in breast cancer. Expert Opin. Ther. Targets 2017, 21, 1001–1016. [Google Scholar] [CrossRef]
  79. Lv, Y.; Wang, X.; Li, X.; Xu, G.; Bai, Y.; Wu, J.; Piao, Y.; Shi, Y.; Xiang, R.; Wang, L. Nucleotide de novo synthesis increases breast cancer stemness and metastasis via CGMP-PKG-MAPK signaling pathway. PLoS Biol. 2020, 18, e3000872. [Google Scholar] [CrossRef]
  80. Stagg, J.; Divisekera, U.; McLaughlin, N.; Sharkey, J.; Pommey, S.; Denoyer, D.; Dwyer, K.M.; Smyth, M.J. Anti-CD73 antibody therapy inhibits breast tumor growth and metastasis. Proc. Natl. Acad. Sci. USA 2010, 107, 1547–1552. [Google Scholar] [CrossRef] [Green Version]
  81. Sadej, R.; Skladanowski, A.C. Dual, Enzymatic and non-enzymatic, function of ecto-5’-nucleotidase (EN, CD73) in migration and invasion of A375 melanoma cells. Acta Biochim. Pol. 2012, 59. [Google Scholar] [CrossRef] [Green Version]
  82. Camici, M.; Tozzi, M.G.; Allegrini, S.; Del Corso, A.; Sanfilippo, O.; Daidone, M.G.; De Marco, C.; Ipata, P.L. Purine salvage enzyme activities in normal and neoplastic human tissues. Cancer Biochem. Biophys. 1990, 11, 201–209. [Google Scholar]
  83. Aghaei, M.; Karami-Tehrani, F.; Salami, S.; Atri, M. Adenosine deaminase activity in the serum and malignant tumors of breast cancer: The assessment of isoenzyme ADA1 and ADA2 activities. Clinical Biochemistry 2005, 38, 887–891. [Google Scholar] [CrossRef] [PubMed]
  84. Mahajan, M.; Tiwari, N.; Sharma, R.; Kaur, S.; Singh, N. Oxidative stress and its relationship with adenosine deaminase activity in various stages of breast cancer. Ind. J. Clin. Biochem. 2013, 28, 51–54. [Google Scholar] [CrossRef] [Green Version]
  85. Mejía-Pedroza, R.A.; Espinal-Enríquez, J.; Hernández-Lemus, E. Pathway-based drug repositioning for breast cancer molecular subtypes. Front. Pharmacol. 2018, 9, 905. [Google Scholar] [CrossRef] [Green Version]
  86. Yang, L.; Wang, Y.; Cai, H.; Wang, S.; Shen, Y.; Ke, C. Application of metabolomics in the diagnosis of breast cancer: A systematic review. J. Cancer 2020, 11, 2540–2551. [Google Scholar] [CrossRef]
  87. Xie, G.; Zhou, B.; Zhao, A.; Qiu, Y.; Zhao, X.; Garmire, L.; Shvetsov, Y.B.; Yu, H.; Yen, Y.; Jia, W. Lowered circulating aspartate is a metabolic feature of human breast cancer. Oncotarget 2015, 6, 33369–33381. [Google Scholar] [CrossRef] [Green Version]
  88. Wang, X.; Zhao, X.; Chou, J.; Yu, J.; Yang, T.; Liu, L.; Zhang, F. Taurine, glutamic acid and ethylmalonic acid as important metabolites for detecting human breast cancer based on the targeted metabolomics. CBM 2018, 23, 255–268. [Google Scholar] [CrossRef]
  89. Redefining Cancer Subtypes Using Multi-Omics and Deep Learning—ARCAS–AI for Life Sciences. Available online: https://arcas.ai/blog/redefining-cancer-subtypes-using-multi-omics-and-deep-learning (accessed on 23 June 2021).
  90. Sathyanarayanan, A.; Gupta, R.; Thompson, E.W.; Nyholt, D.R.; Bauer, D.C.; Nagaraj, S.H. A comparative study of multi-omics integration tools for cancer driver gene identification and tumour subtyping. Brief. Bioinform. 2020, 21, 1920–1936. [Google Scholar] [CrossRef]
  91. Sitter, B.; Bathen, T.F.; Singstad, T.E.; Fjøsne, H.E.; Lundgren, S.; Halgunset, J.; Gribbestad, I.S. Quantification of metabolites in breast cancer patients with different clinical prognosis using HR MAS MR spectroscopy. NMR Biomed. 2010, 23, 424–431. [Google Scholar] [CrossRef]
  92. Tan, J.; Le, A. The heterogeneity of breast cancer metabolism. Adv. Exp. Med. Biol. 2021, 1311, 89–101. [Google Scholar] [CrossRef]
  93. Eghlimi, R.; Shi, X.; Hrovat, J.; Xi, B.; Gu, H. Triple negative breast cancer detection using LC-MS/MS lipidomic profiling. J. Proteome Res. 2020, 19, 2367–2378. [Google Scholar] [CrossRef]
  94. Lampa, M.; Arlt, H.; He, T.; Ospina, B.; Reeves, J.; Zhang, B.; Murtie, J.; Deng, G.; Barberis, C.; Hoffmann, D.; et al. Glutaminase is essential for the growth of triple-negative breast cancer cells with a deregulated glutamine metabolism pathway and its suppression synergizes with MTOR inhibition. PLoS ONE 2017, 12, e0185092. [Google Scholar] [CrossRef]
  95. Timmerman, L.A.; Holton, T.; Yuneva, M.; Louie, R.J.; Padró, M.; Daemen, A.; Hu, M.; Chan, D.A.; Ethier, S.P.; Van’t Veer, L.J.; et al. Glutamine sensitivity analysis identifies the xCT antiporter as a common triple-negative breast tumor therapeutic target. Cancer Cell 2013, 24, 450–465. [Google Scholar] [CrossRef] [Green Version]
  96. Brown, K.K.; Spinelli, J.B.; Asara, J.M.; Toker, A. Adaptive reprogramming of de novo pyrimidine synthesis is a metabolic vulnerability in triple-negative breast cancer. Cancer Discov. 2017, 7, 782. [Google Scholar] [CrossRef] [Green Version]
  97. Suzuki, Y.; Suda, T.; Furuhashi, K.; Suzuki, M.; Fujie, M.; Hahimoto, D.; Nakamura, Y.; Inui, N.; Nakamura, H.; Chida, K. Increased Serum Kynurenine/Tryptophan Ratio Correlates with Disease Progression in Lung Cancer. Lung Cancer 2010, 67, 361–365. [Google Scholar] [CrossRef]
  98. Lyon, D.E.; Walter, J.M.; Starkweather, A.R.; Schubert, C.M.; McCain, N.L. Tryptophan degradation in women with breast cancer: A pilot study. BMC Res. Notes 2011, 4, 156. [Google Scholar] [CrossRef] [Green Version]
  99. Zaugg, K.; Yao, Y.; Reilly, P.T.; Kannan, K.; Kiarash, R.; Mason, J.; Huang, P.; Sawyer, S.K.; Fuerth, B.; Faubert, B.; et al. Carnitine palmitoyltransferase 1c promotes cell survival and tumor growth under conditions of metabolic stress. Genes Dev. 2011, 25, 1041–1051. [Google Scholar] [CrossRef] [Green Version]
  100. Gatza, M.L.; Silva, G.O.; Parker, J.S.; Fan, C.; Perou, C.M. An Integrated genomics approach identifies drivers of proliferation in luminal-subtype human breast cancer. Nat. Genet. 2014, 46, 1051–1059. [Google Scholar] [CrossRef] [Green Version]
  101. Bjørklund, S.S.; Kristensen, V.N.; Seiler, M.; Kumar, S.; Alnæs, G.I.G.; Ming, Y.; Kerrigan, J.; Naume, B.; Sachidanandam, R.; Bhanot, G.; et al. Expression of an estrogen-regulated variant transcript of the peroxisomal branched chain fatty acid oxidase ACOX2 in breast carcinomas. BMC Cancer 2015, 15, 524. [Google Scholar] [CrossRef] [Green Version]
  102. Camarda, R.; Zhou, A.Y.; Kohnz, R.A.; Balakrishnan, S.; Mahieu, C.; Anderton, B.; Eyob, H.; Kajimura, S.; Tward, A.; Krings, G.; et al. Inhibition of fatty acid oxidation as a therapy for MYC-overexpressing triple-negative breast cancer. Nat. Med. 2016, 22, 427–432. [Google Scholar] [CrossRef]
  103. Pucci, S.; Zonetti, M.J.; Fisco, T.; Polidoro, C.; Bocchinfuso, G.; Palleschi, A.; Novelli, G.; Spagnoli, L.G.; Mazzarelli, P. Carnitine palmitoyl transferase-1A (CPT1A): A new tumor specific target in human breast cancer. Oncotarget 2016, 7, 19982–19996. [Google Scholar] [CrossRef] [Green Version]
  104. Patra, K.C.; Hay, N. The pentose phosphate pathway and cancer. Trends Biochem. Sci. 2014, 39, 347–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Ge, T.; Yang, J.; Zhou, S.; Wang, Y.; Li, Y.; Tong, X. The role of the pentose phosphate pathway in diabetes and cancer. Front. Endocrinol. 2020, 11, 365. [Google Scholar] [CrossRef]
  106. Budczies, J.; Brockmöller, S.F.; Müller, B.M.; Barupal, D.K.; Richter-Ehrenstein, C.; Kleine-Tebbe, A.; Griffin, J.L.; Orešič, M.; Dietel, M.; Denkert, C.; et al. Comparative metabolomics of estrogen receptor positive and estrogen receptor negative breast cancer: Alterations in glutamine and beta-alanine metabolism. J. Proteom. 2013, 94, 279–288. [Google Scholar] [CrossRef]
  107. Terunuma, A.; Putluri, N.; Mishra, P.; Mathé, E.A.; Dorsey, T.H.; Yi, M.; Wallace, T.A.; Issaq, H.J.; Zhou, M.; Killian, J.K.; et al. MYC-driven accumulation of 2-hydroxyglutarate is associated with breast cancer prognosis. J. Clin. Investig. 2014, 124, 398–412. [Google Scholar] [CrossRef] [Green Version]
  108. Fares, J.; Fares, M.Y.; Khachfe, H.H.; Salhab, H.A.; Fares, Y. Molecular principles of metastasis: A hallmark of cancer revisited. Signal Transduct. Target. Ther. 2020, 5, 28. [Google Scholar] [CrossRef] [PubMed]
  109. Wei, Q.; Qian, Y.; Yu, J.; Wong, C.C. Metabolic rewiring in the promotion of cancer metastasis: Mechanisms and therapeutic implications. Oncogene 2020, 39, 6139–6156. [Google Scholar] [CrossRef]
  110. Kreuzaler, P.; Panina, Y.; Segal, J.; Yuneva, M. Adapt and conquer: Metabolic flexibility in cancer growth, invasion and evasion. Mol. Metab. 2020, 33, 83–101. [Google Scholar] [CrossRef]
  111. Läsche, M.; Emons, G.; Gründker, C. Shedding new light on cancer metabolism: A metabolic tightrope between life and death. Front. Oncol. 2020, 10, 409. [Google Scholar] [CrossRef]
  112. Li, C.; Li, X.; Li, G.; Sun, L.; Zhang, W.; Jiang, J.; Ge, Q. Identification of a prognosis-associated signature associated with energy metabolism in triple-negative breast cancer. Oncol. Rep. 2020, 44, 819–837. [Google Scholar] [CrossRef]
  113. Feng, J.; Ren, J.; Yang, Q.; Liao, L.; Cui, L.; Gong, Y.; Sun, S. Metabolic gene signature for predicting breast cancer recurrence using transcriptome analysis. Future Oncol. 2021, 17, 71–80. [Google Scholar] [CrossRef]
  114. Oakman, C.; Tenori, L.; Claudino, W.M.; Cappadona, S.; Nepi, S.; Battaglia, A.; Bernini, P.; Zafarana, E.; Saccenti, E.; Fornier, M.; et al. Identification of a serum-detectable metabolomic fingerprint potentially correlated with the presence of micrometastatic disease in early breast cancer patients at varying risks of disease relapse by traditional prognostic methods. Ann. Oncol. 2011, 22, 1295–1301. [Google Scholar] [CrossRef]
  115. Jobard, E.; Pontoizeau, C.; Blaise, B.J.; Bachelot, T.; Elena-Herrmann, B.; Trédan, O. A serum nuclear magnetic resonance-based metabolomic signature of advanced metastatic human breast cancer. Cancer Lett. 2014, 343, 33–41. [Google Scholar] [CrossRef]
  116. Kang, Y.P.; Yoon, J.-H.; Long, N.P.; Koo, G.-B.; Noh, H.-J.; Oh, S.-J.; Lee, S.B.; Kim, H.M.; Hong, J.Y.; Lee, W.J.; et al. Spheroid-induced epithelial-mesenchymal transition provokes global alterations of breast cancer lipidome: A multi-layered omics analysis. Front. Oncol. 2019, 9, 145. [Google Scholar] [CrossRef]
  117. Sigoillot, F.D.; Sigoillot, S.M.; Guy, H.I. Breakdown of the regulatory control of pyrimidine biosynthesis in human breast cancer cells. Int. J. Cancer 2004, 109, 491–498. [Google Scholar] [CrossRef]
  118. Davison, C.; Morelli, R.; Knowlson, C.; McKechnie, M.; Carson, R.; Stachtea, X.; McLaughlin, K.A.; Prise, V.E.; Savage, K.; Wilson, R.H.; et al. Targeting nucleotide metabolism enhances the efficacy of anthracyclines and anti-metabolites in triple-negative breast cancer. NPJ Breast Cancer 2021, 7, 38. [Google Scholar] [CrossRef] [PubMed]
  119. Furuta, E.; Okuda, H.; Kobayashi, A.; Watabe, K. Metabolic genes in cancer: Their roles in tumor progression and clinical implications. Biochim. Biophys. Acta 2010, 1805, 141–152. [Google Scholar] [CrossRef] [Green Version]
  120. Zhang, H.; Liu, X.; Warden, C.D.; Huang, Y.; Loera, S.; Xue, L.; Zhang, S.; Chu, P.; Zheng, S.; Yen, Y. Prognostic and therapeutic significance of ribonucleotide reductase small subunit M2 in estrogen-negative breast cancers. BMC Cancer 2014, 14, 664. [Google Scholar] [CrossRef] [Green Version]
  121. Chen, W.-X.; Yang, L.-G.; Xu, L.-Y.; Cheng, L.; Qian, Q.; Sun, L.; Zhu, Y.-L. bioinformatics analysis revealing prognostic significance of RRM2 gene in breast cancer. Biosci. Rep. 2019, 39. [Google Scholar] [CrossRef] [Green Version]
  122. Mercier, I.; Casimiro, M.C.; Wang, C.; Rosenberg, A.L.; Quong, J.; Minkeu, A.; Allen, K.G.; Danilo, C.; Sotgia, F.; Bonuccelli, G.; et al. Human breast cancer-associated fibroblasts (CAFs) show caveolin-1 downregulation and RB tumor suppressor functional inactivation: Implications for the response to hormonal therapy. Cancer Biol. Ther. 2008, 7, 1212–1225. [Google Scholar] [CrossRef] [Green Version]
  123. Kawano, Y.; Inada, Y.; Sasano, T.; Nishimura, N.; Hata, H.; Matsuoka, M. The purine metabolic enzyme AMPD1 is a novel therapeutic target for multiple myeloma. Blood 2018, 132, 5614. [Google Scholar] [CrossRef]
  124. Wang, N.; Li, Y.; Zhou, J. Downregulation of ribonucleotide reductase subunits M2 induces apoptosis and G1 arrest of cervical cancer cells. Oncol. Lett. 2018, 15, 3719–3725. [Google Scholar] [CrossRef] [Green Version]
  125. Mazzu, Y.Z.; Armenia, J.; Chakraborty, G.; Yoshikawa, Y.; Coggins, S.A.; Nandakumar, S.; Gerke, T.A.; Pomerantz, M.M.; Qiu, X.; Zhao, H.; et al. A novel mechanism driving poor-prognosis prostate cancer: Overexpression of the DNA repair gene, ribonucleotide reductase small subunit M2 (RRM2). Clin. Cancer Res. 2019, 25, 4480–4492. [Google Scholar] [CrossRef] [Green Version]
  126. Zhuang, S.; Li, L.; Zang, Y.; Li, G.; Wang, F. RRM2 elicits the metastatic potential of breast cancer cells by regulating cell invasion, migration and VEGF expression via the PI3K/AKT signaling. Oncol. Lett. 2020, 19, 3349–3355. [Google Scholar] [CrossRef] [Green Version]
  127. Pallante, P.; Terracciano, L.; Carafa, V.; Schneider, S.; Zlobec, I.; Lugli, A.; Bianco, M.; Ferraro, A.; Sacchetti, S.; Troncone, G.; et al. The loss of the CBX7 gene expression represents an adverse prognostic marker for survival of colon carcinoma patients. Eur. J. Cancer 2010, 46, 2304–2313. [Google Scholar] [CrossRef]
  128. Karamitopoulou, E.; Pallante, P.; Zlobec, I.; Tornillo, L.; Carafa, V.; Schaffner, T.; Borner, M.; Diamantis, I.; Esposito, F.; Brunner, T.; et al. Loss of the CBX7 protein expression correlates with a more aggressive phenotype in pancreatic cancer. Eur. J. Cancer 2010, 46, 1438–1444. [Google Scholar] [CrossRef] [Green Version]
  129. Clermont, P.-L.; Sun, L.; Crea, F.; Thu, K.L.; Zhang, A.; Parolia, A.; Lam, W.L.; Helgason, C.D. Genotranscriptomic meta-analysis of the polycomb gene CBX2 in human cancers: Initial evidence of an oncogenic role. Br. J. Cancer 2014, 111, 1663–1672. [Google Scholar] [CrossRef]
  130. Zheng, S.; Lv, P.; Su, J.; Miao, K.; Xu, H.; Li, M. Overexpression of CBX2 in breast cancer promotes tumor progression through the PI3K/AKT signaling pathway. Am. J. Transl. Res. 2019, 11, 1668–1682. [Google Scholar]
  131. Klintman, M.; Buus, R.; Cheang, M.C.U.; Sheri, A.; Smith, I.E.; Dowsett, M. Changes in expression of genes representing key biologic processes after neoadjuvant chemotherapy in breast cancer, and prognostic implications in residual disease. Clin. Cancer Res. 2016, 22, 2405–2416. [Google Scholar] [CrossRef] [Green Version]
  132. Svensson, R.U.; Parker, S.J.; Eichner, L.J.; Kolar, M.J.; Wallace, M.; Brun, S.N.; Lombardo, P.S.; Van Nostrand, J.L.; Hutchins, A.; Vera, L.; et al. Inhibition of acetyl-CoA carboxylase suppresses fatty acid synthesis and tumor growth of non-small-cell lung cancer in preclinical models. Nat. Med. 2016, 22, 1108–1119. [Google Scholar] [CrossRef] [Green Version]
  133. Jones, J.E.C.; Esler, W.P.; Patel, R.; Lanba, A.; Vera, N.B.; Pfefferkorn, J.A.; Vernochet, C. Inhibition of acetyl-CoA carboxylase 1 (ACC1) and 2 (ACC2) reduces proliferation and de novo lipogenesis of EGFRvIII human glioblastoma cells. PLoS ONE 2017, 12, e0169566. [Google Scholar] [CrossRef] [Green Version]
  134. Bai, J.; Zhang, X.; Kang, X.; Jin, L.; Wang, P.; Wang, Z. Screening of core genes and pathways in breast cancer development via comprehensive analysis of multi gene expression datasets. Oncol. Lett. 2019, 18, 5821–5830. [Google Scholar] [CrossRef] [Green Version]
  135. Jeong, D.; Ham, J.; Kim, H.W.; Kim, H.; Ji, H.W.; Yun, S.H.; Park, J.E.; Lee, K.S.; Jo, H.; Han, J.H.; et al. ELOVL2: A novel tumor suppressor attenuating tamoxifen resistance in breast cancer. Am. J. Cancer Res. 2021, 11, 2568–2589. [Google Scholar]
  136. Dang, L.; White, D.W.; Gross, S.; Bennett, B.D.; Bittinger, M.A.; Driggers, E.M.; Fantin, V.R.; Jang, H.G.; Jin, S.; Keenan, M.C.; et al. Cancer-associated IDH1 mutations produce 2-hydroxyglutarate. Nature 2009, 462, 739–744. [Google Scholar] [CrossRef] [Green Version]
  137. Figueroa, M.E.; Abdel-Wahab, O.; Lu, C.; Ward, P.S.; Patel, J.; Shih, A.; Li, Y.; Bhagwat, N.; Vasanthakumar, A.; Fernandez, H.F.; et al. Leukemic IDH1 and IDH2 mutations result in a hypermethylation phenotype, disrupt TET2 function, and impair hematopoietic differentiation. Cancer Cell 2010, 18, 553–567. [Google Scholar] [CrossRef] [Green Version]
  138. Chen, X.; Cao, Q.; Liao, R.; Wu, X.; Xun, S.; Huang, J.; Dong, C. Loss of ABAT-mediated GABAergic system promotes basal-like breast cancer progression by activating Ca2+−NFAT1 axis. Theranostics 2019, 9, 34–47. [Google Scholar] [CrossRef]
  139. Jansen, M.P.H.M.; Sas, L.; Sieuwerts, A.M.; Van Cauwenberghe, C.; Ramirez-Ardila, D.; Look, M.; Ruigrok-Ritstier, K.; Finetti, P.; Bertucci, F.; Timmermans, M.M.; et al. Decreased expression of ABAT and STC2 hallmarks ER-positive inflammatory breast cancer and endocrine therapy resistance in advanced disease. Mol. Oncol. 2015, 9, 1218–1233. [Google Scholar] [CrossRef] [Green Version]
  140. Tenori, L.; Oakman, C.; Claudino, W.M.; Bernini, P.; Cappadona, S.; Nepi, S.; Biganzoli, L.; Arbushites, M.C.; Luchinat, C.; Bertini, I.; et al. Exploration of serum metabolomic profiles and outcomes in women with metastatic breast cancer: A pilot study. Mol. Oncol. 2012, 6, 437–444. [Google Scholar] [CrossRef] [Green Version]
  141. Cao, M.D.; Giskeødegård, G.F.; Bathen, T.F.; Sitter, B.; Bofin, A.; Lønning, P.E.; Lundgren, S.; Gribbestad, I.S. Prognostic value of metabolic response in breast cancer patients receiving neoadjuvant chemotherapy. BMC Cancer 2012, 12, 39. [Google Scholar] [CrossRef] [Green Version]
  142. Wei, S.; Liu, L.; Zhang, J.; Bowers, J.; Gowda, G.A.N.; Seeger, H.; Fehm, T.; Neubauer, H.J.; Vogel, U.; Clare, S.E.; et al. Metabolomics approach for predicting response to neoadjuvant chemotherapy for breast cancer. Mol. Oncol. 2013, 7, 297–307. [Google Scholar] [CrossRef] [Green Version]
  143. Choi, J.S.; Baek, H.-M.; Kim, S.; Kim, M.J.; Youk, J.H.; Moon, H.J.; Kim, E.-K.; Nam, Y.K. Magnetic resonance metabolic profiling of breast cancer tissue obtained with core needle biopsy for predicting pathologic response to neoadjuvant chemotherapy. PLoS ONE 2013, 8, e83866. [Google Scholar] [CrossRef]
  144. Miolo, G.; Muraro, E.; Caruso, D.; Crivellari, D.; Ash, A.; Scalone, S.; Lombardi, D.; Rizzolio, F.; Giordano, A.; Corona, G. Pharmacometabolomics study identifies circulating spermidine and tryptophan as potential biomarkers associated with the complete pathological response to trastuzumab-paclitaxel neoadjuvant therapy in HER-2 positive breast cancer. Oncotarget 2016, 7, 39809–39822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  145. Yoneten, K.K.; Kasap, M.; Akpinar, G.; Gunes, A.; Gurel, B.; Utkan, N.Z. Comparative proteome analysis of breast cancer tissues highlights the importance of glycerol-3-phosphate dehydrogenase 1 and monoacylglycerol lipase in breast cancer metabolism. Cancer Genom. Proteom. 2019, 16, 377–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  146. Shenoy, A.; Belugali Nataraj, N.; Perry, G.; Loayza Puch, F.; Nagel, R.; Marin, I.; Balint, N.; Bossel, N.; Pavlovsky, A.; Barshack, I.; et al. Proteomic patterns associated with response to breast cancer neoadjuvant treatment. Mol. Syst. Biol. 2020, 16. [Google Scholar] [CrossRef]
  147. Yamashita, S.; Hattori, N.; Fujii, S.; Yamaguchi, T.; Takahashi, M.; Hozumi, Y.; Kogawa, T.; El-Omar, O.; Liu, Y.-Y.; Arai, N.; et al. Multi-omics analyses identify HSD17B4 methylation-silencing as a predictive and response marker of HER2-positive breast cancer to HER2-directed therapy. Sci. Rep. 2020, 10, 15530. [Google Scholar] [CrossRef] [PubMed]
  148. Colombié, M.; Jézéquel, P.; Rubeaux, M.; Frenel, J.-S.; Bigot, F.; Seegers, V.; Campone, M. The EPICURE study: A pilot prospective cohort study of heterogeneous and massive data integration in metastatic breast cancer patients. BMC Cancer 2021, 21, 333. [Google Scholar] [CrossRef]
  149. Cava, C.; Sabetian, S.; Castiglioni, I. Patient-specific network for personalized breast cancer therapy with multi-omics data. Entropy 2021, 23, 225. [Google Scholar] [CrossRef]
  150. Kim, H.; Park, J.; Won, H.; Lee, J.; Kong, G. CBX7 Inhibits breast tumorigenicity through DKK-1-mediated suppression of the Wnt/β-catenin pathway. FASEB J. 2015, 29, 300–313. [Google Scholar] [CrossRef]
  151. Piqué, D.G.; Montagna, C.; Greally, J.M.; Mar, J.C. A novel approach to modelling transcriptional heterogeneity identifies the oncogene candidate CBX2 in invasive breast carcinoma. Br. J. Cancer 2019, 120, 746–753. [Google Scholar] [CrossRef] [Green Version]
  152. Clermont, P.-L.; Crea, F.; Chiang, Y.T.; Lin, D.; Zhang, A.; Wang, J.Z.L.; Parolia, A.; Wu, R.; Xue, H.; Wang, Y.; et al. Identification of the epigenetic reader CBX2 as a potential drug target in advanced prostate cancer. Clin. Epigenetics 2016, 8, 16. [Google Scholar] [CrossRef] [Green Version]
  153. Xiao, X.; Huang, X.; Ye, F.; Chen, B.; Song, C.; Wen, J.; Zhang, Z.; Zheng, G.; Tang, H.; Xie, X. The miR-34a-LDHA axis regulates glucose metabolism and tumor growth in breast cancer. Sci. Rep. 2016, 6, 21735. [Google Scholar] [CrossRef] [Green Version]
  154. Li, L.; Kang, L.; Zhao, W.; Feng, Y.; Liu, W.; Wang, T.; Mai, H.; Huang, J.; Chen, S.; Liang, Y.; et al. MiR-30a-5p suppresses breast tumor growth and metastasis through inhibition of LDHA-mediated warburg effect. Cancer Lett. 2017, 400, 89–98. [Google Scholar] [CrossRef]
  155. Wang, Z.-Y.; Loo, T.Y.; Shen, J.-G.; Wang, N.; Wang, D.-M.; Yang, D.-P.; Mo, S.-L.; Guan, X.-Y.; Chen, J.-P. LDH-A silencing suppresses breast cancer tumorigenicity through induction of oxidative stress mediated mitochondrial pathway apoptosis. Breast Cancer Res. Treat. 2012, 131, 791–800. [Google Scholar] [CrossRef] [PubMed]
  156. Marini, C.; Salani, B.; Massollo, M.; Amaro, A.; Esposito, A.I.; Orengo, A.M.; Capitanio, S.; Emionite, L.; Riondato, M.; Bottoni, G.; et al. Direct inhibition of hexokinase activity by metformin at least partially impairs glucose metabolism and tumor growth in experimental breast cancer. Cell Cycle 2013, 12, 3490–3499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  157. Al-Ziaydi, A.G.; Al-Shammari, A.M.; Hamzah, M.I.; Kadhim, H.S.; Jabir, M.S. Hexokinase inhibition using D-mannoheptulose enhances oncolytic newcastle disease virus-mediated killing of breast cancer cells. Cancer Cell Int. 2020, 20, 420. [Google Scholar] [CrossRef]
  158. Shah, K.N.; Mehta, K.R.; Peterson, D.; Evangelista, M.; Livesey, J.C.; Faridi, J.S. AKT-induced tamoxifen resistance is overturned by RRM2 inhibition. Mol. Cancer Res. 2014, 12, 394–407. [Google Scholar] [CrossRef] [Green Version]
  159. Wong, C.C.; Qian, Y.; Yu, J. Interplay between epigenetics and metabolism in oncogenesis: Mechanisms and therapeutic approaches. Oncogene 2017, 36, 3359–3374. [Google Scholar] [CrossRef]
  160. Mazor, T.; Chesnelong, C.; Pankov, A.; Jalbert, L.E.; Hong, C.; Hayes, J.; Smirnov, I.V.; Marshall, R.; Souza, C.F.; Shen, Y.; et al. Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1. Proc. Natl. Acad. Sci. USA 2017, 114, 10743–10748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  161. Park, J.-W.; Turcan, Ş. Epigenetic reprogramming for targeting IDH-mutant malignant gliomas. Cancers 2019, 11, 1616. [Google Scholar] [CrossRef] [Green Version]
  162. Struys, E.A.; Verhoeven, N.M.; Ten Brink, H.J.; Wickenhagen, W.V.; Gibson, K.M.; Jakobs, C. Kinetic characterization of human hydroxyacid-oxoacid transhydrogenase: Relevance to D-2-hydroxyglutaric and γ-hydroxybutyric acidurias. J. Inherit. Metab. Dis. 2005, 28, 921–930. [Google Scholar] [CrossRef]
  163. Mishra, P.; Tang, W.; Putluri, V.; Dorsey, T.H.; Jin, F.; Wang, F.; Zhu, D.; Amable, L.; Deng, T.; Zhang, S.; et al. ADHFE1 is a breast cancer oncogene and induces metabolic reprogramming. J. Clin. Investig. 2018, 128, 323–340. [Google Scholar] [CrossRef] [Green Version]
  164. Mishra, P.; Tang, W.; Ambs, S. ADHFE1 is a MYC-linked oncogene that induces metabolic reprogramming and cellular de-differentiation in breast cancer. Mol. Cell. Oncol. 2018, 5, e1432260. [Google Scholar] [CrossRef] [Green Version]
  165. Chajès, V.; Cambot, M.; Moreau, K.; Lenoir, G.M.; Joulin, V. Acetyl-CoA carboxylase alpha is essential to breast cancer cell survival. Cancer Res. 2006, 66, 5287–5294. [Google Scholar] [CrossRef] [Green Version]
  166. Mauvoisin, D.; Charfi, C.; Lounis, A.M.; Rassart, E.; Mounier, C. Decreasing stearoyl-CoA desaturase-1 expression inhibits β-catenin signaling in breast cancer cells. Cancer Sci. 2013, 104, 36–42. [Google Scholar] [CrossRef] [PubMed]
  167. Bao, J.; Zhu, L.; Zhu, Q.; Su, J.; Liu, M.; Huang, W. SREBP-1 is an independent prognostic marker and promotes invasion and migration in breast cancer. Oncol. Lett. 2016, 12, 2409–2416. [Google Scholar] [CrossRef] [Green Version]
  168. Nie, L.-Y.; Lu, Q.-T.; Li, W.-H.; Yang, N.; Dongol, S.; Zhang, X.; Jiang, J. Sterol regulatory element-binding protein 1 is required for ovarian tumor growth. Oncol. Rep. 2013, 30, 1346–1354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  169. Wen, Y.-A.; Xiong, X.; Zaytseva, Y.Y.; Napier, D.L.; Vallee, E.; Li, A.T.; Wang, C.; Weiss, H.L.; Evers, B.M.; Gao, T. Downregulation of SREBP inhibits tumor growth and initiation by altering cellular metabolism in colon cancer. Cell Death Dis. 2018, 9, 265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  170. Zhang, H.; Liu, S.; Cai, Z.; Dong, W.; Ye, J.; Cai, Z.; Han, Z.; Liang, Y.; Zhuo, Y.; Luo, Y.; et al. Down-regulation of ACACA suppresses the malignant progression of prostate cancer through inhibiting mitochondrial potential. J. Cancer 2021, 12, 232–243. [Google Scholar] [CrossRef]
  171. Yokoyama, C.; Wang, X.; Briggs, M.R.; Admon, A.; Wu, J.; Hua, X.; Goldstein, J.L.; Brown, M.S. SREBP-1, a basic-helix-loop-helix-leucine zipper protein that controls transcription of the low density lipoprotein receptor gene. Cell 1993, 75, 187–197. [Google Scholar] [CrossRef]
  172. Jeon, T.-I.; Osborne, T.F. SREBPs: Metabolic integrators in physiology and metabolism. Trends Endocrinol. Metab. 2012, 23, 65–72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  173. Huang, W.-C.; Li, X.; Liu, J.; Lin, J.; Chung, L.W.K. Activation of androgen receptor, lipogenesis, and oxidative stress converged by SREBP-1 is responsible for regulating growth and progression of prostate cancer cells. Mol. Cancer Res. 2012, 10, 133–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  174. Miyachi, K.; Sawada, Y.; Shida, Y.; Sugawara, A.; Hisatomi, H. Lipogenic gene expression profile in patients with gastric cancer. Mol. Clin. Oncol. 2013, 1, 825–827. [Google Scholar] [CrossRef] [PubMed]
  175. Tan, Z.; Xiao, L.; Tang, M.; Bai, F.; Li, J.; Li, L.; Shi, F.; Li, N.; Li, Y.; Du, Q.; et al. Targeting CPT1A-mediated fatty acid oxidation sensitizes nasopharyngeal carcinoma to radiation therapy. Theranostics 2018, 8, 2329–2347. [Google Scholar] [CrossRef] [PubMed]
  176. Han, S.; Wei, R.; Zhang, X.; Jiang, N.; Fan, M.; Huang, J.H.; Xie, B.; Zhang, L.; Miao, W.; Butler, A.C.-P.; et al. CPT1A/2-mediated FAO enhancement—A metabolic target in radioresistant breast cancer. Front. Oncol. 2019, 9, 1201. [Google Scholar] [CrossRef]
  177. Xiong, X.; Wen, Y.-A.; Fairchild, R.; Zaytseva, Y.Y.; Weiss, H.L.; Evers, B.M.; Gao, T. Upregulation of CPT1A is essential for the tumor-promoting effect of adipocytes in colon cancer. Cell Death Dis. 2020, 11, 736. [Google Scholar] [CrossRef]
  178. Schlaepfer, I.R.; Joshi, M. CPT1A-mediated fat oxidation, mechanisms, and therapeutic potential. Endocrinology 2020, 161. [Google Scholar] [CrossRef]
  179. Pinu, F.R.; Beale, D.J.; Paten, A.M.; Kouremenos, K.; Swarup, S.; Schirra, H.J.; Wishart, D. Systems biology and multi-omics integration: Viewpoints from the metabolomics research community. Metabolites 2019, 9, 76. [Google Scholar] [CrossRef] [Green Version]
  180. Berg, M.; Vanaerschot, M.; Jankevics, A.; Cuypers, B.; Breitling, R.; Dujardin, J.-C. LC-MS metabolomics from study design to data-analysis—Using a versatile pathogen as a test case. Comput. Struct. Biotechnol. J. 2013, 4, e201301002. [Google Scholar] [CrossRef]
  181. Johnson, S.R.; Lange, B.M. Open-access metabolomics databases for natural product research: Present capabilities and future potential. Front. Bioeng. Biotechnol. 2015, 3, 22. [Google Scholar] [CrossRef] [Green Version]
  182. Dias, D.A.; Koal, T. Progress in metabolomics standardisation and its significance in future clinical laboratory medicine. EJIFCC 2016, 27, 331–343. [Google Scholar]
  183. Hasin, Y.; Seldin, M.; Lusis, A. Multi-omics approaches to disease. Genome Biol. 2017, 18, 83. [Google Scholar] [CrossRef] [PubMed]
  184. Menyhárt, O.; Győrffy, B. Multi-omics approaches in cancer research with applications in tumor subtyping, prognosis, and diagnosis. Comput. Struct. Biotechnol. J. 2021, 19, 949–960. [Google Scholar] [CrossRef]
  185. Misra, B.B.; Langefeld, C.D.; Olivier, M.; Cox, L.A. Integrated omics: Tools, advances, and future approaches. J. Mol. Endocrinol. 2018, 62, R21–R45. [Google Scholar] [CrossRef] [Green Version]
  186. Eicher, T.; Kinnebrew, G.; Patt, A.; Spencer, K.; Ying, K.; Ma, Q.; Machiraju, R.; Mathé, E.A. Metabolomics and multi-omics integration: A survey of computational methods and resources. Metabolites 2020, 10, 202. [Google Scholar] [CrossRef]
  187. Lewis, J.E.; Kemp, M.L. Integration of machine learning and genome-scale metabolic modeling identifies multi-omics biomarkers for radiation resistance. Nat. Commun. 2021, 12, 2700. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the most commonly used omic platforms for multi-omics studies.
Figure 1. Schematic representation of the most commonly used omic platforms for multi-omics studies.
Cancers 13 04544 g001
Figure 2. Schematic representation of the impact of metabolic changes on essential molecular processes associated with tumor progression and patient survival. Created with BioRender.com.
Figure 2. Schematic representation of the impact of metabolic changes on essential molecular processes associated with tumor progression and patient survival. Created with BioRender.com.
Cancers 13 04544 g002
Figure 3. Overview of metabolic-related therapeutic targets for the treatment of BC patients identified from multi-omics-based studies. α-KG: alpha-ketoglutarate, ACACA: acetyl-CoA carboxylase alpha, ADHFE1: alcohol dehydrogenase iron containing 1, CBX2: chromobox 2, CBX7: chromobox 7, CPT1: carnitine palmitoyltransferase 1A, CPT2: carnitine palmitoyltransferase 2, ECT: electron transport chain, ELOVL1: ELOVL fatty acid elongase 1, ELOVL2: ELOVL fatty acid elongase 2, FASN: fatty acid synthase, GLUT1: glucose transporter 1, HK: hexokinase, INSIG1: insulin-induced gene 1, LDH: lactate dehydrogenase, OXPHOS: oxidative phosphorylation, PDC: pyruvate dehydrogenase kinase, RRM2: ribonucleotide reductase regulatory subunit M2, SCAP: SREBF chaperone, SCD: stearoyl-CoA desaturase, TCA: tricarboxylic acid, THRSP: thyroid hormone-responsive. Created with BioRender.com.
Figure 3. Overview of metabolic-related therapeutic targets for the treatment of BC patients identified from multi-omics-based studies. α-KG: alpha-ketoglutarate, ACACA: acetyl-CoA carboxylase alpha, ADHFE1: alcohol dehydrogenase iron containing 1, CBX2: chromobox 2, CBX7: chromobox 7, CPT1: carnitine palmitoyltransferase 1A, CPT2: carnitine palmitoyltransferase 2, ECT: electron transport chain, ELOVL1: ELOVL fatty acid elongase 1, ELOVL2: ELOVL fatty acid elongase 2, FASN: fatty acid synthase, GLUT1: glucose transporter 1, HK: hexokinase, INSIG1: insulin-induced gene 1, LDH: lactate dehydrogenase, OXPHOS: oxidative phosphorylation, PDC: pyruvate dehydrogenase kinase, RRM2: ribonucleotide reductase regulatory subunit M2, SCAP: SREBF chaperone, SCD: stearoyl-CoA desaturase, TCA: tricarboxylic acid, THRSP: thyroid hormone-responsive. Created with BioRender.com.
Cancers 13 04544 g003
Table 1. Multi-omics studies focused on the identification of BC metabolic biomarkers.
Table 1. Multi-omics studies focused on the identification of BC metabolic biomarkers.
StudySample Omics DataMajor Findings *
Iqbal et al. [59]TissueM+T↑ glucose-6-phosphate, 3-phosphoglycerate, and fructose-6-phosphate
CBX2 and ↓ CBX7
Hilvo et al. [60]TissueL+TACACA, FASN, INSIG1, and SREBP1
↑ phospholipids, sphingomyelins, and ceramides
Luo et al. [61]Plasma + TissueM+T↓ guanine and hypoxanthine
PNP and HPRT1
Huang et al. [57]Plasma + Serum + TissueM+T↑ hypotaurine, glutamate and
↓ oxoglutarate
GAD1, CSAD and ↓ GPT, GPT2, GLUD1
Dowling et al. [62]SerumM+P↑ glutamate, 12-hydroxyeicosatetraenoic acid, andβ-hydroxybutyrate
↑ coagulation factor V and matrix metalloproteinase 1
ACACA: acetyl-CoA carboxylase alpha, BC: breast cancer, CBX2: chromobox 2, CBX7: chromobox 7, CSAD: cysteine sulfinic acid decarboxylase, GAD1: glutamate decarboxylase 1, GLUD1: glutamate dehydrogenase, GPT: glutamic-pyruvic transaminase, GPT2: glutamic-pyruvic transaminase 2, FASN: fatty acid synthase, HPRT1: hypoxanthine phosphoribosyltransferase, INSIG1: insulin-induced gene 1, L: lipidomics, M: metabolomics, P: proteomics, PNP: purine nucleoside phosphorylase, SREBP1: sterol regulatory element-binding transcription factor 1, T: transcriptomics. * Direction of variation, considering the healthy group as a reference.
Table 2. Multi-omics studies focused on BC tumor metabolic subtyping.
Table 2. Multi-omics studies focused on BC tumor metabolic subtyping.
StudySample Omics DataGroup
Comparison
Major Findings *
Moestue et al. [63]TissueM+TBasal- vs.
luminal-like
↑ GPC/PCho and glycine
CHKA/B
PLA2G4A, PLB1, CHDH and SARDH
Grinde et al. [64]TissueM+TBasal- vs.
luminal-like
↓ PCho/GPC
Putluri et al. [65]Cell lines + TissueM+TBasal- vs.
luminal-like
↓ phenylalanine, tryptophan, tyrosine,
BCAA, lauric acid and oleic acid
↑ guanine, adenine, thymine, uracil, xanthine, and guanosine
Mahendralingam et al. [66]TissueP+TBasal- vs.
luminal-like
↑ glycolysis (PFKM, ALDOC, GAPDH, and PKM) and ↓ OXPHOS
Tang et al. [67]TissueG+M+TER+ vs. ER-↓ carnitine derivates and short- and medium-chain fatty acids
↑ long-chain fatty acids and monoacylglycerols
Barupal et al. [68]TissueM+P+TER+ vs. ER-↑ R5P, adenine, guanosine, guanine, xanthine, and hypoxanthine and
β-alanine
G6PD, PGD, TKT, PGM1, RPIA, DERA
Hilvo et al. [60]TissueL+TER+ vs. ER-↑ palmitate and myristic acid
Haukaas et al. [48]TissueM+T+PPrimary BC subtypingMC1: ↑ GPC, PCho ↓ acetate and
glutamine
CHKAALDH and GLS
MC2: ↑ glucose
MC3: ↑ lactate and alanine
Gong et al. [69]TissueM+TTNBC
subtyping
MPS1: ↑ myristic, palmitoleic, oleic,
and arachidonic acid
ACACA, HMGCR, FASN,
SCD
MPS2: ↑ glucose 1-phosphate,
dihydroxyacetone phosphate,
lactate and
adenosine 3′ 5′ -diphosphate
and ↓ glucose
PFKP, ENO2, TYMS, CTPS1,
SLC2A1, SLC16A1
ACACA: acetyl-CoA carboxylase alpha, ALDH: aldehyde dehydrogenase, ALDOC: aldolase, fructose-bisphosphate C, BC: breast cancer, BCAA: branched-chain amino acid, CHDH: choline dehydrogenase, CHKA: choline kinase alpha, CHKB: choline kinase beta, CTPS1: CPT synthase 1, DERA: deoxyribose-phosphate aldolase, ENO2: enolase 2, ER: estrogen receptor, FASN: fatty acid synthase, G: genomics, GAPDH: glyceraldehyde-3-phosphate dehydrogenase, GLS: glutaminase, GPC: glycerophosphocholine, G6PD: glucose-6-phosphate dehydrogenase, HMGCR: HMG-CoA reductase, L: lipidomics, M: metabolomics, OXPHOS: oxidative phosphorylation, P: proteomics, PCho: phosphocholine, PKFM: phosphofructokinase (muscle), PFKP: phosphofructokinase (platelet), PGD: phosphogluconate dehydrogenase, PGM1: phosphoglucomutase 1, PKM: pyruvate kinase M1/2, PLA2G4A: phospholipase A2 group IVA, PLB1: phospholipase B1, PPP: pentose phosphate pathway, RPIA: ribose 5-phosphate isomerase 1, R5P: ribose-5-phosphate, SARDH: sarcosine dehydrogenase, SCD: stearoyl-CoA desaturase, SLC2A1: solute carrier family 2 member 1, SLC16A1: solute carrier family 16 member 2, T: transcriptomics, TCA: tricarboxylic acid, TKT: transketolase, TNBC: triple negative breast cancer, TYMS: thymidylate synthetase. * Direction of variation, considering luminal or ER+ BC group as a reference.
Table 3. Multi-omics studies focused on identifying metabolic alterations associated with BC prognosis.
Table 3. Multi-omics studies focused on identifying metabolic alterations associated with BC prognosis.
StudySample Omics DataMajor Findings *
Putluri et al. [65]Cell lines + TissueM+TRRM2 (pyrimidine metabolism)
Luo et al. [61]Blood + TissueM+TRRM2 (pyrimidine metabolism) and ↓ AMPD1 (de novo purine metabolism)
Iqbal et al. [59]TissueM+TCBX2 and ↓ CBX7 (glycolysis)
Camarda et al. [102]Cell lines + TissueM+T↓ ACC2 (FAO)
Kang et al. [116]Cell linesL+TELOVL2 (lipid synthesis)
Terunuma et al. [107] TissueM+T+P+E↑ 2HG, SAM and SAH
IDH2 (glutamine metabolism)
Budczies et al. [106]TissueM+TABAT, ↑ β-alanine
(β-alanine metabolism)
ABAT: 4-aminobutyrate aminotransferase, ACC2: acetyl-CoA carboxylase 2, AMPD1: adenosine monophosphate deaminase 1, CBX2: chromobox 2, CBX7: chromobox 7, E: epigenomics, ELOVL2: ELOVL fatty acid elongase 2, FAO: fatty acid oxidation, L: lipidomics, IDH2: isocitrate dehydrogenase (NADP(+)) 2, M: metabolomics, P: proteomics, RRM2: ribonucleotide reductase regulatory subunit M2, SAH: S-adenosylhomocysteine, SAM: S-adenosyl- methionine, T: transcriptomics, 2HG: 2-hydroxyglutarate. * Direction of metabolic alterations directly correlated with worse BC patients’ outcomes.
Table 4. Multi-omics studies focused on developing new therapeutic strategies for the treatment of BC.
Table 4. Multi-omics studies focused on developing new therapeutic strategies for the treatment of BC.
StudyOmics Data BC SubtypePotential Targets
Iqbal et al. [59]M+TTNBC and luminal-likeCBX2 and CBX7
Gong et al. [69]M+TTNBCFASN and LDH
Mahendralingam et al. [66]P+TBasal- and luminal-likeGLUT1, HK, LDH, and PDC
Putluri et al. [65]M+TBasal- and luminal-likeRRM2
Terunuma et al. [107] M+T+P+ETNBC and basal-likeADHFE1
Hilvo et al. [60]L+TTNBC, luminal- and
basal-like
ACACA, ELOVL1, FASN, INSIG1, SCAP, SCD and THRSP
Kang et al. [116]L+TLuminal-likeELOVL2
Camarda et al. [102]M+TTNBC and HER2 +CPT1 and CPT2
ACACA: acetyl-CoA carboxylase alpha, ADHFE1: alcohol dehydrogenase iron containing 1, CBX2: chromobox 2, CBX7: chromobox 7, CPT1: carnitine palmitoyltransferase 1A, CPT2: carnitine palmitoyltransferase 2, E: epigenomics, ELOVL1: ELOVL fatty acid elongase 1, ELOVL2: ELOVL fatty acid elongase 2, FASN: fatty acid synthase, GLUT1: glucose transporter 1, HK: hexokinase, IDH2: isocitrate dehydrogenase (NADP(+)) 2, INSIG1: insulin-induced gene 1, L: lipidomics, LDH: lactate dehydrogenase, M: metabolomics, P: proteomics, PDC: pyruvate dehydrogenase kinase, RRM2: ribonucleotide reductase regulatory subunit M2, SCAP: SREBF chaperone, SCD: stearoyl-CoA desaturase, T: transcriptomics, TNBC: triple negative breast cancer, THRSP: thyroid hormone-responsive.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gómez-Cebrián, N.; Domingo-Ortí, I.; Poveda, J.L.; Vicent, M.J.; Puchades-Carrasco, L.; Pineda-Lucena, A. Multi-Omic Approaches to Breast Cancer Metabolic Phenotyping: Applications in Diagnosis, Prognosis, and the Development of Novel Treatments. Cancers 2021, 13, 4544. https://doi.org/10.3390/cancers13184544

AMA Style

Gómez-Cebrián N, Domingo-Ortí I, Poveda JL, Vicent MJ, Puchades-Carrasco L, Pineda-Lucena A. Multi-Omic Approaches to Breast Cancer Metabolic Phenotyping: Applications in Diagnosis, Prognosis, and the Development of Novel Treatments. Cancers. 2021; 13(18):4544. https://doi.org/10.3390/cancers13184544

Chicago/Turabian Style

Gómez-Cebrián, Nuria, Inés Domingo-Ortí, José Luis Poveda, María J. Vicent, Leonor Puchades-Carrasco, and Antonio Pineda-Lucena. 2021. "Multi-Omic Approaches to Breast Cancer Metabolic Phenotyping: Applications in Diagnosis, Prognosis, and the Development of Novel Treatments" Cancers 13, no. 18: 4544. https://doi.org/10.3390/cancers13184544

APA Style

Gómez-Cebrián, N., Domingo-Ortí, I., Poveda, J. L., Vicent, M. J., Puchades-Carrasco, L., & Pineda-Lucena, A. (2021). Multi-Omic Approaches to Breast Cancer Metabolic Phenotyping: Applications in Diagnosis, Prognosis, and the Development of Novel Treatments. Cancers, 13(18), 4544. https://doi.org/10.3390/cancers13184544

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop