Next Article in Journal
An LC–MS Assay to Measure Superoxide Radicals and Hydrogen Peroxide in the Blood System
Next Article in Special Issue
Analytical Platforms for the Determination of Phospholipid Turnover in Breast Cancer Tissue: Role of Phospholipase Activity in Breast Cancer Development
Previous Article in Journal
Evaluation of MDA-MB-468 Cell Culture Media Analysis in Predicting Triple-Negative Breast Cancer Patient Sera Metabolic Profiles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Volatilomics Reveals Potential Biomarkers for Identification of Renal Cell Carcinoma: An In Vitro Approach

by
Filipa Amaro
1,*,
Joana Pinto
1,*,
Sílvia Rocha
1,2,
Ana Margarida Araújo
1,
Vera Miranda-Gonçalves
3,
Carmen Jerónimo
3,4,
Rui Henrique
3,4,5,
Maria de Lourdes Bastos
1,
Márcia Carvalho
1,6,* and
Paula Guedes de Pinho
1
1
UCIBIO, REQUIMTE, Laboratory of Toxicology, Faculty of Pharmacy, University of Porto, 4050-313 Porto, Portugal
2
Master in Oncology, Institute of Biomedical Sciences Abel Salazar–University of Porto (ICBAS-UP), 4050-313 Porto, Portugal
3
Cancer Biology & Epigenetics Group, Research Centre (CI-IPOP) Portuguese Oncology Institute of Porto (IPO Porto), 4200-072 Porto, Portugal
4
Department of Pathology and Molecular Immunology-Biomedical Sciences Institute (ICBAS), University of Porto, 4050-313 Porto, Portugal
5
Department of Pathology, Portuguese Oncology Institute of Porto (IPO Porto), 4200-072 Porto, Portugal
6
UFP Energy, Environment and Health Research Unit (FP-ENAS), University Fernando Pessoa, 349, 4249-004 Porto, Portugal
*
Authors to whom correspondence should be addressed.
Metabolites 2020, 10(5), 174; https://doi.org/10.3390/metabo10050174
Submission received: 14 April 2020 / Revised: 22 April 2020 / Accepted: 24 April 2020 / Published: 27 April 2020

Abstract

:
The identification of noninvasive biomarkers able to detect renal cell carcinoma (RCC) at an early stage remains an unmet clinical need. The recognition that altered metabolism is a core hallmark of cancer boosted metabolomic studies focused in the search for cancer biomarkers. The present work aims to evaluate the performance of the volatile metabolites present in the extracellular medium to discriminate RCC cell lines with distinct histological subtypes (clear cell and papillary) and metastatic potential from non-tumorigenic renal cells. Hence, volatile organic compounds (VOCs) and volatile carbonyl compounds (VCCs) were extracted by headspace solid-phase microextraction (HS-SPME) and analyzed by gas chromatography–mass spectrometry (GC–MS). Multivariate and univariate analysis unveiled a panel of metabolites responsible for the separation between groups, mostly belonging to ketones, alcohols, alkanes and aldehydes classes. Some metabolites were found similarly altered for all RCC cell lines compared to non-tumorigenic cells, namely 2-ethylhexanol, tetradecane, formaldehyde, acetone (increased) and cyclohexanone and acetaldehyde (decreased). Furthermore, significantly altered levels of cyclohexanol, decanal, decane, dodecane and 4-methylbenzaldehyde were observed in all metastatic RCC cell lines when compared with the non-metastatic ones. Moreover, some alterations in the volatile composition were also observed between RCC histological subtypes. Overall, our results demonstrate the potential of volatile profiling for identification of noninvasive candidate biomarkers for early RCC diagnosis.

Graphical Abstract

1. Introduction

Renal cell carcinoma (RCC) accounts for more than 90% of all kidney malignancies [1] and represents the second most lethal urological cancer [2]. Comprising a heterogeneous group of tumors, RCC can be subdivided into several histological subtypes according to distinct clinical and histological features and different outcomes [3,4]. The main groups include clear cell renal cell carcinoma (ccRCC), papillary renal cell carcinoma (pRCC) and chromophobe renal cell carcinoma (chRCC), among other rare types [3]. ccRCC is the most prevalent histological subtype, representing more than 75% of all RCCs [5], followed by pRCC, which represents 10%–15% of cases [6].
RCC has a favorable prognosis when detected at an early stage [7]. However, this disease prevails asymptomatically until the advanced stages [8] and no satisfactory biomarkers for clinical management are currently used [9,10]. Despite the use of imaging techniques, such as the computed tomography scan, ultrasound exam and/or magnetic resonance, it is indispensable to perform a biopsy of the tumor to obtain a precise diagnosis [11,12]. Hence, the development of accurate and non-invasive diagnostic methods based on detection of specific biomarkers in an early stage is of upmost relevance for RCC prognosis and follow-up [9,10].
The use of metabolomics for investigation of novel biomarkers in oncobiological studies has been rising over the years [13]. Metabolomics, defined as the study of endogenously produced low molecular weight compounds, enables the identification of alterations on metabolite levels that reflect the biological activity of cancer cells [14]. RCC has been considered a metabolic disease [15] due to several alterations in metabolites associated with energy metabolism, particularly those involved in pathways for sustaining cell proliferation [10,16]. The majority of RCC metabolomic studies aims to discriminate RCC patients from controls with diagnostic purposes, but also to investigate the metabolic deviations associated with the onset of RCC, through the use of mass-spectrometry (MS) and nuclear magnetic resonance spectroscopy (NMR) [8,10,17,18,19]. These studies have found increased levels of several types of carnitines and metabolites from fatty acid metabolism in tissue and urine from RCC patients compared to control individuals [17]. The levels of carbohydrate-related metabolites as lactate, pyruvate, glutamate and glucose, among others, have also been found altered in serum/plasma [19], renal tissue [8,18] and urine from RCC patients [19]. Among all biological matrices used in metabolomic studies, urine may be considered the preferential matrix for RCC biomarker discovery since it can offer a non-invasive tool for RCC diagnosis allied with the enormous advantage to be in direct contact with the kidney. However, urine also reflects the metabolic alterations associated factors not related with the disease state, usually called confounding factors, such as age, diet, gender and medication, among others, thus introducing “metabolomic variability” and/or masking valid information [16]. To overcome this issue, we propose herein a bottom-up metabolomic approach (matrices of increasing complexity) focused on the investigation of potential RCC biomarkers in in vitro models to be confirmed and validated in urine from RCC patients in the future.
Volatile organic compounds (VOCs) are end products of cellular metabolism excreted by cells [20,21] and the analysis of those compounds in the exometabolome can provide the metabolic signature directly originated from RCC cells. Scientific evidence has demonstrated that the alterations in the volatile part of exometabolome (VOCs and volatile carbonyl compounds (VCCs) profiles) reflect the genetic mutations, adaptations and modifications in biochemical pathways of tumor cells [21]. Therefore, these compounds can be used as potential biomarkers for RCC detection and diagnosis. The use of in vitro models to study the volatile signatures of cancer cells, thus avoiding the bias introduced by confounding factors, has already been reported by other authors [21,22,23,24,25]. Indeed, deviations in the levels of VOCs have been proven for prostate [21], lung [22,23], bladder [24] and breast [25] cancer cells, whereas the dysregulations occurring in the volatile exometabolome profile of RCC cells are still unexplored. The volatile signature of RCC has only been investigated using more complex matrices, such as urine, for which significant differences in VOC levels have been found between patients and control individuals [26,27].
Hence, this study aimed to investigate, for the first time, the performance of VOCs and VCCs in discriminating RCC cell lines (769-P, 786-O, Caki-1, Caki-2 and ACHN) from one non-tumorigenic renal cell line (HK-2). We further investigated which metabolites might discriminate RCC cell lines according to their histological subtype (ccRCC vs. pRCC) and metastatic potential (metastatic vs. non-metastatic RCC). Prior to analysis by gas chromatography-mass spectrometry (GC–MS), VOCs present in the extracellular medium were directly extracted by headspace solid-phase microextraction (HS-SPME), while VCCs were first derivatized before HS-SPME to enhance the detection of aldehydes (CHO-) and ketones (CO-) [28]. Knowledge on deviations occurring in the levels of volatiles in RCC cells may provide new and specific biomarker candidates to be validated in urines from RCC patients in the future and, therefore, to obtain a powerful tool enabling faster non-invasive RCC diagnosis.

2. Results

