Metabolomic and Lipidomic Profiling Identifies The Role of the RNA Editing Pathway in Endometrial Carcinogenesis
grid.7080.fBiomedical Research Group in Gynecology, Vall Hebron Research Institute (VHIR), Universitat Autònoma de Barcelona, CIBERONC Barcelona, Spain
0000 0001 2186 0438grid.411667.3Department of Oncology, Georgetown University Medical Center, Washington D.C., USA
0000 0001 0675 8654grid.411083.fPathology Department, Vall Hebron University Hospital, Barcelona, Spain
Pathological Oncology Group and Pathology Department, University Hospital Arnau de Vilanova, and University Hospital Bellvitge, IRBLLEIDA and Idibell, University of Lleida, CIBERONC Lleida, Spain
0000 0001 0675 8654grid.411083.fGynecological Oncology Department, Vall Hebron University Hospital, Barcelona, Spain
0000 0001 2325 3084grid.410675.1Basic Sciences Department, International University of Catalonia, CIBERONC Barcelona, Spain
0000 0001 1955 1644grid.213910.8Department of Biochemistry Molecular and Cellular Biology, Georgetown-Lombardi Comprehensive Cancer Center, Washington, D.C., United States
Abstract
Endometrial cancer (EC) remains the most common malignancy of the genital tract among women in developed countries. Although much research has been performed at genomic, transcriptomic and proteomic level, there is still a significant gap in the metabolomic studies of EC. In order to gain insights into altered metabolic pathways in the onset and progression of EC carcinogenesis, we used high resolution mass spectrometry to characterize the metabolomic and lipidomic profile of 39 human EC and 17 healthy endometrial tissue samples. Several pathways including lipids, Kynurenine pathway, endocannabinoids signaling pathway and the RNA editing pathway were found to be dysregulated in EC. The dysregulation of the RNA editing pathway was further investigated in an independent set of 183 human EC tissues and matched controls, using orthogonal approaches. We found that ADAR2 is overexpressed in EC and that the increase in expression positively correlates with the aggressiveness of the tumor. Furthermore, silencing of ADAR2 in three EC cell lines resulted in a decreased proliferation rate, increased apoptosis, and reduced migration capabilities in vitro. Taken together, our results suggest that ADAR2 functions as an oncogene in endometrial carcinogenesis and could be a potential target for improving EC treatment strategies.
Introduction
Endometrial cancer (EC) accounts for 7% of the new cases of female cancers in United States in 2017 and its incidence is increasing1, 2. The most standardized classification divides EC in two different subtypes: type I or endometrioid carcinomas (EEC), which is the most frequent subtype; and type II or non-endometrioid carcinomas (NEEC), which are more aggressive tumors. Among NEEC subtypes, serous is the most prominent histology3, 4. Moreover, EC tumors are classified according to the extent of tumor dissemination (International Federation of Gynecology and Obstetrics or FIGO staging) and histological grade1. About 20% of patients are diagnosed at an advanced stage and/or at a high histological tumor grade and have a low 5-year survival rate associated5, 6. As such, availability of biomarkers for disease stratification and a thorough understanding of biochemical perturbations that underscore EC progression are likely to improve clinical outcomes.
Extensive research performed using high-throughput technologies including genomics, transcriptomics and proteomics has augmented characterization of molecular changes at different levels of cellular expression that are specifically associated with human endometrium7 and onset and progression of EC8, 9. Metabolomics defines the end point of cellular processes and hence provides a readout of current physiological status of the system10. Thus, metabolomics, lipidomics and glycomics have emerged as promising tools for clinical and translational research11. Rapid advancement of metabolomics technologies such as ultra-performance liquid chromatography mass spectrometry (UPLC-MS), enable comprehensive interrogation of the human metabolome and lipidome12–15. Bioinformatics analyses of these data is likely to augment the discovery of new clinical and pharmacological targets16.
It is known that changes at the transcriptomic levels generate a new source of complexity that promotes initiation and progression of cancer and other diseases17–19. The most common RNA editing events are mediated post-transcriptionally by the Adenosine deaminases acting on RNA (ADAR) family of enzymes20. The ADAR gene family catalyzes the deamination of adenosine that is converted to an inosine creating a dysregulation of the adenosine/inosine (A/I) ratio in the cell. Adenosine to inosine (A-to-I) editing can lead to amino acid recoding events21 since inosine is recognized as a guanosine. This family includes three enzymes, ADAR1 (UniProtKB P55265), ADAR2 (UniProtKB P78563) and ADAR3 (UniProtKB Q9NS39); also known as ADAR, ADARB1 and ADARB2, respectively22. ADAR1 and ADAR2 are ubiquitously expressed and, differently that the brain specific ADAR3, they show catalytic activity23.
The overall goal of this study was to identify altered metabolic pathways in EC. Characterization of the metabolome and lipidome of EC tumors was performed using a high resolution mass spectrometry approach in conjunction with UPLC-MS. Pathway validation was performed with an independent set of samples that for the first time, yielded insights into dysregulation of the RNA editing pathway in endometrial tumors. Alterations of the RNA editing pathway correlated with high histological grade and EC serous subtypes that are indicators of poor prognosis in EC. Finally, the expression of ADAR editing enzymes was modulated in vitro to confirm an important oncogenic role of this pathway in EC cell proliferation, apoptosis and migration. These results are instructive of the role of this pathway in EC tumor progression.
Results
Increased expression of ADAR family of enzymes in human EC tumors
Among the metabolites identified in EC tissues, we were particularly interested in the dysregulation of nucleoside inosine since the relative abundance of this metabolite was significantly higher in EC tumors. To our knowledge, the underlying impact of alterations of endogenous levels of inosine has never been investigated in EC. Dysregulated levels of adenosine and inosine (A/I ratio) can be attributed to the modulation of the A-to-I editing pathway24. Consequently, we studied the status of this pathway in EC by analyzing the expression level of members of the A-to-I editing enzymes family (ADAR1 and ADAR2) in three independent sets of samples by immunohistochemistry (IHQ) (Table 4). The first set included the evaluation of 20 EEC samples and their corresponding paired healthy tissues; the second set included 36 EEC tumors diagnosed at different histological grades (low grade, n = 12; intermediate grade, n = 13; high grade, n = 11); and the third set allowed to differentiate between histological subtypes, and so, 78 EEC and 29 NEEC tissue samples were included.
Cohort 1 Cohort 2 Cohort 3 EEC (n = 20) NEEC (n = 0) Control (n = 20) EEC (n = 36) NEEC (n = 0) Control (n = 0) EEC (n = 78) NEEC (n = 29) Control (n = 0) Age > 50 20 — 20 36 — — 78 29 — < 50 — — — —— — — — — — Collection center VHUH 20 — 20 7 — — 78 29 — Lleida — — — 13 — — — — — Santiago — — — 6 — — — — — Virgen Rocío — — — 8 — — — — — MD Anderson — — — 2 — — — — — Uterine condition Pre-menopausal — — — — — — — — — Post-menopausal 20 — 20 36 — — 78 29 — FIGO stage IA 10 — — 13 — — 17 7 — IB 6 — — 10 — — 31 5 — II 2 — — 5 — — 20 4 — III 2 — — 8 — — 7 11 — IV — — — — — — 3 2 — Histologic grade G1 4 — — 12 — — 13 0 G2 7 — — 13 — — 33 2 G3 9 — — 11 — — 32 27 Format TMA — — — 36 — — 78 29 — Individual slide 20 — 20 — — — — — —
Although ADAR1 and ADAR2 staining was clearly localized in the nucleus of epithelial, stromal, and endothelial cells, we specifically analyzed the staining of the epithelial tumor cells (Fig. 2B,D,F and Supplementary Figure 4B). Interestingly, no staining was observed in cells undergoing mitosis (Supplementary Figure 4A). Our results demonstrated that ADAR1 and ADAR2 were both significantly increased in tumor samples compared to healthy endometrial tissue (Fig. 2A), confirming that the RNA editing pathway is activated in EC carcinogenesis. ROC analysis for ADAR1 and ADAR2 expression yielded an AUC of 0.79 and 0.90 respectively, emphasizing the differences observed between control and EC samples (Supplementary Figure 5). Moreover, ADAR expression increased progressively with the presence of poor prognostic factors, such as tumors presenting high grade or a NEEC histology (Fig. 2C and E). These data suggested for the first time an activation of the RNA editing pathway in patients with EC, which positively correlates with aggressive disease, resulting in a dysregulation of the nucleosides/nucleotides balance in the tumor tissue.
Inhibition of ADAR2 reduces viability, increases apoptosis and reduces migration capabilities in human EC cell lines
Next, in order to determine the possible role of the RNA editing pathway in EC progression, we transiently silenced the expression of the ADAR1 and ADAR2 enzymes in three EC cell lines (HEC-1A, Ishikawa and RL95-2). Transfection efficiency was determined by comparing each transfection against the corresponding negative control (NC) (Supplementary Figure 6). Gene silencing for both enzymes was confirmed by western blot analysis (WB) and immunofluorescence (IF) after 96 h of transfection (Supplementary Figure 7), time point in which the modulation of cell growth and cell viability in EC cell lines was assessed. The inhibition of ADAR2, but not ADAR1, resulted in a significant reduction of cell viability and proliferation compared to controls in the three cell lines used in this study (Fig. 3A). Furthermore, the ratio of apoptotic cells significantly increased in the three cell lines upon silencing of ADAR2. Similarly, proliferation assay performed using cells treated with siRNA-ADAR1, did not result in significant changes in apoptosis (Fig. 3B). Finally, the migration capabilities of the EC cell lines were interrogated after inhibiting ADAR1 and ADAR2 in a wound healing assay. Knockdown of ADAR2 expression resulted in significant reduction of the migration rate in HEC-1A and RL95-2 cell lines. Remarkably, similar changes were not observed when the same cell lines were treated with siRNA-ADAR1 (Fig. 3C).
In order to further confirm the functional role of ADAR2 in EC cell lines, we also conducted functional assays silencing ADAR2 expression with two new and different siRNAs (siRNA-ADAR2_B and siRNA-ADAR2_C), independent from the siRNA-ADAR2 used in Fig. 3. We corroborate a significant decrease of cell viability, a significant increase in apoptosis rate and a significant decrease of wound healing rate of the 3 EC cell lines induced by the silencing of ADAR2 expression (Supplementary Figure 8). These results clearly demonstrate that changes in the expression of ADAR2 leads to significant changes in the aggressive behavior of EC cell lines, specifically on cell viability and apoptosis, and cell migration capabilities of EC cell lines.
Discussion
Recently, few studies have reported alterations in the metabolomic or proteomic phenotype of EC in serum25, urine26 or in other sample types27, 28, underscoring the clinical translational relevance of these technologies in furthering the personalized medicine initiative. In this study, we used a global metabolomics profiling approach in order to understand the metabolic changes that take place in EC carcinogenesis and during tumor progression (Fig. 4).
Our study reveals an array of metabolites dysregulated in EC, some of which have been previously described in endometrial carcinogenesis. Glycerophospholipid class of metabolites was found to be upregulated in tumor tissues, including PCs, PEs, and PIs. Lipid biosynthesis and catabolism is known to be altered in several diseases, including cancer28–32. Consistent with our findings, Trousil et al.28 also observed an increase (up to 70%) in PC levels in EC tissues. We also observed a downregulation of the acylamido analogs of endocannabinoids such as palmitamide, stearamide and oleamide in EC. This downregulation has been previously reported33. We also observed a decreased proportion of picolinic acid in tumor samples. Picolinic acid and quinolinic acid are the end products of the Kynurenine pathway and have been shown to have anti-tumoral and pro-tumoral activity respectively34. Moreover, the activity of one of the main enzymes of the pathway, indoleamine 2,3-dioxygenase (IDO) has also been studied in EC35. Decreased levels of Glu Phe Arg Trp and inosine as well as upregulation of 3-Deoxyvitamin D3 and UDP-N-acetyl-D-galactosamine were also found in our EC tumor sample set compared to control tissues. We also analyzed metabolomic changes that underscore tumor progression. Our data suggest significant changes in the lipidome including PC (16:0/20:5), PC (16:0/22:6), PE (16:0/22:6), PE (22:6/P-18:1) and PE (18:1/22:6) at different stages of cancer progression, as well as a significant increase in the levels of UDP-N-acetyl-D-galactosamine and arachidonic acid at advanced stages of EC. Further studies would be needed to fully understand the scope and impact of these alterations in cancer progression.
The dysregulation of inosine was further studied to dissect the functional implications of the RNA editing pathway in EC. The dysregulation of inosine is indicative of a possible imbalance in the I/A ratio, which was also reported by Trousil et al.28 in studies with EC tissue compared to normal endometrium. The A-to-I conversion is the most common type of RNA editing found in mammals mediated by the ADAR enzymes. Although, to date, the RNA editing pathway and the expression and function of the ADAR gene family has not been interrogated in EC, it has been reported that the expression of ADAR enzymes is upregulated in many cancers18, 19, including breast17 and esophageal squamous cell carcinoma36, 37. The ADAR family is comprised of three members: ADAR1 and ADAR2, that are present in most human tissues; and ADAR3, that is brain specific23. Changes in editing frequencies have been described in other diseases including prostate, lung, kidney and testis tumors while reduced RNA levels of ADAR1, ADAR2 and ADAR3 have been observed in brain tumors38, 39. ADAR enzymes are also involved in physiological events such neuronal development, immune response, cell response to viruses and regulation of miRNA expression among others40, 41.
Hence, we asked if the expression of ADAR enzymes had any correlation with EC initiation and progression. Our findings not only elucidate ADAR1 and ADAR2 enzymes to be significantly upregulated in EC tumor tissues compared to healthy endometrium, but also demonstrate a significant correlation between ADARs expression and the malignancy of the tumor. We found that the ADARs expression increased progressively with tumor grade. More importantly, our data showed a significant increase in the expression of ADAR1 and ADAR2 in the most aggressive subtype of EC, the NEEC, that have the worst predicted survival3.
The role of the RNA editing pathway in EC was further investigated by knocking-down the expression of ADAR1 and ADAR2 in HEC-1A, RL95-2 and Ishikawa EC lines. Our results demonstrate the impact of decreased ADAR2 expression on an array of cellular functions in EC cell lines including a significant decrease of cell proliferation and viability, increased apoptosis rate, and reduced migration capabilities in vitro. Similar to our observations, silencing of ADARs in breast cancer cell lines led to less cell proliferation and more apoptosis17 and overexpression of editing enzymes accelerated growth rate and colony formation in esophageal squamous cell carcinoma in vitro 36. Taken together, our results strongly suggest, for the first time, that the RNA editing gene family, specifically ADAR2, may play an important role promoting EC carcinogenesis.
In conclusion, a global molecular profiling approach using high resolution mass spectrometry has been useful, not only to describe changes in the metabolome and lipidome in human endometrial carcinogenesis and EC progression but also led to the discovery of an important alteration of the RNA editing pathway in EC. We further demonstrated the role of this pathway in proliferation and viability, apoptosis, and migration of EC cells, leading us to conclude that the activation of the RNA editing pathway is an oncogenic process in EC. This study opens several avenues for further investigations of ADAR2 as possible target for the development of therapeutic approaches for the treatment of EC patients.
Methods
Patients and tissue collection
Discovery and verification phase
A total of 56 women (39 diagnosed with EC and 17 non-EC patients) were recruited at Vall Hebron University Hospital (VHUH) in Barcelona, Spain. All patients participating signed an informed consent and the study was approved by the Clinical Research Ethics Committee (CREC) at the VHUH (approval number: PR_AMI_50–2012). Endometrial tissues were collected at the histopathology department of VHUH after the surgical intervention. Each tissue piece (about 5 to 10 mg) was placed in a sterile and separate container, properly labeled and stored at −80 °C immediately. A description of the clinical and pathological characteristics of the tumors is detailed in Table 1. Inclusion criteria: post-menopausal women, > 50 years, no previous treatment for pelvic gynecological cancer, negative for HIV and hepatitis viruses.
Validation phase
Patients split in 3 new cohorts were enrolled in VHUH or in University Hospital Arnau de Vilanova of Lleida following the approval of the CREC at each participating institution. Tissue samples were embedded in paraffin block for individual slides or tissue microarray (TMA) construction in VHUH. A description of the clinical and pathological characteristics of the tissues is detailed in Table 4.
Immunohistochemistry
A total of three TMAs were constructed as described previously47. Paraffin-embedded TMA or individual sections of the samples were mounted on slides, deparaffined and rehydrateted. Antigen retrieval was done during 4 min at 115 °C (pH 6) and subsequently blocked with peroxidase (3%) for 5 min and incubated with the primary antibody for 1 h and 30 min at room temperature in a humidified chamber. After blocking for 30 min with EnVision secondary antibody (Ref. K5007-100 ml, Agilent; Santa Clara, CA, USA) slides were treated with DAB (Ref. K3468, Agilent; Santa Clara, CA, USA) reagent and counterstained with hematoxylin for 30 sec. Finally, a dehytratation step and DPX mounting were performed. Pictures were taken using the Olympus BX41 microscope (20X). The following antibodies were tested: ADAR1 (Ref. AMAb90535, dilution 1:100, Sigma-Aldrich; St. Louis, MO, USA) and ADAR2 (Ref. sc-73409, dilution 1:50, Santa Cruz; Heidelberg, Germany, EU). A pathologist evaluated the expression using two criteria: the intensity of staining (0, no staining; 1, weak intensity; 2, moderate intensity; 3, high intensity) and the percentage of endometrial epithelial stained cells (0–100). The product of the two scores yielded final values on a sacale ranging from 0 to 300.
Cell lines
The following EC epithelial cell lines were used in this study. HEC-1A (Ref. HTB112, ATCC; Manassas, VA, USA) cell line was cultured in McCoys 5 A medium (Ref. 80014020, ThermoFisher; Waltham, MA, USA); Ishikawa (Ref. 99040201, Sigma-Aldrich; St. Louis, MO, USA) and RL95-2 (CRL-1671, ATCC; Manassas, VA, USA) were cultured in DMEM/F12 (Ref. 11320-033, ThermoFisher; Waltham, MA, USA). Both supplemented with penicillin/streptomycin (1%) and fetal bovine serum (10%). Incubated at 37 °C in a 5% CO2 humidified chamber. Cell lines were morphologically and genetically authenticated and tested for mycoplasma in accordance with AACR guides.
Western Blot
Protein extraction from cell lines was performed using RIPA buffer (5 nM EDTA, 150 mM NaCl, 1% Triton, 20 nM Tris pH 8 and 1:200 protein inhibitors). After Bradford colorimetric protein quantification samples were run on a SDS/PAGE acrylamide gel for 2 h and transferred to a PVDF membrane. Membranes were blocked for 1 h with 5% milk and incubated with the primary antibodies described above overnight and with the secondary antibodies -Goat anti-rabbit (Ref. P0448, dilution 1:2000, Agilent; Santa Clara, CA, USA) and rabbit anti-mouse (Ref. P0260, dilution 1:2000, Agilent; Santa Clara, CA, USA)- for 1 h at room temperature. Membranes were developed using Immobilon Immobilon Western Chemiluminiscent (Ref. WBKLS0100, Merck Millipore; Billerica, MA, USA). Membranes were finally stained with naphtol blue (Ref. N3393, Sigma-Aldrich; St. Louis, MO, USA) to detect all proteins transferred to the membranes in order to normalize the ADAR protein bands. The following antibodies were tested: ADAR1 (Ref. AMAb90535, dilution 1:100, Sigma-Aldrich; St. Louis, MO, USA) and ADAR2 (Ref. sc-73409, dilution 1:50, Santa Cruz; Heidelberg, Germany, EU).
Immunofluorescence
Cells were fixed for 10 min 96 h after transfection. Cells were blocked for 15 min with 5% BSA solution in 10 ml of 0.5% PBS-T. They were incubated with the primary antibody for 2 h at room temperature and with the secondary antibody for 45 min at room temperature in a dark chamber. Mounting and DAPI staining were performed with ProLong Gold Antifade reagent with DAPI (Ref. P36931, ThermoFisher; Waltham, MA, USA). Pictures were taken under the fluorescence microscope Nikon Eclipse TE2000-S. Primary antibodies: ADAR1 (Ref. AMAb90535, dilution 1:100, Sigma-Aldrich; St. Louis, MO, USA) and ADAR2 (Ref. sc-73409, dilution 1:50, Santa Cruz; Heidelberg, Germany, EU). Secondary antibody: Alexa Fluor 488 (Ref. A11017, Labeling and detection, ThermoFisher; Waltham, MA, USA).
Reverse Transfection of siRNAs
Desired number of cells according to the assay (see proliferation, apoptosis and wound healing assay sections for the exact number of cells) was seeded using media without antibiotic supplementation. Transfection mix was prepared in the following proportion: 25 µl of OptiMEM (Ref. 11058021, Invitrogen-ThermoFisher; Waltham, MA, USA), 0.2 µl of lipofectamine (Ref. 11668019, Invitrogen-ThermoFisher; Waltham, MA, USA), and 0.375 µl of siRNA (50 nM) and added in a 1:3 (v/v) to final volume of the well. After 24 h of transfection, media was changed. Pre-designed siRNAs for ADARs were purchased from Sigma-Aldrich: siRNA-ADAR1: SASI Hs01 00115047- GACUAUCUCUUCAAUGUGU; siRNA-ADAR2: SASI Hs01 00237747-GAGUGAUCGUGGCCUUGCA; siRNA-ADAR2_B: SASI Hs01 00222314- GAGUGAUCGUGGCCUUGCA; siRNA-ADAR2_C: SASI_Hs01_00135204- GAGUGAUCGUGGCCUUGCA; siRNA-NC: negative control, BLOCK-it Fluorescent Oligo, for lipid transfection (Ref. 2013, ThermoFisher; Waltham, MA, USA).
Cell proliferation assay
In order to determine cell proliferation and viability, a total of 4 × 103 cells (HEC-1A); 12 × 103 cells (Ishikawa), and 2 × 104 cells (RL95-2) per well were plated in a 96 well plate (6 replicates per condition). The cell growth rate was measured after 96 h of transfection. Cells were fixed to the plate with glutaraldehide 1% (Ref. G6257, Sigma-Aldrich; St. Louis, MO, USA) for 15 min and stained with Crystal Violet (Ref. C3886, Sigma-Aldrich; St. Louis, MO, USA). After 20 min, cells were washed with H2O and Acetic Acid 15% (Ref. 0641, ThermoFisher; Waltham, MA, USA) was added. After 10 min shaking plates, they were read at 590 nm. Three independent experiments were carried out.
Apoptosis assay
In order to determine cell apoptosis rate, a total of 8 × 104 cells (HEC-1A); 5 × 104 cells (Ishikawa), and 2 × 105 cells (RL95-2) per well were seed in a p24 (8 replicates per condition). After 96 h of transfection, cells were stained with Hoechst (1:1000, Ref. H6024, Sigma-Aldrich; St. Louis, MO, USA) and after 24 h pictures were taken under the fluorescence microscope Nikon Eclipse TE2000-S. The number of nucleus with condensed chromatin and strongly staining of the DNA were quantified as apoptotic cells and divided by the total number of nucleus. Three independent experiments were carried out (4 camps per well were counted).
Wound healing assay
To evaluate cell migration capabilities, a total of 8 × 104 cells (HEC-1A); 5 × 104 cells (Ishikawa), and 2 × 105 cells (RL95-2) per well were seed in a p24 (3 replicates per condition). Wound was generated 72 h post-transfection. Pictures were taken every 8 h using an Olympus FSX100 microscope. Wound healing area was measured using Image J software at the different time points. Three experiments were carried out independently.
All experiments were performed in accordance with the relevant guidelines and regulations.
Statistical analysis
The SPSS statistical package version 23 for Windows® and the GraphPad Prism version 6 were used to perform the statistical analyses and ROC analysis. Each value represents the mean of at least 3 replicates with the corresponding standard deviation. We analyzed the normality of each data set and we used, according to the sample distribution, parametric or no parametric tests. For the IHQ analysis, a Wilcoxon signed rank test was used to compare tumors from controls; and a Kruskal-Wallis and Mann-Whitney tests were applied when comparing protein expression according to grades and histological subtypes, respectively. For the functional analysis, means of the different groups were compared by Kruskal-Wallis followed by Dunn’s multiple comparisons test (in case of significance). We considered significant p-values < 0.05. (*p-v < 0.05; **p-v < 0.01; ***p-v < 0.001; ****p-v: 0.0001).
Electronic supplementary material
Untitled section
Acknowledgements
This work was supported by the Spanish Ministry of Health (RD12/0036/0035), the Spanish Ministry of Economy and Competitivy (PI14/02043), the AECC (Grupos Estables de Investigacion 2011 - AECC- GCB 110333 REVE), the Fundació La Marató TV3 (2/C/2013), the CIRIT Generalitat de Catalunya (2014 SGR 1330) and the European Commission, 7th Framework Program, IRSES (PROTBIOFLUID –269285) – Belgium. The Spanish Ministry of Economy and Competitiveness (IJCI-2015-25000) granted Dr. Colás and and the AGAUR Generalitat de Catalunya (2015FI_B00703) granted Tatiana Altadill. The authors would like to acknowledge the Proteomics and Metabolomics Shared Resource partially supported by Cancer Center Support Grant NIH/NCI grant P30-CA051008. The Institut de Salud Carlos III (FIS (PI13/01701)) also supported this project. Tissue samples were obtained with the support of “Xarxa Catalana de Bancs de Tumors” and “Plataforma de Biobancos” ISCIII (PT13/0010/0014). We would like to acknowledge experimental design guidance provided by Blanca Majem and all clinicians, Irene Campoy, Elena Martinez and Laura Devis for participating in the recruitment of clinical samples. We also acknowledge assistance in the statistical analysis provided by Ms Yiwen Wang. We finally thank volunteers that participated in the study.
Notes
Competing Interests
The authors declare that they have no competing interests.