EV-Associated miRNAs from Peritoneal Lavage are a Source of Biomarkers in Endometrial Cancer

Endometrial cancer (EC) is the sixth most common cancer in women worldwide and is responsible for more than 89,000 deaths every year. Mortality is associated with presence of poor prognostic factors at diagnosis, i.e., diagnosis at an advanced stage, with a high grade and/or an aggressive histology. Development of novel approaches that would permit us to improve the clinical management of EC patients is an unmet need. In this study, we investigate a novel approach to identify highly sensitive and specific biomarkers of EC using extracellular vesicles (EVs) isolated from the peritoneal lavage of EC patients. EVs of peritoneal lavages of 25 EC patients were isolated and their miRNA content was compared with miRNAs of EVs isolated from the ascitic fluid of 25 control patients. Expression of the EV-associated miRNAs was measured using the Taqman OpenArray technology that allowed us to detect 371 miRNAs. The analysis showed that 114 miRNAs were significantly dysregulated in EC patients, among which eight miRNAs, miRNA-383-5p, miRNA-10b-5p, miRNA-34c-3p, miRNA-449b-5p, miRNA-34c-5p, miRNA-200b-3p, miRNA-2110, and miRNA-34b-3p, demonstrated a classification performance at area under the receiver operating characteristic curve (AUC) values above 0.9. This finding opens an avenue for the use of EV-associated miRNAs of peritoneal lavages as an untapped source of biomarkers for EC.


Introduction
Endometrial cancer (EC) is the sixth most common cancer in women worldwide and is responsible for more than 89,000 deaths every year [1]. EC is predominantly a disease that afflicts postmenopausal women, occurs in women older than 50 years in more than 90% of cases, and is detected at a mean age of 65 [2]. Approximately 10% of cases are diagnosed in premenopausal women, 5% of whom are younger than 40 years. Mortality is associated with presence of poor prognostic factors at diagnosis, i.e., diagnosis at an advanced stage, with a high grade and/or an aggressive histology. Patients presenting any of those features are at increased risk of recurrence, and, for them, therapeutic options are limited. Although most ECs are diagnosed early, up to 10% of tumors are diagnosed at a late stage, where the five-year survival drops to 16% compared to 95% in women diagnosed at an early stage [3]. Regarding histology, 80% of EC patients are diagnosed with an endometrioid histology with the average 5-year survival of 75%. Nonetheless, the 20% of patients diagnosed with a non-endometrioid histology account for 47% of EC-related deaths. Grade 3 endometrioid tumors are diagnosed in 15% of all EC patients, although theses tumors are responsible for 27% of EC-related deaths.
The cornerstone treatment of EC is surgery, which is mostly standardized through all hospitals with slight variations. The national comprehensive cancer network (NCCN) recommends both surgical and pathological staging with total hysterectomy, bilateral salpingo-oophorectomy, and peritoneal cytology. Although lymphadenectomy still remains the most reliable way of avowing downstaging and correctly identifying patients who require adjuvant chemotherapy or radiotherapy, recent data have questioned its role in early stage EC due to the high variability in node involvement and the frequently associated comorbidities. Lymphadenectomy might not be performed in early stage EC patients in some centers. Moreover, fertility-sparing surgery (FSS) in reproductive-age patients affected by endometrial cancer has attracted attention in the last decade because the consequences of an approach that is too radical may have a severe impact on a patient's quality of life and psychological well-being [4].
The new era of molecular advances has fostered biomarker research, although the identification of useful biomarkers in liquid biopsies remains a challenge. Among all serum biomarkers, the human epididymis protein 4 (HE4) has been one of the most investigated in EC. HE4 was found to be sufficiently specific but poorly sensitive in patients with EC. The diagnostic performance of HE4 appears to be better than that of the cancer antigen 125 (CA125) in diagnosing EC at an early stage, but its real value and efficacy for management of EC have not been clearly demonstrated in clinical practice [5]. Hence, a clinical challenge in EC is the development of novel molecular approaches to liquid biopsies that permit early diagnosis and recurrence control.
MicroRNAs (miRNAs) are small, non-coding RNA molecules of about 22 nucleotides in length [6] that regulate gene expression at the post-transcriptional level by inhibiting protein translation or destabilizing target transcripts via binding to the 3'-untranslated region (3'UTR), resulting in transcriptional repression or mRNA degradation upon dicer cleavage [7]. miRNAs have been found to play a critical role in almost every physiological process, including differentiation, proliferation, and apoptosis. They have also been described as oncogenes or tumor suppressors in some tumors [8], including EC [6]. Although miRNAs are detected intracellularly, they pass into the extracellular space and can be detected in a broad variety of bodily fluids, either freely in circulation or contained in extracellular vesicles (EVs) [9].
There are different types of EVs, and their size ranges from 20 to 200 nm. The largest vesicles, i.e., microvesicles, are released directly from the budding of the plasma membrane. The smallest vesicles, i.e., exosomes, are formed within intracellular multivesicular bodies and released by their fusion with the cellular membrane. Their function is to mediate intercellular communication, influencing the recipient cell's behavior. Importantly, EVs have attracted the interest of the scientific community as a source of biomarkers, mainly because they carry a broad range of bioactive material (proteins, metabolites, RNA, miRNA, etc.) that is well-protected by the lipid bilayer membrane of EVs, even if they are extracted from circulating or proximal bodily fluids or frozen before any experimental study [10].
Herein, we investigate the use of EVs isolated from the peritoneal lavage, a proximal fluid of EC, as a source of potential EC biomarkers. The peritoneal lavage, at surgery, just before starting the manipulation of the uterus, was used for staging purposes according to the old International Federation of Gynecology and Obstretrics (FIGO) staging rules. In several centers, peritoneal washing is still performed because of the prognostic information that the presence of cancer cells provides by cytologic examination. However, this fluid has not been used for molecular analysis. In this study, we conducted miRNA profiling of EVs isolated from the peritoneal lavages of 25 EC patients and the ascitic fluids of 25 non-cancer patients using the TaqMan OpenArray Human MicroRNA Panel. We identified the most relevant individual miRNAs related to EC and characterized the biological and molecular landscape of the EC milieu. The study was conceived as a proof of concept investigation to demonstrate the feasibility of using the peritoneal lavage as a source of EV-associated miRNA biomarkers of EC.