In this study, the exometabolome volatile profile of several RCC cell lines (769-P, 786-O, Caki-1, Caki-2 and ACHN) and one non-tumor renal cell line (HK-2) was analyzed to investigate whether volatile compounds might be able to discriminate tumorigenic from non-tumorigenic cells. The two most common histological subtypes, ccRCC and pRCC, at different disease stages were selected for this study, thus extending our investigation to the metabolic differences that characterize the histological subtype and tumor stage. The number of peaks detected in the extracellular medium in VOCs and VCCs analyses were 105 and 111, respectively (Figure S1). A very good reproducibility was achieved for both analytical methodologies, since all QCs were clustered in the unsupervised analysis (Figure S2). For VOCs, the principal component analysis (PCA) model including all cell lines and blanks (Figure 1A) unveiled one main cluster including all cell lines and another cluster of blanks (culture medium in the same conditions but without cells). For VCCs, three different clusters were observed in PCA (Figure 1B), namely one including the non-tumor cell line (HK-2) and blanks, a second including the cancer cell lines, and the last composed only by Caki-2. These results suggested that the RCC cell lines exhibited a distinct volatile profile when compared with the non-tumor cell line and blanks.

2.1. Volatile Exometabolome Signature of RCC Cell Lines versus the Non-Tumorigenic Cell Line

Pairwise multivariate supervised analysis (PLS-DA) was performed between each RCC cell line and the non-tumorigenic cell line for both analytical methodologies, VOCs and VCCs (Figure 2A,B, respectively). The models showed a very robust discrimination for all comparisons with high predictive ability (Q2 > 0.5).
The compounds altered between each tumor cell line versus the non-tumor are shown in Table 1 and Table 2. Overall, the compound classes found altered comprised several alcohols, alkanes, alkenes, aldehydes, benzene derivatives and ketones. Interestingly, a higher number of discriminant metabolites was found for ccRCC vs. HK-2 (26 compounds; Table 1) when compared with pRCC vs. HK-2 (15 compounds; Table 2). Some alterations were found in common for all ccRCC and pRCC cell lines compared with HK-2, including increased levels of 2-ethylhexanol, tetradecane, formaldehyde, acetone and an unknown compound (Un RT 10.18, m/z 58), and the decrease in the levels of acetaldehyde and cyclohexanone. These results suggest the existence of a common regulatory mechanism in ccRCC and pRCC cells.
The analysis of the volatile profile of blanks enabled the interpretation of dysregulations occurring in RCC cells in terms of excretion or consumption of culture medium nutrients. The heatmap depicted in Figure 3 shows the average peak area of each altered compound in blanks, HK-2, ccRCC and pRCC cell lines. The hierarchical clustering shown for compounds helps to interpret the patterns of compound concentrations that changed across different groups. Hence, cyclohexanol, decanal, 2,6-dimethyl-7-octen-2-ol, levomenthol, cyclohexanone, acetaldehyde and acetophenone showed a trend for lower levels in ccRCC cell lines compared with HK-2 and blanks, indicating the occurrence of an increased consumption by these cell lines, except for cyclohexanol and decanal in Caki-1 for which the levels were similar to that of HK-2 and blanks. Tetradecane, 2-ethylhexanol and acetone unveiled a tendency for higher levels in ccRCC cell lines compared with HK-2 and blanks, suggesting an increased excretion of those compounds by ccRCC cell lines. 2-Ethoxy-2-methylpropane, 2,4-dimethyl-1-heptene and 4-methylbenzaldehyde were consumed by all cell lines but at a lower extent in 769-P and 786-O compared with HK-2 (Table 1). Benzaldehyde and formaldehyde showed a tendency for consumption in ccRCC and HK-2 cell lines, but a different profile was observed for both metabolites with formaldehyde showing increased levels in ccRCC lines in relation to HK-2, while for benzaldehyde only 786-O unveiled statistically significant results when compared with HK-2 (Table 1). Decane, 3-carene, 3-methylbenzaldehyde, styrene, xylene and ethylbenzene showed a distinct profile in HK-2 and Caki-1, which seem to consume higher levels of those compounds (Table 1 and Figure 3). Regarding the pRCC cell lines, Caki-2 and ACHN showed a larger excretion of 2-ethylhexanol, tetradecane and acetone, whereas only Caki-2 unveiled a higher excretion of cyclohexanol (Table 2 and Figure 3). 3-Methylbenzaldehyde and 4-methylbenzaldehyde showed a different variation pattern for Caki-2 and ACHN (Figure 3). Formaldehyde was consumed by pRCC and HK-2 cell lines but at a lower extent for Caki-2 and ACHN. 2-Pentadecanone, acetophenone and cyclohexanone exhibited a distinct profile in Caki-2 when compared with HK-2, ACHN and blanks, characterized by higher consumption by this cell line. Acetaldehyde revealed an increased consumption in both pRCC cell lines (Table 2).

2.2. Volatile Exometabolome Signature of Metastatic versus Non-Metastatic RCC Cell Lines and ccRCC versus pRCC Cell Lines

To investigate the volatile exometabolome signature of RCC cell lines with different metastatic potential, the same strategy was applied to compare the extracellular medium of metastatic and non-metastatic RCC cell lines (Figure 4). A good separation was obtained for metastatic vs. non-metastatic ccRCC cells in VOCs analysis (Q2 > 0.7; Figure 4a), whereas a poor separation (Q2 = 0.211) was observed in VCCs (Figure 4c). For pRCC, a good separation with high predictive power (Q2 > 0.9) was observed in both analytical methodologies (Figure 4b,d).
The volatile exometabolome signature of metastatic ccRCC cell lines unveiled statistically significant alterations in the levels of 15 compounds when compared with non-metastatic cell lines (Table 3). For pRCC, 17 compounds were found significantly altered between metastatic and non-metastatic cells. Interestingly, most of these alterations were specific of the RCC histological subtype with only a decrease in the levels of decane, dodecane and 4-methylbenzaldehyde found in common for ccRCC and pRCC. In general, metastatic ccRCC cell lines mostly showed dysregulations in the levels of alkanes and alkenes (decreased and increased) and benzene derivatives (decreased) when compared with non-metastatic cells, while metastatic pRCC cell lines mostly showed alterations in the levels of aldehydes (decreased) and ketones (decreased and increased).
Finally, the comparison of the volatile exometabolome between ccRCC and pRCC unveiled poor separation in the PLS-DA models of VOCs and VCCs analysis (Figure S3). In the PLS-DA model of VOCs (Figure S3A), the extracellular media of the Caki-1 cell line was wrongly classified in LV1 positive. Nonetheless, the model showed high predictive power and some statistically significant alterations between ccRCC and pRCC were disclosed, namely in the levels of 2,4-dimethyl-1-heptene, 3-carene, 4-methylbenzaldehyde, 4-methylheptane and an unknown compound (Un RT 8.40 m/z 69; Table 3).

3. Discussion

