Next Article in Journal
RKIP Regulates Differentiation-Related Features in Melanocytic Cells
Next Article in Special Issue
CircEHD2, CircNETO2 and CircEGLN3 as Diagnostic and Prognostic Biomarkers for Patients with Renal Cell Carcinoma
Previous Article in Journal
Current Strategies and Novel Therapeutic Approaches for Metastatic Urothelial Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 12
Issue 6
10.3390/cancers12061450
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extracellular Vesicles Enriched in hsa-miR-301a-3p and hsa-miR-1293 Dynamics in Clear Cell Renal Cell Carcinoma Patients: Potential Biomarkers of Metastatic Disease

by
Francisca Dias
1,2,
Ana Luísa Teixeira
1,*,
Inês Nogueira
1,2,3,
Mariana Morais
1,2,3,
Joana Maia
4,5,
Cristian Bodo
4,
Marta Ferreira
6,
Alexandra Silva
7,
Manuela Vilhena
7,
João Lobo
8,9,
José Pedro Sequeira
9,
Joaquina Maurício
6,
Jorge Oliveira
7,
Klaas Kok
10,
Bruno Costa-Silva
4 and
Rui Medeiros
1,2,3,11,12
1
Molecular Oncology and Viral Pathology Group, IPO-Porto Research Center (CI-IPOP), Portuguese Oncology Institute of Porto (IPO-Porto), Research Center- LAB2, E Bdg 1st floor, Rua Dr António Bernardino de Almeida, 4200-072 Porto, Portugal
2
Institute of Biomedical Sciences Abel Salazar, University of Porto (ICBAS-UP), Rua Jorge Viterbo Ferreira 228, 4050-513 Porto, Portugal
3
Research Department of the Portuguese League Against Cancer Regional Nucleus of the North (LPCC-NRN), Estrada da Circunvalação 6657, 4200-177 Porto, Portugal
4
Systems Oncology Group, Champalimaud Research, Champalimaud Centre for the Unknown, Av. Brasília, 1400-038 Lisbon, Portugal
5
Graduate Program in Areas of Basic and Applied Biology (GABBA), University of Porto, 4200-135 Porto, Portugal
6
Department of Medical Oncology, Portuguese Oncology Institute of Porto (IPO-Porto), Rua Dr António Bernardino de Almeida, 4200-072 Porto, Portugal
7
Department of Urology, Portuguese Oncology Institute of Porto (IPO-Porto), Rua Dr António Bernardino de Almeida, 4200-072 Porto, Portugal
8
Department of Pathology, Portuguese Oncology Institute of Porto (IPO-Porto), Rua Dr António Bernardino de Almeida, 4200-072 Porto, Portugal
9
Cancer Biology and Epigenetics Group, IPO-Porto Research Center (CI-IPOP), Portuguese Oncology Institute of Porto (IPO-Porto), Research Center- LAB3, F Bdg 1st floor, Rua Dr António Bernardino de Almeida, 4200-072 Porto, Portugal
10
Department of Genetics, University Medical Center Groningen (UMCG), University of Groningen, Hanzeplein 1, 9713 GZ Groningen, P.O. Box 30.001, 9700 RB Groningen, The Netherlands
11
Faculty of Medicine, University of Porto (FMUP), Alameda Prof. Hernâni Monteiro, 4200-319 Porto, Portugal
12
Biomedical Research Center (CEBIMED), Faculty of Health Sciences of Fernando Pessoa University (UFP), Praça 9 de Abril 349, 4249-004 Porto, Portugal
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(6), 1450; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061450
Submission received: 19 May 2020 / Accepted: 29 May 2020 / Published: 2 June 2020
(This article belongs to the Special Issue Noncoding RNAs in Renal Cell Carcinoma Landscape)
Browse Figures
Versions Notes

Abstract

:
Clear cell renal cell carcinoma (ccRCC) is the most aggressive subtype of kidney cancer and up to 40% of patients submitted to surgery with a curative intent will relapse. Thus, the aim of this study was to analyze the applicability of an Extracellular vesicle (EV) derived miRNA profile as potential prognosis biomarkers in ccRCC patients. We analyzed a nine-miRNA profile in plasma EVs from 32 ccRCC patients with localized disease (before and after surgery) and in 37 patients with metastatic disease. We observed that the levels of EV-derived hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p, and hsa-miR-301a-3p decreased after surgery, whereas hsa-miR-1293 EV-levels increased. Furthermore, metastatic patients presented higher levels of hsa-miR-301a-3p and lower levels of hsa-miR-1293 when compared to patients with localized disease after surgery. Functional enrichment analysis of the targets of the four miRNAs that decreased after surgery resulted in an enrichment of terms related to cell cycle, proliferation, and metabolism, suggesting that EV-miRNA enrichment in the presence of the tumor could represent an epigenetic mechanism to sustain tumor development. Taken together, these results suggest that EVs content varies depending on the presence or absence of the disease and that an increase of EV-derived hsa-miR-301a-3p, and decrease of EV-derived hsa-miR-1293, may be potential biomarkers of metastatic ccRCC.

1. Introduction

Renal cell carcinoma (RCC) is the most common solid cancer of the adult kidney and one of the most lethal urologic malignancies [1,2,3]. The most common histological subtype is clear cell RCC (ccRCC), which arises from the proximal tubular epithelial cells of the nephron and accounts for approximately 80% of all RCC cases [4]. CcRCC development is associated with two key cancer hallmarks: induction of angiogenesis and metabolic reprogramming, with several studies implying that ccRCC can be considered a metabolic disease [5,6,7,8]. In fact, the majority of known ccRCC related genes interact with cell metabolism pathways and are involved in energy, nutrient, iron and oxygen sensing [9]. One of the most well characterized pathways involved in the development of ccRCC, and the most studied as well, is the VHL/HIF pathway [10,11]. The loss or inactivation of pVHL function in ccRCC leads to a state of “pseudohypoxia” where stabilized HIF-1α and HIF-2α induce the transcription of hypoxia responsive genes, resulting in alterations such as increased glucose and ribose metabolism, pH deregulation, cell proliferation, and angiogenesis, giving ccRCC a high metastatic potential [1,9,12,13,14].
Actually one third of ccRCC patients present metastatic disease at the time of diagnosis, and up to 40% of patients submitted to surgery with a curative intent, will relapse within a five-year period [15]. The therapeutic options for ccRCC patients are limited since the hypoxic microenvironment makes these tumors chemo- and radio- resistant, leaving targeted therapies and immunotherapies as the only options available. However, patients submitted to these therapies tend to develop resistance to therapy within a short period of time [16,17]. Despite the fact that metastatic ccRCC remains incurable, the prognosis for recurrent ccRCC is variable and the detection of early relapse could have an impact on patients’ prognosis [18]. Thus, there is an urgent need for the definition of reliable biomarkers that can help stratify patients according to their metastatic risk. In the recent years, it became clear that one of the key players of tumor microenvironment modulation are the extracellular vesicles (EVs). EVs consist of a mixed population of microvesicles with different sizes, shed by cells, that enable cell-to-cell communication through the transport of active biomolecules from one cell to another [19,20]. Cancers have been found to highjack EV-mediated communication to facilitate several features of the multi-step metastatic process, including cell proliferation, immune suppression, epithelial-to-mesenchymal transition, migration, invasion, angiogenesis, and metastasis [5,21,22]. The release of EVs can be induced by several factors including, hypoxia, pH alterations, injury, platelet activation, irradiation and cellular stress, some of them involved in ccRCC progression [23]. One of the most studied classes of biomolecules carried by EVs are the microRNAs (miRNAs), which consist of small non-coding RNAs that are able to regulate gene expression at a post-transcriptional level, resulting in the attenuated translation of target mRNAs [24,25]. Thousands of protein-coding genes are regulated by miRNAs, and miRNAs are master regulators of diverse biological systems and have an impact in the body physiological responses [26,27]. Several EV-derived miRNAs have been studied and proposed as potential biomarkers in ccRCC [28,29,30,31,32,33,34,35,36,37,38]. However, due to the novelty of the field, most of the studies have been performed in cell lines and only a few used samples from ccRCC patients.
Therefore, the aim of this study is to analyse the impact of EV-derived miRNA profiles of 9 miRNAs (hsa-miR-25-3p, hsa-miR126-5p, hsa-miR-200c-3p, hsa-miR-210-3p, hsa-miR-301a-3p, hsa-miR-519d-3p, has-miR-1233-5p, hsa-miR-1246, and hsa-miR-1293) related to hypoxia and metabolism regulation in plasma EVs from ccRCC patients with localized disease and also in patients with metastatic disease.