Results
We analyzed the miRNA profile of EVs isolated from the ascitic fluids of 25 control patients and the peritoneal lavages of 25 EC patients. Figure 1 illustrates the workflow that was followed in this study. Herein, we investigate the use of EVs isolated from the peritoneal lavage, a proximal fluid of EC, as a source of potential EC biomarkers. The peritoneal lavage, at surgery, just before starting the manipulation of the uterus, was used for staging purposes according to the old International Federation of Gynecology and Obstretrics (FIGO) staging rules. In several centers, peritoneal washing is still performed because of the prognostic information that the presence of cancer cells provides by cytologic examination. However, this fluid has not been used for molecular analysis. In this study, we conducted miRNA profiling of EVs isolated from the peritoneal lavages of 25 EC patients and the ascitic fluids of 25 non-cancer patients using the TaqMan OpenArray Human MicroRNA Panel. We identified the most relevant individual miRNAs related to EC and characterized the biological and molecular landscape of the EC milieu. The study was conceived as a proof of concept investigation to demonstrate the feasibility of using the peritoneal lavage as a source of EV-associated miRNA biomarkers of EC.

Results
We analyzed the miRNA profile of EVs isolated from the ascitic fluids of 25 control patients and the peritoneal lavages of 25 EC patients. Figure 1 illustrates the workflow that was followed in this study. The quality of EVs isolated from the ascitic fluids and peritoneal lavages was measured by size distribution and concentration by Nanoparticle Tracking Analysis. The analysis demonstrated that we analyzed a population that was mostly comprised of small EVs but that also contained a low number of microvesicles. The EVs isolated from EC and non-EC patients did not differ in concentration and size (Supplementary Figure S1). miRNAs were extracted from all EVs and they were used for a systematic miRNA expression analysis using the Taqman OpenArray technology. We detected 371 out of the 754 miRNAs (49.2%) present in the OpenArray. Probes that had a Cycle threshold (Ct) value of 40 in all samples and samples in which more than 80% of the probes had a Ct value greater than 40 were removed, resulting in a study that contained a total of 355 miRNAs from 22 control and 22 EC patients (Table 1). The quality of EVs isolated from the ascitic fluids and peritoneal lavages was measured by size distribution and concentration by Nanoparticle Tracking Analysis. The analysis demonstrated that we analyzed a population that was mostly comprised of small EVs but that also contained a low number of microvesicles. The EVs isolated from EC and non-EC patients did not differ in concentration and size (Supplementary Figure S1). miRNAs were extracted from all EVs and they were used for a systematic miRNA expression analysis using the Taqman OpenArray technology. We detected 371 out of the 754 miRNAs (49.2%) present in the OpenArray. Probes that had a Cycle threshold (Ct) value of 40 in all samples and samples in which more than 80% of the probes had a Ct value greater than 40 were removed, resulting in a study that contained a total of 355 miRNAs from 22 control and 22 EC patients (Table 1). Clinical characteristics of the final cohort of patients included in the study after data normalization.
The differential expression analysis between cancer and control cases yielded a list of 114 miRNAs that were significantly dysregulated (adj. p < 0.05 and abs(logFC) > = 1). Among those, 96 miRNA were found to be downregulated and 18 miRNA were upregulated in EC patients ( Table 2). To evaluate whether these miRNAs can be used as biomarkers, we performed a predictive analysis using the logistic modeling. Eight miRNAs demonstrated predictive performance with area under the receiver operating characteristic curve (AUC) values above 0.90, including miRNA-383-5p, miRNA-10b-5p, miRNA-34c-3p, miRNA-449b-5p, miRNA-34c-5p, miRNA-200b-3p, miRNA-2110, and miRNA-34b-3p (Table 2, rows in bold; Figure 2). All eight miRNAs were significantly downregulated in EC patients (from 3.75-fold to 12.18-fold in the log2_scale).    To further understand the milieu generated by EVs in the context of EC, we performed a bioinformatics study to first identify the predicted transcripts that are regulated by all of the differentially expressed miRNAs and then to assess the biological processes and molecular functions that they participate in. A total of 8074 transcripts were found to be regulated by the 114 differentially expressed miRNAs. To comprehensively integrate the properties of all target transcripts, they were classified with the Gene Ontology (GO) terms shown in Figures 3A and 3B. To further understand the milieu generated by EVs in the context of EC, we performed a bioinformatics study to first identify the predicted transcripts that are regulated by all of the differentially expressed miRNAs and then to assess the biological processes and molecular functions that they participate in. A total of 8074 transcripts were found to be regulated by the 114 differentially expressed miRNAs. To comprehensively integrate the properties of all target transcripts, they were classified with the Gene Ontology (GO) terms shown in Figure 3A  Biological process GO analysis of predicted proteins regulated by differentially expressed miRNA from the ELV of EC and control patients. A total of 8074 genes were predicted to be modulated by 114 miRNA and were included in at least one of the GO Biological process categories as indicated in the pie chart. Others include the following categories that have less than 0.4% of representation: growth; multi-organism process; biological phase; rhythmic process; pigmentation; nitrogen utilization. (B) Molecular Function GO analysis of predicted proteins regulated by differentially expressed miRNA from the ELV of EC and control patients. A total of 8074 genes were predicted to be modulated by 115 miRNA and were included in at least one of the GO Molecular Function categories as indicated in the pie chart. Abbreviations: Reg, Regulation; Pr, Process; Bd, Binding; Act, Activity. Biological process GO analysis of predicted proteins regulated by differentially expressed miRNA from the ELV of EC and control patients. A total of 8074 genes were predicted to be modulated by 114 miRNA and were included in at least one of the GO Biological process categories as indicated in the pie chart. Others include the following categories that have less than 0.4% of representation: growth; multi-organism process; biological phase; rhythmic process; pigmentation; nitrogen utilization. (B) Molecular Function GO analysis of predicted proteins regulated by differentially expressed miRNA from the ELV of EC and control patients. A total of 8074 genes were predicted to be modulated by 115 miRNA and were included in at least one of the GO Molecular Function categories as indicated in the pie chart. Abbreviations: Reg, Regulation; Pr, Process; Bd, Binding; Act, Activity.