The profiling of volatile compounds present in the culture medium (exometabolome) unveiled important alterations associated with the excretion and consumption of nutrients by RCC cells lines when compared with a non-tumor cell line. In addition, different RCC histological subtypes (ccRCC and pRCC), as well as different stages of the disease (metastatic and non-metastatic) disclosed characteristic exometabolome signatures. These metabolic alterations may constitute a panel of potential biomarkers yet to be investigated in urine from RCC patients, compared with cancer-free individuals. Common to all cancer cell lines under study, we found increased levels of 2-ethylhexanol, tetradecane, formaldehyde, acetone and an unknown compound (Un RT 10.18 m/z 58), and decreased levels of acetaldehyde and cyclohexanone. These results suggest that different histological RCC subtypes have, at least in part, a common regulatory mechanism that affects the metabolic pathways in which these compounds are intermediates.
In a general way, RCC cell lines present high levels of hydrocarbons (alkanes and alkenes) in the extracellular medium. The main mechanism that affects the production of hydrocarbons is the oxidative stress, since the production of these volatile metabolites is mainly due to the induction of ROS-mediated polyunsaturated fatty acid peroxidation [30]. In fact, some studies have already shown elevated ROS levels in RCC [31,32] and have also provided evidence that lipid peroxidation is an important mechanism in renal carcinogenesis [33,34].
Another major change observed in this study was a significant increase of some alcohols (namely 2-ethylhexanol) in RCC cell lines when compared to the non-tumorigenic cell line. These results can be explained by a possible up-regulation of cytochrome P450 (CYP450) enzymes, which can hydroxylate several volatiles formed during ROS-mediated lipid peroxidation, leading to the production of the corresponding alcohols [30]. This theory is supported by some studies [33,35,36] whose results demonstrated that some CYP isoforms (such as CYP1B1, CYP1A1, CYP2C19, CYP2D6, CYP2E1, CYP3A4 or CYP3A5) are altered in RCC patients. Furthermore, 2-ethylhexanol has already been found increased in preoperative vs. postoperative urine samples from RCC patients [27], indicating that it may be considered a potential translatable biomarker for non-invasive measurement in urine. However, 2-ethylhexanol has also been reported as a potential urinary biomarker of lung [37] and prostate [38] cancers, which indicates a lack of specificity towards RCC. Nevertheless, this biomarker can be combined with others to define a robust panel of biomarkers for early diagnosis of RCC.
In addition, our metabolomic results also showed that most of the aldehydes (e.g., 4-methylbenzaldehyde, 3-methylbenzaldehyde and formaldehyde) were significantly increased in RCC cell lines. Increased levels of aldehydes may derive from different mechanisms including a direct origin in the reduction of hydroperoxides by CYP450 or from the conversion of some alcohols to their corresponding aldehydes mediated by alcohol dehydrogenases (ADHs) [30]. Since the amount of some lipids in cancer cells membranes is greater than in normal cells, it can be assumed that the increase in aldehyde levels in RCC cell lines may be linked to changes in the membrane lipid composition subsequent to increased oxidative stress [30,31,32]. Furthermore, it is also known that RCC cells have a significantly higher activity of total ADH and ADH class I that is exacerbated in patients with a more advanced stage [39,40], meaning that tumoral cells have an increased ability to produce aldehydes, which can intensify carcinogenesis. In agreement, a significant enhancement of aldehydes has been previously shown in urine from RCC patients compared to healthy controls [27].
Remarkably, acetaldehyde, another well-known aldehyde, appears significantly decreased in all RCC cell lines when compared to the non-tumorigenic cell line. This finding is opposite to what would be expected based on aforementioned higher activity of ADH in renal cancer cells. The ability of acetaldehyde to bind to DNA and cellular proteins to form adducts may explain the results verified in the exometabolome [40]. Moreover, acetaldehyde can lead to activation of proto-oncogenes, inactivation of tumor suppressor genes in replicating cells and inhibition of many important enzymes of DNA synthesis pathways, playing an important role in carcinogenesis [40]. Further studies are needed to explore whether RCC cells possess/induce detoxifying mechanisms for this toxic and reactive aldehyde.
Ketones were also found altered in the exometabolome of RCC cells compared with the non-tumor cell line, with special emphasis for acetone whose levels are significantly increased. Several studies [15,23,41] have associated this class of volatile compounds with altered pathways in cancer, namely the propensity of tumor cells to favor their metabolism via glycolysis (Warburg effect) with a consequent increase in lactate secretion. These will lead to an accumulation of acetyl-CoA, which in turn results in an increased production of ketone bodies, namely acetone [30].
In contrast, cyclohexanone was found decreased in the exometabolome of all RCC cells. Of note, this finding is not in agreement with the previous study by Wang et al. [27] where they showed significantly higher levels of cyclohexanone in preoperative urines from RCC patients. Although these contradictory results may be a result of sources of residual variation as diet, medication use, gender, age and comorbidities, among others, which may affect urine, further studies should be carried out in the future. Noticeably, this volatile compound is highly consumed by the pRCC cell line Caki-2 compared to other cancer cell lines. Caki-2 revealed a different profile since the levels of 2-pentadecanone and an unknown compound (Un RT 9.80, m/z 59) were found significantly altered only in this pRCC cell line. Contrary to other RCC histological subtypes, ccRCC and pRCC arise from adjacent parts of the proximal tubule [42], which may justify the poor discrimination found between both subtypes. Nevertheless, some volatile compounds were found to be discriminative between ccRCC and pRCC subtypes, mainly belonging to hydrocarbons, aldehydes and ketone classes. Although differences in the metabolic profile of ccRCC and pRCC subtypes have already been reported in terms of amino acid and fatty acid metabolism, whose levels were found to be increased in ccRCC [43], thus being associated with the higher aggressiveness of ccRCC comparing to the pRCC subtype [42,44], our study reports for the first time differences in the volatile exometabolome of these two RCC subtypes.
Regarding the comparison metastatic vs. non-metastatic stage in both the clear cell and the papillary RCC cell lines, the results suggest a variation of the volatile compounds in a stage-dependent manner since statistically significant differences in the levels of several metabolites were found. Curiously, the separation between metastatic and non-metastatic cell lines was more robust in the pRCC subtype, which justifies the markedly altered profile of the pRCC cell line, Caki-2. These results are consistent with the findings reported by Schaeffeler et al. [42], which demonstrated that ccRCC-derived metastases showed similar metabolite levels when compared with primary ccRCC. Our results revealed common alterations in the levels of several alkanes (decanal, decane and dodecane), one aldehyde (4-methylbenzaldehyde) and one alcohol (cyclohexanol) in all metastatic cell lines compared with the non-metastatic. Indeed, several ketones and aldehydes were only found altered in the pRCC metastatic cell line, ACHN. However, our results disclosed decreased levels of aldehydes in the metastatic cell lines, which may be used in amino acid synthesis, as proposed by some other authors [23].
This study demonstrates for the first time, to our knowledge, the potential of volatile compounds in discriminating RCC cells from nonmalignant cells, according to the histological subtype and tumor stage. However, the metabolic pathways that lead to VOCs and VCCs production are not yet clearly understood, hindering an accurate interpretation of the results. Nevertheless, our findings support the hypothesis that altered energy production in RCC cells enables them to survive under oxidative stress conditions and to migrate to other tissues and form metastases [45]. This study also reinforces that the investigation of the volatile signature of the exometabolome using cancer cell lines can provide candidate biomarkers with potential to be translatable for urine of RCC patients, which may ultimately lead to the development of volatile sensor-based approach for non-invasive diagnosis of RCC.

4. Materials and Methods

4.1. Chemicals

All chemicals were of analytical grade and dissolved in deionized water unless otherwise indicated. RPMI-1640 medium, sodium chloride (NaCl, ≥ 99.5%), O-(2,3,4,5,6-pentafluorobenzyl) hydroxylamine (PFBHA, ≥ 99%) and thymol (≥ 98.5%) were obtained from Sigma-Aldrich (St. Louis, Missouri, USA). The antibiotic mixture penicillin/streptomycin (10.000 U/mL/10.000 mg/mL), heat-inactivated fetal bovine serum (FBS) and trypsin-EDTA (0.25%) were purchased from GIBCO Invitrogen (Barcelona, Spain). Hydrochloric acid (HCl) and sodium hydrogencarbonate were obtained from Merck (Darmstadt, Germany).

4.2. Cell lines and Culture Conditions

Five human RCC immortalized cell lines (769-P, 786-O, ACHN, Caki-1 and Caki-2) and one non-tumorigenic renal cell line (HK-2) were purchased from the American Type Culture Collection (ATCC; Manassas, VA, USA) and their characteristics are summarized in Table S1. The selection of tumorigenic cell lines was carried out to allow the study of the most common RCC subtypes (ccRCC and pRCC) at different disease stages. Thus, 769-P and 786-O are recognized as primary ccRCC while Caki-1 is a metastatic ccRCC cell line [46]. Furthermore, and despite some controversies in histological classification, Caki-2 has been recognized as a primary pRCC cell line expressing papillary characteristics and ACHN is considered a metastatic pRCC cell line [46]. HK-2 is an immortalized proximal tubule cell line derived from normal human kidney [47]. All cell lines were grown in RPMI-1640 supplemented with 10% of FBS and 1% penicillin/streptomycin and incubated at 37 °C and 5% CO2. All cell lines were routinely tested for Mycoplasma spp. contamination (PCR Mycoplasma Detection Set, TaKaRa).
For each cell line, two different cell stock vials were grown in parallel and the experiments carried out up to 7 consecutive passages in plastic T25 culture flasks, after an adaptation stage of at least three passages for all cell lines. After 24 h, a confluence of 50%–60% was observed, the culture medium was discarded and replaced with 7 mL of fresh culture medium (which was prepared in enough amount to maintain the same composition in all collected flasks, thus reducing possible bias). The culture flasks were maintained in the incubator for another 48 h reaching near confluence, together with blanks (T25 culture flasks with the same culture medium but without cells). After 48 h, the extracellular medium from each culture flask with cells and blanks was collected and centrifuged (5 min, 2000g, 4 °C). Then, the supernatant was divided into two 2.5 mL aliquots and immediately stored at −80 °C until analysis. A pool of all samples and blanks, called quality control (QCs), was also prepared under the same collection and storage conditions.

4.3. Sample Preparation and Volatile Extraction by HS-SPME

Sample preparation was performed as described in previous studies published by our group [21,24,28,48]. Briefly, stored samples were first thawed slowly on ice to minimize the loss of volatiles. For VOCs analysis, 2 mL of each extracellular culture medium was transferred to a 10 mL glass vial with 30 μL of the internal standard (thymol, 2 mg/L) and 0.43 g of NaCl (salting out effect). Then, VOCs were extracted through HS-SPME using a divinylbenzene/carboxen/polydimethylsiloxane (DVB/CAR/PDMS) fiber with an incubation time of 11 min and extraction of 30 min at 44 °C, under continuous stirring (250 rpm) [21,24,48]. For VCCs analysis, 2 mL of each culture medium was placed into a 10 mL glass vial with 47 μL of the derivatizing agent (O-(2,3,4,5,6-pentafluorobenzyl) hydroxylamine, 40 g/L) and the compounds were extracted using a polydimethylsiloxane/divinylbenzene (PDMS/DVB) fiber with an incubation time of 6 min at 62 °C and extraction using the same temperature during 51 min, under continuous stirring (250 rpm) [28,48]. All samples were randomly injected, and the QCs were injected at the same conditions, on every 6 samples, to evaluate the reproducibility in both analyses.