2. Results

2.1. EVs Characterization

The EVs were characterized according to size, shape and purity. The NTA analysis indicated that the vast majority of isolated EVs presented a size range between 50 and 200 nm, which is consistent with the size of exosomes and small microvesicles (Figure 1A). A transmission electron microscopy (TEM) image (Figure 1B) shows the variability of sizes and morphology present in EVs from purified PFP. We also utilized EVs flow cytometry to confirm the purity of our EV isolates [39]. In all cases, more than 80% of particles present in our isolates corresponded to CFSE+ vesicular structures (Figure 1C).

2.2. EV-Derived miRNA Levels in ccRCC Patients’ Plasma Samples

Two patients from Group A relapsed during the follow-up and were removed from the analysis. The remaining 30 patients from group A are currently alive with no evidence of disease. In addition to that, two EV-derived miRNAs were excluded from the analysis. Hsa-miR-519d-3p was only detected in a small number of samples and hsa-1233-5p was not detected at all, which did not allow a statistical analysis.
The expression levels of the remaining eight miRNAs are represented in Figure 2. Focusing on the patients from Group A, we observed that the levels of EV-derived hsa-miR-25-3p (p = 0.003), hsa-miR-126-5p (p < 0.001), hsa-miR-200c-3p (p < 0.001) and hsa-miR-301a-3p (p = 0.006) decreased after surgery (Figure 2A–C,E). Hsa-miR-210-3p also decreased after tumor removal but the decrease was lower and only statistically significant when the localized disease samples were compared with the follow-up samples (p = 0.010) (Figure 2D). On the other hand, we observed that hsa-miR-1293 (p = 0.002) EV-levels increased after tumor removal (Figure 2G). Hsa-miR-1246 EV-levels also increased after tumor removal, but at a slower rate since the levels were only statistical significantly higher when we compared the localized and follow-up samples (p = 0.044) (Figure 2F).
When we compare the patients from Group A with the patients from Group B, we observe two scenarios. First, we see that some miRNAs have differences in their expression levels when comparing the samples from localized disease (Group A) with the metastatic disease (Group B). This is the case for hsa-miR-126-5p and hsa-miR-200c-3p that are downregulated in the metastatic patients (p < 0.001 and p < 0.001, respectively), (Figure 2B,C). On the other hand, we observed differences in miRNA levels when we compare the follow-up samples with the metastatic samples, with the metastatic samples presenting higher levels of hsa-miR-301a-3p (p = 0.026) and lower levels of miR-1293 (p = 0.004) (Figure 2E,G). In addition to that, hsa-miR-210-3p presented a tendency for being upregulated (p = 0.053) and hsa-miR-1246 presented a tendency for being downregulated (p = 0.053), both in the metastatic samples, when compared with the follow-up group (Figure 2D,F).
Regarding the other clinical pathological characteristics, we observed that patients with localized disease with tumors larger than 7 cm presented higher levels of hsa-miR-126-5p (p = 0.013) (Table S1). Additionally, patients who smoke presented higher levels of hsa-miR-1293 (p = 0.006) and hsa-miR-210-3p (p = 0.034) when compared to non-smokers and ex-smokers, respectively. We did not observe any statistical differences between the EV-derived miRNAs and hypertension and diabetes mellitus.

2.3. EV-Derived miRNAs Impact on Overall Survival of Metastatic ccRCC Patients

Since the patients from group A were alive with no evidence of disease on their last observation, we carried out the survival analysis on the patients from group B (Figure 3). The patients were divided in tertiles according the EV-derived miRNA levels using the –ΔCq values of each miRNA (high, intermediate and low levels). Regarding the 8 EV-derived miRNAs analyzed, only the hsa-miR-200c-3p presented statistical significant differences. Patients with a high or low expression level had a lower overall survival compared to patients with intermediate levels of hsa-miR-200c-3p (log rank Mantel Cox test, p = 0.025) (Figure 3C). In addition to that, patients with lower levels of has-miR-25-3p presented a tendency for a lower overall survival (log rank Mantel Cox test, p = 0.054) (Figure 3A).

2.4. Hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p, hsa-miR-301a-3p Overlaping Target Genes

Hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p, and hsa-miR-301a-3p significantly decreased after tumor removal, suggesting that their expression may be related to the presence of the tumor. As such, we used miRTarBase (version 8.0), the largest known online database of validated miRNA:mRNA target interactions, to establish a network of the target mRNAs of this miRNA profile and further evaluate its impact on ccRCC [40]. According to miRTarBase there are a total of 1381 validated target genes for these four miRNAs, but we only focused on the the 135 whose miRNA:mRNA interaction was validated according to strong validation methods (such as luciferase reporter assay, western blot and qPCR) (Table S2). From the 135 validated target genes, PTEN, VEGFA, TIMP2, and BCL2L11 were common to more than one of the miRNAs (Figure 3B). PTEN was the most common target, being regulated by hsa-miR-25-3p, hsa-miR-200c-3p, and hsa-miR-126-5p (Figure 4).

2.5. Functional Annotation and Pathway Enrichment Analysis of has-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p Targets Network

In order to explore the biological impact of the four-miRNA signature in patients with localized disease, we analyzed their 135 target genes with the STRINGapp from Cytoscape software (v3.7.2). A total of 134 of the 135 protein coding genes were filtered into a protein-protein interaction (PPI) network with 134 nodes and 889 edges that presented significant enrichment (p = 1.0 × 10−16). For a deeper understanding of the protein interactions we performed a Markov clustering (MCL), which resulted in the clustering of the proteins into 26 clusters according to their STRING interaction score (Figure 5, Table S3).
For the functional enrichment analysis, we focused on the largest cluster, which reduced the initial network of 134 proteins to 53 proteins. This cluster contained 53 nodes, with PTEN and BCL2L11 among them, and 328 edges and retained the significant PPI enrichment (p = 1.0 × 10−16). The functional enrichment analysis was performed with an FDR threshold of 1%, and the redundant terms were eliminated using a redundancy cutoff of 0.5, which resulted in a total of 549 functional enriched terms among the Reactome, KEGG, and GO categories (Tables S4–S8). The top 20 enriched terms for each category are represented in Figure 6.
Among the functionally enriched terms in the Reactome and KEGG pathways we can find AKT, BCL-2, PTEN and p53 signaling. Regarding the GO terms, we can observe that terms related with the Biological Processes were mostly about cellular metabolic processes and also regulation of cell cycle and cell proliferation.

3. Discussion

