MUC1 Tissue Expression and Its Soluble Form CA15-3 Identify a Clear Cell Renal Cell Carcinoma with Distinct Metabolic Profile and Poor Clinical Outcome

An altered metabolism is involved in the development of clear cell renal carcinoma (ccRCC). MUC1 overexpression has been found to be associated with advanced disease and poor prognosis. In this study, we evaluated the metabolomic profile of human ccRCC, according to MUC1 expression, and integrated it with transcriptomic data. Moreover, we analyzed the role of MUC1 in sustaining ccRCC aggressiveness and the prognostic value of its soluble form CA15-3. Integrated metabolomic and transcriptomic analysis showed that MUC1-expressing ccRCC was characterized by metabolic reprogramming involving the glucose and lipid metabolism pathway. In addition, primary renal cancer cells treated with a small interfering RNA targeting MUC1 (siMUC1) migrated and proliferated at a slower rate than untreated cancer cells. After cisplatin treatment, the death rate of cancer cells treated with siMUC1 was significantly greater than that of untreated cells. Kaplan–Meier curves showed significant differences in CSS and PFS among groups of patients with high versus low levels of CA15-3. In a multivariate analysis, CA15-3 was an independent adverse prognostic factor for cancer-specific and progression-free survival. In conclusion, MUC1 expressing ccRCC is characterized by a particular metabolic reprogramming. The inhibition of MUC1 expression decreases cell motility and viability and improves cisplatin susceptibility, suggesting that this pathway can regulate de novo chemotherapy resistance in ccRCC.


Introduction
Renal cell carcinoma (RCC) represents around 3-5% of all cancers, and recent estimates have calculated that in 2022 in the United States, 79,000 new cases of kidney cancer (50,290 in men and 28,710 in women) will be diagnosed and about 13,920 people (8960 men and 4960 women) will die from this disease [1].
Although the pathogenesis of RCC is still not fully understood, recent evidence has confirmed that RCC is fundamentally a metabolic disease. In fact, many studies have showed that alterations in metabolic pathways are involved in the development of RCC, and that in this cancer many mutated genes have a fundamental role in controlling cell metabolic activities [2][3][4][5]. Furthermore, the recent re-evaluation of cancer as a disease of cell metabolism has led to the discovery of specific oncometabolites that play a key role in tumor growth and progression, especially in other urologic tumors such as prostate and bladder cancer [6][7][8][9]. The data of the cancer genome atlas (TCGA) research network have showed that a reprogrammed metabolism represents a major feature in clear cell RCC (ccRCC) [10,11], and as a result of the introduction of high-throughput technology, researchers have been 2. Results 2.1. Global Metabolic Profile Distinguishes MUC1H from MUC1L ccRCC MUC1 expression was evaluated by immunohistochemistry on normal and pathological tissues, stratifying samples according to staining intensity as tumor with high MUC1 expression (MUC1 H ) and tumor with low MUC1 expression (MUC1 L ) ( Figure 1A).
Next, untargeted metabolomic analysis was performed on 56 kidney-derived tissues, including 20 normal tissues, 12 MUC1 H and 24 MUC1 L ccRCC (Supplementary Table S1), using LC-MS and GC-MS platforms. In total, 516 metabolites were identified, and 116 were found to be differentially expressed in MUC1 H compared to MUC1 L tumors ( Figure 1B). The application of principal component analysis (PCA) to differentiate normal and pathological samples as a function of the global tissue metabolome showed that they cluster into three separate groups and that they can be metabolically differentiated ( Figure 1C,D). cluster into three separate groups and that they can be metabolically differentiated ( Figure  1C,D). To obtain a global overview of altered biochemical processes, we performed a metabolite set enrichment analysis (MSEA) using MetaboAnalyst 4.0 [35]. This functional approach showed that alterations in glucose utilization, including PPP, one-carbon and lipid metabolism (linoleic acid and glycerolipids metabolism, and unsaturated fatty acid biosynthesis) had the highest impact in MUC1 H cancer cell metabolism compared to MUC1 L tumors ( Figure 2).