4.4. GC–MS Analysis: Equipment and Conditions

VOCs detection was accomplished using a 436-GC system (Bruker Daltonics, Fremont, CA) coupled to a SCION Single Quadrupole (SQ) mass detector and a Bruker Daltonics MS workstation software (version 8.2). A fused silica capillary column Rxi-5Sil MS (30 m × 0.25 mm × 0.25 μm; RESTEK Corporation, U.S., Bellefonte, Pennsylvania) was used and helium C-60 (Gasin, Portugal) was the carrier gas (flow rate 1 mL/min). The oven temperature was fixed at 40 °C for 1 min, increasing to 250 °C (rate 5 °C/min), held for 5 min, followed by increasing to 300 °C (rate 5 °C/min) and held for 1 min. The MS detector was operated in electron impact mode (70 eV) at 260 °C, the transfer line temperature was 250 °C and the manifold temperature was 40 °C. Data acquisition was performed in full scan mode with a mass range between 40 and 350 m/z at a scan rate of 6 scans/s. VCCs were analyzed using a 436-GC system (Bruker Daltonics, Fremont, CA) coupled to a EVOQ Triple Quadrupole (TQ) mass detector and a Bruker MS workstation software version 8.2. The column, carrier gas and oven temperature conditions were the same described above for VOCs analysis. The MS detector was operated in the electron impact mode (70 eV) at 270 °C, the temperature of transfer line and manifold were 260 °C and 40 °C, respectively. Data acquisition was performed in the full scan mode considering a 50–600 m/z mass range.

4.5. Compound Identification and GC–MS Data Pre-Processing

Compounds detected in both VOCs and VCCs analytical methodologies were identified by comparison of the MS fragmentation with the mass spectra present in the National Institute of Standards and Technology (NIST 14) database and by comparing the experimental Kovats retention index (RI) with literature. When possible, the identification was confirmed by comparison of the retention time and MS spectra of the samples with the commercially available standard compounds analyzed using the same conditions.
Before statistical analysis, the data obtained from GC–MS were converted into a manageable format (netcdf) and preprocessed using MZmine 2.53 [49]. The preprocessing steps consisted of filtering, peak detection, chromatogram deconvolution and alignment. The parameters used for preprocessing of VOCs profile were: RT range 2–32 minutes; m/z range 50–250; MS data noise level 1 × 104; m/z tolerance 0.3; chromatogram baseline level 5 × 104 and peak duration range 0.05–1 min. For VCCs analysis, the parameters used were RT range 9–45 minutes; m/z range 50–500; MS data noise level 1 × 104; m/z tolerance 0.3; chromatogram baseline level 1 × 104; peak duration range 0.03–0.2 min. GC contaminants from the chromatographic column or fiber (among others) were manually removed from the final data matrix.

4.6. Statistical Analysis

The statistical analysis pipeline used for VOCs and VCCs final matrix was the same. First, data were normalized by the total area (TA) of the chromatogram and scaled to Pareto. Then, PCA (unsupervised analysis) was applied to detect trends, outliers and to evaluate the distribution of the QCs, followed by PLS-DA (supervised analysis) to discriminate classes and to identify the specific metabolic signature of each cell line, in SIMCA 15.0.2 (Umetrics Umea, Sweden). A default sevenfold internal cross-validation and permutation testing were used to confirm the robustness of PLS-DA models. All volatile compounds with VIP (variable importance in the projection) higher than one were considered important for group discrimination and were subsequently submitted to the Mann–Whitney test. The false discovery rate (FDR) method was used to adjust p-values for multiple comparison using the MetaboAnalyst 4.0 [50], and compounds with p-value < 0.05 were considered statistically significant. In addition, the AUC (area under the curve) was computed for each discriminant compound using the same analytical platform. The effect size and the standard error were also computed as described by Berben et al. [29]. Finally, the average normalized area of each statistically significant compound in each sample group was represented in a heatmap in MetaboAnalyst 4.0.

5. Conclusions

In this study, we used an in vitro metabolomics approach to obtain the volatile signatures of RCC cells with distinct histological subtypes and stages. The results showed that alterations in the volatile profile are capable of discriminating RCC from non-tumorigenic renal cells. This unique volatile profile of RCC cells can be useful to unveil promising noninvasive biomarkers for early diagnosis in more complex biological matrices, such as urine. A multi-biomarker panel can provide a more robust and accurate tool for RCC diagnosis than the use of only one metabolite as a biomarker. Indeed, 2-ethylhexanol, tetradecane, formaldehyde, acetone, cyclohexanone and acetaldehyde were able to discriminate RCC cell lines from the non-tumor cell line and may constitute a potential panel of candidate biomarkers for RCC diagnosis. This study also demonstrated that the volatile profile of metastatic cell lines markedly differed from the non-metastatic in each RCC subtype, indicating some promising biomarkers applicable to the categorization of tumor stage, pending further investigation and validation.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/5/174/s1, Table S1: Characteristics of RCC cell lines (769-P, 786-O, Caki-1, Caki-2 and ACHN) and non-tumoral cell line (HK-2); Figure S1: HS-SPME-GC–MS chromatograms of (A) VOCs and (B) VCCs of QC culture medium. 1: 2-ethyl-hexanol; 2: acetophenone; 3: dodecane; 4: decanal; 5: dodecanal; 6: 2-pentadecanone; 7: formaldehyde; 8: acetaldehyde; 9: acetone; 10: tetradecane; 11: 2-octanone; 12: 3-methyl-benzaldehyde; Figure S2: PCA scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of QCs (dark blue, n = 15) vs. all cell lines and blanks (green, n = 88); Figure S3: PLS-DA scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of ccRCC cell lines (769-P, 786-O and Caki-1, light blue, n = 38) vs. pRCC cell lines (Caki-2 and ACHN, red, n = 24); In vitro RCC data.xlsx file containing all data obtained in this study.

Author Contributions

