Identification of non-invasive miRNAs biomarkers for prostate cancer by deep sequencing analysis of urinary exosomes
Molecular Cancer volume 16, Article number: 156 (2017)
The aim of this study was to identify microRNAs in urinary exosomes that are differently expressed in prostate cancer patients and healthy donors. For this purpose, RNA was extracted from urinary exosomes from 20 prostate cancer patients and 9 healthy males and the microRNAs were analyzed by next generation sequencing. Interestingly, 5 microRNAs – miR-196a-5p, miR-34a-5p, miR-143-3p, miR-501-3p and miR-92a-1-5p – were significantly downregulated in exosomes from prostate cancer patients. Furthermore, RT-qPCR analysis of an independent cohort of 28 prostate cancer patients and 19 healthy males confirmed that miR-196a-5p and miR-501-3p were downregulated in prostate cancer samples. These results suggest that specific microRNAs in urinary exosomes might serve as non-invasive biomarkers for prostate cancer. In particular, miR-196a-5p and miR-501-3p are promising biomarkers that need to be further studied in large patient cohorts.
Prostate cancer (PCa) accounts for 12% of all cancer cases worldwide, and it is the second most commonly diagnosed cancer in men . There is an urgent need for early diagnosis and prognosis biomarkers for the disease . Biomarkers that can be measured in biofluids, often referred to as liquid biopsies , avoid the discomforts and may reflect the heterogeneity of the tumor better than prostate needle biopsies. Remarkably, exosomes released by cancer cells into biological fluids contain molecules that reflect the disease status, and are considered today as a new type of liquid biopsies [4, 5]. In particular, exosomes contain microRNAs (miRNAs) , which are small noncoding RNAs involved in the regulation of gene expression. miRNAs have been show to play several functions in cancer, and specific miRNA signatures have been identified in several cancer types [7, 8]. Interestingly, cancer-specific miRNAs have been identified in exosomes , and several exosomal miRNAs may serve as PCa biomarkers [10,11,12,13]. We have here investigated the use of next generation sequencing to identify novel miRNA-based PCa biomarkers in urinary exosomes and used RT-qPCR to validate the results.
Isolation of exosomal RNA from urine samples
Urinary exosomes were isolated by sequential centrifugation as previously described . We have previously characterized these vesicles by several methods such as Western blot, electron microscopy, and mass spectrometry-based proteomics . Based on these analyses we concluded that our samples are enriched in exosomes, but we cannot exclude that other extracellular vesicles co-isolate with exosomes to some extent. A material corresponding to 2.5–10 μg exosomal proteins (as measured with the BCA assay) was used for RNA isolation. Before RNA was isolated, the samples were treated with proteinase K and RNAse A to degrade potential RNA/protein complexes. Several methods have been described to isolate RNA from urinary exosomes , and in this study mirCury, miRNeasy and trizol were tested. After RNA isolation, the samples were treated with DNAse to degrade potential DNA. The RNA concentrations obtained were too low to be measured by Nanodrop, and therefore the samples were analyzed in Agilent RNA 6000 Pico chips. The miRNeasy Micro RNA isolation kit gave the highest RNA yield, 1–3 ng RNA from a material corresponding to 4–5 μg of exosomal protein. As shown in a typical Bioanalyzer RNA profile (Additional file 1: Figure S1), no evidence of 18S and 28S RNA was observed.
Next generation sequencing of urinary exosomes from PCa patients and healthy controls
In this study, RNA isolated from urinary exosomes from 20 PCa patients and 9 healthy male donors was used (Table 1). In order to optimize sequencing performance, the libraries were built using specific adaptors designed for small amounts of starting material (lower limit 1 ng total RNA). The libraries were PCR amplified, pooled and size selected and NGS was performed using Illumina high-throughput RNA sequencing technology. In average 5.8 million raw reads per sample were obtained and approximately 60% of the reads were mapped to the genome (Additional file 1: Table S1), mainly to miRNA (24%) and mRNA (29%) (Fig. 1a). The rest of the reads mapped to scRNA, miscRNA, and lincRNA, and rRNA, snoRNA, snRNA and tRNA together made 0.6% of the reads (Fig. 1a). The miRNA read counts normalized to reads per million (RPM) were used in further analyses.
As shown in Fig. 1b, the most abundant miRNAs was miR-10b-5p, followed by let-7b-5p, miR-30a-5p, miR-10a-5p and let-7a-5p. miR-10b-5p has also been shown to be very abundant in other studies of urinary exosomes [15, 16]. miRNAs with less than 10 reads per million in less than 5 samples were filtered out, leaving 254 miRNAs (Additional file 1: Table S2), or 217 miRNAs if miRNAs from different precursors were considered as one. In total 80 miRNA with 10 or more reads were found in all the samples. Hierarchical clustering analysis showed that the samples did not cluster into specific groups (Additional file 1: Figure S2). However, the level of 5 miRNAs was significantly different in control versus PCa patients: miR-196a-5p, miR-34a-5p, miR-143-3p, miR-501-3p and miR-92a-1-5p (Table 2, Fig. 1c, d, Additional file 1: Figure S3). All of these miRNAs were shown to be down-regulated in urinary exosomes from PCa patients. The most promising miRNA, miR-196a, was able to distinguish the two groups with 89% specificity and 100% sensitivity. Receiver operating characteristic (ROC) curve analysis showed that the area under the curve (AUC) for this miRNA was 0.92 (95% CI 0.79–1.06) and for miR143-3p was 0.72 (95% CI 0.48–0.97) (Additional file 1: Figure S4), thus showing the diagnostic potential of these miRNAs for PCa. Finally, no relevant differences were observed when PCa patients were subdivided by GS (GS 7a and GS 7b) or as intermediate and aggressive based on the D’Amico classification (see as an example miR-196a in Additional file 1: Figure S5).
All of the 5 miRNAs shown to be significantly altered in PCa samples versus healthy controls in this study were downregulated in urinary exosomes of PCa patients. It is not clear whether the downregulation of these miRNAs is due to miRNA downregulation in the tumor or to nonspecific responses to the presence of the tumor. It has been suggested that, at least at early tumor stages, a reduction in the level of a circulating miRNA is probably not the result of a down regulation of the miRNA in the tumor itself, but to a negative effect of the tumor in the expression of the miRNA in other cells . In a previous NGS of urinary exosomes, the majority of the miRNAs significantly altered between the control and the PCa group were also found to be downregulated, but the entire list of altered miRNA was not provided . It is not clear why more miRNAs were altered between the control and the PCa group in that study compared to our study, but the different patient cohorts used or different NGS pipelines could be possible explanations. Interestingly, in the above mentioned study, PCR validation studies of 3 specific miRNAs (miR-21, miR-204 and miR-375) showed that only isomiRs (miRNA isoforms) with 3′ end modifications, but not the mature miRNAs, were able to discriminate between controls and PCa patients . Since we only used the mature forms in our experimental PCR settings, it would be interesting to analyze potential differences in miRNAs isoforms in future studies. Finally, it should be mentioned that urine contains several population of exosomes, and that it is not clear to which extent prostate cells contribute to the total population. It is possible that if prostate-derived exosomes are isolated and analyzed separately, both upregulated and downregulated miRNAs are found in PCa patients.
Validation of the RNA-seq data by quantitative PCR analysis
To validate the miRNAs identified by NGS, they were analyzed by quantitative RT-PCR on a new patient set including urine samples from 28 PCa patients and 19 healthy male donors (Table 1). For normalization we used the average of three miRNAs/small RNAs; miR-10b-5p, let-7b-5p and U6 snRNA. miR-10b-5p and let-7b-5p were chosen because, as shown in the NGS study, they were similarly expressed in the control and in the patient group (Additional file 1: Figure S6), and U6 snRNA was chosen as it has previously been used for normalization in other studies . Importantly, the level of miR-196a-5p and miR-501-3p was significantly reduced in the PCa group compared to the healthy control group, thus validating the sequencing data by an additional analytical method in an independent sample set (Fig. 2a, b). ROC curves for these miRNAs showed that the AUC for miR-196-5p and miR-501-3p was 0.73 (95% CI 0.56–0.89) and 0.69 (95% CI 0.52–0.85) respectively (Fig. 2c, d). For comparison, it has been reported that the AUC for PSA to discriminate between any PCa and cancer-free controls in the Prostate Cancer Prevention Trial was 0.68 (95% CI 0.67–0.69) . A combination (sum) of the two miRNAs did not significantly improve the diagnostic potential (AUC 0.74; 95% CI 0.59 to 0.89). Finally, we could not find a significant reduction of miR-34a-5p and miR-92a-1-5p level in the PCa group by RT-qPCR (Additional file 1: Figure S7), and miR-143-3p was not detected in the experimental conditions that were used.
Based on both NGS and PCR studies, miR-196a-5p appears as the most promising biomarker. The diagnostic performance calculated from NGS results (0.92; 95% CI 0.79–1.06) was better than from PCR results (0.73; 95% CI 0.56–0.89). This is often observed when biomarkers are measured by different methodologies, but the fact that both methods show a downregulation of miR-196a-5p opens for additional studies and optimizations. miR-196a has been found to be dysregulated, mainly upregulated, in several cancer types, and to act as an oncogene by affecting cell proliferation, migration and invasion [20, 21]. In fact, not much is known about miR-196a in PCa, though it has been shown that miR-196a is upregulated in PCa tissue compared to normal tissue [22, 23], while miR-196b has been shown to be downregulated . However, it should be mentioned that exosomal miRNAs do not necessarily represent the cellular levels since miRNA sorting mechanism may affect the incorporation of miRNA into exosomes.
Among the numerous targets of miR-196a we find Radixin, ANXA1, several HOX genes and ETS-related gene (ERG). Remarkably, ERG is the most frequently overexpressed oncogene in PCa, and high expression of ERG is associated with advanced tumour stage, high Gleason score, metastasis and shorter survival times . The oncogene function of ERG is related to the regulation of a plethora of genes implicated in several cellular processes such as differentiation, growth, motility and invasion . It is possible then that the high levels of miR-196a-5p in normal prostate cells help to maintain the levels of ERG low. In terms of PCa, it is also interesting to find that analysis of miR-196a-5p in DIANA-miRPath v3.0 identified 16 target genes belonging to the KEGG pathway Prostate cancer (hsa05215) (Additional file 1: Figure S8) . However, when looking for potential targets of the downregulated miRNAs that may be relevant for PCa we should keep in mind that, as previously mentioned, the expression of the miRNAs is not necessarily changed in the tumor.
miR-501-3p was also shown to be downregulated in urinary exosomes from PCa patients both by NGS and by PCR. It has recently been shown that miR-501-3p promoted the invasiveness of pancreatic ductal adenocarcinoma cells possibly by suppressing E-cadherin . Downregulation of miR-501-3p in serum has recently been found to be useful for prediction and prognosis of lymph node metastasis in gastric cancer , but miR-501-3p has been found to be upregulated in serum samples of PCa patients compared to benign prostatic hyperplasia .
The discovery of biomarkers that can supplement PSA is a main goal of PCa research. In particular, biomarkers for early diagnosis and risk stratification are urgently needed. This study confirms the feasibility of NGS analysis of small amount of exosomal RNA and shows that NGS is a valid method to identify novel RNA-based biomarkers for PCa. We have shown that miR-196a-5p and miR-501-3p have diagnostic potential for PCa. These miRNAs will have to be the subject of future studies in order to determine their specific clinical utility in additional patient cohorts. Furthermore, the functional role of these miRNA in relation to PCa is also an interesting subject that requires further study.
Area under the curve
Next Generation Sequencing
Prostate specific antigen
Receiver operating curve
Globocan 2012 - http://globocan.iarc.fr.
Silberstein JL, Pal SK, Lewis B, Sartor O. Current clinical challenges in prostate cancer. Transl Androl Urol. 2013;2:122–36.
Crowley E, Di Nicolantonio F, Loupakis F, Bardelli A. Liquid biopsy: monitoring cancer-genetics in the blood. Nat Rev Clin Oncol. 2013;10:472–84.
Torrano V, Royo F, Peinado H, Loizaga-Iriarte A, Unda M, Falcón-Perez JM, et al. Vesicle-MaNiA: extracellular vesicles in liquid biopsy and cancer. Curr Opin Pharmacol. 2016;29:47–53.
Fais S, O’Driscoll L, Borras FE, Buzas E, Camussi G, Cappello F, et al. Evidence-Based Clinical Use of Nanoscale Extracellular Vesicles in Nanomedicine. ACS Nano. 2016;10:3886–99.
Valadi H, Ekström K, Bossios A, Sjöstrand M, Lee JJ, Lötvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9:654–9.
Price C, Chen J. MicroRNAs in Cancer Biology and Therapy: Current Status and Perspectives. Genes Dis. 2014;1:53–63.
Hayes J, Peruzzi PP, Lawler S. MicroRNAs in cancer: biomarkers, functions and therapy. Trends Mol Med. 2014;20:460–9.
Kinoshita T, Yip KW, Spence T, Liu F-F. MicroRNAs in extracellular vesicles: potential cancer biomarkers. J Hum Genet. 2017;62:67–74.
Hessvik NP, Sandvig K, Llorente A. Exosomal miRNAs as Biomarkers for Prostate Cancer. Front Genet. 2013;4:36.
Endzeliņš E, Melne V, Kalniņa Z, Lietuvietis V, Riekstiņa U, Llorente A, et al. Diagnostic, prognostic and predictive value of cell-free miRNAs in prostate cancer: a systematic review. Mol Cancer. 2016;15:41.
Fendler A, Stephan C, Yousef GM, Kristiansen G, Jung K. The translational potential of microRNAs as biofluid markers of urological tumors. Nat Rev Urol. 2016;13:734–52.
Valentino A, Reclusa P, Sirera R, Giallombardo M, Camps C, Pauwels P, et al. Exosomal microRNAs in liquid biopsies: future biomarkers for prostate cancer. Clin Transl Oncol. 2017;19:651–7.
Øverbye A, Skotland T, Koehler CJ, Thiede B, Seierstad T, Berge V, et al. Identification of prostate cancer biomarkers in urinary exosomes. Oncotarget. 2015;6:30357–76.
Cheng L, Sun X, Scicluna BJ, Coleman BM, Hill AF. Characterization and deep sequencing analysis of exosomal and non-exosomal miRNA in human urine. Kidney Int. 2014;86:433–44.
Koppers-Lalic D, Hackenberg M, de Menezes R, Misovic B, Wachalska M, Geldof A, et al. Non-invasive prostate cancer detection by measuring miRNA variants (isomiRs) in urine extracellular vesicles. Oncotarget. 2016;7:22566–78.
Witwer KW. Circulating microRNA biomarker studies: pitfalls and potential solutions. Clin Chem. 2015;61:56–63.
Pesta M, Klecka J, Kulda V, Topolcan O, Hora M, Eret V, et al. Importance of miR-20a expression in prostate cancer tissue. Anticancer Res. 2010;30:3579–83.
Thompson IM, Ankerst DP, Chi C, Lucia MS, Goodman PJ, Crowley JJ, et al. Operating characteristics of prostate-specific antigen in men with an initial PSA level of 3.0 ng/ml or lower. JAMA. 2005;294:66–70.
Chen Z-Y, Chen X, Wang Z-X. The role of microRNA-196a in tumorigenesis, tumor progression, and prognosis. Tumor Biol. 2016;37:15457–66.
Chen C, Zhang Y, Zhang L, Weakley SM, Yao Q. MicroRNA-196: critical roles and clinical applications in development and cancer. J Cell Mol Med. 2011;15:14–23.
Volinia S, Calin GA, Liu C-G, Ambs S, Cimmino A, Petrocca F, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci. 2006;103:2257–61.
Ambs S, Prueitt RL, Yi M, Hudson RS, Howe TM, Petrocca F, et al. Genomic profiling of microRNA and messenger RNA reveals deregulated microRNA expression in prostate cancer. Cancer Res. 2008;68:6162–70.
Fuse M, Kojima S, Enokida H, Chiyomaru T, Yoshino H, Nohata N, et al. Tumor suppressive microRNAs (miR-222 and miR-31) regulate molecular pathways based on microRNA expression signature in prostate cancer. J Hum Genet. 2012;57:691–9.
Adamo P, Ladomery MR. The oncogene ERG: a key factor in prostate cancer. Oncogene. 2016;35:403–14.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43:W460–6.
Ling Q, Xu X, Ye P, Xie H, Gao F, Hu Q, et al. The prognostic relevance of primary tumor location in patients undergoing resection for pancreatic ductal adenocarcinoma. Oncotarget. 2017;8:15159–67.
Jiang X, Wang W, Yang Y, Du L, Yang X, Wang L, et al. Identification of circulating microRNA signatures as potential noninvasive biomarkers for prediction and prognosis of lymph node metastasis in gastric cancer. Oncotarget. 2017;8:65132–42.
Haldrup C, Kosaka N, Ochiya T, Borre M, Høyer S, Orntoft TF, et al. Profiling of circulating microRNAs for prostate cancer biomarker discovery. Drug Deliv Transl Res. 2014;4:19–30.
We would like to acknowledge Ivana Malovic for excellent technical assistance, Solveig Mjelstad Olafsrud and Erik Høye for help with statistics and data analysis and Edgars Endzeliņš for literature review. We are grateful for the sequencing services provided by the Helse Sør-Øst Genomics Core Facility at Oslo University Hospital. We would also like to thank the Prostate Biobank at Oslo University Hospital and to the persons that donated urine samples for this study.
This work was supported by The Norwegian Financial Mechanism 2009–2014 under Project Contract number NFI/R/2014/045 and The Norwegian Cancer Society. B.F. was supported by the South-Eastern Norway Regional Health Authority (grant number 2014041).
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Ethics approval and consent to participate
The collection of urine samples was approved by the Norwegian Regional Committees for Medical and Health Research Ethics and the participants gave informed written consent.
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Rodríguez, M., Bajo-Santos, C., Hessvik, N.P. et al. Identification of non-invasive miRNAs biomarkers for prostate cancer by deep sequencing analysis of urinary exosomes. Mol Cancer 16, 156 (2017). https://doi.org/10.1186/s12943-017-0726-4