Discussion
In this study, we investigated for the first time the miRNA content of EVs isolated from peritoneal lavages and ascitic fluids of EC and control patients, respectively. Our study shows that the EV-associated miRNAs can be consistently extracted from those proximal bodily fluids and that the miRNA expression profiles can indicate and represent the status of EC patients. The EV-associated miRNAs were analyzed using the Taqman OpenArray technology. The differential expression analysis yielded 114 miRNAs that were significantly dysregulated in EC patients.
An abundance of scientific research has been published regarding the role of miRNAs in EC [11]. Torres et al. published the first study focused on miRNA expression both in tissue and plasma samples of EC patients. They investigated the expression of miRNA-99a, miRNA-100, and miRNA-199b, which target the mTOR kinase. A combined signature of miRNA-99a and miRNA-199b in plasma samples resulted in 88% sensitivity and 93% specificity, indicating a good diagnostic potential [12]. Despite these findings, they were not applied in the clinical setting [13]. In this respect, EVs arise as a source of biomarkers with an unexploited potential. They can be isolated from bodily fluids, such as saliva, blood, urine, malignant pleural effusion, and ascitis [14]. In EC, miRNAs isolated from EVs have been scarcely studied. Akhil et al. evaluated the potential of the miRNA content of urine-derived EVs as a diagnostic biomarker in EC patients [15], and Hanzi Xu et al. isolated EVs from serum samples and identified 209 upregulated and 66 downregulated circular RNA (circRNAs) in EVs from serum of patients with EC compared with those from serum of healthy controls [16].
Although plasma, serum, and urine biopsies are the most common liquid biopsies, the use of proximal bodily fluids as a source of biomarkers has attracted the attention of the biomarker research community. Proximal bodily fluids, such as uterine fluid for EC, offer an improved representation of the molecular alterations that take place in the tumor [17]. The peritoneal fluid is another proximal fluid of EC; however, this type of proximal fluid has not been yet exploited to investigate EC-related biomarkers or any other cancer originating within the peritoneal cavity.
To the best of our knowledge, our study is the first to report the value of this proximal fluid for the identification of miRNAs associated with EVs in EC. Importantly, this study identified the dysregulation of 114 miRNAs, among which miRNA-383-5p, miRNA-10b-5p, miRNA-34c-3p, miRNA-449b-5p, miRNA-34c-5p, miRNA-200b-3p, miRNA-2110, and miRNA-34b-3p are of special interest, as they demonstrated a high classification potential. Interestingly, some of these miRNAs were found in previous EC studies. In concordance with our study, miRNA-10b and miRNA-34b were found to be downregulated in endometrial serous adenocarcinoma versus normal endometrial tissue, indicating that these miRNAs might also be associated with the aggressive subtype of the serous EC. In fact, in that study, reduced miRNA-10b expression was found to be significantly correlated with shorter overall survival [18]. MiRNA-34 has been described as a fundamental regulator of tumor suppression; it controls multiple protein targets involved in the cell cycle and apoptosis, and was associated with metastasis and chemoresistance [19]. In contrast to our study, miRNA-200b was found to be upregulated in endometrial serous adenocarcinoma versus normal endometrial tissue [18]. However, it has been reported to be downregulated in various human malignancies, and its function has been postulated to be oncogenic (i.e., involved in proliferation, motility, apoptosis, stemness, and the epithelial-to-mesenchymal transition) [20]. Interestingly, accumulating evidence in the field of endometriosis suggests that apoptosis that occurs in the peritoneal cavity may play a pivotal role in addressing the immune homeostasis in the peritoneal microenvironment [21]. This causes scavenging mechanisms to fail, allowing for the survival of endometriotic cells in patients with endometriosis [22]. In the context of endometrial cancer, we speculate that EVs derived from endometrial cancer cells might reach the peritoneal cavity and target the mesothelial liner cells of the peritoneum and ovaries to modulate the immune homeostasis and create a more favorable milieu to metastasize. In fact, most of the metastasis associated with endometrial cancer occurs either in the vagina, in the lymph nodes, or within the peritoneal cavity. Several of the miRNAs identified in our study are related to tumor progression in endometrial cancer or in other tumor types. It is possible that endometrial cancer cells may spread into the peritoneum and the ovaries through the Fallopian tube, and develop peritoneal metastasis and ovarian metastasis in the absence of lympho-vascular space invasion. It is important to note that normal endometrium may spread into the peritoneum and the ovaries through retrograde menstruation, and give rise to endometriosis. Metastasis that originates through this pathway may be associated with indolent behavior. A subset of ECs from ovarian metastasis are associated with such a good prognosis that they were interpreted as synchronous tumors [23]. Next-generation sequence analysis has confirmed that they are metastatic. Assessing transtubal exosomal release from endometrial cancer may help us to understand the mechanisms involved in this type of indolent metastasis.
The study leaves some questions open. This study enabled us to identify a large number of dysregulated miRNAs associated with EVs in the peritoneal lavage of EC patients. We think that the EVs that we have identified come from endometrial and mesothelial cells in the EC group, and mostly from mesothelial cells in the control group. The comparative analysis of EC patients with non-cancer patients with ascites suggests that these selected miRNAs come from EVs of EC tissue. However, there is obviously the probability that a subset of EVs came from an inflammatory reaction associated with EC and some other factors distinct from the EC pathology, such as the source of EVs (peritoneal lavage vs. ascitic fluids), the surgery (after induction of general anesthesia), and the control samples, which were obtained by means of paracentesis (with local anesthesia). Nevertheless, these promising biomarkers should be further validated as well as combined in order to increase the already excellent accuracy of each individual miRNA. This should be done in an independent study involving a larger cohort of EC patients versus a control group with a higher biologic variability, including, for example, patients with leiomyomas or women requiring tubal ligation for definitive contraception. Although we tested whether or not differentially expressed miRNAs were dependent on gender, further studies should include only female controls. Moreover, further research should be directed to an evaluation of the prognostic potential of each specific dysregulated miRNA, as this might help to guide the surgical treatment of EC patients.