F.A. was responsible for the execution of the experimental work and data analysis. S.R. and A.M.A. supported cell culture and GC–MS analysis. J.P. helped with the statistical analysis of the data. V.M.-G., C.J. and R.H. were responsible for testing for Mycoplasma spp. contamination and provided conceptual advice. P.G.d.P., M.C. and M.d.L.B. contributed to the conception, design and management of the study. F.A. wrote the manuscript with input from J.P. and M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financed by FEDER—Fundo Europeu de Desenvolvimento Regional funds through the COMPETE 2020—Operacional Programme for Competitiveness and Internationalisation (POCI), by Portuguese funds through FCT—Fundação para a Ciência e a Tecnologia in the framework of the project POCI-01-0145-FEDER-030388 - PTDC/SAU-SER/30388/2017 and by national funds from FCT/MCTES (UIDB/MULTI/04378/2020). A.M.A. has a PhD fellowship funded by SFRH/BD/107708/2015 and V.M-G is a junior researcher funded by POCI-01-0145-FEDER-29043.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ljungberg, B.; Albiges, L.; Abu-Ghanem, Y.; Bensalah, K.; Dabestani, S.; Fernandez-Pello, S.; Giles, R.H.; Hofmann, F.; Hora, M.; Kuczyk, M.A.; et al. European Association of Urology Guidelines on Renal Cell Carcinoma: The 2019 Update. Eur. Urol. 2019, 75, 799–810. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, X.; Zhang, M.; Liu, X.; Sun, H.; Guo, Z.; Tang, X.; Wang, Z.; Li, J.; Li, H.; Sun, W.; et al. Urine Metabolomics for Renal Cell Carcinoma (RCC) Prediction: Tryptophan Metabolism as an Important Pathway in RCC. Front. Oncol. 2019, 9, 663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Muglia, V.F.; Prando, A. Renal cell carcinoma: histological classification and correlation with imaging findings*. Radiol. Bras. 2015, 48, 166–174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Cheville, J.C.; Lohse, C.M.; Zincke, H.; Weaver, A.L.; Blute, M.L. Comparisons of outcome and prognostic features among histologic subtypes of renal cell carcinoma. Am. J. Surg. Pathol. 2003, 27, 612–624. [Google Scholar] [CrossRef]
  5. Monteiro, M.; Carvalho, M.; Bastos, M.L.; De Pinho, P.G. Biomarkers in renal cell carcinoma: A metabolomics approach. Metabolomics 2014, 10, 1210–1222. [Google Scholar] [CrossRef]
  6. Deng, J.; Li, L.; Xia, H.; Guo, J.; Wu, X.; Yang, X.; Hong, Y.; Chen, Q.; Hu, J. A comparison of the prognosis of papillary and clear cell renal cell carcinoma: Evidence from a meta-analysis. Medicine 2019, 98, e16309. [Google Scholar] [CrossRef]
  7. Zheng, H.; Ji, J.; Zhao, L.; Chen, M.; Shi, A.; Pan, L.; Huang, Y.; Zhang, H.; Dong, B.; Gao, H. Prediction and diagnosis of renal cell carcinoma using nuclear magnetic resonance-based serum metabolomics and self-organizing maps. Oncotarget 2016, 7, 59189–59198. [Google Scholar] [CrossRef] [Green Version]
  8. Popławski, P.; Tohge, T.; Bogusławska, J.; Rybicka, B.; Tański, Z.; Trevino, V.; Fernie, A.R.; Piekiełko-Witkowska, A. Integrated transcriptomic and metabolomic analysis shows that disturbances in metabolism of tumor cells contribute to poor survival of RCC patients. Biochim. et Biophys. Acta (BBA) - Mol. Basis Dis. 2017, 1863, 744–752. [Google Scholar] [CrossRef] [PubMed]
  9. Farber, N.J.; Kim, C.J.; Modi, P.K.; Hon, J.D.; Sadimin, E.; Singer, E.A. Renal cell carcinoma: the search for a reliable biomarker. Transl. Cancer Res. 2017, 6, 620–632. [Google Scholar] [CrossRef] [Green Version]
  10. Gupta, A.; Nath, K.; Bansal, N.; Kumar, M. Role of metabolomics-derived biomarkers to identify renal cell carcinoma: a comprehensive perspective of the past ten years and advancements. Expert Rev. Mol. Diagn. 2019, 20, 5–18. [Google Scholar] [CrossRef]
  11. Weiss, R.; Lin, P.-Y. Kidney cancer: Identification of novel targets for therapy. Kidney Int. 2006, 69, 224–232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hutson, T.E. Renal Cell Carcinoma: Diagnosis and Treatment, 1994–2003. Bayl. Univ. Med Cent. Proc. 2005, 18, 337–340. [Google Scholar] [CrossRef] [Green Version]
  13. Armitage, E.G.; Ciborowski, M. Applications of Metabolomics in Cancer Studies. Adv. Exp. Med. Biol. 2017, 965, 209–234. [Google Scholar] [CrossRef] [PubMed]
  14. Verma, M.; Banerjee, H.N. Metabolomic Approaches in Cancer Epidemiology. Diseases 2015, 3, 167–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wettersten, H.I.; Abu Aboud, O.; Lara, P.N.; Weiss, R.H. Metabolic reprogramming in clear cell renal cell carcinoma. Nat. Rev. Nephrol. 2017, 13, 410–419. [Google Scholar] [CrossRef]
  16. Rodrigues, D.; Monteiro, M.; Jerónimo, C.; Henrique, R.; Belo, L.; Bastos, M.L.; De Pinho, P.G.; Carvalho, M. Renal cell carcinoma: a critical analysis of metabolomic biomarkers emerging from current model systems. Transl. Res. 2017, 180, 1–11. [Google Scholar] [CrossRef]
  17. Nizioł, J.; Bonifay, V.; Ossoliński, K.; Ossoliński, T.; Ossolińska, A.; Sunner, J.; Beech, I.; Arendowski, A.; Ruman, T. Metabolomic study of human tissue and urine in clear cell renal carcinoma by LC-HRMS and PLS-DA. Anal. Bioanal. Chem. 2018, 410, 3859–3869. [Google Scholar] [CrossRef] [Green Version]
  18. Gao, H.; Dong, B.; Jia, J.; Zhu, H.; Diao, C.; Yan, Z.-H.; Huang, Y.; Li, X. Application of ex vivo 1H NMR metabonomics to the characterization and possible detection of renal cell carcinoma metastases. J. Cancer Res. Clin. Oncol. 2012, 138, 753–761. [Google Scholar] [CrossRef]
  19. Falegan, O.S.; Ball, M.; Shaykhutdinov, R.A.; Pieroraio, P.M.; Farshidfar, F.; Vogel, H.J.; Allaf, M.E.; Hyndman, M.E. Urine and Serum Metabolomics Analyses May Distinguish between Stages of Renal Cell Carcinoma. Metabolites 2017, 7, 6. [Google Scholar] [CrossRef] [Green Version]
  20. Schmidt, K.; Podmore, I. Current Challenges in Volatile Organic Compounds Analysis as Potential Biomarkers of Cancer. J. Biomark. 2015, 2015, 981458. [Google Scholar] [CrossRef] [Green Version]
  21. Lima, A.R.C.M.; Araujo, A.M.; Pinto, J.; Jerónimo, C.; Henrique, R.; Bastos, M.L.; Carvalho, M.; De Pinho, P.G. Discrimination between the human prostate normal and cancer cell exometabolome by GC-MS. Sci. Rep. 2018, 8, 5539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jia, Z.; Patra, A.; Kutty, V.K.; Venkatesan, V. Critical Review of Volatile Organic Compound Analysis in Breath and In Vitro Cell Culture for Detection of Lung Cancer. Metabolites 2019, 9, 52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Jia, Z.; Zhang, H.; Ong, C.N.; Patra, A.; Lu, Y.; Lim, C.T.; Venkatesan, T. Detection of Lung Cancer: Concomitant Volatile Organic Compounds and Metabolomic Profiling of Six Cancer Cell Lines of Different Histological Origins. ACS Omega 2018, 3, 5131–5140. [Google Scholar] [CrossRef] [PubMed]
  24. Rodrigues, D.; Pinto, J.; Araujo, A.M.; Reis, S.; Jerónimo, C.; Henrique, R.; Bastos, M.L.; De Pinho, P.G.; Carvalho, M. Volatile metabolomic signature of bladder cancer cell lines based on gas chromatography–mass spectrometry. Metabolomics 2018, 14, 62. [Google Scholar] [CrossRef]
  25. Silva, C.; Perestrelo, R.; Silva, P.; Tomás, H.; Câmara, J.S. Volatile metabolomic signature of human breast cancer cell lines. Sci. Rep. 2017, 7, 43969. [Google Scholar] [CrossRef] [Green Version]
  26. Monteiro, M.; Moreira, N.; Pinto, J.; Pires-Luís, A.S.; Henrique, R.; Jerónimo, C.; Bastos, M.L.; Gil, A.M.; Carvalho, M.; De Pinho, P.G. GC-MS metabolomics-based approach for the identification of a potential VOC-biomarker panel in the urine of renal cell carcinoma patients. J. Cell. Mol. Med. 2017, 21, 2092–2105. [Google Scholar] [CrossRef]
  27. Wang, D.; Wang, C.; Pi, X.; Guo, L.; Wang, Y.; Li, M.; Feng, Y.; Lin, Z.; Hou, W.; Li, E. Urinary volatile organic compounds as potential biomarkers for renal cell carcinoma. Biomed. Rep. 2016, 5, 68–72. [Google Scholar] [CrossRef] [Green Version]
  28. Lima, A.R.C.M.; Pinto, J.; Azevedo, A.I.; Barros-Silva, D.; Jerónimo, C.; Henrique, R.; Bastos, M.L.; De Pinho, P.G.; Carvalho, M. Identification of a biomarker panel for improvement of prostate cancer diagnosis by volatile metabolic profiling of urine. Br. J. Cancer 2019, 121, 857–868. [Google Scholar] [CrossRef]
  29. Berben, L.; Sereika, S.M.; Engberg, S. Effect size estimation: Methods and examples. Int. J. Nurs. Stud. 2012, 49, 1039–1047. [Google Scholar] [CrossRef]
  30. Janfaza, S.; Khorsand, B.; Nikkhah, M.; Zahiri, J. Digging deeper into volatile organic compounds associated with cancer. Boil. Methods Protoc. 2019, 4, bpz014. [Google Scholar] [CrossRef]
  31. Ganesamoni, R.; Bhattacharyya, S.; Kumar, S.; Chauhan, A.; Mete, U.K.; Agarwal, M.M.; Mavuduru, R.; Kaushik, G.; Mandal, A.K.; Singh, S.K. Status of Oxidative Stress in Patients With Renal Cell Carcinoma. J. Urol. 2012, 187, 1172–1176. [Google Scholar] [CrossRef] [PubMed]
  32. Shanmugasundaram, K.; Block, K. Renal Carcinogenesis, Tumor Heterogeneity, and Reactive Oxygen Species: Tactics Evolved. Antioxid. Redox Signal. 2016, 25, 685–701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Šverko, A.; Sobočanec, S.; Kušić, B.; Mačak-Šafranko, Ž.; Šarić, A.; Leniček, T.; Kraus, O.; Andrisic, L.; Korolija, M.; Balog, T.; et al. Superoxide dismutase and cytochrome P450 isoenzymes might be associated with higher risk of renal cell carcinoma in male patients. Int. Immunopharmacol. 2011, 11, 639–645. [Google Scholar] [CrossRef] [PubMed]
  34. Gago-Dominguez, M.; Castelao, J.E. Lipid peroxidation and renal cell carcinoma: further supportive evidence and new mechanistic insights. Free. Radic. Boil. Med. 2006, 40, 721–733. [Google Scholar] [CrossRef]
  35. Li, H.; Li, W.-X.; Dai, S.; Guo, Y.; Zheng, J.; Liu, J.; Wang, Q.; Chen, B.; Li, G.-H.; Huang, J. Identification of metabolism-associated genes and pathways involved in different stages of clear cell renal cell carcinoma. Oncol. Lett. 2017, 15, 2316–2322. [Google Scholar] [CrossRef]
  36. McFadyen, M.; Melvin, W.T.; I Murray, G. Cytochrome P450 CYP1B1 activity in renal cell carcinoma. Br. J. Cancer 2004, 91, 966–971. [Google Scholar] [CrossRef] [Green Version]
  37. Hanai, Y.; Shimono, K.; Matsumura, K.; Vachani, A.; Albelda, S.; Yamazaki, K.; Beauchamp, G.K.; Oka, H. Urinary Volatile Compounds as Biomarkers for Lung Cancer. Biosci. Biotechnol. Biochem. 2012, 76, 679–684. [Google Scholar] [CrossRef] [Green Version]
  38. Jiménez-Pacheco, A.; Salinero-Bachiller, M.; Iribar, M.C.; López-Luque, A.; Miján-Ortiz, J.L.; Peinado, J.M. Furan and p-xylene as candidate biomarkers for prostate cancer. Urol. Oncol. Semin. Orig. Investig. 2018, 36, 243.e21–243.e27. [Google Scholar] [CrossRef]
  39. Orywal, K.; Jelski, W.; Werel, T.; Szmitkowski, M. The activity of class I, II, III and IV alcohol dehydrogenase isoenzymes and aldehyde dehydrogenase in renal cell carcinoma. Exp. Mol. Pathol. 2015, 98, 403–406. [Google Scholar] [CrossRef]
  40. Orywal, K.; Szmitkowski, M. Alcohol dehydrogenase and aldehyde dehydrogenase in malignant neoplasms. Clin. Exp. Med. 2016, 17, 131–139. [Google Scholar] [CrossRef] [Green Version]
  41. Weiss, R.H. Metabolomics and Metabolic Reprogramming in Kidney Cancer. Semin. Nephrol. 2018, 38, 175–182. [Google Scholar] [CrossRef] [PubMed]
  42. Schaeffeler, E.; Büttner, F.; Reustle, A.; Klumpp, V.; Winter, S.; Rausch, S.; Fisel, P.; Hennenlotter, J.; Kruck, S.; Stenzl, A.; et al. Metabolic and Lipidomic Reprogramming in Renal Cell Carcinoma Subtypes Reflects Regions of Tumor Origin. Eur. Urol. Focus 2019, 5, 608–618. [Google Scholar] [CrossRef] [Green Version]
  43. Jing, L.; Guigonis, J.-M.; Borchiellini, D.; Durand, M.; Pourcher, T.; Ambrosetti, D. LC-MS based metabolomic profiling for renal cell carcinoma histologic subtypes. Sci. Rep. 2019, 9, 15635. [Google Scholar] [CrossRef] [PubMed]
  44. Hakimi, A.A.; Reznik, E.; Lee, C.-H.; Creighton, C.J.; Brannon, A.R.; Luna, A.; Aksoy, B.A.; Liu, E.M.; Shen, R.; Lee, W.; et al. An Integrated Metabolic Atlas of Clear Cell Renal Cell Carcinoma. Cancer Cell 2016, 29, 104–116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hanahan, D.; A Weinberg, R. Hallmarks of Cancer: The Next Generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  46. Brodaczewska, K.K.; Szczylik, C.; Fiedorowicz, M.; Porta, C.; Czarnecka, A.M. Choosing the right cell line for renal cell cancer research. Mol. Cancer 2016, 15, 83. [Google Scholar] [CrossRef] [Green Version]
  47. Ryan, M.J.; Johnson, G.; Kirk, J.; Fuerstenberg, S.M.; Zager, R.A.; Torok-Storb, B. HK-2: An immortalized proximal tubule epithelial cell line from normal adult human kidney. Kidney Int. 1994, 45, 48–57. [Google Scholar] [CrossRef] [Green Version]
  48. Araujo, A.M.; Moreira, N.; Lima, A.R.; Bastos, M.L.; Carvalho, F.; Carvalho, M.; De Pinho, P.G. Analysis of extracellular metabolome by HS-SPME/GC–MS: Optimization and application in a pilot study to evaluate galactosamine-induced hepatotoxicity. Toxicol. Lett. 2018, 295, 22–31. [Google Scholar] [CrossRef]
  49. Pluskal, T.; Castillo, S.; Villar-Briones, A.; Orešič, M. MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinform. 2010, 11, 395. [Google Scholar] [CrossRef] [Green Version]
  50. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Principal component analysis (PCA) scores scatter plots obtained for the headspace solid-phase microextraction (HS-SPME)/gas chromatography–mass spectrometry (GC–MS) chromatogram data of (A) volatile organic compounds (VOCs) and (B) volatile carbonyl compounds (VCCs) of culture media of HK-2 (green, n = 12), 769-P (dark blue, n = 13), 786-O (red, n = 14), Caki-1 (orange, n = 11), Caki-2 (purple, n = 12), ACHN (light blue, n = 12) and blanks (pink, n = 14).
