Quantitative proteomic analysis reveals potential diagnostic markers and pathways involved in pathogenesis of renal cell carcinoma

There are no serum biomarkers for the accurate diagnosis of clear cell renal cell carcinoma (ccRCC). Diagnosis and decision of nephrectomy rely on imaging which is not always accurate. Non-invasive diagnostic biomarkers are urgently required. In this study, we preformed quantitative proteomics analysis on a total of 199 patients including 30 matched pairs of normal kidney and ccRCC using isobaric tags for relative and absolute quantitation (iTRAQ) labeling and LC-MS/MS analysis to identify differentially expressed proteins. We found 55 proteins significantly dysregulated in ccRCC compared to normal kidney tissue. 54 were previously reported to play a role in carcinogenesis, and 39 are secreted proteins. Dysregulation of alpha-enolase (ENO1), L-lactate dehydrogenase A chain (LDHA), heat shock protein beta-1 (HSPB1/Hsp27), and 10 kDa heat shock protein, mitochondrial (HSPE1) was confirmed in two independent sets of patients by western blot and immunohistochemistry. Pathway analysis, validated by PCR, showed glucose metabolism is altered in ccRCC compared to normal kidney tissue. In addition, we examined the utility of Hsp27 as biomarker in serum and urine. In ccRCC patients, Hsp27 was elevated in the urine and serum and high serum Hsp27 was associated with high grade (Grade 3–4) tumors. These data together identify potential diagnostic biomarkers for ccRCC and shed new light on the molecular mechanisms that are dysregulated and contribute to the pathogenesis of ccRCC. Hsp27 is a promising diagnostic marker for ccRCC although further large-scale studies are required. Also, molecular profiling may help pave the road to the discovery of new therapies.


INTRODUCTION
Renal cell carcinoma (RCC) is the most common neoplasm in the adult kidney, with an increasing incidence over the past 20 years [1]. About 80% of RCCs are of the clear-cell type (ccRCC). Early diagnosis of RCC is associated with favorable prognosis (5-year survival rate ~ 85%). Unfortunately, RCC is often asymptotic, with about 30% of patients diagnosed at the metastatic stage when the prospects for cure are dismal (5-year survival rate ~9%) [2]. Traditional manifestations of pain, mass and hematuria are ineffective for early diagnosis [3]. The diagnosis of RCC, and the subsequent resection of the kidney are based on imaging findings, which are not always accurate. There are currently no serum biomarkers available to confirm the identity of renal masses, whether benign or malignant. A non-invasive test in serum or urine will have a significant impact on patient management.
Proteomics combined with mass spectrometry (MS) offers great promise for unveiling the complex molecular events of tumorigenesis and identification of cancer biomarkers. Tissue proteomics is a promising alternative strategy for the discovery and identification of tumor markers. A major advantage of tissues proteomics is that relevant proteins are much more abundant. Additionally, the link between proteins that are differentially expressed in diseased tissue with the disease itself is stronger than serum.
In this study, we performed quantitative proteomic analysis using isobaric tags for relative and absolute quantitation (iTRAQ) labeling and LC-MS to identify proteins that are dysregulated in ccRCC compared to normal kidney. We identified a number of proteins that can distinguish between tumor and normal tissues. We verified the dysregulated expression of four of the most interesting proteins on two independent sets of tissue using both western blot and immunohistochemistry. Also, these proteins showed to have the ability to distinguish between normal and cancer tissue with accuracy. Identified secreted proteins may serve as diagnostic makers. Furthermore, using pathway analysis and protein-protein analysis, we elucidated the potential involvement of these proteins in RCC pathogenesis and identified key pathways that are dysregulated in RCC. We also examined the use of Hsp27 as a biomarker in the serum and urine and its prognostic utility.