It is now well established that the formation and maintenance of a cancer niche deeply relies in the surrounding microenvironment that is composed by different cell types such as endothelial cells, stem cells, fibroblasts and immune cells [41]. Moreover, EVs play an important role in the complex network of communication that occurs inside the tumor microenvironment and also between the tumor microenvironment and the rest of the human body. This is due to the bioactive molecules they transfer between cells that are able to make a biological impact through the alteration of the recipient cell phenotype.
In our study, we observed that patients with localized ccRCC presented an alteration of the pattern of plasmatic EV-derived miRNAs after tumor removal, with a decrease of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p, and hsa-miR-301a-3p, and an increase of hsa-miR-1293 approximately four months after surgery, suggesting that the presence of the tumor has a key role in the EV network that is established in the patients’ body. In fact, hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p, and hsa-miR-301a-3p have already been reported in EVs, with an impact on cell proliferation, pre-metastatic niche formation, invasion, and metastization in other tumor models [42,43,44,45]. When we looked into the validated target genes of the four miRNAs that decreased after tumor removal, we observed that PTEN was the most common target, being regulated by hsa-miR-25-3p, hsa-miR-200c-3p and miR-301a-3p. Since one miRNA is capable of regulating several targets, and the same mRNA can be regulated by multiple miRNAs, some authors defend that the repressing capability of a single miRNA is quite modest and not capable of inducing significant changes on protein expression [46]. However, the fact that PTEN is repressed by three of the studied miRNAs suggests that this gene might be under a stronger repressing environment in the patients with localized disease. PTEN is a tumor suppressor gene whose major function relies on the inhibition of the PI3K/Akt pathway and consequent regulation of pathways related to cell growth and proliferation [47,48]. This gene is frequently altered in ccRCC and recent data from The Cancer Genome Atlas (TCGA) showed that PTEN loss of function in ccRCC patients was associated with a more vigorous cell metabolism and cell growth and also with a worse prognosis in survival and disease recurrence [49,50]. Moreover, a study performed by Zhang and colleagues demonstrated that astrocyte-derived EVs transferred PTEN-targeting miRNAs to metastatic cells, with an impact on mRNA and protein downregulation [51]. This epigenetic silencing of PTEN in metastatic cells leads to CCL2 secretion and myeloid cell recruitment, which promoted metastatic cell expansion through reduced apoptosis and enhanced proliferation through PI3K/Akt pathway activation [51].
Our functional annotation analysis on the validated target genes of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p showed that the most enriched terms were related to cell growth, proliferation and metabolism, with PTEN, AKT, BCL-2 and p53 signalling among the most enriched terms. Since these EV-derived miRNAs decrease after tumor removal, our hypothesis is that their EV-enrichment in the presence of the tumor could represent an epigenetic silencing mechanism used by ccRCC to sustain tumor development and growth through activation of the PI3K/Akt pathway. In fact, despite the fact that the overall mutation rate of PI3K/Akt pathway in ccRCC is relatively low, the overall activation of PI3K/Akt in ccRCC is high compared to other cancers, suggesting that the dyregulation of the PI3K/Akt pathway in ccRCC could be a consequence of epigenetic mechanisms mediated by EVs [52,53,54]. However, additional studies are necessary in order to validate this hypothesis.
Regarding the metastatic patients, EV-derived hsa-miR-301a-3p and hsa-miR-1293 were the most promising miRNAs to differentiate patients in follow-up with no evidence of disease from patients with metastatic disease. Hsa-mir-301a-3p presented a decreasing expression pattern after tumor removal and its levels kept decreasing until follow-up. However, the expression levels were significantly increased in the metastatic group, suggesting that this miRNA could play an important role in the metastization process and might have the potential to be used as prognostic biomarker. To support this hypothesis, a study performed by Wang and colleagues showed that, under a hypoxic microenvironment, pancreatic cancer cells generate EVs enriched in hsa-miR-301a-3p that induce M2 macrophages polarization through the PTEN/PI3K/Akt signalling pathway and promote metastization [43]. In addition to that, Yan and colleagues described an oncogenic role for hsa-miR-301a-3p in Laryngeal Squamous Cell Carcinoma development through inhibition of Smad4 and participation in the EMT process [55]. Regarding hsa-miR-1293, we observed the opposite effect since its EV-derived levels started to gradually increase after tumor removal until follow-up, but were significantly decreased in the metastatic group, suggesting a tumor suppresser role for this miRNA. In fact, Takagawa reported that hsa-miR-1293 was able to suppress in vivo tumor growth in a xenograft mouse model through the inhibition of BRD4 and DNA repair genes, and hypothesized that this miRNA could be a candidate for the development of miRNA-based therapeutics [56]. In addition to that, TIMP-1 is also a validated target of hsa-miR-1293 and its increased expression is associated with a worse prognosis in ccRCC patients, which reinforces the tumor suppressor role of hsa-miR-1293. [57,58].
Regarding the impact of the EV-derived miRNA profile in the metastatic patients’ overall survival, only hsa-miR-200c-3p presented significant, but contradictory, results since both its high and low levels were associated with a worse survival. Although based on a limited number of cases, our results may highlight the importance of a balance in EV-derived miRNA abundance and how disruption of that balance can have an impact on the patients’ prognosis. To the best of our knowledge this is the first study monitoring the differences of EV-derived miRNAs in several time points during the course of a one-year follow-up period of ccRCC patients with localized disease. This type of study is useful for a better understanding of the disease impact on the host EV dynamics and search for potential new prognosis biomarkers, such as hsa-miR-301a-3p and hsa-miR-1293. However, the small samples size in our study may limit the ability to distinguish meaningful differences, especially in terms of the survival analysis in the metastatic group, being essential to the replication of the associations reported in a larger sample size. In addition to that, the patients included in this study should continue to be monitored in order to get a deeper understanding of EV-miRNA enrichment during the course of the disease.

4. Materials and Methods

4.1. Ethics Statment

This study was conducted according to the principles of the Helsinki Declaration, having been approved by the ethics committee of the Portuguese Oncology Institute of Porto (IPO-Porto) (project reference: 251/015). All individuals signed a written informed consent in order to participate in the study.

4.2. Study Population

The analysis of the EV-derived miRNA profile was conducted through a hospital-based study involving a total of 69 ccRCC patients. All individuals were Caucasian from the north of Portugal, with histopathological diagnosis of ccRCC, and admitted and treated at the IPO-Porto between November 2015 and June 2019. The patients were divided into two groups: 32 patients diagnosed with localized disease that underwent surgical intervention composed Group A; and Group B was composed of 37 patients with metastatic disease. Clinical characteristics of patients were obtained from their medical records (Table 1). Blood from Group A patients was collected three times during the study: before undergoing surgery, approximately 4 months after surgery, and approximately 1 year after surgery, while from Group B, patients’ blood was collected once. All blood collections were performed during the morning period and stored at 4 °C immediately until sample processing.
Tumor classification and staging was established according to the tumor-node-metastasis (TNM) classification system of the American Joint Committee on Cancer (AJCC) 8th edition (2018) and the International Society of Urological Pathology (ISUP) Classification of Renal Neoplasia [59].

4.3. EVs Isolation

EVs were isolated from the plasma fraction using the Total Exosome Isolation from Plasma Kit (InvtrogenTM, Waltham, MA, USA) with additional protocol optimizations. Firstly, 8 mL of peripheral blood was collected from the patients in EDTA tubes and centrifuged 5 min at 1800× g to obtain the plasma fraction. The plasma fraction was then centrifuged 3 additional times at increasing speeds (300× g, 2100× g and 10,000× g) for a period of 15 min each in order to obtain platelet-free-plasma (PFP). After centrifugation, the supernatant was recovered and filtered through a 0.22 μM filter (GE Healthcare WhatmanTM, Chicago, IL, USA). After a 10-min treatment with proteinase K, the Total Exosome Isolation (TEI) reagent was added to 200 μL of PFP and the solution was incubated for 30 min at 4 °C. The precipitated EVs were recovered by 5 min centrifugation at 10,000× g at room temperature. The pellet containing the pre-enriched EVs was resuspended in filtered PBS (0.22 µm membrane filters) and stored at −80 °C until further analysis.