Figure 1. Principal component analysis (PCA) scores scatter plots obtained for the headspace solid-phase microextraction (HS-SPME)/gas chromatography–mass spectrometry (GC–MS) chromatogram data of (A) volatile organic compounds (VOCs) and (B) volatile carbonyl compounds (VCCs) of culture media of HK-2 (green, n = 12), 769-P (dark blue, n = 13), 786-O (red, n = 14), Caki-1 (orange, n = 11), Caki-2 (purple, n = 12), ACHN (light blue, n = 12) and blanks (pink, n = 14).
Metabolites 10 00174 g001
Figure 2. Pairwise partial least squares-discriminant analysis (PLS-DA) scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of each tumor cell line vs. the non-tumor cell line. (a,f) 769-P (dark blue, n = 13) vs. HK-2 (green, n = 12), (b,g) 786-O (red, n = 14) vs. HK-2 (green, n = 12), (c,h) Caki-1 (yellow, n = 11) vs. HK-2 (green, n = 12), (d,i) Caki-2 (purple, n = 12) vs. HK-2 (green, n = 12) and (e,j) ACHN (light blue, n = 12) vs. HK-2 (green, n = 12).
Figure 2. Pairwise partial least squares-discriminant analysis (PLS-DA) scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of each tumor cell line vs. the non-tumor cell line. (a,f) 769-P (dark blue, n = 13) vs. HK-2 (green, n = 12), (b,g) 786-O (red, n = 14) vs. HK-2 (green, n = 12), (c,h) Caki-1 (yellow, n = 11) vs. HK-2 (green, n = 12), (d,i) Caki-2 (purple, n = 12) vs. HK-2 (green, n = 12) and (e,j) ACHN (light blue, n = 12) vs. HK-2 (green, n = 12).
Metabolites 10 00174 g002
Figure 3. Heatmaps illustrating the mean levels (normalized peak areas) of metabolites changing in ccRCC and pRCC cell lines when compared with the non-tumor cell line (HK-2). The levels of those compounds in blanks (culture medium without cells) are also represented to aid in the interpretation of the results in terms of excretion and consumption. Rows correspond to the mean normalized peak area of each metabolite, while the columns represent each sample group.
Figure 3. Heatmaps illustrating the mean levels (normalized peak areas) of metabolites changing in ccRCC and pRCC cell lines when compared with the non-tumor cell line (HK-2). The levels of those compounds in blanks (culture medium without cells) are also represented to aid in the interpretation of the results in terms of excretion and consumption. Rows correspond to the mean normalized peak area of each metabolite, while the columns represent each sample group.
Metabolites 10 00174 g003
Figure 4. Partial least squares-discriminant analysis (PLS-DA) scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of the metastatic cell lines vs. non-metastatic cell lines of each RCC histological subtype. (a,c) metastatic ccRCC (Caki-1, red, n = 11) vs. non-metastatic ccRCC (769-P and 786-O, dark blue, n = 27) and (b,d) metastatic pRCC (ACHN, green, n = 12) vs. non-metastatic pRCC (Caki-2, orange, n = 12).
Figure 4. Partial least squares-discriminant analysis (PLS-DA) scores scatter plots obtained for the HS-SPME/GC–MS chromatogram data of (A) VOCs and (B) VCCs of culture media of the metastatic cell lines vs. non-metastatic cell lines of each RCC histological subtype. (a,c) metastatic ccRCC (Caki-1, red, n = 11) vs. non-metastatic ccRCC (769-P and 786-O, dark blue, n = 27) and (b,d) metastatic pRCC (ACHN, green, n = 12) vs. non-metastatic pRCC (Caki-2, orange, n = 12).
Metabolites 10 00174 g004
Table 1. List of volatile compounds (VOCs and VCCs) significantly altered between each clear cell renal cell carcinoma (ccRCC; 769-P, 786-O and Caki-1) cell line and the non-tumoral cell line (HK-2).
Table 1. List of volatile compounds (VOCs and VCCs) significantly altered between each clear cell renal cell carcinoma (ccRCC; 769-P, 786-O and Caki-1) cell line and the non-tumoral cell line (HK-2).
Metabolite769-P (n = 13) vs. HK-2 (n = 12) 786-O (n = 14) vs. HK-2 (n = 12) Caki-1 (n = 11) vs. HK-2 (n = 12)
ES ± SEp-Value aAUC ES ± SEp-Value aAUC ES ± SEp-Value aAUC
Alcohols
2-Ethylhexanol b, L2, *↑1.42 ± 0.86 E1.26 × 10−30.885 ↑0.94 ± 0.79 E4.88 × 10−20.738 ↑2.39 ± 1.05 E6.21 × 10−40.955
2,6-Dimethyl-7-octen-2-ol b, L2 ↓−1.79 ± 0.94 C2.02 × 10−30.916
Cyclohexanol b, L2 ↓−1.29 ± 0.82 C8.06 × 10−30.821
Levomenthol b, L2↓−0.84 ± 0.794.90 × 10−20.734 ↓−1.30 ± 0.87 C2.96 × 10−20.818
Alkanes
2-Ethoxy-2-methyl-propane b, L2 ↑1.01 ± 0.80 C3.00 × 10−20.762
Decane b, L2↑0.84 ± 0.794.90 ×10−20.734 ↑1.12 ± 0.81 E2.23 × 10−20.780
Tetradecane c, L2, *↑3.68 ± 1.25 E2.90 × 10−71.000 ↑2.81 ± 1.07 E2.90 × 10−71.000 ↑1.88 ± 0.94 E5.61 × 10−40.909
Dodecane b, L2 ↓−1.29 ± 0.87 C2.96 × 10−20.811
Alkenes
2,4-Dimethyl-1-heptene b, L2↑0.80 ± 0.79 C3.00 × 10−20.769 ↑0.89 ± 0.78 C2.23 × 10−20.780
3-Carene b, L1↑2.32 ± 0.994.33 × 10−50.962 ↑1.84 ± 0.902.59 × 10−40.940
Aldehydes
4-Methylbenzaldehyde b, L2↑2.00 ± 0.94 C1.60 × 10−40.936 ↑1.69 ± 0.88 C1.10 × 10−30.899
Benzaldehyde b, L1 ↓−1.59 ± 0.86 C2.74 × 10−30.863
3-Methylbenzaldehyde c, L2↑1.75 ± 0.89 E2.34 × 10−50.946 ↑2.03 ± 0.93 E9.32 × 10−60.959
Acetaldehyde c, L1↓−10.93 ± 3.06 C2.90 × 10−71.000 ↓−4.48 ± 1.43 C2.90 × 10−71.000 ↓−4.77 ± 1.55 C5.18 × 10−60.992
Formaldehyde c, L2↑3.34 ± 1.18 C2.90 × 10−71.000 ↑2.58 ± 1.02 C2.90 × 10−71.000 ↑3.13 ± 1.18 C5.18 × 10−61.000
Decanal b, L1, *↓−1.72 ± 0.90 C1.61 × 10−30.872 ↓−1.87 ± 0.90 C3.00 × 10−30.857
Benzene Derivatives
Ethylbenzene b, L2↑1.46 ± 0.863.58 × 10−40.917 ↑1.00 ± 0.79 E2.42 × 10−30.875
Styrene b, L2, *↑1.32 ± 0.84 E1.61 × 10−30.871
Xylene b, L2↑1.61 ± 0.88 E4.33 × 10−50.962 ↑1.11 ± 0.81 E6.27 × 10−40.917
Ketones
Acetophenone b, L2↓− 0.96 ± 0.80C4.31 × 10−20.750
Cyclohexanone b, L2, *↓− 3.29 ± 1.19 C2.88 × 10−61.000 ↓−1.53 ± 0.85 C2.74 × 10−30.863 ↓−1.54 ± 0.91 C6.04 × 10−30.879
Acetone c, L2↑2.12 ± 0.94 E5.63 × 10−50.929 ↑2.22 ± 0.96 E7.25 × 10−60.964 ↑1.33 ± 0.86 E6.56 × 10−30.833
Unknowns
Un (RT 8.40, m/z 69) b, L4↑1.44 ± 0.86 C1.43 × 10−20.801 ↑1.53 ± 0.85 C3.31 × 10−30.851
Un (RT 10.18, m/z 58) c, L4 ↑3.64 ± 1.24 E2.90 × 10−71.000 ↑3.07 ± 1.12 E2.90 × 10−71.000 ↑2.04 ± 0.97 E5.18 × 10−60.992
Un (RT 12.82. m/z 69) b, L4 ↑0.98 ± 0.792.40 × 10−20.774
Un (RT 16.64, m/z 61) b, L4↓− 4.16 ± 1.382.88 × 10−61.000 ↓−4.39 ± 1.41 E3.73 × 10−61.000
ES: Effect size, SE: standard error, AUC: area under the curve. ES and SE were determined as described in reference Berben et al. (2012) [29]. a False discovery rate adjusted p-value; b, c Compounds detected through VOCs and VCCs analytical methods, respectively. ↑ and ↓: volatile compound found increased and decreased, respectively, in the extracellular medium of RCC compared with the non-tumoral cell line. E, C Metabolites excreted and consumed, respectively, when compared with blanks. No indication was added when the levels of the metabolites in tumoral cell lines were similar to that of blanks. L1 Identified metabolites (GC–MS analysis of the metabolite of interest and a chemical reference standard under identical analytical conditions within the same laboratory).L2 Putatively annotated compounds (spectral MS similarity with the NIST database).L4 Unknown compounds. * Metabolites previously reported as potential biomarkers of RCC in urine [27].
Table 2. List of volatile compounds (VOCs and VCCs) significantly altered between papillary renal cell carcinoma (pRCC; Caki-2 and ACHN) cell line and the non-tumoral cell line (HK-2).
Table 2. List of volatile compounds (VOCs and VCCs) significantly altered between papillary renal cell carcinoma (pRCC; Caki-2 and ACHN) cell line and the non-tumoral cell line (HK-2).
MetaboliteCaki-2 (n = 12) vs. HK-2 (n = 12) ACHN (n = 12) vs. HK-2 (n = 12)
ES ± SEp-Value aAUC ES ± SEp-Value aAUC
Alcohols
2-Ethylhexanol b, L2, *↑2.51 ± 1.05 E2.22 × 10−50.972 ↑4.29 ± 1.44 E7.40 × 10−61.000
Cyclohexanol b, L2↑1.45 ± 0.87 E2.87 × 10−30.868
Alkanes
Tetradecane c, L2, *↑3.15 ± 1.18 E8.63 × 10−71.000 ↑2.50 ± 1.05 E1.73 × 10−61.000
Aldehydes
4-Methylbenzaldehyde b, L2 ↓−1.84 ± 0.93 C4.78 × 10−30.875
3-Methylbenzaldehyde c, L2↑0.81 ± 0.81 C1.21 × 10−20.799
Acetaldehyde c, L1↓−9.59 ± 2.82 C8.63 × 10−71.000 ↓−10.95 ± 3.19 C1.73 × 10−61.000
Formaldehyde c, L2↑4.39 ± 1.46C8.63 × 10−71.000 ↑2.34 ± 1.02 C7.25 × 10−60.979
Ketones
2-Pentadecanone b, L1↓−2.64 ± 1.079.86 × 10−60.986
Acetophenone b, L2↓−4.85 ± 1.57 C3.70 × 10−61.000
Cyclohexanone b, L2, *↓−7.09 ± 2.15 C3.70 × 10−61.000 ↓−1.40 ± 0.87 C4.78 × 10−30.868
Acetone c, L2↑4.63 ± 1.52 E8.63 × 10−61.000 ↑2.89 ± 1.13 E1.04 × 10−50.972
Unknowns
Un (RT 9.80, m/z 59) b, L4↑1.55 ± 0.89 C3.05 × 10−30.861
Un (RT 10.18, m/z 58) c, L4 ↑4.65 ± 1.52 E8.63 × 10−71.000 ↑3.38 ± 1.23 E1.73 × 10−61.000
Un (RT 12.82. m/z 69) b, L4↓−0.99 ± 0.82 C2.47 × 10−60.785
Un (RT 16.64, m/z 61) b, L4 ↓−1.23 ± 0.85 E1.13 × 10−20.833
ES: effect size, SE: standard error, AUC: area under the curve. ES and SE were determined as described in reference Berben et al. (2012) [29]. a False discovery rate adjusted p-value; b, c Compounds detected through VOCs and VCCs analytical methods, respectively. ↑ and ↓: volatile compound found increased and decreased, respectively, in the extracellular medium of RCC compared with non-tumoral cell line. E, C Metabolites excreted and consumed, respectively, when compared with blanks. No indication was added when the levels of the metabolites in tumoral cell lines were similar to that of blanks. L1 Identified metabolites (GC–MS analysis of the metabolite of interest and a chemical reference standard under identical analytical conditions within the same laboratory). L2 Putatively annotated compounds (spectral MS similarity with the NIST database). L4 Unknown compounds. * Metabolites previously reported as potential biomarkers of RCC in urine [27].
Table 3. List of volatile compounds (VOCs and VCCs) significantly altered between metastatic cell lines compared with non-metastatic cell lines and between RCC histological subtype (ccRCC and pRCC).
Table 3. List of volatile compounds (VOCs and VCCs) significantly altered between metastatic cell lines compared with non-metastatic cell lines and between RCC histological subtype (ccRCC and pRCC).
MetaboliteMetastatic ccRCC (n = 11) vs.
Non-Metastatic ccRCC (n = 27)
Metastatic pRCC (n = 12) vs.
Non-Metastatic pRCC (n = 12)
ccRCC (n = 38) vs. pRCC (n = 24)
ES ± SEp-Value aAUC ES ± SEp-Value aAUC ES ± SEp-Value aAUC
Alcohols
2-Ethylhexanol b, L2 ↑1.20 ± 0.84 E1.19 × 10−20.828 ↓−0.80 ± 0.52 E1.88 × 10−30.730
Cyclohexanol b, L2↑1.83 ± 0.80 E1.71 × 10−40.882 ↓−2.61 ± 1.077.40 × 10−60.993 ↑1.15 ± 0.54 C2.04 × 10−50.813
Alkanes
2-Ethoxy-2-methyl-propane b, L2↓−0.99 ± 0.72 C7.42 × 10−30.781 ↑0.95 ± 0.531.13 × 10−40.763
Decane b, L2↓−2.54 ± 0.89 C1.05 × 10−70.989 ↓−0.96 ± 0.82 C1.67 × 10−20.806 ↓−0.98 ± 0.532.89 × 10−40.728
Dodecane b, L2↓−2.33 ± 0.861.35 × 10−60.959 ↓−1.28 ± 0.851.81 × 10−20.799 ↑1.03 ± 0.549.68 × 10−50.782
4-Methylheptane b, L2↓−1.40 ± 0.75 C2.60 × 10−40.872 ↑1.28 ± 0.55 C1.31 × 10−60.839
Tetradecane c, L2 ↓−1.64 ± 0.90 E5.78 × 10−40.896
Alkenes
2,4-Dimethyl-1-heptene b, L2↓−1.60 ± 0.77 C5.02 × 10−50.906 ↑1.29 ± 0.55 C1.31 × 10−60.842
3-Carene b, L1↓−2.29 ± 0.86 C2.99 × 10−81.000 ↑1.36 ± 0.566.99 × 10−70.850
Aldehydes
4-Methylbenzaldehyde b, L2↓−2.26 ±0.85 C1.89 × 10−70.983 ↓−1.23 ± 0.85 C2.79 × 10−20.778 ↑1.28 ± 0.55 C2.22 × 10−60.827
3-Methylbenzaldehyde c, L2 ↓−0.97 ± 0.82 C3.84 × 10−40.910
Acetaldehyde c, L1 ↓−1.90 ± 0.94 C3.84 × 10−40.917
Formaldehyde c, L2 ↓−1.29 ± 0.85 C5.56 × 10−30.826
Decanal b, L1↑2.58 ± 0.90 E7.25 × 10−70.969 ↓−1.07 ± 0.83 E1.13 ×10−20.833 ↑0.54 ± 0.51 E3.29 × 10−20.654
Benzene Derivatives
Styrene b, L2↓−1.46 ± 0.766.20 × 10−40.852
Xylene b, L2↓−0.67 ± 0.70 E3.38 × 10−30.808
Ethylbenzene b, L2↓−0.51± 0.70 C7.12 × 10−30.784
Ketones
2-Pentadecanone b, L1 ↑2.80 ± 1.11 E1.11 × 10−50.986 ↑0.84 ± 0.53 E8.89 × 10−40.734
4-Methyl-2-hexanone b, L2 ↑0.97 ± 0.539.68 × 10−50.768
Acetone c, L2 ↓−1.97 ± 0.95 E3.84 × 10−40.910
Acetophenone b, L2 ↑5.83 ± 1.82 C5.55 × 10−61.000 ↓−1.40 ± 0.56 C2.03 × 10−40.761
Cyclohexanone b, L2 ↑8.14 ± 2.44 C7.40 × 10−61.000
Unknowns
Un (RT 8.40, m/z 69) b, L4↓−2.02 ± 0.82 C2.86 × 10−60.949 ↓−1.31 ± 0.55 E6.89 × 10−60.795
Un (RT 9.45 m/z 71) b, L4 ↑1.79 ± 0.921.56 × 10−20.813
Un (RT 9.80, m/z 59) b, L4 ↓−1.79 ± 0.924.33 × 10−40.924
Un (RT 10.18, m/z 58) c, L4 ↓−1.69 ± 0.91 E3.84 × 10−4 0.917
Un (RT 12.82. m/z 69) b, L4↑1.70 ± 0.79 C2.00 × 10−50.918
Un (RT 16.64, m/z 61) b, L4↑3.58 ± 1.06 E8.00 × 10−70.966 ↑1.40 ± 0.56 C6.99 × 10−70.844
ES: effect size, SE: standard error, AUC: area under the curve. ES and SE determined as described in reference Berben et al. (2012) [29]. a False discovery rate adjusted p-value; b, c Compounds detected through VOCs and VCCs analytical methods, respectively. ↑ and ↓: Volatile compound found increased and decreased, respectively, in the extracellular medium of metastatic cell lines compared with non-metastatic cell lines and between RCC histological subtype (ccRCC and pRCC). E, C Metabolites excreted and consumed, respectively, when compared with blanks. No indication was added when the levels of the metabolites in tumoral cell lines were similar to that of blanks. L1 Identified metabolites (GC–MS analysis of the metabolite of interest and a chemical reference standard under identical analytical conditions within the same laboratory). L2 Putatively annotated compounds (spectral MS similarity with the NIST database). L4 Unknown compounds.