MUC1 Expression Distinguishes ccRCC with Different Perturbations of Glucose Metabolism
The Warburg effect is defined as the increase in the rate of glucose uptake and preferential production of lactate, even in the presence of oxygen. We previously showed that in ccRCC, oncogenic signaling pathways promote cancer through rerouting the sugar metabolism [14]. MUC1 H renal cancer cells showed a signature of an increased glucose uptake (glucose transporter GLUT1 was overexpressed in this subset compared to MUC1 L To obtain a global overview of altered biochemical processes, we performed a metabolite set enrichment analysis (MSEA) using MetaboAnalyst 4.0 [35]. This functional approach showed that alterations in glucose utilization, including PPP, one-carbon and lipid metabolism (linoleic acid and glycerolipids metabolism, and unsaturated fatty acid biosynthesis) had the highest impact in MUC1 H cancer cell metabolism compared to MUC1 L tumors ( Figure 2).  Glycogenolysis was partially inhibited in MUC1 H ccRCC, and significant reduction in glycogen degradation products were observed. MUC1 H compared to MUC1 L tumor and normal tissues showed increased glucose levels associated with elevations in upstream glycolytic intermediates (glucose 6-phosphate and fructose 6-phosphate), reductions in downstream intermediates (fructose 1,6-bisphosphate, 3-phosphoglycerate, 2-phosphoglycerate and phosphoenolpyruvate) and increased lactate production (MUC1 H vs. MUC1 L vs. normal) (Figure 3). Analysis of the glycolytic enzymes glucose-6-phosphate

MUC1 Expression Distinguishes ccRCC with Different Perturbations of Glucose Metabolism
The Warburg effect is defined as the increase in the rate of glucose uptake and preferential production of lactate, even in the presence of oxygen. We previously showed that in ccRCC, oncogenic signaling pathways promote cancer through rerouting the sugar metabolism [14]. MUC1 H renal cancer cells showed a signature of an increased glucose uptake (glucose transporter GLUT1 was overexpressed in this subset compared to MUC1 L cells, Figure 1) and utilization, with glucose levels significantly elevated, along with higher levels of other sugars (such as fructose) and their phospho-derivatives ( Figure 3). M2 (PKM2) demonstrated an increased expression of these proteins in MUC1 H ccRCC as compared with MUC1 L tumors and non-neoplastic tissue. Conversely, the expression of L-lactate dehydrogenase B chain (LDH-B) was significantly higher in normal kidney tissue, and the LDH-A/LDH-B ratio between the three groups was in accordance with the higher-efficiency lactate production observed in MUC1 H ccRCC ( Figure 4). These results had an expression profile compatible with that obtained by real-time PCR, as shown in Supplementary Figure 2. The increased lactate production found in MUC1 H ccRCC is likely related to the particular LDH isoform expressed in this tumor. LDH exhibits five isoforms assembled in tetramers of either of the two types of subunits, namely LDH-A and LDH-B. MUC1 H tumors showed a reduced expression of LDH-B and higher levels of LDH-A, thus increasing the levels of LDH-5 (LDH-A4) which is the most effective isoenzyme in producing lactate.
Taken together, these findings are consistent with the observation of significantly higher hypoxia-inducible factor 1-α (HIF-1α) levels in MUC1H ccRCC (p < 0.0001; Figure  4), indicating that MUC1 and HIF-1α transcriptionally activate these metabolic genes. In fact, it has been shown that MUC1 physically occupies the promoter of glycolytic genes such as GLUT1 and LDHA, and also MUC1 interacts with HIF-1α, regulating its stability and activity [36]. This signature, characterized by the higher levels of glucose 6-phosphate and fructose 6-phosphate, associated with the increased tissue level of ribulose 5-phosphate, elevated levels of G6PDH and an increased activity of this enzyme (Figure 4), suggested an increased shunting of the metabolites of the upper part of glycolysis into the pentose phosphate pathway in MUC1 H tumors. Glycogenolysis was partially inhibited in MUC1 H ccRCC, and significant reduction in glycogen degradation products were observed. MUC1 H compared to MUC1 L tumor and normal tissues showed increased glucose levels associated with elevations in upstream glycolytic intermediates (glucose 6-phosphate and fructose 6-phosphate), reductions in downstream intermediates (fructose 1,6-bisphosphate, 3-phosphoglycerate, 2-phosphoglycerate and phosphoenolpyruvate) and increased lactate production (MUC1 H vs. MUC1 L vs. normal) ( Figure 3). Analysis of the glycolytic enzymes glucose-6-phosphate isomerase (G6PI), L-lactate dehydrogenase A chain (LDH-A) and pyruvate kinase isoform M2 (PKM2) demonstrated an increased expression of these proteins in MUC1 H ccRCC as compared with MUC1 L tumors and non-neoplastic tissue. Conversely, the expression of L-lactate dehydrogenase B chain (LDH-B) was significantly higher in normal kidney tissue, and the LDH-A/LDH-B ratio between the three groups was in accordance with the higher-efficiency lactate production observed in MUC1 H ccRCC ( Figure 4). These results had an expression profile compatible with that obtained by real-time PCR, as shown in Supplementary Figure S2.
The increased lactate production found in MUC1 H ccRCC is likely related to the particular LDH isoform expressed in this tumor. LDH exhibits five isoforms assembled in tetramers of either of the two types of subunits, namely LDH-A and LDH-B. MUC1 H tumors showed a reduced expression of LDH-B and higher levels of LDH-A, thus increasing the levels of LDH-5 (LDH-A4) which is the most effective isoenzyme in producing lactate.

Alterations in Pentose Phosphate Pathway (PPP) Activity in MUC1 H Versus MUC1 L ccRCC
In addition to the described alterations in glycolytic intermediates, metabolites related to PPP were also perturbed in the tumor samples. The PPP plays an important role in cancer cell biology, generating precursors for the nucleotides' biosynthesis and producing NADPH for oxidoreductive and anabolic reactions. An increased PPP activity was observed mainly in MUC1 H ccRCC, as shown by higher levels of the intermediates sedoheptulose 7-phosphate, ribose 5-phosphate and ribulose 5-phosphate/xylulose 5phosphate (isobaric compounds) ( Figure 5).
These findings, in association with the increased enzymatic activity of G6PDH and overexpression of transketolase (TKT), are suggestive of an increased activation of this pathway for anabolic reactions, especially in MUC1 H tumors ( Figure 4). Taken together, these findings suggest that the metabolic flux through glycolysis is differentially partitioned, especially in MUC1 H ccRCC. In fact, while the sugars produced in the upper part of glycolysis are shifted to the PPP to promote anabolic reactions and redox homeostasis, the triose phosphates generated in the lower part are diverted towards the TCA cycle and one-carbon metabolism. Taken together, these findings are consistent with the observation of significantly higher hypoxia-inducible factor 1-α (HIF-1α) levels in MUC1H ccRCC (p < 0.0001; Figure 4), indicating that MUC1 and HIF-1α transcriptionally activate these metabolic genes. In fact, it has been shown that MUC1 physically occupies the promoter of glycolytic genes such as GLUT1 and LDHA, and also MUC1 interacts with HIF-1α, regulating its stability and activity [36]. This signature, characterized by the higher levels of glucose 6-phosphate and fructose 6-phosphate, associated with the increased tissue level of ribulose 5-phosphate, elevated levels of G6PDH and an increased activity of this enzyme (Figure 4), suggested an increased shunting of the metabolites of the upper part of glycolysis into the pentose phosphate pathway in MUC1 H tumors.

Alterations in Pentose Phosphate Pathway (PPP) Activity in MUC1 H Versus MUC1 L ccRCC
In addition to the described alterations in glycolytic intermediates, metabolites related to PPP were also perturbed in the tumor samples. The PPP plays an important role in cancer cell biology, generating precursors for the nucleotides' biosynthesis and producing NADPH for oxidoreductive and anabolic reactions. An increased PPP activity was observed mainly in MUC1 H ccRCC, as shown by higher levels of the intermediates sedoheptulose 7-phosphate, ribose 5-phosphate and ribulose 5-phosphate/xylulose 5-phosphate (isobaric compounds) ( Figure 5).

Changes in Tricarboxylic Acid (TCA) Cycle Intermediates and Related Metabolites
In MUC1 H tumors, citrate and succinate were markedly increased when compared to MUC1 L ccRCC and normal tissue, whereas fumarate and malate were significantly reduced in cancer tissue ( Figure 6).  These findings, in association with the increased enzymatic activity of G6PDH and overexpression of transketolase (TKT), are suggestive of an increased activation of this pathway for anabolic reactions, especially in MUC1 H tumors ( Figure 4). Taken together, these findings suggest that the metabolic flux through glycolysis is differentially partitioned, especially in MUC1 H ccRCC. In fact, while the sugars produced in the upper part of glycolysis are shifted to the PPP to promote anabolic reactions and redox homeostasis, the triose phosphates generated in the lower part are diverted towards the TCA cycle and one-carbon metabolism.

Changes in Tricarboxylic Acid (TCA) Cycle Intermediates and Related Metabolites
In MUC1 H tumors, citrate and succinate were markedly increased when compared to MUC1 L ccRCC and normal tissue, whereas fumarate and malate were significantly reduced in cancer tissue ( Figure 6).
The altered levels of TCA cycle intermediates, in association with a reduced reactive oxygen species (ROS) production ( Figure 7A) and a decreased membrane potential ( Figure 7B), suggests that the rate of oxidative phosphorylation is impaired in MUC1 H tumors, which is consistent with NDUFA4L2 overexpression observed MUC1 H samples ( Figure 7C,D) [17].
The observed changes in the levels of Krebs cycle metabolites are suggestive of a metabolic signature that may be consistent with an increased glutaminolysis. In fact, nearly all the amino acids were significantly reduced in tumor tissue, with the exception of glutamine, glutamate and cysteine, which were significantly elevated in MUC1 H tumors. In this scenario it has been shown that cancer cells with defective mitochondria use glutaminedependent reductive carboxylation rather than oxidative metabolism as the major pathway of citrate formation [37].

Changes in Tricarboxylic Acid (TCA) Cycle Intermediates and Related Metabolites
In MUC1 H tumors, citrate and succinate were markedly increased when compared to MUC1 L ccRCC and normal tissue, whereas fumarate and malate were significantly reduced in cancer tissue ( Figure 6).  The altered levels of TCA cycle intermediates, in association with a reduced reactive oxygen species (ROS) production ( Figure 7A) and a decreased membrane potential (Figure 7B), suggests that the rate of oxidative phosphorylation is impaired in MUC1 H tumors, which is consistent with NDUFA4L2 overexpression observed MUC1 H samples ( Figure  7C,D) [17]. The observed changes in the levels of Krebs cycle metabolites are suggestive of a metabolic signature that may be consistent with an increased glutaminolysis. In fact, nearly all the amino acids were significantly reduced in tumor tissue, with the exception of glutamine, glutamate and cysteine, which were significantly elevated in MUC1 H tumors. In this scenario it has been shown that cancer cells with defective mitochondria use glutamine-dependent reductive carboxylation rather than oxidative metabolism as the major pathway of citrate formation [37].

Alterations in Redox Homeostasis
Cancer cells are exposed to high amounts of oxidative stress as the result of increased metabolic activity. In MUC1 H tumor samples, the accumulation of several powerful antioxidants, including alpha-tocopherol (p < 0.0001), beta-tocopherol (p < 0.0001), gammatocopherol (p < 0.0001), ascorbate (vitamin C) (p = 0.0005) and ergothioneine (p < 0.0001), was observed. The levels of the powerful cellular antioxidant glutathione (p = 0.01) were increased as well. Increased levels of glutathione precursors and maintenance of cysteine levels despite significant reductions in nearly all the other amino acids (except glutamine and glutamate) are suggestive of an increased glutathione synthesis in an attempt to maintain redox homeostasis ( Figure 8). Conversely, the levels of the glutathione degradation product 5-oxoproline were significantly reduced in MUC1 H cancer tissues (p < 0.001). Moreover, the polyamines putrescine and spermidine, which bind to and stabilize DNA to promote cellular proliferation, were significantly elevated in MUC1 H versus MUC1 L tumor tissues ( Figure 8).

Changes in Amino Acids and Lipid Metabolism
In addition to the alteration in the biochemical pathways described above, modifications in amino acid metabolism were prevalent in MUC1 H cancer tissue, as evidenced by significant reductions in nearly all the free amino acids measured with the exception of cysteine, glutamate and glutamine ( Figure 9). Taken together, these changes are indicative of a dynamic remodeling that includes protein degradation and the utilization of amino acids for different metabolic activities in cancer cells. Moreover, kynurenine, a metabolite of the amino acid tryptophan [38], was elevated in MUC1 H ccRCC compared to MUC1L tumors (p = 0.0004; Figure 9).
Alteration of lipid metabolism represents another fundamental feature of ccRCC [15,39]. In a recent study [16], we showed that ccRCC was characterized by increased levels of very long-chain and polyunsaturated fatty acids and that this accumulation was sustained by overexpression of stearoyl-CoA desaturase and fatty acid elongase 2 and 5. Alterations in nearly all the free fatty acids measured were observed, and an accumulation of many essential and long-chain fatty acids was found in MUC1 H samples ( Figure 10). These changes may be associated with an increased fatty acid uptake and/or synthesis for subsequent membrane biosynthesis and cell signaling to promote cellular growth and proliferation. Additional information on metabolites evaluated in each lipid class has been provided in the Supplementary Table S2.

Changes in Amino Acids and Lipid Metabolism
In addition to the alteration in the biochemical pathways described above, modifications in amino acid metabolism were prevalent in MUC1 H cancer tissue, as evidenced by significant reductions in nearly all the free amino acids measured with the exception of cysteine, glutamate and glutamine ( Figure 9). Taken together, these changes are indicative of a dynamic remodeling that includes protein degradation and the utilization of amino acids for different metabolic activities in cancer cells. Moreover, kynurenine, a metabolite of the amino acid tryptophan [38], was elevated in MUC1 H ccRCC compared to MUC1L tumors (p = 0.0004; Figure 9).  Alteration of lipid metabolism represents another fundamental feature of ccRCC [15,39]. In a recent study [16], we showed that ccRCC was characterized by increased levels of very long-chain and polyunsaturated fatty acids and that this accumulation was sustained by overexpression of stearoyl-CoA desaturase and fatty acid elongase 2 and 5. Alterations in nearly all the free fatty acids measured were observed, and an accumulation of many essential and long-chain fatty acids was found in MUC1 H samples ( Figure 10). These changes may be associated with an increased fatty acid uptake and/or synthesis for subsequent membrane biosynthesis and cell signaling to promote cellular growth and proliferation. Additional information on metabolites evaluated in each lipid class has been provided in the Supplementary Table S2.

Integrated Metabolomics/Transcriptomic Signature
To compare the relative changes in gene expression and metabolite abundance in MUC1 H versus MUC1 L ccRCC, we integrated the metabolomics data with gene expression data from the Jones cohort (GSE15641) including 23 normal kidney samples and 32 ccRCC samples, stratifying the patients according to MUC1 expression. The combined analysis identified 16 significantly enriched biochemical pathways (p < 0.05), including those of glycolysis, PPP, glutaminolysis and unsaturated fatty acid metabolism (Supplementary Figure S3).

MUC1 Expression Is Associated with Increased Cancer Cell Proliferation and Migration and Has a Role in Chemotherapy Resistance
MUC1 has been implicated in various aspects of tumor biology including cell proliferation, migration and cancer drug resistance [41][42][43]. Therefore, to study the role of MUC1 in these processes, in vitro assays were performed.
The scratch wound-healing assay showed that primary MUC1 H ccRCC cells treated with siMUC1 had a decreased migratory ability compared with MUC1 L cancer cells ( Figure 12).

MUC1 Expression Is Associated with Increased Cancer Cell Proliferation and Migration and Has a Role in Chemotherapy Resistance
MUC1 has been implicated in various aspects of tumor biology including cell proliferation, migration and cancer drug resistance [41][42][43]. Therefore, to study the role of MUC1 in these processes, in vitro assays were performed.
The scratch wound-healing assay showed that primary MUC1 H ccRCC cells treated with siMUC1 had a decreased migratory ability compared with MUC1 L cancer cells (Figure 12).   In addition, MUC1-silenced renal cancer cells proliferated at a slower rate than nonsilenced cancer cells. After cisplatin treatment, the death rate of tumor cells treated with siMUC1 was significantly greater than that of untreated cancer cells (p = 0.008, Figure 13). The MTT assay confirmed these findings, demonstrating a decreased cell viability when tumor cells were pre-treated with siMUC1 before cisplatin incubation ( Figure 13). In addition, MUC1-silenced renal cancer cells proliferated at a slower rate than nonsilenced cancer cells. After cisplatin treatment, the death rate of tumor cells treated with siMUC1 was significantly greater than that of untreated cancer cells (p = 0.008, Figure 13). The MTT assay confirmed these findings, demonstrating a decreased cell viability when tumor cells were pre-treated with siMUC1 before cisplatin incubation ( Figure 13).

MUC1 Soluble Form (Serum CA15-3) Is a Risk Factor for ccRCC Progression and Mortality
MUC1 soluble form, also known as CA15-3, is widely used as a serum marker for breast cancer [44]. Spearman's test showed a positive correlation between MUC1 tissue expression and CA15-3 serum levels (rs = 0.63; p = 0.001). To evaluate its diagnostic and prognostic role in ccRCC, serum CA 15-3 (normal range: 0-25 U/mL) was preoperatively measured in a cohort of 428 consecutive patients who underwent radical or partial nephrectomy for ccRCC at our institution and in 140 healthy adult volunteers with no evi-

MUC1 Soluble Form (Serum CA15-3) Is a Risk Factor for ccRCC Progression and Mortality
MUC1 soluble form, also known as CA15-3, is widely used as a serum marker for breast cancer [44]. Spearman's test showed a positive correlation between MUC1 tissue expression and CA15-3 serum levels (rs = 0.63; p = 0.001). To evaluate its diagnostic and prognostic role in ccRCC, serum CA15-3 (normal range: 0-25 U/mL) was preoperatively measured in a cohort of 428 consecutive patients who underwent radical or partial nephrectomy for ccRCC at our institution and in 140 healthy adult volunteers with no evidence of malignancy. At the time of nephrectomy, CA15-3 serum levels were significantly higher in patients with ccRCC as compared with non-neoplastic patients (p < 0.0001; Figure 14A). One month after surgery, CA15-3 serum levels were significantly reduced (p = 0.0002; Figure 14A). Moreover, at the time of diagnosis significantly higher CA15-3 serum levels were observed in patients with higher nuclear grade (G3-4 vs. G1-2, p < 0.0001; Figure 14B), with lymph node involvement (N1 vs. N0, p = 0.0001; Figure 14C) and with visceral metastases (M1 vs. M0, p < 0.0001; Figure 14D). Kaplan-Meier survival curves for cancer-specific survival (CSS) and progression-free survival (PFS), stratified by CA15-3 serum levels, are shown in Figure 15. Both CSS and PFS were significantly decreased in patients with high serum levels of CA15-3. In addition, subgroup analyses in localized (pT1-2, N0/M0) disease confirmed the prognostic value of CA15-3 in this subset of patients ( Figure 15).
Univariate analysis for the predefined variables showed that the pathological stage, presence of nodal and visceral metastases, nuclear grade, presence of necrosis, tumor size, and high levels of CA15-3 were significantly associated with the risk of death (Table 1) and progression (Table 2). In a multivariate analysis by Cox regression modeling, the   Univariate analysis for the predefined variables showed that the pathological stage, presence of nodal and visceral metastases, nuclear grade, presence of necrosis, tumor size, and high levels of CA15-3 were significantly associated with the risk of death (Table 1) and progression ( Table 2). In a multivariate analysis by Cox regression modeling, the pathological stage, presence of nodal and visceral metastases, nuclear grade and increased CA15-3 serum level were independent adverse prognostic factors for CSS and PFS (Tables 1 and 2).

Discussion
In this study, we delineated the metabolic alterations that characterize MUC1-expressing ccRCC using a combination of chromatography-coupled tandem mass spectrometry, bioinformatics and in vitro cell-based assays. MUC1 is a transmembrane glycoprotein associated with an aggressive phenotype in different epithelial cancers [26,27]. In addition to its role in regulating tumor cell survival, invasion and migration, recent findings indicate that MUC1 induces transcriptional modifications that result in metabolic reprogramming. Shukla et al. showed that MUC1 and HIF-1α crosstalk mediated an anabolic glucose metabolism reprogramming that increased pyrimidine pools in pancreatic cancer [43]. The increased nucleotide biosynthesis led to accumulation of deoxycytidine and reduced the effectiveness of gemcitabine by competitive inhibition.
Moreover, Gunda et al. demonstrated that MUC1-associated metabolic alterations had a role in reducing radiation-induced cytotoxicity and DNA damage in pancreatic cancer, by increasing PPP and nucleotide biosynthesis [45].
MUC1 is also aberrantly expressed in triple-negative breast cancer, where it contributes to epigenetic reprogramming, chromatin remodeling and chemoresistance [46]. In addition, MUC1 regulates glutamine metabolism, enhancing its uptake and carbon flux to glutaminolysis, amino acid metabolism and the TCA cycle [47].
In ccRCC, MUC1 expression is induced by the HIF/hypoxia pathway and promotes the migration and invasion of cancer cells [34]; however, the metabolic reprogramming associated with its overexpression is poorly understood.
In our study, the analysis of metabolic alterations and integration with transcriptomic data showed that MUC1 H ccRCC is characterized by metabolic alterations involving glucose and the lipid metabolism pathway. In particular, we found that in MUC1 H tumors, high glucose levels were accompanied by an increase in the upstream glycolytic intermediates and a decrease in the downstream metabolites. These findings in association with high expression levels of G6PDH, TKT and PPP intermediates suggest an increased rerouting of the sugar metabolism toward this pathway, with the aim of promoting both anabolic reactions and redox homeostasis, especially in MUC1 H ccRCC.
Alteration of lipid metabolism represents another important hallmark of ccRCC. In addition to their role in membrane formation (phospholipids, glycolipids and sterols), lipids represent a fundamental energetic substrate and produce second messengers through their phospholipase-dependent hydrolysis. The production of many bioactive second messengers, such as phosphatidic acid, diacylglycerol and arachidonic acid, can trigger the activation of several signaling pathways that promote cancer cell proliferation, survival, migration, invasion and metastatization [48,49].
It has been shown that ccRCC is characterized by a significant accumulation of polyunsaturated fatty acids (PUFAs), in association with overexpression of stearoyl-CoA desaturase-1 (∆-9-desaturase; SCD1) and fatty acid elongase 2 and 5 (ELOVL2 and ELOVL5) [16]. In the present study, we found that these perturbations were more evident in MUC1 H ccRCC compared to MUC1 L tumors and normal tissues. In addition, MUC1 H tumors were characterized by accumulation of several powerful antioxidants, including glutathione, in association with increased expression of NDUFA4L2 and reduced mitochondrial membrane potential and ROS production, as shown by reduced signals of the fluorescence probe TMRE and MitoSOX.
Overexpression of MUC1 provided additional advantages in different biological processes such as cell proliferation, cancer cell migration and resistance to chemotherapy.
RCC is typically a chemoresistant tumor, and although numerous mechanisms of chemoresistance have been described, the overexpression of P-glycoprotein is one of the most significant [50][51][52]. In addition to approaches targeting P-glycoprotein, other tactics to overcome drug resistance include MUC-1 targeting approaches.
In accordance with previous studies performed in different tumors, we found that MUC1 knockdown decreased cell migration and viability and improved cisplatin susceptibility, suggesting that this protein can be involved in de novo chemotherapy resistance in ccRCC.
Finally, we analyzed the prognostic role of CA15-3, a serum marker derived from MUC1 and widely used for breast cancer clinical management [44].
Previous studies suggested that high serum levels of CA15-3 could be associated with an advanced disease in RCC patients. Grankvist et al. [53] observed increased levels of CA15-3 in up to 30% of patients with high-grade and high-stage RCC. In a later study, Briasoulis et al. identified abnormal levels of this serum biomarker in 23% of RCC patients with metastatic disease [54]. In accordance with these results, in a more recent study we showed that CSS and PFS were significantly shorter for clear cell and non-clear cell RCC patients with elevated CA15-3 levels [20].
In the present study, multivariate analysis showed that high levels of CA15-3, together with the presence of nodal and visceral metastases, advanced stage and nuclear grade, were significantly predictive of risk of death. Moreover, this biomarker remained an independent prognosticator of outcome for PFS. Kaplan-Meier curves showed clear differences in CSS and PFS between patients with CA15-3 values below and above the normal values. Moreover, subgroup analyses for localized disease confirmed significant differences in CSS and PFS for CA15-3 values ≤ 25 or >25 U/mL in this clinical setting.
In conclusion, this study suggests that MUC1 acts as a metabolic regulator in ccRCC, inducing a particular metabolic reprogramming. The inhibition of MUC1 expression decreases cell motility and viability and improves cisplatin susceptibility, suggesting that this pathway can regulate de novo chemotherapy resistance in ccRCC. Finally, we showed that CA15-3 could serve as a promising biomarker to identify ccRCC patients with poor prognosis.

Study Population and Tissue Collection
Investigation has been conducted in accordance with the ethical standards and according to the Declaration of Helsinki and according to national and international guidelines and has been approved by the authors' institutional review board. Written informed consent to take part was given by all participants.
Primary renal tumor (n = 36) and non-neoplastic tissues (n = 20) were collected from patients who underwent radical or partial nephrectomy for ccRCC (Supplementary Table S1). All specimens were immediately used to obtain primary cell cultures after frozen-section confirmation of the diagnosis. Two pathologists confirmed the presence of clear cell RCC in the neoplastic tissues and excluded tumor cells in the healthy specimens. Patients with diabetes mellitus and/or an estimated glomerular filtration rate (MDRD equation) < 60 mL/min/1.73 m 2 were excluded from the study.
For CA15-3 level evaluation, serum CA15-3 (normal range: 0-25 U/mL) was preoperatively measured in a cohort of 428 consecutive patients who underwent radical or partial nephrectomy for ccRCC at our institution and in 140 healthy adult volunteers with no evidence of malignancy. Patients with an estimated glomerular filtration rate (eGFR calculated using MDRD equation) <60 mL/min/1.73 m 2 and with other known causes of elevated levels of this biomarker were excluded from the study (namely, patients with other malignant tumors, breast and liver diseases, endometriosis, sarcoidosis, tuberculosis and systemic lupus erythematosus). Detailed clinical and pathological characteristics of these patients are summarized in Supplementary Table S3. All patients were preoperatively staged by thoraco-abdominal Computed Tomography or Magnetic Resonance Imaging. Tumor staging was reassigned according to the seventh edition of the AJCC-UICC TNM classification. The 2016 World Health Organization classification was used to attribute histological type and nuclear grade.

Data Quality: Instrument and Process Variability
Metabolic profiling of RCC samples was carried out at Metabolon Inc. Instrument variability was determined by calculating the median relative standard deviation (RSD) for the internal standards that were added to each sample prior to injection in the mass spectrometers. Overall process variability was determined by calculating the median RSD for all endogenous metabolites (i.e., non-instrument standards) present in 100% of the Client Matrix samples, which are technical replicates of pooled client samples.

Sample Preparation
All tissue samples were maintained at −80 • C until processed. At the time of analysis, samples were thawed, and extracts prepared according to Metabolon's standard protocol, which is designed to remove protein, dislodge small molecules bound to protein or physically trapped in the precipitated protein matrix and recover a wide range of chemically diverse metabolites.
In particular, the sample preparation process was carried out using the Hamilton Company automated MicroLab STAR ® system. Recovery standards were added prior to the first step in the extraction process for QC purposes. Sample preparation was conducted using a proprietary series of organic and aqueous extractions to remove the protein fraction, while allowing maximum recovery of small molecules. The resulting extract was divided into two fractions: one for analysis by LC and one for analysis by GC. Samples were placed briefly on a TurboVap ® (Zymark, Hopkinton, MA, USA) to remove the organic solvent. Each sample was then frozen and dried under vacuum. Samples were then prepared for the appropriate instrument, either LC/MS or GC/MS. For QA/QC purposes, a number of additional samples were included in each day's analyses. Furthermore, a selection of QC compounds was added to every sample, including those under test. These compounds were carefully chosen so as not to interfere with the measurement of the endogenous compounds. The Supplementary Tables S4 and S5 describe the QC samples and compounds. These QC samples were primarily used to evaluate the control process for each study as well as assisting in data curation.

Liquid Chromatography/Mass Spectrometry (LC/MS, LC/MS)
The LC/MS portion of the platform was based on a Waters ACQUITY UPLC and a Thermo-Finnigan LTQ mass spectrometer, which consisted of an electrospray ionization (ESI) source and linear ion-trap (LIT) mass analyzer. The sample extract was split into two aliquots, dried, then reconstituted in acidic or basic LC-compatible solvents, each of which contained 11 or more injection standards at fixed concentrations. One aliquot was analyzed using acidic positive ion optimized conditions and the other using basic negative ion optimized conditions in two independent injections using separate dedicated columns. Extracts reconstituted in acidic conditions were gradient-eluted using water and methanol, both containing 0.1% formic acid, while the basic extracts, also reconstituted using water/methanol, contained 6.5 mM ammonium bicarbonate. The MS analysis alternated between MS and data-dependent MS2 scans using dynamic exclusion.

Gas Chromatography/Mass Spectrometry (GC/MS)
The samples destined for GC/MS analysis were re-dried under vacuum desiccation for a minimum of 24 h prior to being derivatized under dried nitrogen using bistrimethyl-silyltrifluoroacetamide (BSTFA). The GC column was 5% phenyl and the temperature ramped from 40 to 300 • C in a 16 min period. Samples were analyzed on a Thermo-Finnigan Trace DSQ fast-scanning single-quadrupole mass spectrometer using electron impact ionization. The instrument was tuned and calibrated for mass resolution and mass accuracy on a daily basis. The information output from the raw data files was automatically extracted as discussed below.

Accurate Mass Determination and MS/MS Fragmentation (LC/MS, LC/MS/MS)
The LC/MS portion of the platform was based on a Waters ACQUITY UPLC and a Thermo-Finnigan LTQ-FT mass spectrometer, which had a linear ion-trap (LIT) front end and a Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer backend. For ions with counts exceeding 2 million, an accurate mass measurement could be performed. Accurate mass measurements could be made on the parent ion as well as fragments. The typical mass error was less than 5 ppm. Characterizing ions with less than two million counts is a more laborious process. Fragmentation spectra (MS/MS) were typically generated in a data-dependent manner, but if necessary, targeted MS/MS could be employed, as in the case of lower level signals.

Compound Identification
Compounds were identified by comparison to library entries of purified standards or recurrent unknown entities. The identification of known chemical entities was based on comparison to metabolomic library entries of purified standards.

Bioinformatics and Statistical Analyses
Significance tests were performed with MedCalc 9.2.0.1 (MedCalc software, Mariakerke, Belgium) and "R" (http://cran.r-project.org (accessed on 1 March 2022)). Global biochemical profiles were determined in human kidney tissue/tumor samples and compared across groups stratified by renal tissue pathology status. An estimate of the false discovery rate (q-value) was calculated to consider the multiple comparisons that normally occur in metabolomic-based studies. Comparisons of metabolite median values between different groups were evaluated by Mann-Whitney U test. Spearman's correlation was applied to evaluate the association between MUC1 tissue expression and CA15-3 serum levels. In the cancer-specific survival (CSS) analysis, patients still alive or lost to follow-up were censored, as well as patients who died of RCC-unrelated causes. Progression-free survival (PFS) was calculated from the date of surgery to the date of disease recurrence. Disease progression was assessed radiographically (using CT scan or MRI) with a surveillance schedule based on EAU guidelines. Estimates of CSS and PFS were calculated according to the Kaplan-Meier method and compared with the log-rank test. Univariate and multivariate analyses were performed using the Cox proportional hazards regression model to identify the most significant variables for predicting CSS and PFS. A backward selection procedure was performed with removal criterion p > 0.10 based on likelihood ratio tests. A p-value of < 0.05 was considered statistically significant.

Integration of Metabolomic and Transcriptomic Data
To identify genes associated with high versus low MUC1 expression, gene expression data from the Jones cohort (GSE15641) including 23 normal kidney samples and 32 ccRCC samples were used. We first stratified the patients by MUC1 expression, and next we generated a rank file for each expressed gene by the log2 fold change in high MUC1 samples over low MUC1 samples. Next, we ran a gene set enrichment analysis (GSEA) to determine which pathways were statistically enriched across the renal cancer datasets [40]. The normalized enrichment score (NES) was used to evaluate the extent and direction of enrichment of each pathway.
MetaboAnalyst 4.0 (https://www.metaboanalyst.ca (accessed on 1 March 2022)) was used for the metabolite set enrichment (MSEA) and to integrate the data from the transcriptomics and metabolomics experiments [35]. The main function of this module is to pinpoint the pathways involved in the underlying biological processes by combining the evidence based on alterations in both gene expression and metabolite concentrations. In addition, a topology analysis was used to evaluate the relative importance of the gene/compounds based on their relative locations within a pathway. Over-representation analysis of the pathways was based on the hypergeometric test.

Real-Time PCR
Total RNA of normal and tumor tissues was reverse transcribed with the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems Foster City, CA, USA), following the manufacturer's instructions. Quantitative real-time polymerase chain reactions (PCR) were performed using the iQTM SYBR Green Supermix buffer (6mM-MgCl2, dNTPs, iTaq DNA polymerase, SYBR Green I, fluorescein and stabilizers) (BIO-RAD Laboratories, Hercules, CA, USA). The primers used for rea-time PCR are listed in Supplementary Table S6.
Quantification of the mRNA levels was performed on a MiniOpticon Real-Time PCR detection system (BIO-RAD Laboratories). In the PCR reactions, the following protocol was used: polymerase activation at 95 • C for 3 min, followed by 45 cycles at 95 • C for 10 s, 60 • C for 30 s. Melting curves were generated through 60 additional cycles (65 • C for 5 s with an increment of 0.5 • C/cycle). Gene expression results were obtained as mean Ct (threshold cycle) values of triplicate samples. Expression was determined using the 2-∆∆Ct method. Expression values were normalized to β-Actin.
For the immunoassay procedures, 25 µL of each dilute lysate sample in assay buffer (5 µg total protein/well) and HeLa cell lysate (positive control) were added into wells in duplicate, according to the manufacturer's instructions. To each well, 25 µL of the Mixed Beads were added and the plate was incubated for 2 h at room temperature. Human glycolysis pathway detection biotinylated antibodies were added for 1 h; each captured a specific bead. After that, the reaction mixture was incubated for 30 min with Streptavidin-PE conjugate to complete the reaction on the surface of each microsphere. Finally, the MILLIPLEX ® MAP was analyzed by Luminex xMAP ® technology. The immunoassay on the surface of each fluorescent-coded magnetic bead, MagPlex-C microsphere, was identified and quantified based on fluorescent signals. The median fluorescence intensity (MFI) was read with the Luminex 200TM instrument and measured with xPONENT ® software.

Determination of Glucose-6-Phosphate Dehydrogenase Activity
Glucose-6-phosphate dehydrogenase (G6PDH) activity was assessed in tissue lysates using the Glucose-6-Phosphate Dehydrogenase Assay Kit (Abcam, Cambridge, UK). The assay was performed by adding 50 µL of dilute samples (1:2 with assay buffer) into duplicate wells of 96-well plates. Then, 50 µL of reaction mix (G6PDH substrate and G6PDH developer) was added to each well containing samples or positive control. At two (T1) and thirty (T2) minutes after addition of the reagents and incubation at 37 • C, the plate was read in an ELISA microplate reader at 450 nm. G6PDH activity was spectrophotometrically measured by monitoring the production of NADH, and calculated as a change in optical density (∆A450 nm) at T1 and T2 per sample volume added into the reaction well, expressed in nmol/mg.

Primary Cell Cultures from Renal Tissues
Tumor (ccRCC) and normal (N) kidney tissue specimens were immediately placed in a Petri dish with phosphate-buffered saline (PBS) 1× and cut into small pieces of about 1 mm 3 . Each small piece of tissue was placed on the surface of the Petri dish with Dulbecco's Modified Eagle Medium (DMEM, Invitrogen, Life Technologies, Monza, Italy) supplemented with 10% fetal bovine serum (FBS, Sigma-Aldrich, Milan, Italy) and 1% penicillin-streptomycin-L-glutamine (Sigma-Aldrich, Milan, Italy). Cell proliferation was obtained around the kidney specimens as previously described [55]. Kidney epithelial tubular and neoplastic cells were isolated with EpCAM (CD326) Ab-conjugated magnetic microbeads (Miltenyi Biotec, Bergisch Gladbach, Germany) under the effect of a magnetic field generated by the Mini MACS Separation Unit (Miltenyi Biotec), as previously described [55]. These cells were then characterized for EpCAM, CA IX and MUC1 by immunocytochemistry. Primary ccRCC cells were used for proliferation studies at the second passage.

Small Interfering (siRNA) Transfection
Isolated normal and tumor renal cells were cultured at 2 × 105 cells per well in a 12-well plate with Keratinocyte Serum-Free Medium (KSFM), supplemented with 5 ng/mL recombinant epidermal growth factor (rEGF), 50 µg/mL bovine pituitary extract (BPE) (Gibco, Life Technologies, Monza, Italy) and 30 ng/mL cholera toxin (Sigma-Aldrich, Milan, Italy). The transfection of siRNA was carried out using Lipofectamine 3000 (Life Technologies, Monza, Italy) in accordance with the manufacturer's procedure. For each transfection, 50 nM of small interfering RNA targeting MUC1 (siMUC1) (Qiagen, Hilden, Germany) was used. In transfection experiments, a mock transfection control was performed by putting cells through the transfection procedure without adding siRNA. The validated non-silencing siRNA sequence AllStars Negative Control siRNA (50 nM, Qiagen, Hilden, Germany) was used as negative control. Each transfection experiment was performed in triplicate. After transfection, normal and tumor renal cells were incubated for 12 h and 24 h at 37 • C in 5% CO 2 and used for total RNA extraction, wound-healing and cell viability assays.

Immunohistochemistry
Immunohistochemical evaluation of MUC1, GLUT1 and NDUFA4L2 protein expression was conducted on paraffin-embedded tissue sections. These sections (3 mm) were deparaffinized and rehydrated through xylenes and graded alcohol series. Slides were subjected to specific epitope demasking by microwave treatment at 700 W in citrate buffer (0.01 M, pH 6.0). After antigen retrieval, the tissue samples were incubated for 10 min with 3% H 2 O 2 to block endogenous peroxidase activity. Sections were blocked with Protein Block Serum-Free (Dako) at room temperature for 10 min and then incubated with anti-MUC1 antibody (1:200, Novus Biologicals, Littleton, CO, USA), anti-NDUFA4L2 antibody (1:100 dilution, Proteintech, Chicago, IL, USA) and anti-GLUT1 antibody (1:200, Novus Biologicals, Littleton, CO, USA) at 4 • C overnight. Binding of the secondary biotinylated antibody was detected using the Dako Real EnVision Detection System, Peroxidase/DAB kit (Dako, Agilent, Santa Clara, CA, USA), according to the manufacturer's instructions. Sections were counterstained with Mayer's hematoxylin (blue) and mounted with glycerol (Dako Cytomation). Negative controls were obtained by incubating serial sections with the blocking solution and then omitting the primary antibodies. Digital images were obtained using the Aperio ScanScope CS2 device (Aperio Technologies, Vista, CA, USA) and further analyses of the scanned images were performed with the ImageScope V12.1.0.5029 (Aperio). Specific staining was quantified by applying the Positive Pixel Count v9_v10.0.0.1805 algorithm (Aperio) and expressed as percentage of positive pixels in the analyzed area.

Wound-Healing Assay
The 2 × 10 5 normal and primary tumor renal cells were seeded onto a six-well plate with 2 mL of Keratinocyte Serum-Free Medium (KSFM) to create a confluent monolayer, supplemented with 5 ng/mL recombinant epidermal growth factor (rEGF), 50 µg/mL bovine pituitary extract (BPE) (Gibco) and 30 ng/mL cholera toxin (Sigma-Aldrich, Milan, Italy). Cells were incubated overnight at 37 • C, 5% CO 2 , and then exposed to 50 nM of siMUC1 (Quiagen) or incubated in medium without siRNA for 72 h. After 24 h from the treatment, a wound was manually created by scraping the cell monolayer with a P200 pipette tip. A reference mark was created on the dish and a time 0 image was acquired. After 12 and 24 h, additional images were taken in the matched region, and the woundhealing area was quantified with ImageJ software (http://rsbweb.nih.gov/ij/ (accessed on 1 March 2022)). Each experimental condition was performed in triplicate.
For the dye exclusion test and MTT assay, normal and primary tumor cells were seeded at a density of 1.5 × 105 and 1.5 × 104 cells in 6-well and 96-well plates, respectively, (Sigma Aldrich, Milan, Italy) and incubated overnight at 37 • C, 5% CO 2 , 1 day before exposure to siMUC1 or to medium alone. The cells were exposed to 50 nM of siMUC1 for 24 h and then treated with 10 µM cis-Diamminedichloroplatinum(II) (cisplatin) for 1 h and 2 h. After several washes to remove cisplatin, the cells were again incubated in medium with 50 nM of siMUC1 for 48 h. Lastly, the cells were trypsinized and the viable cells were counted using trypan blue or were exposed to MTT. Each experimental condition was performed in triplicate.

Immunofluorescence Microscopy
Primary tumor cells were seeded at a density of 2 × 10 5 on glass coverslips and left to adhere overnight at 37 • C in 5% CO 2 . These preparations were stained for NDUFA4L2 (Proteintech, Chicago, IL, USA) and MitoSOX red (Molecular Probes, Life Technologies, Monza, Italy). The expression and localization of proteins was evaluated by indirect im-munofluorescence and confocal microscopy analysis. In particular, the cells were fixed using ice-cold 4% paraformaldehyde for 10 min at room temperature. Then, the preparations were blocked with 1% BSA in PBS for 1 h at room temperature and incubated overnight at 4 • C with a primary antibody against NDUFA4L2 (1:25 in blocking), followed by incubation for 1 h at 37 • C with the secondary antibody goat anti-rabbit IgG FITC (1:200; Novus Biologicals). To study mitochondrial superoxide generation, live cells were stained with 5 µM MitoSOX red for 10 min at 37 • C. All preparations were counterstained with TO-PRO-3 (Molecular Probes). Negative controls were performed by omitting the primary antibodies. Specific fluorescence was acquired by a Leica TCS SP2 (Leica, Wetzlar, Germany) confocal laser-scanning microscope using a ×63 objective lens.

Mitochondrial Membrane Potential
Mitochondrial membrane potential was determined by incubating primary tumor cells with the fluorescent dye tetramethylrhodamine ethyl ester (TMRE; Sigma-Aldrich, Milan, Italy). Cells were incubated with 2.5 µM TMRE for 1 h at 37 • C and subsequently analyzed by immunofluorescence.

Informed Consent Statement:
Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement:
The datasets generated and/or analyzed during the current study are available in the GEO repository: accession number GSE15641.