Identification of dysregulated proteins between ccRCC and normal tissues
A schematic of the work flow is shown in Supplementary Figure 1. A total of 40 samples were analyzed in the discovery phase. Briefly, equal amounts of protein from each tissue type were digested with trypsin, labeled with iTRAQ and combined. Samples were then separated by off-line SCX liquid chromatography and analyzed by RP LC-MS. We identified a total of 1591 nonredundant proteins with local false discovery rate (FDR) < 5% (Supplementary Table 5); 345 of these proteins were reliably quantified (Supplementary Table 6). Fifty five proteins fulfilled our criteria for dysregulation between ccRCC and normal (see materials section, Supplementary Table 7): 15 were upregulated (iTRAQ ratios of ≥1.5) and 40 were downregulated (iTRAQ ratios of ≤0.67). Table  1 shows a heat map of the 55 significantly dysregulated proteins.
Unsupervised clustering was performed on the 55 significantly dysregulated proteins. The samples clustered into two main groups: one that contained 9 of 10 ccRCC samples; and second included all normal samples and one cancer sample ( Figure 1). The difference in expressions of dysregulated proteins between ccRCC and normal Clustering analysis was performed based on 345 proteins for which quantitative information was available. The samples clustered into two main groups: one that contained the 9 of 10 ccRCC samples (C1-C5, C7-C10); and second contained all normal samples (N1-N10) plus one cancer sample; C6. The difference in expressions of dysregulated proteins between ccRCC and normal kidney samples was statistically significant (p<0.001). kidney was statistically significant (p<0.001). To verify our analysis, we performed control clustering based on 345 quantified proteins excluding the 55 significantly dysregulated proteins. As expected, cancer and normal samples didn't cluster into distinct groups (p>0.5; Supplementary Figure 2).

Identification of dysregulated proteins that can serve as serum biomarkers
Of 55 dysregulated proteins 39 (70.9%) satisfied one or more of the four criteria to be classified as "secretory" proteins (Supplementary Table 7); 5 proteins (9.1%) were found to be extracellular; 10 (18.2%) were membrane-bound; 8 (14.5%) could be released from cells

Figure 2: Verification of protein dysregulation in ccRCC by Western blot and immunohistochemical analyses. A:
Representative blots showing the expression of proteins in normal kidney tissues (N1, N2) and ccRCC (C1, C2). For AHNAK, ENO1, and HSP27, expression was significantly higher and for HSPE1 it was significantly lower in cancer compared to normal kidney tissue. β-actin was used as a loading control. B: Graphical representation of the average fold change in expression of the proteins between ten ccRCCs and matched normal specimens as determined by densitometry (C/N). Expressions of protein in normal samples were normalized. C-J: Representative photomicrographs showing differential expression of ENO1, HSPB1, HSPE1 and LDHA in ccRCC compared to normal kidney tissue by immunohistochemistry (Original magnification × 200). proteins are enzymes which catalyze the reactions of glycolysis, citric acid cycle, metabolism and catabolism of Acetyl-CoA. ETFB serves as a specific electron acceptor for several dehydrogenases, including five acytyl-CoA dehydrogenases (AD), glutaryl-CoA and sarcosine dehydrogenase; and HSP27 forms a complex with G6PDH that increased its activity. Upregulated proteins are shown in red and downregulated ones in green. B: Visualization of protein-protein interactions for dysregulated proteins in ccRCC using STRING analysis. Dysregulated proteins were used as input for STRING and are represented as spheres of distinct colors. Blue lines represent interactions between proteins and the thickness of the lines display the level of confidence associated with each interaction. Our dysregulated proteins formed one main cluster. Within this cluster, we identified three mini-clusters which contained proteins related to cell metabolism. C: PCR validation of the down-regulation of genes involved in the TCA cycle. We found that 18 of the 28 (64%) genes examined had significant decreased expression in ccRCC tissues when compared to normal kidney tissue. SUCLG1, p=0.003; IDH3A, p=0.026; PDHB, p<0.001; SUCLG2, p=0.006; SDHD, p=0.012; DLD, p=0.020; SDHB, p=0.007; PCK2, p=0.037; MDH2, p=0.043; MDH1, p<0.001; IDH2, p=0.030; DLST, p=0.037; FH, p=0.006; SUCLA2, p=0.001; DLAT, p<0.001; SDHC, p<0.001; IDH3B, p=0.006; IDH3G, p<0.001. D: Bar graph showing differential expression of Hsp27 in urine. We assayed the expression of Hsp27 by ELISA in 21 pre-operative RCC patients and 9 individuals with no malignancy. The average expression in RCC patients was significantly higher than those with no malignancy (1.822ng/ mL vs. 0.3365ng/mL, respectively, p<0.05). www.impactjournals.com/oncotarget via the exosome pathway; 1 (1.8%) is predicted to be classically secreted according to the SignalP analysis; and SecretomeP analyses predicted that 32 (58.2%) proteins are likely to be non-classically secreted.