Patients and Ascitic Fluid and Peritoneal Lavage Collection
All subjects provided informed consent before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by The Clinical Research Ethics Committee of Hospital Arnau de Vilanova in Lleida, Spain (Approval number: CEIC-1630). Samples were obtained with support from the IRBLleida Biobank (B.0000682) and Plataforma biobancos PT17/0015/0027. Ascitic fluids and peritoneal lavages were extracted from a cohort of 50 patients, corresponding to 25 control patients with decompensated cirrhosis and 25 patients with EC who underwent curative surgery. In the control patients, the collection of ascitic fluids was performed as follows: Ascitic fluids were aspirated using 18 or 21G needles (for diagnostic paracentesis) or an over-the-needle catheter device (for therapeutic paracentesis). The procedure was performed under sterile conditions, the needle insertion site was selected by ultrasound guidance, and the skin and parietal peritoneum were previously anesthetized with 2% mepivacine. A total of 100 mL of ascitic fluid was gently aspirated, collected into a 50 mL tube, and stored at −80 • C until use. In EC patients, the collection of peritoneal lavage was performed during surgery, once the abdominal cavity was opened and prior to any manipulation of the uterus. A total of 100 mL of physiological saline was instilled into the abdominal cavity with a 50 mL syringe, mobilizing patients for the correct distribution of saline, which was then extracted with a 50 mL syringe connected to a 14-gauge aspiration needle. The peritoneal lavage was gently aspirated. A volume ranging from 50 to 100 mL was collected and stored at −80 • C until use. The clinical features of each patient are listed in Supplementary Table S1.