4.4. EVs NTA Analysis

Samples were analyzed for particle concentration and size distribution by the NS300 Nanoparticle Tracking Analysis (NTA) system (NanoSight–Malvern Panalytical, Malvern, UK). Samples were pre-diluted in filtered PBS to achieve a concentration within the range for optimal NTA analysis. Video acquisitions were performed using a camera level of 16 and a threshold between 5 and 7. Five to nine videos of 30 s were captured per sample. Analysis of particle concentration per mL and size distribution was performed with NTA software v3.4.

4.5. Quantification of Vesicular Structures by EVs Flow Cytometry

We employed EVs flow cytometry procedure for quantification of vesicular structures in our EV isolates, as recently described by Maia and colleagues [39]. Briefly, 2 × 109 particles of purified EVs were mixed with 40 µL of PBS containing Carboxyfluorescein Diacetate Succinimidyl Ester (CFSE – Thermo Fisher Scientific-LTI C34554, Waltham, MA, USA) in a final concentration of 40 µM and incubated for 90 min at 37 °C. For removal of unbound CFSE, size exclusion chromatography (SEC) columns (iZON - qEV original columns SP1, Oxford, UK) were used. Samples containing unstained or stained EVs, and appropriate controls, were diluted up to 500 µL of PBS and processed by qEV following manufacturer’s instructions. EVs-enriched fractions #7, #8 and #9 were then compiled and retrieved for analysis with the Flow Cytometer Apogee A60-Micro-Plus (Apogee Flow Systems, London, UK). The A60-Micro-Plus machine is equipped with three spatially separated lasers (488 nm – Position C, 405 nm – Position A and 638 nm–Position B), 7 fluorescence color detectors (525/50, LWP590, 530/30, 574/26, 590/40, 695/40, 676/36), and 3 light scatter detectors (SALS, MALS and LALS). For internal control across assays, before each FC experiment, we used two mixes of beads (Apogee–1493 and Apogee-1517). Before being loaded, samples were diluted in filtered PBS (0.22 µm membrane filters) to bring their concentration within the operational range of the equipment (maximum of 3000 events/second). All samples were run at a flow rate of 1.5 µL/minute using a 405 nm LALS threshold of 70. The 405 nm LALS PMT noise level was monitored and always maintained below 0.35. For the experiments depicted, the stopping criteria utilized was the number of events acquired, so samples were run until a minimum of 250,000 events was reached. The acquired data were exported and analyzed with FlowJo software v10.4.2 (FlowJo LLC, Ashland, OR, USA).

4.6. miRNA Extraction and cDNA Synthesis

MiRNA isolation and purification of cell-derived EVs was done using the Plasma/Serum RNA Purification Mini Kit from NORGEN (Norgen Biotek Corporation, Therold, ON, Canada) according to the manufacturer supplementary protocol for EV RNA purification from EVs already isolated from precipitation methods. RNA concentration and purity were measured using the NanoDrop Lite spectrophotometer (Thermo Scientific®, Waltham, MA, USA) and served as template for cDNA synthesis using a TaqManTM Advanced miRNA cDNA Synthesis Kit (Applied Biosystems®, Foster City, CA, USA) according to the manufacturer protocol.

4.7. Quantitative Real Time PCR

MiRNA expression levels were analyzed by quantitative real-time PCR. The reactions were carried out in a StepOnePlusTM qPCR Real-Time PCR machine, in a volume of 10 µL containing 1× TaqManTM Fast Advanced Master mix (Applied Biosystems), with 1X TaqManTM Advanced miRNA Assays probes (hsa-miR-25-3p 477994_mir; hsa-miR-126-5p 477888_mir; hsa-miR-200c-3p 478351_mir; hsa-miR210-3p 477970_mir; hsa-miR-301a-3p 477815_mir; hsa-miR-519d-3p 478986_mir; hsa-miR-1233-5p 479549_mir; hsa-miR-1246 477881_mir, and hsa-miR-1293 478692_mir-Applied Biosystems), and 2.5 µL cDNA. For miRNA expression normalization two housekeeping controls were used: hsa-let7a-5p (478575_mir-Applied Biosystems) and hsa-miR-16-5p (477860_mir-Applied Biosystems). These housekeeping miRNAs were chosen based on the fact that they are reported as typical EVs cargo [60]. The amplification conditions were as follows: holding stage 95 °C for 20 s, followed by 45 cycles of 95 °C for 1 s and 60 °C for 20 s. Three technical replicates were made for each sample. Data analysis was done using StepOneTM Sofware v2.2 (Applied Biosystems) with the same baseline and threshold set for each plate, in order to generate quantification cycle (Cq) values for all the miRNAs in each sample.

4.8. Statistical Analysis

Statistical analyses were done using IBM SPSS Statistics software for Windows (Version 22.0). According the mRNAs levels distribution, the Student t´test or Mann–Whitney U test were used in order to evaluate any statistical differences in the normalized expression of the EV-derived miRNAs. The quality of the housekeeping miRNAs was tested using the BestKeeper software [61]. Only hsa-let-7a-5p presented a stable behavior among all samples and was used to data normalization. The Kaplen Meier method and Log Rank test were used to establish the association of the EV-derived mRNA levels (low, intermediate, and high) to the overall survival.

4.9. Protein-Protein Interaction (PPI) Network and Cluster Analysis

The Search Tool for the Retrieval of Interacting Genes (STRING) database is an online tool that is used to develop protein-protein interaction (PPI) networks [62]. We used the STRINGapp in the Cytoscape software (v3.7.2) to construct the protein interaction relationship of the selected target genes encoding proteins. Those with a combined score of >0.4 were selected as significant. In order to get a better understanding of the resulting PPI network, we grouped the proteins according to the strength of their STRING interaction score using the clusterMaker app for Cytoscape to perform Markov clustering (MCL) [63,64].

4.10. Functional Annotation and Pathway Enrichment Analysis

Gene ontology analysis (GO) is a common useful method for annotating genes and gene products, allowing the identification of characteristic biological attributes for transcriptome data. The functional enrichment analysis of GO, Reactome and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was performed with the enrichment anaylsis tool of the STRINGapp, with a False Discovery Rate (FDR) of p < 0.01 as significance threshold. The enrichment results were filtered in order to remove redundant terms. The redundancy filtering takes the list of enriched terms sorted by FDR value and removes the terms that are too similar to any of the previous, better scoring terms that were not themselves removed. The similarity between two terms is measured by the Jaccard índex of the sets of genes annotated by the two terms. A term is added to the filtered list only if it has Jaccard similarity less than the used specified redundancy cutoff to any other term already in the filtered list.

5. Conclusions

Taken together, these results suggest that EVs content varies depending on the presence or absence of the disease and an increase of EVs enriched in hsa-miR-301a-3p, and decrease of EVs enriched in hsa-miR-1293, may be potential biomarkers of metastatic disease in ccRCC patients.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/12/6/1450/s1, Table S1: Association of EV-derived miRNA profile expression with Clinical-pathological characteristics of patients with localized disease (Group A)., Table S2: Hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p mRNA targets validated with strong evidence according miRTarBase V8.0, Table S3: Protein-protein-interaction (PPI) network clusters according to Marckov Clusterring (MCL) analysis, Table S4: Reactome pathway enrichment analysis for terms with FDR p value <0.01, Table S5: KEGG pathway enrichment analysis, Table S6: Gene Ontology (GO) Cellular Components enrichment analysis, Table S7: Gene Ontology (GO) Biological data enrichment analysis, Table S8: Gene Ontology (GO) Molecular Function data enrichment analysis.