Validation of the protein dysregulation
Of 39 "secretory" proteins we selected the 5 most promising for further verification as potential ccRCC biomarkers. The upregulated proteins ENO1 and Hsp27 were previously documented to be involved in tumorigenesis. AHNAK was a promising candidate with no record of its involvement in other malignancies. The downregulated HSPE1 was also documented to be involved in the pathogenesis of different malignancies.
In the first step, we verified differential expressions of these four proteins by Western blot analysis in samples of the discovery cohort ( Figure 2A). The expressions of AHNAK, ENO1, and Hsp27 were found to be significantly elevated [1.68 fold (p<0.002), 1.62 fold (p<0.01), and 1.47 fold (p<0.01), respectively] in ccRCC compared to matched normal tissues. HSPE1 was significantly downregulated in ccRCC [0.47 fold (p<0.002)], in agreement with our MS findings ( Figure 2B). Overall, these results confirm our MS analysis.
We additionally verified the dysregulation of ENO1, HSPB1, HSPE1, and LDHA in an independent cohort of patients by IHC using TMAs consisting of 85 cases of ccRCC and matched normal kidney tissue from the same patient ( Figure 2C-J). We found that ENO1 immunoexpression was increased in 56 (70%) cases of ccRCC when compared to normal kidney tissue from the same patient (Supplementary Table 8), which was in agreement with both the MS and Western blot results.
When we examined the expression of HSPB1/ Hsp27, 53 (69%) cases showed increased expression in ccRCC compared to normal matched tissues. 11 (14%) cases showed no significant change in expression between ccRCC and normal, and 13 (17%) cases showed decreased expression in ccRCC. The downregulation of HSPE1 was also seen immunohistochemically as 71 (91%) cases of ccRCC showed decreased expression compared to matched normal. Only 6 (8%) cases showed similar expression levels between normal and cancer. Comparable findings were seen for LDHA where there was increased expression in 76 (96%) ccRCC tissues examined. Only one case showed decreased expression in cancer and two cases (3%) showed no change in expression. The results of the IHC analyses are summarized in Supplementary Table  8.
We also performed protein-protein interaction analysis and found significant protein-protein interactions among our dysregulated proteins. These proteins formed one main cluster ( Figure 3B), suggesting the presence of many interactions within this protein set. We also identified three smaller clusters within the main cluster. All three mini-clusters were related to glucose metabolism. In particular, we found the largest of the mini-clusters contained proteins that are involved in the TCA cycle.

Experimental validation of metabolic dysregulation in ccRCC
We confirmed downregulation of genes involved in the citric acid cycle by PCR ( Figure 3C). Using PCR arrays, we found 18 of the 28 genes (64%) involved in the citric acid cycle showed significant decreased expression in ccRCC tissues when compared to normal kidney tissue. The most downregulated genes were succinate-CoA ligase (SUCLG1, -13.0 times fold change, p=0.003) and isocitrate dehydrogenase 3 (NAD+) alpha (IDH3A, -7.31 times fold change, p=0.026). Another gene that showed a high degree of downregulation in ccRCC tissues was pyruvate dehydrogenase (PHDB, , -7.0 times fold change, p<0.001).