EV Isolation
EVs were isolated with a differential centrifugation method as previously described [24] with slight modifications. Briefly, ascitic fluids and peritoneal lavages were centrifuged at 300× g for 10 min, followed by centrifugation at 2500× g for 20 min at the moment that the sample was collected, and frozen at −80 • C. Then, samples were centrifuged at 10,000 g for 30 min (Thermo Scientific Heraeus MultifugeX3R Centrifuge (FiberLite rotor F15-8x-50c)). The supernatant was filtered through 0.22 µm filters (Merck Millipore), and the obtained sample was transferred to ultracentrifuge tubes (Beckman Coulter), which were filled with phosphate-buffered saline (PBS), to perform two consecutive ultracentrifugation steps at 100,000 g for 2 hours each on a Thermo Scientific Sorvall WX UltraSeries Centrifuge with an AH-629 rotor. The pellet containing the EVs was resuspended in 50 µL of PBS. From those, 5 µL were isolated for nanoparticle tracking analysis (NTA) and quantification, and the rest was frozen at −80 • C with 500 µL of Qiazol for RNA extraction.

Nanoparticle Tracking Analysis
The size and number of EVs were determined using a Nanosight LM10 instrument equipped with a 405 nm laser and a Hamamatsu C11440 ORCA-Flash 2.8 camera (Hamamatsu) with Nanoparticle Tracking Analysis (NTA, Malvern Instruments, UK). Each sample was diluted appropriately with Milli-Q water (Milli-Q Synthesis, Merck Millipore, MA, USA) to give counts in the linear range of the instrument. The particles in the laser beam underwent Brownian motion, and a video was recorded for 60 s in triplicate. The analysis was performed by following the manufacturer's instructions, and data were analyzed using version 2.3 of the NTA software.