Author Contributions

Conceptualization, F.D. and A.L.T.; investigation and methodology, F.D., I.N., M.M., J.M. (Joana Maia) and C.B.; software, F.D.; formal analysis, F.D., A.L.T. and B.C.-S.; resources, A.S., M.V., J.O., B.C.-S. and R.M.; data curation, F.D., M.F., J.L., J.P.S., and J.M. (Joaquina Maurício); writing—original draft preparation, F.D. and A.L.T.; writing—review and editing, K.K., B.C.-S. and R.M.; supervision, R.M.; project administration, A.L.T. and R.M.; funding acquisition, R.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund (ERDF) [NORTE-01-0247-FEDER-033399] and by the IPO-Porto Research Center [CI-IPOP-21-2015].

Acknowledgments

The authors would like to thank the nursing team of the Urology Department from the Portuguese Oncology Institute of Porto (IPO-Porto) for their support regarding the patients’ recruitment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ricketts, C.J.; De Cubas, A.A.; Fan, H.; Smith, C.C.; Lang, M.; Reznik, E.; Bowlby, R.; Gibb, E.A.; Akbani, R.; Beroukhim, R.; et al. The Cancer Genome Atlas Comprehensive Molecular Characterization of Renal Cell Carcinoma. Cell Rep. 2018, 23, 313–326. [Google Scholar] [CrossRef] [Green Version]
  2. Capitanio, U.; Montorsi, F. Renal cancer. Lancet 2016, 387, 894–906. [Google Scholar] [CrossRef]
  3. Moch, H.; Cubilla, A.L.; Humphrey, P.A.; Reuter, V.E.; Ulbright, T.M. The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs-Part A: Renal, Penile, and Testicular Tumours. Eur. Urol. 2016, 70, 93–105. [Google Scholar] [CrossRef]
  4. Escudier, B.; Porta, C.; Schmidinger, M.; Rioux-Leclercq, N.; Bex, A.; Khoo, V.; Grunwald, V.; Gillessen, S.; Horwich, A. Renal cell carcinoma: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-updagger. Ann. Oncol. 2019, 30, 706–720. [Google Scholar] [CrossRef] [Green Version]
  5. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  6. Ward, P.S.; Thompson, C.B. Metabolic reprogramming: A cancer hallmark even warburg did not anticipate. Cancer Cell 2012, 21, 297–308. [Google Scholar] [CrossRef] [Green Version]
  7. Courtney, K.D.; Bezwada, D.; Mashimo, T.; Pichumani, K.; Vemireddy, V.; Funk, A.M.; Wimberly, J.; McNeil, S.S.; Kapur, P.; Lotan, Y.; et al. Isotope Tracing of Human Clear Cell Renal Cell Carcinomas Demonstrates Suppressed Glucose Oxidation In Vivo. Cell Metab. 2018, 28, 793–800. [Google Scholar] [CrossRef] [Green Version]
  8. Wettersten, H.I.; Aboud, O.A.; Lara, P.N., Jr.; Weiss, R.H. Metabolic reprogramming in clear cell renal cell carcinoma. Nat. Rev. Nephrol. 2017, 13, 410–419. [Google Scholar] [CrossRef]
  9. Linehan, W.M.; Srinivasan, R.; Schmidt, L.S. The genetic basis of kidney cancer: A metabolic disease. Nat. Rev. Urol. 2010, 7, 277–285. [Google Scholar] [CrossRef]
  10. Baldewijns, M.M.; van Vlodrop, I.J.H.; Vermeulen, P.B.; Soetekouw, P.M.M.B.; van Engeland, M.; de Bruine, A.P. VHL and HIF signalling in renal cell carcinogenesis. J. Pathol. 2010, 221, 125–138. [Google Scholar] [CrossRef]
  11. Haase, V.H. The VHL tumor suppressor: Master regulator of HIF. Curr. Pharm. Des. 2009, 15, 3895–3903. [Google Scholar] [CrossRef] [Green Version]
  12. Linehan, W.M.; Schmidt, L.S.; Crooks, D.R.; Wei, D.; Srinivasan, R.; Lang, M.; Ricketts, C.J. The Metabolic Basis of Kidney Cancer. Cancer Discov. 2019. [Google Scholar] [CrossRef] [Green Version]
  13. Morais, M.; Dias, F.; Teixeira, A.L.; Medeiros, R. MicroRNAs and altered metabolism of clear cell renal cell carcinoma: Potential role as aerobic glycolysis biomarkers. Biochim. Biophys. Acta 2017, 1861, 2175–2185. [Google Scholar] [CrossRef]
  14. Dias, F.; Teixeira, A.L.; Santos, J.I.; Gomes, M.; Nogueira, A.; Assis, J.; Medeiros, R. Renal cell carcinoma development and miRNAs: A possible link to the EGFR pathway. Pharmacogenomics 2013, 14, 1793–1803. [Google Scholar] [CrossRef]
  15. Dabestani, S.; Thorstenson, A.; Lindblad, P.; Harmenberg, U.; Ljungberg, B.; Lundstam, S. Renal cell carcinoma recurrences and metastases in primary non-metastatic patients: A population-based study. World J. Urol. 2016, 34, 1081–1086. [Google Scholar] [CrossRef]
  16. Ravaud, A.; Gross-Goupil, M. Overcoming resistance to tyrosine kinase inhibitors in renal cell carcinoma. Cancer Treat. Rev. 2012, 38, 996–1003. [Google Scholar] [CrossRef]
  17. Auvray, M.; Auclin, E.; Barthelemy, P.; Bono, P.; Kellokumpu-Lehtinen, P.; Gross-Goupil, M.; de Velasco, G.; Powles, T.; Mouillet, G.; Vano, Y.-A.; et al. Second-line targeted therapies after nivolumab-ipilimumab failure in metastatic renal cell carcinoma. Eur. J. Cancer 2019, 108, 33–40. [Google Scholar] [CrossRef]
  18. Dabestani, S.; Beisland, C.; Stewart, G.D.; Bensalah, K.; Gudmundsson, E.; Lam, T.B.; Gietzmann, W.; Zakikhani, P.; Marconi, L.; Fernández-Pello, S.; et al. Long-term Outcomes of Follow-up for Initially Localised Clear Cell Renal Cell Carcinoma: RECUR Database Analysis. Eur. Urol. Focus 2019, 5, 857–866. [Google Scholar] [CrossRef] [Green Version]
  19. Minciacchi, V.R.; Freeman, M.R.; Di Vizio, D. Extracellular vesicles in cancer: Exosomes, microvesicles and the emerging role of large oncosomes. Semin. Cell Dev. Biol. 2015, 40, 41–51. [Google Scholar] [CrossRef] [Green Version]
  20. Nogues, L.; Benito-Martin, A.; Hergueta-Redondo, M.; Peinado, H. The influence of tumour-derived extracellular vesicles on local and distal metastatic dissemination. Mol. Asp. Med. 2018, 60, 15–26. [Google Scholar] [CrossRef]
  21. Zhao, H.; Achreja, A.; Iessi, E.; Logozzi, M.; Mizzoni, D.; Di Raimo, R.; Nagrath, D.; Fais, S. The key role of extracellular vesicles in the metastatic process. Biochim. Biophys. Acta Rev. Cancer 2018, 1869, 64–77. [Google Scholar] [CrossRef]
  22. Jurj, A.; Zanoaga, O.; Braicu, C.; Lazar, V.; Tomuleasa, C.; Irimie, A.; Berindan-Neagoe, I. A Comprehensive Picture of Extracellular Vesicles and Their Contents. Molecular Transfer to Cancer Cells. Cancers 2020, 12, 298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kucharzewska, P.; Belting, M. Emerging roles of extracellular vesicles in the adaptive response of tumour cells to microenvironmental stress. J. Extracell. Vesicles 2013, 2. [Google Scholar] [CrossRef]
  24. Gambari, R.; Fabbri, E.; Borgatti, M.; Lampronti, I.; Finotti, A.; Brognara, E.; Bianchi, N.; Manicardi, A.; Marchelli, R.; Corradini, R. Targeting microRNAs involved in human diseases: A novel approach for modification of gene expression and drug development. Biochem. Pharm. 2011, 82, 1416–1429. [Google Scholar] [CrossRef]
  25. Mohr, A.M.; Mott, J.L. Overview of microRNA biology. Semin. Liver Dis. 2015, 35, 3–11. [Google Scholar] [CrossRef] [Green Version]
  26. Cai, Y.; Yu, X.; Hu, S.; Yu, J. A brief review on the mechanisms of miRNA regulation. Genom. Proteom. Bioinform. 2009, 7, 147–154. [Google Scholar] [CrossRef] [Green Version]
  27. Lai, X.; Wolkenhauer, O.; Vera, J. Understanding microRNA-mediated gene regulatory networks through mathematical modelling. Nucleic Acids Res. 2016, 44, 6019–6035. [Google Scholar] [CrossRef] [Green Version]
  28. Grange, C.; Tapparo, M.; Collino, F.; Vitillo, L.; Damasco, C.; Deregibus, M.C.; Tetta, C.; Bussolati, B.; Camussi, G. Microvesicles released from human renal cancer stem cells stimulate angiogenesis and formation of lung premetastatic niche. Cancer Res. 2011, 71, 5346–5356. [Google Scholar] [CrossRef] [Green Version]
  29. Teixeira, A.L.; Dias, F.; Gomes, M.; Fernandes, M.; Medeiros, R. Circulating biomarkers in renal cell carcinoma: The link between microRNAs and extracellular vesicles, where are we now? J. Kidney Cancer VHL 2014, 1, 84. [Google Scholar] [CrossRef]
  30. Grange, C.; Brossa, A.; Bussolati, B. Extracellular Vesicles and Carried miRNAs in the Progression of Renal Cell Carcinoma. Int. J. Mol. Sci. 2019, 20, 1832. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, L.; Yang, G.; Zhao, D.; Wang, J.; Bai, Y.; Peng, Q.; Wang, H.; Fang, R.; Chen, G.; Wang, Z.; et al. CD103-positive CSC exosome promotes EMT of clear cell renal cell carcinoma: Role of remote MiR-19b-3p. Mol. Cancer 2019, 18, 86. [Google Scholar] [CrossRef] [PubMed]
  32. Song, S.; Long, M.; Yu, G.; Cheng, Y.; Yang, Q.; Liu, J.; Wang, Y.; Sheng, J.; Wang, L.; Wang, Z.; et al. Urinary exosome miR-30c-5p as a biomarker of clear cell renal cell carcinoma that inhibits progression by targeting HSPA5. J. Cell Mol. Med. 2019, 23, 6755–6765. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, W.; Ni, M.; Su, Y.; Wang, H.; Zhu, S.; Zhao, A.; Li, G. MicroRNAs in Serum Exosomes as Potential Biomarkers in Clear-cell Renal Cell Carcinoma. Eur. Urol. Focus 2018, 4, 412–419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wang, X.; Wang, T.; Chen, C.; Wu, Z.; Bai, P.; Li, S.; Chen, B.; Liu, R.; Zhang, K.; Li, W.; et al. Serum exosomal miR-210 as a potential biomarker for clear cell renal cell carcinoma. J. Cell Biochem. 2018. [Google Scholar] [CrossRef] [PubMed]
  35. Fujii, N.; Hirata, H.; Ueno, K.; Mori, J.; Oka, S.; Shimizu, K.; Kawai, Y.; Inoue, R.; Yamamoto, Y.; Matsumoto, H.; et al. Extracellular miR-224 as a prognostic marker for clear cell renal cell carcinoma. Oncotarget 2017, 8, 109877–109888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. He, J.; He, J.; Min, L.; He, Y.; Guan, H.; Wang, J.; Peng, X. Extracellular vesicles transmitted miR-31-5p promotes sorafenib resistance by targeting MLH1 in renal cell carcinoma. Int. J. Cancer 2019. [Google Scholar] [CrossRef]
  37. Zhang, J.; Wang, W.; Chen, S.; Wang, Y.; Xia, S. Distinct non-invasive evaluation values of tumor-derived cell-free microRNAs, circulating microvesicles and exosomal microRNAs after renal carcinoma cryoablation. Biochem. Biophys. Res. Commun. 2018, 503, 1278–1283. [Google Scholar] [CrossRef]
  38. Butz, H.; Nofech-Mozes, R.; Ding, Q.; Khella, H.W.Z.; Szabó, P.M.; Jewett, M.; Finelli, A.; Lee, J.; Ordon, M.; Stewart, R.; et al. Exosomal MicroRNAs Are Diagnostic Biomarkers and Can Mediate Cell–Cell Communication in Renal Cell Carcinoma. Eur. Urol. Focus 2016, 2, 210–218. [Google Scholar] [CrossRef]
  39. Maia, J.; Batista, S.; Couto, N.; Gregorio, A.; Bodo, C.; Elzanowska, J.; Moraes, M.; Costa-Silva, B. Population Analysis of Extracellular Vesicles in Microvolumes of Biofluids. bioRxiv 2020. [Google Scholar] [CrossRef]
  40. Chou, C.-H.; Shrestha, S.; Yang, C.-D.; Chang, N.-W.; Lin, Y.-L.; Liao, K.-W.; Huang, W.-C.; Sun, T.-H.; Tu, S.-J.; Lee, W.-H.; et al. miRTarBase update 2018: A resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018, 46, D296–D302. [Google Scholar] [CrossRef]
  41. Brücher, B.L.D.M.; Jamall, I.S. Cell-cell communication in the tumor microenvironment, carcinogenesis, and anticancer treatment. Cell Physiol. Biochem. 2014, 34, 213–243. [Google Scholar] [CrossRef] [PubMed]
  42. Zeng, Z.; Li, Y.; Pan, Y.; Lan, X.; Song, F.; Sun, J.; Zhou, K.; Liu, X.; Ren, X.; Wang, F.; et al. Cancer-derived exosomal miR-25-3p promotes pre-metastatic niche formation by inducing vascular permeability and angiogenesis. Nat. Commun. 2018, 9, 5395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, X.; Luo, G.; Zhang, K.; Cao, J.; Huang, C.; Jiang, T.; Liu, B.; Su, L.; Qiu, Z. Hypoxic Tumor-Derived Exosomal miR-301a Mediates M2 Macrophage Polarization via PTEN/PI3Kγ to Promote Pancreatic Cancer Metastasis. Cancer Res. 2018, 78, 4586–4598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kawakubo-Yasukochi, T.; Morioka, M.; Hazekawa, M.; Yasukochi, A.; Nishinakagawa, T.; Ono, K.; Kawano, S.; Nakamura, S.; Nakashima, M. miR-200c-3p spreads invasive capacity in human oral squamous cell carcinoma microenvironment. Mol. Carcinog. 2018, 57, 295–302. [Google Scholar] [CrossRef]
  45. Mormile, R. Prevention and management of allograft rejection in heart transplantation: Are miR-126-5p and miR-126-5p-expressing mesenchymal stem cell exosomes the Ariadne’s thread? Cardiovasc. Pathol. 2020, 47, 107208. [Google Scholar] [CrossRef]
  46. Baek, D.; Villén, J.; Shin, C.; Camargo, F.D.; Gygi, S.P.; Bartel, D.P. The impact of microRNAs on protein output. Nature 2008, 455, 64–71. [Google Scholar] [CrossRef] [Green Version]
  47. Maehama, T.; Taylor, G.S.; Dixon, J.E. PTEN and myotubularin: Novel phosphoinositide phosphatases. Annu. Rev. Biochem. 2001, 70, 247–279. [Google Scholar] [CrossRef]
  48. Geva, G.A.; Gielchinsky, I.; Aviv, N.; Max, K.E.A.; Gofrit, O.N.; Gur-Wahnon, D.; Ben-Dov, I.Z. Urine cell-free microRNA as biomarkers for transitional cell carcinoma. BMC Res. Notes 2017, 10, 641. [Google Scholar] [CrossRef] [Green Version]
  49. Fan, C.; Zhao, C.; Wang, F.; Li, S.; Wang, J. Significance of PTEN Mutation in Cellular Process, Prognosis, and Drug Selection in Clear Cell Renal Cell Carcinoma. Front. Oncol. 2019, 9, 357. [Google Scholar] [CrossRef]
  50. Guo, S.-P.; Zhai, Y.-Q.; Wang, W.-L.; Ma, F.-C. Expression and significance of a new tumor suppression gene PTEN in primary renal cell carcinoma. Ai Zheng 2002, 21, 582–587. [Google Scholar]
  51. Zhang, L.; Zhang, S.; Yao, J.; Lowery, F.J.; Zhang, Q.; Huang, W.-C.; Li, P.; Li, M.; Wang, X.; Zhang, C.; et al. Microenvironment-induced PTEN loss by exosomal microRNA primes brain metastasis outgrowth. Nature 2015, 527, 100–104. [Google Scholar] [CrossRef] [PubMed]
  52. Akbani, R.; Ng, P.K.S.; Werner, H.M.J.; Shahmoradgoli, M.; Zhang, F.; Ju, Z.; Liu, W.; Yang, J.-Y.; Yoshihara, K.; Li, J.; et al. A pan-cancer proteomic perspective on The Cancer Genome Atlas. Nat. Commun. 2014, 5, 3887. [Google Scholar] [CrossRef] [PubMed]
  53. Guo, H.; German, P.; Bai, S.; Barnes, S.; Guo, W.; Qi, X.; Lou, H.; Liang, J.; Jonasch, E.; Mills, G.B.; et al. The PI3K/AKT Pathway and Renal Cell Carcinoma. J. Genet. Genom. 2015, 42, 343–353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Li, J.; Lu, Y.; Akbani, R.; Ju, Z.; Roebuck, P.L.; Liu, W.; Yang, J.-Y.; Broom, B.M.; Verhaak, R.G.W.; Kane, D.W.; et al. TCPA: A resource for cancer functional proteomics data. Nat. Methods 2013, 10, 1046–1047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lu, Y.; Gao, W.; Zhang, C.; Wen, S.; Huangfu, H.; Kang, J.; Wang, B. Hsa-miR-301a-3p Acts as an Oncogene in Laryngeal Squamous Cell Carcinoma via Target Regulation of Smad4. J. Cancer 2015, 6, 1260–1275. [Google Scholar] [CrossRef] [Green Version]
  56. Takagawa, Y.; Gen, Y.; Muramatsu, T.; Tanimoto, K.; Inoue, J.; Harada, H.; Inazawa, J. miR-1293, a Candidate for miRNA-Based Cancer Therapeutics, Simultaneously Targets BRD4 and the DNA Repair Pathway. Mol. Ther. 2020, in press. [Google Scholar] [CrossRef]
  57. Kawata, N.; Yamaguchi, K.; Igarashi, T.; Takahashi, S. TIMP-1 as well as Microvessel Invasion and High Nuclear Grade Is a Significant Determinant Factor for Extension of Tumor Diameter in Localized RCC. J. Oncol. 2016, 2016, 5035127. [Google Scholar] [CrossRef] [Green Version]
  58. Li, P.; Ma, Y.; Wang, Y.; Chen, T.; Wang, H.; Chu, H.; Zhao, G.; Zhang, G. Identification of miR-1293 potential target gene: TIMP-1. Mol. Cell Biochem. 2013, 384, 1–6. [Google Scholar] [CrossRef]
  59. Srigley, J.R.; Delahunt, B.; Eble, J.N.; Egevad, L.; Epstein, J.I.; Grignon, D.; Hes, O.; Moch, H.; Montironi, R.; Tickoo, S.K.; et al. The International Society of Urological Pathology (ISUP) Vancouver Classification of Renal Neoplasia. Am. J. Surg. Pathol. 2013, 37, 1469–1489. [Google Scholar] [CrossRef]
  60. Zeringer, E.; Li, M.; Barta, T.; Schageman, J.; Pedersen, K.W.; Neurauter, A.; Magdaleno, S.; Setterquist, R.; Vlassov, A.V. Methods for the extraction and RNA profiling of exosomes. World J. Methodol. 2013, 3, 11–18. [Google Scholar] [CrossRef]
  61. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol. Lett. 2004, 26, 509–515. [Google Scholar] [CrossRef] [PubMed]
  62. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Morris, J.H.; Apeltsin, L.; Newman, A.M.; Baumbach, J.; Wittkop, T.; Su, G.; Bader, G.D.; Ferrin, T.E. clusterMaker: A multi-algorithm clustering plugin for Cytoscape. BMC Bioinform. 2011, 12, 436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Enright, A.J.; Van Dongen, S.; Ouzounis, C.A. An efficient algorithm for large-scale detection of protein families. Nucleic Acids Res. 2002, 30, 1575–1584. [Google Scholar] [CrossRef]