Hsp27 is a potential biomarker in RCC
We performed a preliminary analysis to evaluate the potential utility of Hsp27 and ENO1 as serum and urine markers in ccRCC. We found significant differences in the levels of urinary Hsp27 between patients with ccRCC compared to normal individuals (average 1.822ng/mL vs. 0.3365ng/mL, respectively (p<0.05; Figure 3D). These data suggest Hsp27 may be a useful diagnostic marker for RCC patients.
Serum analysis showed that average Hsp27 expression in patients with ccRCC was 2.334ng/mL compared to 1.233ng/mL in non-cancerous individuals. Using a cut-off of 2.529ng/mL (average Hsp27 expression of non-cancerous patient's ± 2 SD) patients were classified as having high or low Hsp27 expression. Patients with high Hsp27 were more likely to have ccRCC (p=0.008; Table 3). Also, higher tumor grades (grade 3-4) were associated with higher Hsp27 expression (p=0.013). We also examined the expression levels of ENO1 in serum and urine and found that expression was higher in ccRCC compared to normal but the differences did not reach statistical significance. The average expression of ENO1 in the serum of RCC patients and patients with no malignancy was 0.0061ng/mL and 0.0057ng/mL, respectively. Urine samples had 0.007187 ng/mL vs. 0.006842 ng/mL, respectively (data not shown).

DISCUSSION
A significant challenge in ccRCC is the lack of tools that can distinguish between benign and malignant kidney masses. Imaging has limited accuracy and cannot be used reliably to confirm the nature of the lesion. This results in some patients having their kidneys unnecessarily removed for benign lesions. Biopsy is an invasive procedure with limited success and a number of side effects. A non-invasive urine or blood-based test will represent a revolutionary step in RCC diagnosis. We performed comprehensive quantitative proteomics analysis and identified promising biomarkers that can be further investigated as non-invasive tests for the accurate diagnosis of kidney cancer. We identified 55 proteins that are significantly dysregulated in ccRCC relative to normal kidney tissues (Table 1). Hierarchical clustering showed the ability of these proteins to distinguish between RCC and normal kidney tissues (Figure 1). Our results are in agreement with previous reports [9;18;19;22-26].
There are a number of unique aspects in our study. The two-dimensional liquid chromatographic approach employed is considered superior to 2D-PAGE. Proteins with molecular masses higher than 150 kDa and lower than 15 kDa as well as proteins with isoelectric points outside the range of pH 3 -10 cannot be identified by 2D-gels and hydrophobic membrane proteins are underrepresented on 2D gels. Advanced features of our mass spectrometric analyses include simultaneous profiling of multiple admixed specimens which helps to eliminate artifactual differences due to differential protein losses during purification or separation. In addition, the use of four isotopic labels allows quantitative comparison among four different tissue samples.
Of 55 dysregulated proteins, 39 (70.9%) were classified as "secretory" proteins and thus have the potential to serve as serum/urine diagnostic biomarker for ccRCC. Interestingly, ENO1, LDHA, and Hsp27, AHNAK have been reported in other malignancies to be upregulated in serum of cancer patients, compared to healthy controls: serum ENO1 was elevated in patients with small-cell lung carcinoma; and high serum LDH was linked with significantly poor survival of colorectal cancer [27]; high level of Hsp27 was present in serum from patients with breast cancer [28]; and serum AHNAK level was elevated in patients with ovarian cancer [29]. That suggests that upregulation of these proteins can be detected in serum of ccRCC patients as well. Our pilot study shows a potential of Hsp27 as a serum/urine biomarker. This needs to be validated in a larger cohort of patients.
Our results also strongly suggest that ccRCCs have a protein "signature" that is required for carcinogenesis. All 55 dysregulated proteins had previously been reported to be involved in tumorigenesis processes. Interestingly, the majority of these proteins are involved in "metabolic" processes ( Figure 3A and Supplementary Table 7). Dysregulated cellular metabolism is adapted by cancer cells to meet the requirements of rapid cell proliferation, growth, negative regulation of apoptosis, survival under hypoxia etc.
Our findings are comparable to earlier reports on functional analyses in ccRCC [30;31]. The link between carbohydrate metabolism and RCC is not surprising. In contrast to normal proliferating cells, tumor cells have to survive in environments with varying oxygen and nutrient supplies [32]. The increase in lactate dehydrogenase and the activation of the pyruvate kinase pathways indicate active anaerobic glycolysis which is a reflection of the hypoxic conditions known to be an integral component of the pathogenesis of RCC [33;34]. Also the "clear cell" morphology of RCC is known to result from the accumulation of glycogen as a result of disturbed carbohydrate metabolism.
Global analysis through molecular profiling and pathway analysis can have significant clinical applications. In addition to identified new therapeutic targets, pathwayderived metabolic products can be also used as diagnostic, predictive and prognostic markers [35].
In conclusion, through quantitative proteomic analysis, we identified differential protein expressions that can distinguish between ccRCC and normal kidney tissues. Most of these proteins are involved in biological pathways pertinent to carcinogenesis. About 70% of these proteins can be potentially shed into serum. If the upregulation of the most promising proteins would be confirmed in serum/urine of ccRCC patients, then that would suggest their potential as noninvasive biomarkers for confirming the diagnosis of RCC, which may greatly improve patient management and increase overall survival.

Specimen preparation and protein extraction
This study included a total of 199 patients. For the discovery set, 40 samples (normal and cancerous tissues) were used. For the first validation set, we analyzed 85 cases with matching normal and cancerous tissue from the same patient. Serum analysis was done using 54 preoperative samples of RCC and 36 normal individuals. Urine analysis was done on 9 normal persons compared to 21 pre-operative urines from ccRCC patients.
ccRCC and matching normal kidney tissues from the same patient were obtained from nephrectomy specimens at St. Michael's Hospital, Toronto, Canada within 15 min post-surgery and snap frozen in liquid nitrogen. As ccRCC is known to arise from the proximal tubules [8], the kidney cortex is considered a suitable representation of normal kidney [9]. All specimens were histologically confirmed.  Table 1.
Tissues were prepared as described previously [10;11]. Briefly, tissues were homogenized in a protease-inhibitor cocktail (Roche). Cell debris was separated and clear supernatant was used for analysis. A reference sample was prepared from a pool of 30 normal kidney tissues. Protein concentrations were determined using the Bradford assay (Sigma-Aldrich, St. Louis, USA) [11;12].

iTRAQ sample labeling and strong cation exchange (SCX) chromatography
For iTRAQ LC-MS analysis, 100 µg of sample were denatured, disulfide bonds were reduced, and the cysteine residues were blocked as per the iTRAQ protocol (Applied Biosystems, Foster City, CA) and digested with trypsin. Samples were then labeled with the iTRAQ tags (Supplementary Table 2). Labeling of the reference sample was randomized for each set to eliminate any potential bias associated with a particular iTRAQ reporter tag. Samples were then dried using a vacuum centrifuge.
The iTRAQ sets were dissolved in 1.7 mL of Buffer A (Supplementary Table 3) and filtered using a 0.45-µm syringe filter (Millipore, Cambridge, ON, Canada). Each set was then separated by off-line SCX chromatography using an HP1050 HPLC instrument (Agilent, Palo Alto, CA). Separation was performed using a linear binary gradient over 1 h. Buffer C was used to strip the column after the run. A total of 30 SCX fractions were collected per iTRAQ set. These fractions were dried using a vacuum centrifuge. www.impactjournals.com/oncotarget

Reverse phase (RP) LC-MS
The SCX fractions were analyzed in triplicate using a Nanobore LC system (LC Packings, Amsterdam, Netherlands) and a QSTAR Pulsar mass spectrometer (Applied Biosystems/MDS SCIEX, Foster City, CA) in positive ion mode. Fractions were re-dissolved in 16 µL of eluant A [94.9% deionized water, 5.0% methanol, and 0.1% formic acid (pH 3)]. For subsequent fractions, the amount of eluant A was incremented by 2 µL over the preceding fraction to accommodate the increase in the amount of KCl. A 1-µL aliquot of the sample (~1 µg of total peptides) was loaded onto a C18 RP pre-column (LC Packings: 300 µm x 5 mm) and desalted before separation on an RP analytical column (75-µm x 150-mm packed inhouse with 3-µm Kromasil C18 beads with 100 Å pores, The Nest Group, Southborough, USA). Eluant A was used to load the sample onto the C18 pre-column at a flow rate of 25 µL min -1 . After 4 min, the C18 pre-column was switched in-line with the RP analytical column. Separation was performed at 100 nL min -1 using a nonlinear binary gradient (Supplementary Table 4) starting with eluant A and transitioning to eluant B (5.0% deionized water, 94.9% methanol, and 0.1% formic acid).
MS data were acquired in information-dependent acquisition (IDA) mode using the Analyst QS 1.1 software (Applied Biosystems/MDS SCIEX). The LC-MS analysis was performed using a 1-s TOF-MS survey scan from 400 to 1500 Da, followed by four, 2-s product-ion scans, from 80 to 2000 Da, of the four most-abundant ion peaks in the survey scan. The collision energy (CE) was automatically controlled by the IDA CE parameter script. Switching criteria were set for ions with m/z ≥ 400 and <1500, charge states of +2 to +4, and abundances of ≥10 counts. Using Analyst QS 1.1 controlled dynamic exclusion, former target ions were excluded for 30 s, and ions within a 100ppm window were ignored. Precursor ion exclusion lists were used to minimize redundancy.

Protein identification by Protein Pilot
MS data of each fraction was used to identify proteins by searching a concatenated Swissprot/Panther database of 66082 distinct human protein entries (version June 2, 2010). The database was searched using Protein Pilot software, version 2.0.1 (AB SCIEX, Foster City, USA), which uses the Paragon algorithm [13]. Protein identification was performed at a confidence threshold of 95% (Protein Pilot Unused score ≥1.3) with MMTS selected as cysteine modification, with the search option 'emphasis on biological modifications' checked, and with one of missed and/or non-specific cleavages permitted. Peptide and protein summaries were generated.
To minimize redundancy in subsequent iteration, a precursor ion exclusion list, generated in-house, was added to the acquisition method after each iteration as described previously [14]. Tolerance windows for exclusion were set at 100 ppm for m/z and 360 s for elution time.
iTRAQ ratio re-calculation and identification of dysregulated proteins To identify non-redundant proteins, data acquired for all 25 fractions from each iTRAQ set injected in triplicate were searched against a database that was created by concatenating the Swissprot human protein database and its reverse (as of June 2, 2010). Only proteins identified with local false discovery rate (FDR) ≤ 5% were considered for further analysis [15].
Proteins identified in seven iTRAQ sets were compiled and matched by accession numbers. Redundant proteins and peptides, and proteins identified in reverse sequence were removed from the list. To improve the confidence of protein quantitation, the mean expression iTRAQ ratios of the proteins were re-calculated, using a script written in Matlab (version 7.7.0.471), based on the criteria that the protein must be identified by a minimum of three peptides, with ≥95% confidence, and with an expression ratio error factor (EF) <11.1%. To enhance confidence in the protein quantitation even more, we included only 95% of all quantified proteins with the lowest computed EF (which corresponds to a confidence > 0.05 in Supplementary Table 5) for further consideration. Proteins were considered to be dysregulated if iTRAQ ratios were ≥1.5 or ≤0.67 in ≥50% in ccRCC relative to reference samples.

Clustering analysis of dysregulated proteins in ccRCC samples
Proteins were included in the analysis if quantification was available in at least 50% of the samples. The average iTRAQ ratios were logarithmically transformed, for hierarchical clustering was used the City-block distance method. As a control, the samples were hierarchically clustered based on quantified proteins without dysregulated proteins. Hierarchical clustering analysis was performed using the Cluster 3.0 software and the result was visualized with the TreeView software [16]. To assess the statistical significance of the difference between proteins expressions in ccRCC and normal kidney samples were used the Welch's two-tailed t-test. Value of p<0.05 was considered as significant. www.impactjournals.com/oncotarget

Selection of candidate ccRCC markers
Dysregulated proteins were selected for further verification if they had a documented role in tumorigenesis and a potential to be identified in serum. We classified a dysregulated protein as "secreted" if it satisfied at least one of the following four criteria: (1) its subcellular location is extracellular or membrane-bound; (2) it is non-classically secreted by the exosome pathway, (3) it is classically secreted, according to SignalP 4.0 analysis; or (4) it is nonclassically secreted, according to SecretomeP 2.0 analysis.

Pathway and protein-protein interaction analysis
Dysregulated proteins were subjected to pathway analysis using Reactome and KEGG software. Predicted protein-protein interactions were generated and visualized using STRING 9.0 software [17]. Parameters used for species and confidence were, "homo sapiens" and "medium confidence (0.400)," respectively. Clusters of proteins were determined through confidence levels of surrounding protein-protein interactions.

Tissue microarray construction and IHC
Tissue microarray (TMA) blocks were constructed as described previously [14] and contained 85 cases of ccRCC and matched normal kidney tissue from the same patient. Sections were deparaffinized, hydrated in ethanol, pre-treated in a microwave oven for 20 min at 800 W in 1 L of citrate buffer (0.01 M, pH 6.0), and incubated with hydrogen peroxide (0.3% v/v) in PBS for 15 min. Sections were blocked with 10% fetal bovine serum and incubated with primary antibodies for Hsp27, ENO1, LDHA, and HSPE1 overnight at 4°C. Protein expression was detected with Dako LSAB+ kit (Dako Cytomation) and counterstained with Mayer's hematoxylin.

ELISA
Human ELISA kits were purchased for Hsp27 and ENO1 (AbCam) and procedures were performed as recommended by the manufacturer. Urine samples were obtained from 11 patients with primary ccRCC and six patients with no malignancy. Serum samples were analyzed for a total of 21 patients with ccRCC and 9 patients with no malignancy.

RT2 Glucose Metabolism Profiler Assay
We used the Human Glucose Metabolism RT² Profiler™ PCR Array (SA Biosceinces) to examine the expression of 84 key genes involved in the regulation and enzymatic pathways of glucose and glycogen metabolism. Total RNA was extracted using the RNeasy Kit (Qiagen) and RNA integrity was assessed using the Agilent 2100 Bioanalyzer (Agilent). Four matched pairs of tumor and normal matched kidney tissue from the same patient were analyzed (a total of eight samples). Briefly, cDNA was synthesized from 1μg total RNA using the RT 2 First Strand Kit (SA Biosceinces) according to the manufacture's protocol. Reaction mixtures were incubated at 42ºC for 15 minutes, 95ºC for 5 minutes and held on ice until the PCR reaction. PCR cycling conditions were as follows: 95ºC for 10 min, followed by 40 cycles of 95ºC for 15 sec and 60ºC for one minute, on the Step One Plus (Life Technologies). Relative quantification Changes in expression were measured by obtaining the threshold cycle and normalizing to the average of three housekeeping genes including Hypoxanthine Phosphoribosyltransferase 1 (HPRT1), Ribosomal Protein L13a (RPL13A), and β-actin (ACTB). Fold change was calculated by 2 (-ΔΔCt) where ΔΔC t = [CAKI-1 cells transfected with miR-215 (C t target-C t control)] -[untransfected CAKI-1 cells (C t target-C t control)].