Total RNA Extraction and OpenArray Analysis
Total RNA, including miRNAs and other RNAs, was isolated from the EV samples using the miRNeasy MiniKit (Qiagen, Hilden, Germany) according to the manufacturer's protocol. RNA from EVs was eluted with 30 µL of Nuclease-free water (Qiagen, Hilden, Germany). MiRNA expression was determined using a TaqMan OpenArray Human MicroRNA Panel, QuantStudio 12K Flex (Catalog number: 4470187, Thermo Fisher Scientific, Waltham, MA, USA), a fixed-content panel containing 754 well-characterized human miRNA sequences from the Sanger miRBase v14, and according to the manufacturer's instructions. Reverse transcription (RT) was performed on 2 µL RNA using Megaplex™ Primer Pools A and B and the supporting TaqMan ® MicroRNA Reverse Transcription Kit as follows: 15 min at 42 • C and 5 min at 85 • C. Then, 5 µL of the resulting cDNA was preamplified prior to real-time PCR analysis using Megaplex™ PreAmp Pools and the TaqMan ® PreAmp Master Mix using the following conditions: one single step at 95 • C for 5 min, 20 cycles of a two-step program (3 sec, 95 • C and 30 sec, 60 • C) followed by a single cycle of 10 min at 99 • C to inactivate the enzyme. The preamplified products were diluted to 1:20 in 0.1× TE buffer pH 8.0, and mixed in a 1:1 ratio with TaqMan ® OpenArray ® Real-Time PCR Master Mix in the 384-well OpenArray ® Sample Loading Plate. TaqMan ® OpenArray ® MicroRNA Panels were automatically loaded using the AccuFill™ System.

Preprocessing and Differential Expression Analysis
All bioinformatics analyses were performed with BioConductor (version 3.7) [25] in the R statistical environment (version 3.5.0) [26]. For the data preprocessing, the HTqPCR (version 1.34) R package [27] was used. Probes that had a "Cycle threshold" (Ct) value of 40 in all samples were removed. Samples in which more than 80% of the probes had a Ct value greater than 40 were retained. To assure comparability across samples, the Ct values were delta-normalized. The average Ct values of the probes hsa−miR−150−5p, hsa−let−7g-5p, hsa−miR−598−3p, and hsa−miR−361−3p were used for normalization. These probes had Ct values of 40 in a maximum of three samples and the lowest interquartile range across all samples. The differential expression analysis was carried out with the empirical Bayes approach on linear models using the limma (version 3.36) R Package [28]. Results were corrected for multiple testing using the False Discovery Rate (FDR) [29].

Development of Predictors
For predictive analysis, the whole patient cohort was randomly divided into training and validation sets with a ratio of 3:2. Calculated (with the limma R Package) relative miRNA expression values were used as input variables into a logistic regression model between groups. Each miRNA (adjusted p-value <0.05) was fitted in the logistic regression model to differentiate the EC and the control patient groups in the training set, and its classification ability was evaluated using the area under the ROC curve (AUC), accuracy, sensitivity, and specificity values on the validation set. The procedures for the division of the patient cohort into training and validation sets and fitting the logistic model were repeated 500 times and statistics were collected.

miRNA Target Gene Prediction and Bioinformatics Analysis
miRNA target genes were predicted using the Predictive Target Module of the miRWalk2.0 online software [30] (https://goo.gl/ajG9ja). Only genes for which miRNAs recognize a minimum 7 bp seed length, seed start at position 1, and sequence localize at 3´UTR were considered as valid targets. Moreover, to improve target gene prediction accuracy, we considered only those transcripts that were predicted in at least eight out of the 12 databases (miRWalk, miRanda, MicroT4, miRDB, miRMap, miRBridge, miRNAMap, PICTAR2, RNA22, PITA, TargetScan, and RNAhybrid) presented in the miRWalk2.0 tool.
The online Panther software [31] (http://www.pantherdb.org/) was used for the Gene Ontology (GO) functional analysis to analyze the potential functions of the predicted target genes. Biological process (BP) and molecular function (MF) GO terms were analyzed and plotted.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2072-6694/11/6/839/s1, Figure S1: EV characterization. (A) Box-plot representing the average mode of EVs isolated from the peritoneal lavage and the ascitic fluid of EC and control patients, respectively (Mean ± SD), measured by Nanoparticle Tracking Analysis. (B) Size distribution and concentration of isolated EVs of a peritoneal lavage of an EC patient (left) and an ascitic fluid of a control patient (right), measured by Nanoparticle Tracking Analysis. Table S1: Clinicopathological characteristics of all patients. Clinical characteristics of the total cohort of patients recruited in the study. Acknowledgments: The authors would like to acknowledge the work of all who have participated in the recruitment of clinical samples. We also thank the patients for their willingness to participate in the study.

Conflicts of Interest:
The authors declare no conflicts of interest.