Cancers 12 01450 g001 550
Figure 1. (A) Nanoparticle tracking analysis (NTA) of EVs derived from plasma of ccRCC patients. The red error bars indicate ± 1 standard error of the mean. (B) Transmission electron microscopy (TEM) of EVs from purified platelet-free plasma (PFP). The TEM image was acquired in the Histology and Electron Microscopy platform from I3S Porto using a Transmission Electron Microscope Jeol JEM 1400. (C) Proportion of fluorescent EVs as observed by flow cytometry from a sample EV isolate derived from a patient, either previously stained with CFSE (left panel) or without previous staining (right panel).
Figure 1. (A) Nanoparticle tracking analysis (NTA) of EVs derived from plasma of ccRCC patients. The red error bars indicate ± 1 standard error of the mean. (B) Transmission electron microscopy (TEM) of EVs from purified platelet-free plasma (PFP). The TEM image was acquired in the Histology and Electron Microscopy platform from I3S Porto using a Transmission Electron Microscope Jeol JEM 1400. (C) Proportion of fluorescent EVs as observed by flow cytometry from a sample EV isolate derived from a patient, either previously stained with CFSE (left panel) or without previous staining (right panel).
Cancers 12 01450 g001
Cancers 12 01450 g002 550
Figure 2. EV-derived miRNA expression (-ΔCq) in pre-surgery, post-surgery, and follow-up samples (Group A) and metastatic disease samples (Group B) of ccRCC patients. (A) hsa-miR-25-3p plasmatic EV-expression; (B) hsa-miR-126-5p plasmatic EV-expression; (C) hsa-miR-200c-3p plasmatic EV-expression; (D) hsa-miR-210-3p plasmatic EV-expression; (E) hsa-miR-301a-3p plasmatic EV-expression; (F) hsa-miR-1246 plasmatic EV-expression; (G) hsa-miR-1293 plasmatic EV-expression; * p < 0.05, ** p < 0.001.
Figure 2. EV-derived miRNA expression (-ΔCq) in pre-surgery, post-surgery, and follow-up samples (Group A) and metastatic disease samples (Group B) of ccRCC patients. (A) hsa-miR-25-3p plasmatic EV-expression; (B) hsa-miR-126-5p plasmatic EV-expression; (C) hsa-miR-200c-3p plasmatic EV-expression; (D) hsa-miR-210-3p plasmatic EV-expression; (E) hsa-miR-301a-3p plasmatic EV-expression; (F) hsa-miR-1246 plasmatic EV-expression; (G) hsa-miR-1293 plasmatic EV-expression; * p < 0.05, ** p < 0.001.
Cancers 12 01450 g002
Cancers 12 01450 g003 550
Figure 3. Overall survival analysis of metastatic ccRCC patients according to EV-derived miRNAs expression. Comparison of survival curves was made with Log-rank (Mantel-Cox) test. (A) hsa-miR-25-3p (low levels–12 patients; intermediate levels–11 patients and high levels–12 patients). (B) hsa-miR-125-5p (low levels–12 patients; intermediate levels–13 patients and high levels–11 patients). (C) hsa-miR-200-3p (low levels–12 patients; intermediate levels–12 patients and high levels–12 patients). (D) hsa-miR-210-3p (low levels–11 patients; intermediate levels–13 patients and high levels–12 patients). (E) hsa-miR-301a-3p (low levels–11 patients; intermediate levels–10 patients and high levels–11 patients). (F) hsa-miR-1246 (low levels–11 patients; intermediate levels–11 patients and high levels–12 patients). (G) hsa-miR-1293 (low levels–12 patients; intermediate levels–11 patients and high levels–12 patients).
Figure 3. Overall survival analysis of metastatic ccRCC patients according to EV-derived miRNAs expression. Comparison of survival curves was made with Log-rank (Mantel-Cox) test. (A) hsa-miR-25-3p (low levels–12 patients; intermediate levels–11 patients and high levels–12 patients). (B) hsa-miR-125-5p (low levels–12 patients; intermediate levels–13 patients and high levels–11 patients). (C) hsa-miR-200-3p (low levels–12 patients; intermediate levels–12 patients and high levels–12 patients). (D) hsa-miR-210-3p (low levels–11 patients; intermediate levels–13 patients and high levels–12 patients). (E) hsa-miR-301a-3p (low levels–11 patients; intermediate levels–10 patients and high levels–11 patients). (F) hsa-miR-1246 (low levels–11 patients; intermediate levels–11 patients and high levels–12 patients). (G) hsa-miR-1293 (low levels–12 patients; intermediate levels–11 patients and high levels–12 patients).
Cancers 12 01450 g003
Cancers 12 01450 g004 550
Figure 4. (A) Venn diagram of the validated target genes of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p obtained using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/). (B) Representation of the overlaping target genes of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p.
Figure 4. (A) Venn diagram of the validated target genes of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p obtained using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/). (B) Representation of the overlaping target genes of hsa-miR-25-3p, hsa-miR-126-5p, hsa-miR-200c-3p and hsa-miR-301a-3p.
Cancers 12 01450 g004
Cancers 12 01450 g005 550
Figure 5. Clustered PPI network of the 134 target genes. Clustering was performed using the Markov clustering (MCL) algorithm in the clusterMaker2 Cytoscape app using an inflation value of 4.0. Proteins without any interaction partners within the network (singletons) are omitted from the visualization.
Figure 5. Clustered PPI network of the 134 target genes. Clustering was performed using the Markov clustering (MCL) algorithm in the clusterMaker2 Cytoscape app using an inflation value of 4.0. Proteins without any interaction partners within the network (singletons) are omitted from the visualization.
Cancers 12 01450 g005
Cancers 12 01450 g006 550
Figure 6. Reactome, KEGG and GO analysis of the 53 selected genes. The functional enrichment anaylis was made with the STRINGapp from Cytoscape.
Figure 6. Reactome, KEGG and GO analysis of the 53 selected genes. The functional enrichment anaylis was made with the STRINGapp from Cytoscape.
Cancers 12 01450 g006
Table
Table 1. Clinical-pathological characteristics of the study population. Group A is constituted by the patients with localized disease and Group B is constituted by the patients with metastatic disease.
Table 1. Clinical-pathological characteristics of the study population. Group A is constituted by the patients with localized disease and Group B is constituted by the patients with metastatic disease.
Clinical-Pathological CharacteristicsGroup A (N = 32)Group B (N = 37)
N%N%
Gender
Male2475.02670.3
Female825.01129.7
Age
Average ± SD61.9 ± 12.462.4 ± 9.9
Type of surgery
Partial nephrectomy1856.300.0
Radical nephrectomy1443.83491.9
No surgery0038.1
Tumor size
<7 cm2475.0924.32
≥7 cm825.02259.46
Undetermined00616.22
T
T11753.11027.03
T213.1616.22
T31340.61643.24
T40038.12
Tx13.125.40
N
N0001540.54
N1–N200513.51
Nx001745.95
M
M0271002259.46
M1001540.54
Mx5000.0
Clinical Stage
I1753.11027.03
II13.1616.22
III1340.61643.24
IV0038.12
Not available13.135.40
ISUP classification
126.312.70
21856.31437.84
31134.3821.62
413.1924.32
Not available00513.51
Smoking status
Smoker618.838.11
Ex-smoker1134.4718.92
Non-smoker1443.82772,97
Not available13.100.0
Hypertension
Yes2062.51540.54
No1134.42259.45
Not available13.100.0
Diabetes mellitus
DM I00.000.0
DM II928.11335.14
No2268.82362.16
Not available13.112.70