Share and Cite

MDPI and ACS Style

Amaro, F.; Pinto, J.; Rocha, S.; Araújo, A.M.; Miranda-Gonçalves, V.; Jerónimo, C.; Henrique, R.; Bastos, M.d.L.; Carvalho, M.; Guedes de Pinho, P. Volatilomics Reveals Potential Biomarkers for Identification of Renal Cell Carcinoma: An In Vitro Approach. Metabolites 2020, 10, 174. https://doi.org/10.3390/metabo10050174

AMA Style

Amaro F, Pinto J, Rocha S, Araújo AM, Miranda-Gonçalves V, Jerónimo C, Henrique R, Bastos MdL, Carvalho M, Guedes de Pinho P. Volatilomics Reveals Potential Biomarkers for Identification of Renal Cell Carcinoma: An In Vitro Approach. Metabolites. 2020; 10(5):174. https://doi.org/10.3390/metabo10050174

Chicago/Turabian Style

Amaro, Filipa, Joana Pinto, Sílvia Rocha, Ana Margarida Araújo, Vera Miranda-Gonçalves, Carmen Jerónimo, Rui Henrique, Maria de Lourdes Bastos, Márcia Carvalho, and Paula Guedes de Pinho. 2020. "Volatilomics Reveals Potential Biomarkers for Identification of Renal Cell Carcinoma: An In Vitro Approach" Metabolites 10, no. 5: 174. https://doi.org/10.3390/metabo10050174

APA Style

Amaro, F., Pinto, J., Rocha, S., Araújo, A. M., Miranda-Gonçalves, V., Jerónimo, C., Henrique, R., Bastos, M. d. L., Carvalho, M., & Guedes de Pinho, P. (2020). Volatilomics Reveals Potential Biomarkers for Identification of Renal Cell Carcinoma: An In Vitro Approach. Metabolites, 10(5), 174. https://doi.org/10.3390/metabo10050174

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