Abstract
Yi-Fan Xu1, Bethany N. Hannafon1, Yan D. Zhao2, Russell G. Postier3 and Wei-Qun Ding1
1Department of Pathology, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, OK 73104, USA
2Department of Biostatistics and Epidemiology, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, OK 73104, USA
3Department of Surgery, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, OK 73104, USA
Correspondence to:
Wei-Qun Ding, email: [email protected]
Keywords: exosome, microRNAs, pancreatic cancer, plasma biomarkers, early detection
Received: April 27, 2017 Accepted: July 28, 2017 Published: August 18, 2017
ABSTRACT
Patients with localized pancreatic cancer (stage I and stage IIA) have a much higher survival rate than those presenting at later stages, yet early detection remains a challenge to this malignancy. The aim of this study was to evaluate whether exosome miRNA signatures are indicative of localized pancreatic cancer. Exosomes were collected from the conditioned media of pancreatic cancer cell lines and plasma samples of localized pancreatic cancer patients (Stage I-IIA, n=15), and healthy subjects (n=15). Cellular and exosome miRNAs from pancreatic cancer cell lines were profiled by next-generation small RNA sequencing. Plasma exosome miRNA expression was analyzed by qRT-PCR. We found that certain miRNAs, such as miR-196a and miR-1246, are highly enriched in pancreatic cancer exosomes. Consistently, plasma exosome miR-196a and miR-1246 levels were significantly elevated in pancreatic cancer patients as compared to healthy subjects. An analysis of the cancer subtypes indicated that plasma exosome miR-196a is a better indicator of pancreatic ductal adenocarcinoma (PDAC), whereas plasma exosome miR-1246 is significantly elevated in patients with intraductal papillary mucinous neoplasms (IPMN). In contrast, there were no differences in the plasma exosome miR-196a and miR-1246 levels between patients with pancreatic neuroendocrine tumors (NET) and healthy subjects. In conclusion, we demonstrate that certain miRNA species, such as miR-196a and miR-1246, are highly enriched in pancreatic cancer exosomes and elevated in plasma exosomes of patients with localized pancreatic cancer.
INTRODUCTION
Pancreatic cancer is the fourth-leading cause of cancer deaths in the United States [1]. The 5-year survival rate for pancreatic cancer patients was 7.8% from 2006 to 2012 [2], contrasting to a 66.9% 5-year survival rate for cancer at all sites [3]. The high mortality rate is primarily attributed to late stage diagnoses, when treatment options are often limited, and to the aggressive nature of pancreatic cancer with KRAS being a master regulator [4, 5]. Statistics have shown that localized pancreatic cancer (stage I-IIA) has a much greater 5-year survival rate (29.3%) than pancreatic cancer at later stages (11% for regionally invasive and 2.6% for distant metastasis) [2, 6]. Therefore new strategies for early detection, especially for detecting pancreatic cancer at the localized stage, are critical to improving pancreatic cancer outcomes.
There are several challenges for early detection of pancreatic cancer. First, the clinical signs and symptoms of pancreatic cancer, such as weight loss, abdominal pain, and jaundice, often manifest at a late stage and are not generally specific to pancreatic cancer [7]. Second, screening by radiological imaging of pancreatic cancer in the general population is impractical due to the high cost (e.g. computed tomography) or invasive procedures (e.g. endoscopic ultrasound or endoscopic cholangiopancreatography). Third, a sensitive and specific circulating biomarker for effective screening of pancreatic cancer is not currently available. Diagnostic evaluation often relies on the patient symptoms and appropriate imaging examination, with computed tomography (CT) being the most common initial test for identifying a pancreatic mass [8]. Cancer antigen 19-9 (CA 19-9) has been utilized for pancreatic cancer detection and the plasma level of CA 19-9 has been shown to be elevated in 75-85% of pancreatic cancer patients [9, 10]. However, circulating CA 19-9 is not specific to pancreatic cancer [11]. Furthermore it has reported that only 40% of stage I pancreatic cancer patients are CA 19-9 positive [12]. New approaches for the development of circulating biomarkers to detect localized pancreatic cancer are urgently needed in order to better manage this devastating disease.
Deregulated microRNA (miRNA) expression in pancreatic cancer cells and tissues has been reported by different groups [13–17]. Because pancreatic cancer cells display different miRNA profiles compared to normal pancreatic ductal cells, circulating miRNAs have been explored as biomarkers for pancreatic cancer [18–27]. However, a miRNA signature has not been established for clinical diagnosis or screening of pancreatic cancer. A major obstacle in the process of developing miRNAs as biomarkers for pancreatic cancer is presumably the heterogeneous nature of the circulating miRNA populations. The miRNAs isolated from the circulation are undoubtedly derived from a variety of cell types and molecular complexes. This intrinsic heterogeneity likely compromises the sensitivity and specificity of circulating miRNA detection.
Exosomes are endosome-derived membrane-bound vesicles that are secreted from different cell types and present in almost all biological fluids [28, 29]. Exosomes are 40-100 nm in size, and contain proteins, lipids and miRNAs derived from their parental cell [30–32]. Studies have shown that cancer cells secrete more exosomes than normal cells, and that cancer exosomes are released into the tumor microenvironment and circulation [33–38], suggesting that exosome contents, such as miRNAs, are potential biomarkers to improve the sensitivity and specificity of cancer detection [32, 36, 39–41]. Specific to pancreatic cancer, a recent study reported that cancer specific exosomes are detectable in patient serum that can indicate the presence of pancreatic cancer [42]. However, circulating exosome miRNAs have not been well-evaluated as biomarkers for localized pancreatic cancer largely due to the highly limited availability of the clinical specimens. To determine whether exosome miRNA signatures are indicative of localized pancreatic cancer, we characterized the miRNA content of exosomes derived from pancreatic cancer cells and analyzed candidate miRNAs in plasma exosomes from patients with localized pancreatic cancer. We report here that two highly enriched miRNAs in pancreatic cancer exosomes, miR-196a and miR-1246, are potential indicators of localized pancreatic cancer and may serve as circulating biomarkers for early detection of this malignancy.
RESULTS
Isolation and characterization of pancreatic cancer cell-derived exosomes
PANC-1, MIA-PaCa-2, BxPC-3 and hTERT-HPNE cells were cultured for 3 days, and whole exosomes were isolated from the conditioned medium by ultracentrifugation. The isolated exosomes were subjected to nanoparticle tracking analysis. As shown in Figure 1, the size of the exosomes from these four cell lines were comparable, ranging around 100 nm in diameter, similar to our recent reports [43, 44]. The isolated exosomes had detectable CD63, a well-established surface marker for exosomes, as evidenced by western blot analysis (Figure 1B). The PANC-1-derived exosomes were also verified by electron microscopy (EM) (Figure 1A, inset) as we recently described [44].
Figure 1: Characterization of exosomes isolated from pancreatic cancer cell lines. Exosomes were isolated from the conditioned media of PANC-1, MIA-PaCa-2, BxPC-3 and hTERT-HPNE cells. (A) Nanoparticle tracking analysis of 50X diluted exosomes (n=3). Inset shows an exosome image observed under EM (30,000x). (B) Western blot of CD63 (30-75 kDa) under non-reducing conditions.
miRNA profiling of PANC-1- and HPNE-derived exosomes
Total RNA was isolated from PANC-1- and hTERT-HPNE-derived exosomes and a small RNA library was generated. The cDNA library was sequenced by Illumina sequencers, and the sequences were analyzed against the miRBase database and quantified as reads per million for each individual miRNA species. There were around 150 miRNAs detected in exosomes derived from both cell lines (Supplementary Table 1). By comparing miRNA levels in PANC-1 and hTERT-HPNE exosomes, we identified the top 30 miRNAs that are highly expressed in PANC-1 exosomes compared to hTERT-HPNE exosomes (Figure 2A). We also identified top 30 highly abundant miRNAs within PANC-1 exosomes (Figure 2B). We selected miR-196a, miR-196b and miR-1246 as candidate exosome miRNAs for biomarker development, due to their selective expression and high abundance in PANC-1 exosomes (Table 1). Expression of exosome miR-196a, miR-196b, and miR-1246 was verified by qRT-PCR analysis using exosomes derived from hTERT-HPNE, PANC-1, MIA-PaCa-2, and BxPC-3 cells. Significantly higher levels of these three miRNAs were evident in exosomes derived from cancer cells as compared to exosomes derived from hTERT-HPNE cells (Figure 3), consistent with the RNA sequencing results. We therefore focused on exosome miR-196a, miR-196b and miR-1246 in the following experiments.
Table 1: Expression of miR-1246, miR-196a and miR-196b in PANC-1 and hTERT-HPNE exosomes (reads per million)
miRNA | Group Mean | Group Mean | PANC-1 Exo/ |
---|---|---|---|
miR-1246 | 815.942 | 10481.01 | 12.58 |
miR-196a | 0.48568 | 1179.84 | 2429.25 |
miR-196b | 0.48568 | 944.4833 | 1944.66 |
Figure 2: Selectively expressed and highly enriched miRNAs in pancreatic cancer exosomes. Cellular and exosome RNA was isolated from hTERT-HPNE and PANC-1 cells. microRNA expression was analyzed using Small RNA next generation sequencing (n=3). (A) Top 30 miRNAs highly expressed in PANC-1 exosomes compared to hTERT-HPNE exosomes (fold of changes). (B) Top 30 miRNAs highly abundant in PANC-1 exosomes (reads per million).
Figure 3: qRT-PCR analysis of miR-1246, miR-196a and miR-196b expression in PANC-1, MIA-PaCa-2, and BxPC-3 exosomes versus hTERT-HPNE exosomes. Total RNA was extracted from the exosomes and miR-196a, miR196b and miR-1246 expression was measured by qRT-PCR. Fold-change in expression is shown for miR-1246, miR-196a and miR-196b relative to their levels in hTERT-HPNE exosomes and normalized to the spike-in cel-miR-54 control (n=3, means ± SEM). *p<0.05; **p<0.01; Student’s t-test.
Plasma exosome miR-196a and miR-1246 are indicative of localized pancreatic cancer
To determine whether plasma exosome miRNAs are indicative of localized pancreatic cancer, 15 plasma samples from patents with localized pancreatic cancer (stage I or stage IIA; see Table 2 for patient information) and matched plasma samples from healthy controls were collected. Plasma exosomes were isolated and total RNA was extracted from the isolated plasma exosomes as we recently reported [44]. The expression levels of plasma exosome miR-1246, miR-196a and miR-196b were analyzed by qRT-PCR. As shown in Figure 4, the expression of plasma exosome miR-196a (p=0.0105) and miR-1246 (p = 0.0217) was significantly higher in the group of patients with localized pancreatic cancer as compared to the control population. No such difference was detected for miR-196b. There was no significant difference in plasma exosome miR-1246 and miR-196a expression levels between stage I and stage IIA patients and no correlation between tumor size and the plasma exosome miRNA expression levels (data not shown). Receiver operating characteristic curves (ROC) were constructed for miR-196a, miR-196b and miR-1246 to compare the diagnostic value of their plasma exosome levels in predicting stage I or stage IIA pancreatic cancer (Figure 5). The area under curve (AUC) for miR-196a was 0.81 (95 % CI 0.64, 0.97; p<0.001) and for miR-1246 was 0.73 (95 % CI 0.54, 0.92; p = 0.019), indicating their fair predictive power. The AUC for miR-196b was 0.71 (95 % CI 0.52, 0.91; p=0.033). These data show that plasma exosome miRNAs, especially miR-196a and miR-1246, are indicative of localized pancreatic cancer.
Table 2: Clinicopathological features of pancreatic cancer patients included in the study
Plasma Sample | Age | Gender | Diagnosis | Tumor size | Stagea |
---|---|---|---|---|---|
PL1 | 49 | male | PDAC | 2.5×2.8×3 cm | Stage IIA |
PL2 | 85 | female | PDAC | 1.8×0.7×0.4 cm | Stage IIA |
PL3 | 34 | female | NET | 1.4×1.2×1.8 cm | Stage I |
PL4 | 74 | female | NET | 1.6×1.5×1.8 cm | Stage I |
PL5 | 74 | male | IPMN | 1.4×0.5×1.6cm | Stage I |
PL6 | 70 | male | PDAC | 0.6×0.4×1 cm | Stage IIA |
PL7 | 65 | female | NET | 5.8×5.1×7.2 cm | Stage IIA |
PL8 | 47 | male | NET | 2.5×1.7×2.7 cm | Stage IIA |
PL9 | 84 | male | PDAC | 2.5×4.5 cm | Stage IIA |
PL10 | 64 | female | PDAC | 0.7×0.5×0.8 cm | Stage I |
PL11 | 60 | female | PDAC | NA | Stage I |
PL12 | 73 | male | IPMN | 2×1.5×1.8 cm | Stage I |
PL13 | 67 | male | IPMN | 5×5×1 cm | Stage I |
PL14 | 70 | male | PDAC | 1.5×1.0×2.0 cm | Stage IIA |
PL15 | 83 | female | IPMN | 9×2.5×3.5 cm | Stage IIA |
Mean Age | 66.66 |
aPathologic TNM tumor staging, (American Joint Committee on Cancer Care 7th Edition).
Figure 4: qRT-PCR analysis of plasma exosome miRNA expression in normal and pancreatic cancer patients. Plasma samples were collected from healthy subjects without a history of cancer (n=15) and patients with localized pancreatic cancer (n=15; see Table 2). The Exoquick reagent was used to isolate exosomes from plasma samples. Total RNA was extracted from the plasma exosomes and miR-196a, miR196b and miR-1246 expression was measured by qRT-PCR. Cel-miR-54 was spiked-in to serve as experimental control. Data are expressed as relative expression levels. *p<0.05, Student’s t-test.
Figure 5: ROC analysis of microRNA expression in normal and pancreatic cancer patient plasma exosomes. Receiver Operator Characteristic (ROC) curves were constructed using each microRNA expression value for miR-1246 (top left), miR-196a (top right) and miR-196b (low) respectively. The average of triplicate expression values was computed for individual patient sample and microRNA. The Area Under the Curve (AUC) with 95% CI were computed and shown for each ROC curve. The Wilcoxon-Mann-Whitney test was used to test the null hypothesis that the AUC is equal to .5 (i.e., no predictive power) and the P values for each test were shown.
Plasma exosome miR-196a and miR-1246 expression in subtypes of localized pancreatic cancer
Different histologic types of pancreatic cancer have different clinical outcomes [45]. The 15 patients with localized pancreatic cancer included 7 pancreatic ductal adenocarcinomas (PDAC), 4 intraductal papillary mucinous neoplasms (IPMN), and 4 well differentiated neuroendocrine tumors (NET). Figure 6 shows that plasma exosome miR-196a is a better indicator for PDAC (Figure 6A, P=0.0053) and that plasma exosome miR-1246 levels are significantly elevated in patients with IPMN (Figure 6B, P<0.0001). Plasma exosome miR-196a levels were also significantly higher in IPMN patients (P=0.0246). There was no significant difference in plasma exosome miR-196a and miR-1246 levels in patients with NET as compared to healthy controls (Figure 6C). Similar levels of plasma exosome miR-196b were observed in all comparisons.
Figure 6: qRT-PCR analysis of plasma exosome miRNA expression in PADC, IPMN and NET patients relative to healthy subjects. Expression of miR-1246, miR-196a and miR-196b was analyzed using plasma exosomes from PDAC patients (A, n=7), IPMN patients (B, n=4), NET patients (C, n=4), and healthy subjects (n=15). Cel-miR-54 was spiked-in to serve as experimental control. Data are expressed as relative expression levels. *p<0.05; **p<0.01; ***p<0.001; Student’s t-test.
Immunoaffinity isolation of plasma exosomes using a GPC1 antibody does not improve the diagnostic value of exosome miR-1246 for pancreatic adenocarcinoma
A recent report demonstrated that GPC1 positive plasma exosomes are indicative of pancreatic adenocarcinoma in a large patient cohort study [42]. To explore whether miRNAs in the GPC1 positive exosomes are better indicators of localized pancreatic adenocarcinoma, we utilized GPC1 antibody-conjugated beads to isolate GPC1 positive exosomes from the total plasma exosome populations. The miR-1246 expression in GPC1 positive exosomes from 3 patients with localized pancreatic adenocarcinoma and 3 healthy donors were compared (Figure 7). miR-1246 was analyzed due to its high abundance in pancreatic cancer exosomes (Table 1) and would therefore be detectable in the limited amount of RNA isolated from the immunoaffinity isolated exosomes. There were more GPC1 positive exosomes in the patient plasma than in the control population as suggested by western blot for GPC1 after immunoaffinity isolation (Figure 7A). This was further indicated by the quantity of GPC1 positive exosome protein (Figure 7B left) in the patient plasma versus normal plasma. The GPC1 positive exosome RNA concentration was also slightly higher in patient plasma than in normal plasma even though statistical significance was not evident (Figure 7B center). However, plasma exosome miR-1246 expression levels in patients with localized pancreatic adenocarcinoma were not increased after GPC1 immunoaffinity isolation (Figure 7B right). The same is true when cell line-derived exosomes were analyzed using GPC1 immunoaffinity isolation (data not shown), suggesting that GPC1 positive exosomes do not provide an advantage for detection of localized pancreatic cancer in plasma exosome miRNA. Further verification of this conclusion using a larger set of plasma exosome samples from localized pancreatic cancer patients is warranted.
Figure 7: miR-1246 expression in GPC1 positive plasma exosomes. Total plasma exosomes from three patients with localized pancreatic cancer and three healthy controls were subjected to immunoaffinity isolation using the GPC1 antibody complexed magnetic beads. Exosome protein and RNA were extracted from the antibody-beads complexes using TRIzol reagent. (A) Western blot analysis of the immunoaffinity isolated exosomes from patients (PDAC) and healthy controls (HC), using an antibody against GPC1 (equal amount of exosome protein was loaded in each lane). Note that GPC1 was detected at significantly higher levels in PDAC GPC1 positive exosomes (bands range from 70-130 kDa). (B) Protein concentrations in GPC1 positive plasma exosomes (left), **p<0.01, Student’s t-test; RNA concentrations in GPC1 positive plasma exosomes (center), and miR-1246 expression in GPC1 positive plasma exosomes (right, n=3, means ± SEM).
DISCUSSION
The most interesting finding from the present study is that miR-196a and miR-1246 are highly enriched in pancreatic cancer exosomes, and that these two miRNAs are detected at significantly higher levels in plasma exosome samples from patients with localized pancreatic cancer as compared to those from matched healthy controls. These results indicate that plasma exosome miRNAs, such as miR-196a and miR-1246, are potential biomarkers for early detection of pancreatic cancer. Subtype analysis further indicated that plasma exosome miR-196a is a better indicator for localized PDAC, and plasma exosome miR-1246 is an excellent indicator for localized IPMN. Interestingly, both plasma exosome miR-196a and miR-1246 were expressed at similar levels among patients with NET and the healthy controls. Correlation of histological subtypes and miRNA expression has been described in various malignancies, including lung cancer [46], gastric cancer [47], malignant mesothelioma [48], melanoma [49], cholangiocarcinoma [50], and periampullary carcinoma [51].
miR-1246 contained in microvesicles of colon cancer cells was recently found to promote tumor angiogenesis via targeting the promyelocytic leukemia protein and the Smad signaling pathway in endothelial cells [52]. In hepatocellular carcinoma, miR-1246 was reported to enhance cell migration and invasion through reducing CADM1 expression [53], and increase cancer stemness via targeting Octamer 4 [54]. The expression levels of miR-1246 were reported to be associated with CCNG2-mediated chemoresistance and stemness in pancreatic cancer cells [55]. miR-1246 was also found to promote tumor progression in cervical [56, 57] and lung cancer [58–61]. A recent study indicated that circulating exosome miR-1246 levels are elevated in pancreatic cancer patients [62]. Similar to miR-1246, miR-196a seems to be involved in tumorigenesis and tumor progression in different types of cancer [63]. Recent reports suggested that miR-196a promotes pancreatic cancer progression by targeting the NF-κB inhibitor α [64] and high expression of miR-196a may predict poor outcomes of patients with pancreatic cancer [65]. In addition, studies have shown that circulating miR-196a, in combination with other miRNAs, may serve as biomarkers for pancreatic cancer [65, 66]. These previous reports are consistent with our current findings, suggesting that both miR-1246 and miR-196a are likely oncomiRs and may serve as diagnostic or prognostic markers for pancreatic cancer.
An important feature of this study is that we focused on plasma exosome miRNAs in patients with localized pancreatic cancer, which is at the early stage when local lymph nodes are not involved and there are no symptoms or signs manifested [2]. These patients, if treated, have a significantly higher 5 year survival rate [2]. Serum or plasma miRNAs have been explored for their potential as biomarkers for pancreatic cancer [18–27]. Recent studies have shown that exosome miRNAs are potential indicators for pancreatic cancer [41, 67]. Serum miR-1246 levels were also demonstrated to be significantly higher in pancreatic cancer patients [62]. However, most of these studies have included a general population of pancreatic cancer patients of mixed stages. To our knowledge, there have been no reports focusing exclusively on the potential of plasma or serum exosome miRNAs as biomarkers for localized pancreatic cancer. We have demonstrated in this study, in a pilot cohort of 15 patients with localized pancreatic cancer, that plasma exosome miR-196a and miR-1246 levels are significantly higher in early stage pancreatic cancer as compared to matched healthy controls. When this group of patients was divided into histologic subtypes, miR-196a was found to be a better and significant indicator of PDAC, and miR-1246 an excellent indicator of IPMN. These observations are generally in line with previous reports showing that serum levels of miR-1246 is elevated in pancreatic cancer patients [62], and that miR-196a is deregulated in pancreatic cancer tissues, including those at early stages [66]. The findings from this pilot study signify new plasma exosome miRNA signatures that are potential biomarkers for early detection of pancreatic cancer. Given the relative small patient sample size of this pilot study, as limited by the availability of localized pancreatic cancer specimens, further evaluation of these two plasma exosome miRNAs in a larger cohort is required to establish these new biomarkers for detection of localized pancreatic cancer.
Another feature of the current study, which is different from many other biomarker studies, is that we began the investigation with small RNA sequencing to characterize pancreatic cancer exosome miRNA signatures using human pancreatic cell line-derived exosomes. We have recently applied and validated this strategy for biomarker discovery in breast cancer [44]. The RNA sequencing results provided us with an indication as to which exosome miRNAs we might focus on for biomarker development. Two criteria were used to select candidate exosome miRNAs for clinical study: first, these miRNAs are selectively encapsulated in pancreatic cancer exosomes versus normal pancreatic ductal epithelial cell exosomes; and second, these miRNAs are highly abundant in pancreatic cancer exosomes. These criteria are based on the consideration that at the localized pancreatic cancer stages (stage I-IIA), there would likely have fewer cancer exosomes available in the circulation as the tumor is still localized to the pancreas without lymph node involvement and distal metastasis. We selected miR-1246, miR-196a and miR-196b for further evaluation using pancreatic cancer patient plasmas. Interestingly, while plasma exosome miR-196b expression was similar in pancreatic cancer patients and the healthy controls, which is inconsistent with cell line studies, plasma exosome miR-196a and miR-1246 expression was significantly higher in pancreatic cancer patients than in the matched healthy controls. These results indicate that starting with a characterization of human pancreatic cancer cell line-derived exosomes seems to be a valid approach for plasma exosome biomarker discovery; but one needs to be cautious as the results from cell lines may not always be applicable to patient plasma exosomes.
The ROC curve analysis for our patient study indicates that while the difference of plasma exosome miR-196a and miR-1246 expression between pancreatic cancer patients and the controls is statistically significant, the sensitivity and specificity of these plasma exosome miRNAs as indictors for localized pancreatic cancer may be further improved (see AUC for each ROC curve in Figure 5). In this context, we attempted to evaluate pancreatic cancer exosomes from whole plasma exosomes using an immunoaffinity isolation method, in order to obtain more specific and sensitive analysis of plasma exosome miRNA expression. As GPC1 was recently reported as a pancreatic adenocarcinoma exosome surface marker [42], we performed immunoaffinity isolation from whole plasma exosomes using the reported GPC1 antibody. Whereas immunoaffinity isolation with the GPC1 antibody captured more exosomes from pancreatic adenocarcinoma patient plasma samples, miR-1246 expression in GPC1 positive exosomes is not elevated in patients with localized PDAC compared to the matched controls. This unexpected result suggests that immunoaffinity isolation using GPC1 antibodies for plasma exosome miRNA analysis may not be an ideal approach to capture differentially expressed exosome miRNAs in patients with localized PDAC. New strategies are needed to furthering this avenue of investigation.
In conclusion, we have demonstrated that certain miRNAs, such as miR-196a and miR-1246, are highly and selectively enriched in pancreatic cancer cell-derived exosomes and that the plasma exosome levels of these two miRNAs are significantly elevated in patients with localized pancreatic cancer. Plasma exosome miR-196a and miR-1246 may serve as biomarkers for different subtypes of localized pancreatic cancer, with miR-196a for PDAC and miR-1246 for IPMN. As early diagnosis remains the key to improving pancreatic cancer outcomes, further studies are warranted to firmly establish plasma exosome miRNA signatures indicative of localized pancreatic cancer.
MATERIALS AND METHODS
Cell culture
The human pancreatic epithelial cell line hTERT-HPNE, and pancreatic cancer cell lines PANC-1, MIA-PaCa-2 and BxPC-3 were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). Cells were cultured according to ATCC’s instructions with the exception of exosome-depleted fetal bovine serum (FBS) and horse serum being used wherever appropriate. Exosome-depleted FBS and horse serum were prepared by pelleting the serum exosomes by ultracentrifugation at 100,000 × g for 2 h at 4°C, and the resulting supernatant was filtered through a 0.2-μm pore filter. Cells were routinely maintained in a humidified chamber at 37°C and 5% CO2.
Patient plasma samples
Plasma samples (n=15) were collected from donors with no history of pancreatic cancer (mean age = 48, 11 females and 4 males) from the Oklahoma Blood Institute, Oklahoma City, OK, USA under an approved IRB protocol. Localized pancreatic cancer plasma samples (n=15, Table 2) were collected from patients who were seen at the Stephenson Cancer Center at the University of Oklahoma Health Science Center (OUHSC), and underwent primary tumor biopsy or resection. The study was approved by the OUHSC IRB (#5535) and written informed consent was obtained from all participants. Blood was collected before biopsy or resection as described [44]. Plasma was separated by centrifugation and stored at -80 °C. Clinical stages were classified using the tumor, node, metastasis (TNM) system (AJCC 7th edition). Patient information including age, tumor size, clinical stages and histological grades were obtained from the clinical and pathological records.
Exosome isolation
Exosomes were isolated from the culture medium utilizing a combination of centrifugation, ultracentrifugation, and filtration as we previously described [43, 44]. Plasma exosomes were isolated using the Exoquick reagent (System Biosciences, Mountain View, CA, USA) following the manufacturer’s protocol (500μl plasma per sample).
Western blot analysis
Total exosome protein was prepared by resuspending the exosomes in RIPA Buffer (50mM Tris-HCl pH 7.4, 150mM NaCl, 0.5% sodium deoxycholate, 1% NP-40, and 0.1% sodium dodecyl sulfate) containing 1 mM phenlymethylsulfonyl fluoride, 5μg/ml leupeptin, 2μg/ml aprotinin, and 1μg/ml pepstatin A [44]. Approximately 30–40μg of protein from each sample was separated under non-reducing conditions (for CD63) or reducing conditions (for glypican-1 (GPC1)) on a 10% SDS-PAGE gel, transferred to a polyvinylidene fluoride membrane, blotted with a primary antibody against CD63 (sc-5275, Santa Cruz Biotechnology, Santa Cruz, CA, USA), GPC1 (PA-24972, Thermo Scientific, Waltham, MA USA), and a secondary IgG-HRP (Santa Cruz Bio Technology Inc., CA, USA).
Electron microscopy
Exosomes were suspended in 1 × PBS and fixed with 2% paraformaldehyde. The fixed sample was absorbed onto formvar-coated copper grids for 20 min in a dry environment and fixed in 1% glutaraldehyde for 5 min. After being rinsed in distilled water, samples were stained with uranyl oxalate for 5 min followed by methyl cellulose uranyl acetate for 10 min on ice. Excess liquid was wicked off the grid, and grids were stored at room temperature until imaging. Imaging was performed using a Hitachi H7600 electron microscope.
Nanoparticle tracking analysis
Isolated exosomes were diluted in PBS and analyzed using the Nanosight NS300 System (Malvern Instruments, UK) which is equipped with a blue laser (405 nm) [44]. Nanoparticles illuminated by the laser were captured for 60 seconds. The Nanosight Tracking Analysis software was used to provide particle concentrations and size distribution profiles.
RNA extraction
Total RNA was extracted from exosome pellets using the TRIzol reagent (Invitrogen/Life Technologies, Carlsbad, California) following the manufacturer’s protocol [44]. RNA concentration was quantitated using the NanoDrop ND-100 Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).
Small RNA library preparation and next generation sequencing
Small RNA libraries were constructed using the New England Biolabs (NEB) NEBNext Multiplex Small RNA Library Prep Set for Illumina sequencers and the NEB standard protocol, as we recently described [44]. Each library was indexed to multiplex four samples per sequencing run on the Illumina MiSeq platform using MiSeq 50 cycle Reagent Kits v2. A minimum of 17 million 50-bp sequencing reads were collected from each sample and data were analyzed using the Genesifter software (formerly Geospiza, PerkinElmer, Santa Clara, CA, USA). Raw data for each sample were aligned to the most recent miRBASE database (miRBase.org [68]). Pairwise comparisons of the alignment results for identification of miRNAs that are differentially expressed at a significant level were conducted using the Genesifter software.
Quantitative real-time reverse transcription PCR
For miRNA expression analysis, complementary DNA from 80 ng of total RNA was synthesized using the qScript miRNA cDNA Synthesis Kit (Quanta BioSciences, Inc., Beverly, MA, USA). A cDNA aliquot, equivalent of 4 ng of the original RNA, was mixed with Perfecta SYBR Green SuperMix, Universal PCR Primer and specific forward primers for miR-1246 (5’-GCGCGATGGATTTTTGGAGCAG-3’), miR-196a (5’-GCGTAGGTAGTTTCCTGTTGTTGGG-3’), or miR-196b (5’-GCGTAGGTAGTTTCCTGTTGTTGGG-3’) in 20 μl qPCR reactions. A synthetic Caenorhabditis elegans miR-54 (cel-miR-54) RNA oligonucleotide (Integrated DNA Technologies, Coralville, IA, USA) was spiked-into RNA samples as a control [69] (5’-UACCCGUAAUCUUCAUAAUCCGAG-3’). PCR reactions were run on a Bio-Rad CFX 96 Real-Time PCR (Bio-Rad, Hercules, CA, USA) instrument under following conditions: 95 °C for 2 min, 40 cycles of 95 °C for 5 s, 60 °C for 15 s and 70°C for 15 s. Changes in gene expression were calculated using the ΔΔCT method as follows: ΔCT = CT(target miRNA) − CT (cel-miR-54); ΔΔCT = ΔCT (cancer exosome miRNA)- ΔCT (normal exosome miRNA); the fold change = 2−ΔΔCT.
Immunoaffinity magnetic bead-based exosome isolation
Protein G Dynabeads were purchased from Thermo Fisher (Waltham, MA, USA). The process of preparing GPC1 antibody (PA-24972, Thermo Scientific, Waltham, MA USA) coated beads and immunoaffinity capture of the exosomes were described in our recent report [44]. Proteins and RNA was extracted from the exosome-bound beads using the TRIzol reagent.
Statistical analysis
Statistical analyses were completed using GraphPad Prism software (GraphPad Software, Inc. La Jolla, CA, USA) and SAS software (SAS institute, Cary, North Carolina, USA). Student’s t test was used to determine significant differences among control and experimental groups. Receiver operating characteristic (ROC) curves were constructed using each miRNA expression value. The area under the curve (AUC) with 95% CI was calculated for each ROC curve. The Wilcoxon-Mann-Whitney test was used to test the null hypothesis that the AUC is equal to 0.5 (i.e., no predictive power).
Abbreviations
Cancer antigen 19-9 (CA 19-9); Area under the curve (AUC); fetal bovine serum (FBS); glypican-1 (GPC1); microRNA (miRNA); intraductal papillary mucinous neoplasms (IPMN); neuroendocrine tumors (NET); New England Biolabs (NEB); pancreatic ductal adenocarcinomas (PDAC), polyvinylidene fluoride (PVDF); Quantitative Real-Time PCR (qRT-PCR); Receiver operating characteristic (ROC).
Author contributions
YFX carried out the experiments and drafted the manuscript. BNH participated in the design and coordination of the experiments, and helped to draft the manuscript. YDZ conducted the statistical analysis on patient plasma data. RGP participated in the design and coordination of the study. WQD conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
ACKNOWLEDGMENTS
We greatly appreciate the services and support provided by the Functional Genomics Core facility at the Peggy and Charles Stephenson Cancer Center, the DNA Sequencing and Genomics facility in the Laboratory for Molecular Biology and Cytometry Research at the University of Oklahoma Health Sciences Center and the Imaging Core Facility at the Oklahoma Medical Research Foundation.
CONFLICTS OF INTEREST
There are no conflicts to disclose.
FUNDING
This study was supported in part by grants from the National Institute of General Medical Sciences of the National Institutes of Health (U54GM104938), the Oklahoma Center for the Advancement of Science and Technology (HR14-147), the Presbyterian Health Foundation, and the Peggy and Charles Stephenson Cancer Center.
REFERENCES
1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016; 66:7-30.
2. NATIONAL CANCER INSTITUTE-SEER Stat Fact Sheets: Pancreas Cancer. http://seer.cancer.gov/statfacts/html/pancreas.html. Accessed on 8 March, 2017.
3. NATIONAL CANCER INSTITUTE--SEER Stat Fact Sheets: All Cancer Sites. http://seer.cancer.gov/statfacts/html/all.html. Accessed on 8 March, 2017.
4. Hanada K, Okazaki A, Hirano N, Izumi Y, Minami T, Ikemoto J, Kanemitsu K, Hino F. Effective screening for early diagnosis of pancreatic cancer. Best practice & research Clinical gastroenterology. 2015; 29:929-939.
5. Donahue TR, Dawson DW. Leveraging Mechanisms Governing Pancreatic Tumorigenesis To Reduce Pancreatic Cancer Mortality. Trends in endocrinology and metabolism: TEM. 2016; 27:770-781.
6. Kamarajah SK, Burns WR, Frankel TL, Cho CS, Nathan H. Validation of the American Joint Commission on Cancer (AJCC) 8th Edition Staging System for Patients with Pancreatic Adenocarcinoma: A Surveillance, Epidemiology and End Results (SEER) Analysis. Annals of surgical oncology. 2017.
7. Porta M, Fabregat X, Malats N, Guarner L, Carrato A, de Miguel A, Ruiz L, Jariod M, Costafreda S, Coll S, Alguacil J, Corominas JM, Sola R, et al. Exocrine pancreatic cancer: symptoms at presentation and their relation to tumour site and stage. Clinical & translational oncology. 2005; 7:189-197.
8. Yabar CS, Winter JM. Pancreatic Cancer: A Review. Gastroenterology clinics of North America. 2016; 45:429-445.
9. Ballehaninna UK, Chamberlain RS. The clinical utility of serum CA 19-9 in the diagnosis, prognosis and management of pancreatic adenocarcinoma: An evidence based appraisal. Journal of gastrointestinal oncology. 2012; 3:105-119.
10. Medscape-Pancreatic Cancer. http://emedicine.medscape.com/article/280605-overview. Accessed on 1 March, 2017.
11. Winter JM, Yeo CJ, Brody JR. Diagnostic, prognostic, and predictive biomarkers in pancreatic cancer. Journal of surgical oncology. 2013; 107:15-22.
12. Frebourg T, Bercoff E, Manchon N, Senant J, Basuyau JP, Breton P, Janvresse A, Brunelle P, Bourreille J. The evaluation of CA 19-9 antigen level in the early detection of pancreatic cancer. A prospective study of 866 patients. Cancer. 1988; 62:2287-2290.
13. Frampton AE, Krell J, Jamieson NB, Gall TM, Giovannetti E, Funel N, Mato Prado M, Krell D, Habib NA, Castellano L, Jiao LR, Stebbing J. microRNAs with prognostic significance in pancreatic ductal adenocarcinoma: A meta-analysis. European journal of cancer. 2015; 51:1389-1404.
14. Hernandez YG, Lucas AL. MicroRNA in pancreatic ductal adenocarcinoma and its precursor lesions. World journal of gastrointestinal oncology. 2016; 8:18-29.
15. Sun L, Chua CY, Tian W, Zhang Z, Chiao PJ, Zhang W. MicroRNA Signaling Pathway Network in Pancreatic Ductal Adenocarcinoma. Journal of genetics and genomics. 2015; 42:563-577.
16. Yao J, Li Z, Wang X, Xu P, Zhao L, Qian J. MiR-125a regulates chemo-sensitivity to gemcitabine in human pancreatic cancer cells through targeting A20. Acta biochimica et biophysica Sinica. 2016; 48:202-208.
17. Li H, Xiang H, Ge W, Wang H, Wang T, Xiong M. Expression and functional perspectives of miR-184 in pancreatic ductal adenocarcinoma. International journal of clinical and experimental pathology. 2015; 8:12313-12318.
18. Le Large TY, Meijer LL, Mato Prado M, Kazemier G, Frampton AE, Giovannetti E. Circulating microRNAs as diagnostic biomarkers for pancreatic cancer. Expert review of molecular diagnostics. 2015; 15:1525-1529.
19. Permuth-Wey J, Chen DT, Fulp WJ, Yoder SJ, Zhang Y, Georgeades C, Husain K, Centeno BA, Magliocco AM, Coppola D, Malafa M. Plasma MicroRNAs as Novel Biomarkers for Patients with Intraductal Papillary Mucinous Neoplasms of the Pancreas. Cancer prevention research. 2015; 8:826-834.
20. Cote GA, Gore AJ, McElyea SD, Heathers LE, Xu H, Sherman S, Korc M. A pilot study to develop a diagnostic test for pancreatic ductal adenocarcinoma based on differential expression of select miRNA in plasma and bile. The American journal of gastroenterology. 2014; 109:1942-1952.
21. Liu R, Chen X, Du Y, Yao W, Shen L, Wang C, Hu Z, Zhuang R, Ning G, Zhang C, Yuan Y, Li Z, Zen K, et al. Serum microRNA expression profile as a biomarker in the diagnosis and prognosis of pancreatic cancer. Clinical chemistry. 2012; 58:610-618.
22. Komatsu S, Ichikawa D, Takeshita H, Morimura R, Hirajima S, Tsujiura M, Kawaguchi T, Miyamae M, Nagata H, Konishi H, Shiozaki A, Otsuji E. Circulating miR-18a: a sensitive cancer screening biomarker in human cancer. In Vivo. 2014; 28:293-297.
23. Morimura R, Komatsu S, Ichikawa D, Takeshita H, Tsujiura M, Nagata H, Konishi H, Shiozaki A, Ikoma H, Okamoto K, Ochiai T, Taniguchi H, Otsuji E. Novel diagnostic value of circulating miR-18a in plasma of patients with pancreatic cancer. British journal of cancer. 2011; 105:1733-1740.
24. Liffers ST, Munding JB, Vogt M, Kuhlmann JD, Verdoodt B, Nambiar S, Maghnouj A, Mirmohammadsadegh A, Hahn SA, Tannapfel A. MicroRNA-148a is down-regulated in human pancreatic ductal adenocarcinomas and regulates cell survival by targeting CDC25B. Laboratory investigation. 2011; 91:1472-1479.
25. Munding JB, Adai AT, Maghnouj A, Urbanik A, Zollner H, Liffers ST, Chromik AM, Uhl W, Szafranska-Schwarzbach AE, Tannapfel A, Hahn SA. Global microRNA expression profiling of microdissected tissues identifies miR-135b as a novel biomarker for pancreatic ductal adenocarcinoma. International journal of cancer. 2012; 131:E86-95.
26. Zhang J, Zhao CY, Zhang SH, Yu DH, Chen Y, Liu QH, Shi M, Ni CR, Zhu MH. Upregulation of miR-194 contributes to tumor growth and progression in pancreatic ductal adenocarcinoma. Oncology reports. 2014; 31:1157-1164.
27. Zhao C, Zhang J, Zhang S, Yu D, Chen Y, Liu Q, Shi M, Ni C, Zhu M. Diagnostic and biological significance of microRNA-192 in pancreatic ductal adenocarcinoma. Oncology reports. 2013; 30:276-284.
28. Kowal J, Tkach M, Thery C. Biogenesis and secretion of exosomes. Current opinion in cell biology. 2014; 29:116-125.
29. S ELA, Mager I, Breakefield XO, Wood MJ. Extracellular vesicles: biology and emerging therapeutic opportunities. Nature reviews Drug discovery. 2013; 12:347-357.
30. Keller S, Sanderson MP, Stoeck A, Altevogt P. Exosomes: from biogenesis and secretion to biological function. Immunology letters. 2006; 107:102-108.
31. Goetzl EJ, Goetzl L, Karliner JS, Tang N, Pulliam L. Human plasma platelet-derived exosomes: effects of aspirin. FASEB journal. 2016.
32. Hannafon BN, Ding WQ. Intercellular Communication by Exosome-Derived microRNAs in Cancer. International journal of molecular sciences. 2013; 14:14240-14269.
33. Kosaka N. Decoding the Secret of Cancer by Means of Extracellular Vesicles. Journal of clinical medicine. 2016; 5.
34. Ciardiello C, Cavallini L, Spinelli C, Yang J, Reis-Sobreiro M, de Candia P, Minciacchi VR, Di Vizio D. Focus on Extracellular Vesicles: New Frontiers of Cell-to-Cell Communication in Cancer. International journal of molecular sciences. 2016; 17.
35. Falasca M, Kim M, Casari I. Pancreatic cancer: Current research and future directions. Biochimica et biophysica acta. 2016; 1865:123-132.
36. Shao H, Chung J, Issadore D. Diagnostic technologies for circulating tumour cells and exosomes. Bioscience reports. 2015; 36.
37. De Toro J, Herschlik L, Waldner C, Mongini C. Emerging roles of exosomes in normal and pathological conditions: new insights for diagnosis and therapeutic applications. Frontiers in immunology. 2015; 6:203.
38. Robbins PD, Morelli AE. Regulation of immune responses by extracellular vesicles. Nature reviews Immunology. 2014; 14:195-208.
39. Sadovska L, Eglitis J, Line A. Extracellular Vesicles as Biomarkers and Therapeutic Targets in Breast Cancer. Anticancer research. 2015; 35:6379-6390.
40. Thind A, Wilson C. Exosomal miRNAs as cancer biomarkers and therapeutic targets. Journal of extracellular vesicles. 2016; 5:31292.
41. Lai X, Wang M, McElyea SD, Sherman S, House M, Korc M. A microRNA signature in circulating exosomes is superior to exosomal glypican-1 levels for diagnosing pancreatic cancer. Cancer Lett. 2017.
42. Melo SA, Luecke LB, Kahlert C, Fernandez AF, Gammon ST, Kaye J, LeBleu VS, Mittendorf EA, Weitz J, Rahbari N, Reissfelder C, Pilarsky C, Fraga MF, et al. Glypican-1 identifies cancer exosomes and detects early pancreatic cancer. Nature. 2015; 523:177-182.
43. Hannafon BN, Carpenter KJ, Berry WL, Janknecht R, Dooley WC, Ding WQ. Exosome-mediated microRNA signaling from breast cancer cells is altered by the anti-angiogenesis agent docosahexaenoic acid (DHA). Mol Cancer. 2015; 14:133.
44. Hannafon BN, Trigoso YD, Calloway CL, Zhao YD, Lum DH, Welm AL, Zhao ZJ, Blick KE, Dooley WC, Ding WQ. Plasma exosome microRNAs are indicative of breast cancer. Breast Cancer Res. 2016; 18:90.
45. Hosoda W, Wood LD. Molecular Genetics of Pancreatic Neoplasms. Surg Pathol Clin. 2016; 9:685-703.
46. Li H, Jiang Z, Leng Q, Bai F, Wang J, Ding X, Li Y, Zhang X, Fang H, Yfantis HG, Xing L, Jiang F. A prediction model for distinguishing lung squamous cell carcinoma from adenocarcinoma. Oncotarget. 2017; 8:50704-50714. doi: 10.18632/oncotarget.17038.
47. Yepes S, Lopez R, Andrade RE, Rodriguez-Urrego PA, Lopez-Kleine L, Torres MM. Co-expressed miRNAs in gastric adenocarcinoma. Genomics. 2016; 108:93-101.
48. Lamberti M, Capasso R, Lombardi A, Di Domenico M, Fiorelli A, Feola A, Perna AF, Santini M, Caraglia M, Ingrosso D. Two Different Serum MiRNA Signatures Correlate with the Clinical Outcome and Histological Subtype in Pleural Malignant Mesothelioma Patients. PloS one. 2015; 10:e0135331.
49. Poliseno L, Haimovic A, Segura MF, Hanniford D, Christos PJ, Darvishian F, Wang J, Shapiro RL, Pavlick AC, Berman RS, Hernando E, Zavadil J, Osman I. Histology-specific microRNA alterations in melanoma. The Journal of investigative dermatology. 2012; 132:1860-1868.
50. Plieskatt JL, Rinaldi G, Feng Y, Peng J, Yonglitthipagon P, Easley S, Laha T, Pairojkul C, Bhudhisawasdi V, Sripa B, Brindley PJ, Mulvenna JP, Bethony JM. Distinct miRNA signatures associate with subtypes of cholangiocarcinoma from infection with the tumourigenic liver fluke Opisthorchis viverrini. Journal of hepatology. 2014; 61:850-858.
51. Kalluri Sai Shiva UM, Kuruva MM, Mitnala S, Rupjyoti T, Guduru Venkat R, Botlagunta S, Kandagaddala R, Siddapuram SP, Sekaran A, Chemalakonda R, Rebala P, Duvvuru NR. MicroRNA profiling in periampullary carcinoma. Pancreatology. 2014; 14:36-47.
52. Yamada N, Tsujimura N, Kumazaki M, Shinohara H, Taniguchi K, Nakagawa Y, Naoe T, Akao Y. Colorectal cancer cell-derived microvesicles containing microRNA-1246 promote angiogenesis by activating Smad 1/5/8 signaling elicited by PML down-regulation in endothelial cells. Biochim Biophys Acta. 2014; 1839:1256-1272.
53. Sun Z, Meng C, Wang S, Zhou N, Guan M, Bai C, Lu S, Han Q, Zhao RC. MicroRNA-1246 enhances migration and invasion through CADM1 in hepatocellular carcinoma. BMC Cancer. 2014; 14:616.
54. Chai S, Ng KY, Tong M, Lau EY, Lee TK, Chan KW, Yuan YF, Cheung TT, Cheung ST, Wang XQ, Wong N, Lo CM, Man K, et al. Octamer 4/microRNA-1246 signaling axis drives Wnt/beta-catenin activation in liver cancer stem cells. Hepatology. 2016; 64:2062-2076.
55. Hasegawa S, Eguchi H, Nagano H, Konno M, Tomimaru Y, Wada H, Hama N, Kawamoto K, Kobayashi S, Nishida N, Koseki J, Nishimura T, Gotoh N, et al. MicroRNA-1246 expression associated with CCNG2-mediated chemoresistance and stemness in pancreatic cancer. Br J Cancer. 2014; 111:1572-1580.
56. Zhang B, Chen J, Ren Z, Chen Y, Li J, Miao X, Song Y, Zhao T, Li Y, Shi Y, Ren D, Liu J. A specific miRNA signature promotes radioresistance of human cervical cancer cells. Cancer Cell Int. 2013; 13:118.
57. Chen J, Yao D, Zhao S, He C, Ding N, Li L, Long F. MiR-1246 promotes SiHa cervical cancer cell proliferation, invasion, and migration through suppression of its target gene thrombospondin 2. Arch Gynecol Obstet. 2014; 290:725-732.
58. Huang W, Li H, Luo R. The microRNA-1246 promotes metastasis in non-small cell lung cancer by targeting cytoplasmic polyadenylation element-binding protein 4. Diagn Pathol. 2015; 10:127.
59. Kim G, An HJ, Lee MJ, Song JY, Jeong JY, Lee JH, Jeong HC. Hsa-miR-1246 and hsa-miR-1290 are associated with stemness and invasiveness of non-small cell lung cancer. Lung Cancer. 2016; 91:15-22.
60. Yuan D, Xu J, Wang J, Pan Y, Fu J, Bai Y, Zhang J, Shao C. Extracellular miR-1246 promotes lung cancer cell proliferation and enhances radioresistance by directly targeting DR5. Oncotarget. 2016; 7:32707-32722. doi: 10.18632/oncotarget.9017.
61. Zhang WC, Chin TM, Yang H, Nga ME, Lunny DP, Lim EK, Sun LL, Pang YH, Leow YN, Malusay SR, Lim PX, Lee JZ, Tan BJ, et al. Tumour-initiating cell-specific miR-1246 and miR-1290 expression converge to promote non-small cell lung cancer progression. Nature communications. 2016; 7:11702.
62. Madhavan B, Yue S, Galli U, Rana S, Gross W, Muller M, Giese NA, Kalthoff H, Becker T, Buchler MW, Zoller M. Combined evaluation of a panel of protein and miRNA serum-exosome biomarkers for pancreatic cancer diagnosis increases sensitivity and specificity. Int J Cancer. 2015; 136:2616-2627.
63. Chen ZY, Chen X, Wang ZX. The role of microRNA-196a in tumorigenesis, tumor progression, and prognosis. Tumour biology. 2016.
64. Huang F, Tang J, Zhuang X, Zhuang Y, Cheng W, Chen W, Yao H, Zhang S. MiR-196a promotes pancreatic cancer progression by targeting nuclear factor kappa-B-inhibitor alpha. PloS one. 2014; 9:e87897.
65. Bloomston M, Frankel WL, Petrocca F, Volinia S, Alder H, Hagan JP, Liu CG, Bhatt D, Taccioli C, Croce CM. MicroRNA expression patterns to differentiate pancreatic adenocarcinoma from normal pancreas and chronic pancreatitis. Jama. 2007; 297:1901-1908.
66. Xue Y, Abou Tayoun AN, Abo KM, Pipas JM, Gordon SR, Gardner TB, Barth RJ Jr, Suriawinata AA, Tsongalis GJ. MicroRNAs as diagnostic markers for pancreatic ductal adenocarcinoma and its precursor, pancreatic intraepithelial neoplasm. Cancer genetics. 2013; 206:217-221.
67. Nuzhat Z, Kinhal V, Sharma S, Rice GE, Joshi V, Salomon C. Tumour-derived exosomes as a signature of pancreatic cancer - liquid biopsies as indicators of tumour progression. Oncotarget. 2017; 8:17279-17291. doi: 10.18632/oncotarget.13973.
68. Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011; 39:D152-157.
69. Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, Sloan CA, Rosenbloom KR, Roe G, Rhead B, Raney BJ, Pohl A, Malladi VS, et al. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013; 41:D64-69.