Share and Cite

MDPI and ACS Style

Dias, F.; Teixeira, A.L.; Nogueira, I.; Morais, M.; Maia, J.; Bodo, C.; Ferreira, M.; Silva, A.; Vilhena, M.; Lobo, J.; et al. Extracellular Vesicles Enriched in hsa-miR-301a-3p and hsa-miR-1293 Dynamics in Clear Cell Renal Cell Carcinoma Patients: Potential Biomarkers of Metastatic Disease. Cancers 2020, 12, 1450. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061450

AMA Style

Dias F, Teixeira AL, Nogueira I, Morais M, Maia J, Bodo C, Ferreira M, Silva A, Vilhena M, Lobo J, et al. Extracellular Vesicles Enriched in hsa-miR-301a-3p and hsa-miR-1293 Dynamics in Clear Cell Renal Cell Carcinoma Patients: Potential Biomarkers of Metastatic Disease. Cancers. 2020; 12(6):1450. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061450

Chicago/Turabian Style

Dias, Francisca, Ana Luísa Teixeira, Inês Nogueira, Mariana Morais, Joana Maia, Cristian Bodo, Marta Ferreira, Alexandra Silva, Manuela Vilhena, João Lobo, and et al. 2020. "Extracellular Vesicles Enriched in hsa-miR-301a-3p and hsa-miR-1293 Dynamics in Clear Cell Renal Cell Carcinoma Patients: Potential Biomarkers of Metastatic Disease" Cancers 12, no. 6: 1450. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061450

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