Candidate microRNA biomarkers of pancreatic ductal adenocarcinoma: meta-analysis, experimental validation and clinical significance
Journal of Experimental & Clinical Cancer Research volume 32, Article number: 71 (2013)
The diagnostic and prognostic value of microRNA (miRNA) expression aberrations in pancreatic ductal adenocarcinoma (PDAC) has been studied extensively in recent years. However, differences in measurement platforms and lab protocols as well as small sample sizes can render gene expression levels incomparable.
A comprehensive meta-review of published studies in PDAC that compared the miRNA expression profiles of PDAC tissues and paired neighbouring noncancerous pancreatic tissues was performed to determine candidate miRNA biomarkers for PDAC. Both a miRNA vote-counting strategy and a recently published Robust Rank Aggregation method were employed. In this review, a total of 538 tumour and 206 noncancerous control samples were included.
We identified a statistically significant miRNA meta-signature of seven up- and three down-regulated miRNAs. The experimental validation results showed that the miRNA expression levels were in accordance with the meta-signature. The results from the vote-counting strategy were consistent with those from the Robust Rank Aggregation method. The experimental validation confirmed that the statistically unique profiles identified by the meta-review approach could discriminate PDAC tissues from paired nonmalignant pancreatic tissues. In a cohort of 70 patients, the high expression of miR-21 (p=0.018, HR=2.610; 95% CI=1.179-5.777) and miR-31 (p=0.039, HR=2.735; 95% CI=1.317-6.426), the low expression of miR-375 (p=0.022, HR=2.337; 95% CI=1.431-5.066) were associated with poor overall survival following resection, independent of clinical covariates.
The identified miRNAs may be used to develop a panel of diagnostic and prognostic biomarkers for PDAC with sufficient sensitivity and specificity for use in a clinical setting.
Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive of all malignancies . Less than 20% of PDAC patients present with localised, potentially curable tumours. The overall 5-year survival rate is <5%. Because the chemotherapeutic options prolong life only minimally, the current PDAC mortality is nearly identical to its incidence . Extensive studies have been performed to identify biomarkers for this disease. At the messenger RNA (mRNA) level, quite a few, including some very specific molecular variations have been found in cancerous tissues . MicroRNAs (miRNAs), a class of short non-coding RNA molecules that range in size from 19 to 25 nucleotides, have been proposed as promising biomarkers of early cancer detection and accurate prognosis as well as targets for more efficient treatment [4, 5]. MiRNAs play important roles in regulating the translation of many genes and the degradation of mRNAs through base pairing to partially complementary sites, predominately in the 3′ untranslated region [6, 7]. Several studies have implicated miRNAs in the regulation of tumour biology [8–10]. Model biomarkers should be easily quantifiable and associate strongly with clinical outcome, and miRNAs may match these criteria.
High-throughput technologies have been employed to identify differences in miRNA expression levels between normal and cancerous tissues. These studies have the potential to identify dozens or hundreds of differentially expressed miRNAs, although only a small fraction of them may be of actual clinical utility as diagnostic/prognostic biomarkers. Finding a meaningful way in which to combine different data sources is often a non-trivial task. Differences in measurement platforms and lab protocols as well as small sample sizes can render gene expression levels incomparable. Hence, it may be better to analyse datasets separately and then aggregate the resulting gene lists. This strategy has been applied to identify gene co-expression networks  and to define more robust sets of cancer-related genes [12, 13] and miRNAs [14, 15].
In the meta-review approach, the results of several individual studies are combined to increase statistical power and subsequently resolve any inconsistencies or discrepancies among different profiling studies. In this study, we applied two meta-review approaches: the well-known vote-counting strategy [12, 13], which is based on the number of studies reporting a gene as being consistently expressed and then further ranking these genes with respect to total sample size and average fold-change, and the recently published Robust Rank Aggregation method [16, 17]. Pathway analysis was then performed to identify the physiological impact of miRNA deregulation in PDAC progression. Moreover, we further validated the most up-regulated and down-regulated miRNAs from the meta-review in a clinical setting. The expression levels of a subset of candidate miRNAs were assessed by quantitative real-time polymerase chain reaction (qRT-PCR). With the validation of candidate miRNAs, we selected the most promising miRNAs based on factors such as fold-change to explore their potential effects on the survival of PDAC patients after surgical resection.
Materials and methods
Selection of studies and datasets
The Scopus database (http://www.scopus.com) was searched for PADC miRNA expression profiling studies using search term TITLE-ABS-KEY [(mirna* OR microrna* OR mir-* OR mir) AND profil* AND (pancreas* cancer OR pancreatic carcinoma OR pancreas* neoplas* OR pancreatic tumo* OR carcinoma of pancreas* OR cancer of pancreas*)]. The same strategy was also applied to searches of the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/), ArrayExpress (http://www.ebi.ac.uk/arrayexpress/), and PubMed (http://www.ncbi.nlm.nih.gov/pubmed). The last search was performed on May 11, 2013. The titles and abstracts of the articles were screened, and the full text of the articles of interest was evaluated. We included only original experimental articles that were published in English and that compared the expression of miRNAs in PDAC tissue and noncancerous pancreatic tissue in humans. Articles were excluded based on the following criteria: (i) review articles, case reports or letters; (ii) non-English articles; (iii) studies of individual pre-selected candidate miRNAs; (iv) studies that used RT-PCR for initial selection (the reasons for this exclusion criterion are explained in the Discussion section); (v) studies using cell lines or serum from PDAC patients; (vi) studies that did not use a miRNA microarray platform; (vii) studies profiling different histological subtypes; (viii) studies that did not include noncancerous tissue.
Two investigators (MM and XK) independently evaluated and extracted the data using standard protocols, and all discrepancies were resolved by a third investigator (MW). From the full text and corresponding supplemental information, the following eligibility items were collected and recorded for each study: author, region, period, selection and characteristics of the recruited PDAC patients, platform of miRNA expression profiling, and the list of up- and down-regulated miRNAs and their corresponding fold-changes. When the gene list was not available, the authors were contacted directly. All miRNA names were standardised according to miRBase version 20.
The miRNAs were ranked according to their importance as follows: (i) number of comparisons in agreement (i.e., listing the same miRNAs as having a consistent direction of change and being differentially expressed, respectively); (ii) total number of samples for comparison in agreement; (iii) average fold-changes reported for comparisons in agreement. Total sample size was considered more important than average fold-change because many studies did not report a fold-change. Furthermore, the average fold-change was based solely on the subset of studies for which a fold change value was available.
Robust rank aggregation method
The list of extracted miRNAs was ranked based on their associated p-values (less than 0.05 was considered significant) when their fold-changes were not reported. All of the protocols for the Robust Rank Aggregation method are free to download at the comprehensive R Archive Network website (http://cran.r-project.org/). Details can be found in the package documentation (http://cran.r-project.org/web/packages/RobustRankAggreg/RobustRankAggreg.pdf). This method assigns a p- value to each element in the aggregated list, which indicates how much better it is ranked compared with a null model, expecting random ordering. To assess the stability of the acquired p- values, leave-one-out cross-validation was applied in the Robust Rank Aggregation algorithm. This analysis was repeated 10,000 times, and each time, one random gene list was left out of the analysis. The p- values acquired from each round for each miRNA were then averaged.
MiRNA target prediction and enrichment analysis
The mRNA targets of the miRNA genes were predicted using TargetScan (http://www.targetscan.org/), miRDB (http://mirdb.org/miRDB/), and miRANDA (http://www.microrna.org/microrna/getGeneForm.do), as each algorithm determines target binding differently. We selected targets from the miRANDA/miSVR search with scores less than −1.25 for further analysis. Enrichment analyses for KEGG and Panther pathways and Gene Ontology terms were performed with the GeneCodis tool (http://genecodis.dacya.ucm.es/). The potential targets of each miRNA were used as input.
Ethical approval for this study was obtained from the Department of General Surgery of Ruijin Hospital at Shanghai Jiaotong University (Shanghai, China). All patients provided informed written consent for their tissues to be used for scientific research and to publish their case details.
Seventy-eight PDAC tissue samples and neighbouring noncancerous pancreatic tissue samples (collected postoperatively from September 2010 to August 2011) used in this study were obtained from the Department of General Surgery of Ruijin Hospital at Shanghai Jiaotong University (Shanghai, China). The specimens were obtained from patients undergoing PDAC resection with curative intent. All diagnoses were based on pathological and/or cytological evidence. The histological features of the specimens were evaluated by a senior pathologist according to the WHO (World Health Organization) classification criteria. The tissues were obtained before chemotherapy and radiation therapy. Upon removal of the surgical specimen, research personnel immediately transported the tissue to the surgical pathology lab. Pathology faculty performed a gross analysis of the specimen and selected pancreatic tissues that appeared to be cancerous and pancreatic tissues that appeared to be normal for analysis. Each sample was immediately frozen in liquid nitrogen and stored at −80°C prior to RNA isolation and qRT-PCR analysis. A second level of quality control was performed on the adjacent benign tissues. Histological slides were prepared from the section of frozen tissue that was directly adjacent to the tissue from which the RNA was isolated. These slides were examined by experienced pathologists to determine if the benign tissues contained any pancreatic tumour cells. Benign tissues that contained residual tumour tissues were not included in the study. Complete clinicopathological follow-up data of the PDAC patients from which the specimens were collected were available.
Validation of the most up-regulated or down-regulated miRNAs using qRT-PCR
Total RNA was isolated from the frozen tissue sample with TRIzol (Invitrogen) according to the manufacturer’s instructions. First-strand complementary DNA (cDNA) was synthesised from 2 μg of the total RNA using an oligo-dT primer and superscript II reverse transcriptase (Invitrogen). Then, quantification of the most up-regulated or down-regulated miRNAs was performed by qRT-PCR using SYBRR Premix Ex Taq (TakaRa). The U6 primers were obtained from TakaRa. PCR was performed in a real-time PCR system (Bio-Rad) as follows: 95°C for 3 min, followed by 35 cycles of 95°C for 5 sec, 60°C for 20 sec and 72°C for 30 sec, and then 94°C for 1 min and 60°C for 1 min, with an increase of 0.5°C per cycle. The expression level values were normalised to those of the small nuclear RNA U6 as a control. Relative fold-changes of miRNA expression were calculated using the △△CT method, and the values were expressed as 2-△△CT. The primer sequences were as follows: U6, 5′-CTCGCTTCGGCAGCACA-3′ (forward), 5′-AACGCTTCACGAATTTGCGT-3′ (reverse); miR-155, 5′-cgGCGGTTAATGCTAATCGTG-3′ (forward), 5′-GTGCAGGGTCCGAGGT-3′ (reverse); miR-100, 5′-GAATTCCCATACTGGTTGGCTCCCGC-3′ (forward), 5′-CTCGAGACGAATTCAATCGAAATATTC-3′ (reverse); miR-21, 5′-ACACTCCAGCTGGGTAGCTTATCAGACTGA-3′ (forward), 5′-TGGTGTCGTGGAGTCG-3′ (reverse); miR-221, 5′-CCCAGCATTTCTGACTGTTG-3′ (forward), 5′-TGTGAGACCATTTGGGTGAA-3′ (reverse); miR-31, 5′-ACGCGGCAAGATGCTGGCA-3′ (forward), 5′-CAGTGCTGGGTCCGAGTGA-3′ (reverse); miR-143, 5′-CCTGGCCTGAGATGAAGCAC-3′ (forward), 5′-CAGTGCTGGGTCCGAGTGA-3′ (reverse); miR-23a, 5′-CTTGAACTCCTGGCCTGAAG-3′ (forward), 5′-GCCAAAGAAACACTCACAGCT-3′ (reverse); miR-217, 5′-GCGTACTGCATCAGGAACTGATTGGA-3′ (forward), 5′-GGGCACACAAAGGCAACTTTTGT-3′ (reverse); miR-148a, 5′-TCAGTGCACTACAGAACTTTGT-3′ (forward), 5′-GCTGTCAACGATACGCTACGT-3′ (reverse); miR-375, 5′-GAAGATCTTGAGGTACATCGCAGAGGCCAG-3′ (forward), 5′-CATGCCATGGGGGCCGGAGCGGAAGACCC-3′ (reverse).
Kaplan-Meier survival analysis was used to analyse the association between postoperative survival and the miRNA expression level measured by qRT-PCR, and the resulting curves were divided into two classes (high and low expression in comparison to the mean level of miRNA expression as the threshold). Survival analysis was performed for each clinical covariate to assess their influence on outcome using a log-rank test. A multivariate Cox regression model was used to adjust for competing risk factors, and the hazard ratio (HR) with a 95% confidence interval (CI) was reported as an estimate of overall survival risk. The variables that were found to be significant in univariate analysis at p<0.05 were included in the final multivariate analysis in a backwards stepwise fashion. The statistical analyses were performed using the SPSS 18.0 for Windows software package (SPSS Inc.). Differences were considered to be statistically significant when the p-value was <0.05.
Five hundred and ninety-eight relevant publications were indexed in the databases mentioned above (Scopus, GEO, PubMed and ArrayExpress). According to the inclusion criteria and the identification of duplicate publications, only fourteen independent studies [18–31] were included. However, one article was excluded for the unavailability of a ranked gene list both publically and in response to a request from the corresponding author . The selection process is shown in Figure 1. Among the analysed studies, some of the studies employed patient samples as low as 5  or 3 , which was too small to provide any reliable data. Not surprisingly, these two studies [19, 20] were the basis for excluding numerous candidates that were consistently reported as either up- or down-regulated in other studies. The most glaring example of the strategic error of including these two studies in our meta-analysis is miRNA-100, which, despite being reported to be up-regulated in 7 studies [21, 23, 24, 26, 28, 29, 31], was considered to be down-regulated in one of the aforementioned studies , which only employed 5 tumour samples. Therefore, if Ref 19 was included, miR-100 would be listed as a miRNA with an inconsistent direction and would be subsequently excluded from the list of most consistently reported miRNAs. In addition, the fold-change in this study  was very low (less than 2) and may not have been significant if a large sample size was analysed. Other examples include miR-145, miR-141, miR-379, miR-200c, and miR-125b, which were reported in an opposite direction solely in these two studies. To avoid these deviations, these two small-sample-size studies were excluded from our meta-analysis. A brief description of the eleven included studies [21–31] and the acronyms by which the studies are referred to in the following text are provided in Table 1.
The number of patients with PDAC that were investigated in these eleven studies ranged from 8 to 160 (median 47). The studies employed a diversity of microarray platforms (either commercial or custom), and the average number of miRNAs assayed was 778 (ranging from 377 to 1200; data were missing in three papers [23, 24, 29]). Only five studies [21–23, 26, 27] provided the whole list of differentially expressed miRNAs, while the others presented only a portion of their data. Our pooled dataset included a total of 538 tumour samples and 206 noncancerous control samples (at least), as in some studies, the number of noncancerous control samples was not specified .
A total of 439 differentially expressed miRNAs were reported in the eleven miRNA expression profiling studies; 254 were up-regulated and 185 were down-regulated in at least one study. Among the 439 miRNAs, 98 were reported in at least two studies; 77 (78.57%) with a consistent direction (Tables 2 and 3) and 21 with an inconsistent direction (Table 4) among the studies. Among the 77 miRNAs with a consistent direction, 50 were reported to be up-regulated (Table 2) and 27 were reported to be down-regulated (Table 3). One miRNA (miR-155) was reported in eight studies, three miRNAs (miR-21, miR-100 and miR-221) were reported in seven studies and twelve miRNAs were reported in at least five studies, with a consistent direction in all reports (Table 5). The miRNAs that were consistently reported in at least five studies are shown in Table 5. Although there were no strong disagreements between the individual miRNA profiling studies, the top lists varied considerably from study to study.
Using the Robust Rank Aggregation method, we identified a statistically significant meta-signature of 7 up- and 3 down-regulated miRNAs in PDAC samples compared to noncancerous pancreatic tissues (Table 6). All meta-signature miRNAs that reached statistical significance after Bonferroni correction were reported by at least 5 datasets. Majority of the meta-signature miRNAs belong to the broadly conserved seed family (conserved across most vertebrates and bony fish).
Our results from the vote-counting strategy were almost the same with those from the Robust Rank Aggregation method. The ten GO processes and pathways that were most strongly enriched with respect to the meta-signature miRNA candidates (miR-155, miR-100, miR-21, miR-221, miR-31, miR-143, miR-23a, miR-217, miR-148a and miR-375) are shown in Table 7.
Experimental validation of the expression levels of the most deregulated miRNAs in patients with PDAC
To determine if the ten most deregulated miRNAs from the meta-analysis (miR-155, miR-100, miR-21, miR-221, miR-31, miR-143, miR-23a, miR-217, miR-148a and miR-375) could be used as diagnostic biomarkers of PDAC, the expression levels of these miRNAs were compared between PDAC tissues and neighbouring noncancerous tissues by qRT-PCR analysis. The results showed that the expression levels of miR-155, miR-100, miR-21, miR-221, miR-31, miR-143 and miR-23a were increased, whereas the levels of miR-217, miR-148a and miR-375 were decreased in the PDAC tissues (all p<0.05). Detailed data are available in Table 8.
Determination of prognostic significance of the candidate miRNAs in PDAC
The clinicopathological characteristics of 78 PDAC patients are shown in Table 9. The expression levels of individual miRNAs along with other well-known potential prognostic clinicopathological factors, such as histology, T category, lymph node metastasis, tumour size, perineural invasion, venous invasion and margin were included in a univariate analysis. With respect to the miRNA expression levels, for the up-regulated miRNAs, a fold-change of ≥2 was defined as high expression, and a fold-change of <2 was defined as low expression; for the down-regulated miRNAs, a fold-change of ≥2 was defined as low expression, and a fold-change of <2 was defined as high expression. Patients with advanced disease (UICC stage IV and concomitance of distant metastases) were excluded because we assumed that the prognosis of these patients (n=8) is determined by the occurrence of relapse or metastasis rather than other biological characteristics, such as miRNA expression levels.
Kaplan-Meier survival analysis was used to analyse the association between postoperative survival and the miRNA expression level, and the resulting curves were divided into two classes (high and low expression in comparison with the mean level of miRNA expression as the threshold), as shown in Figure 2.
A univariate analysis using the Cox hazard regression model demonstrated that a high expression level of miR-21 (p=0.018, HR=2.610; 95% CI=1.179-5.777) and miR-155 (p=0.035, HR=2.414; 95% CI=1.064-5.478), a low expression level of miR-375 (p=0.022, HR=2.337; 95% CI=1.431-5.066), T category (p=0.039, HR=2.282; 95% CI=1.043-4.994) and margin involvement (p=0.026, HR=2.550; 95% CI=1.120-5.805) are associated with poor patient survival.
Identification of two prognostic miRNAs by multivariate analysis
In a multivariate analysis using the Cox hazard regression model, a high expression level of miR-21 (p=0.021, HR=2.599; 95% CI=1.151-5.867), a low expression level of miR-375 (p=0.034, HR=2.451; 95% CI=1.429-5.135) and margin involvement (p=0.030, HR=2.543; 95% CI=1.093-5.918) were identified as significant unfavourable prognostic factors (Table 10).
The common drawback of miRNA expression profiling studies is the lack of agreement among several studies. Differences in measurement platforms and lab protocols as well as small sample sizes can render gene expression levels incomparable. Sato et al.  and Wang et al.  systematically analysed representative miRNA profiling platforms and revealed that each platform is relatively stable in terms of its own intra-reproducibility; however, the inter-platform reproducibility among different platforms is low. Although the ideal method involves the analysis the raw miRNA expression datasets that are pooled together, such a rigorous approach is often impossible due to the unavailability of raw data and the low inter-platform concordance of results among different studies would bring difficulties to the analysis. To overcome these limitations, it might be better to analyse datasets separately and then aggregate the resulting gene lists. In this study, we used a meta-analysis approach to analyse PDAC-specific miRNAs derived from independent profiling experiments. The well-known vote-counting strategy [12, 13] and the recently published Robust Rank Aggregation method [16, 17] were employed. The core elements of the two methods were searches for the most recognised miRNAs among the included studies.
Two principal methods are used to measure miRNA expression levels: qRT-PCR and microarray hybridisation. The technological merits and drawbacks of qRT-PCR and microarrays for miRNA analysis are similar to those for RNA or genomic DNA quantification . RT-PCR, a semiquantitative method, is labour intensive and provides data for only one, or very few, miRNA(s) per assay. However, the rapid increase in the number of known miRNAs renders this method inefficient on a genomic scale, and it is most likely better used as a tool for validation rather than discovery. Microarrays are the best option for a standardised genome-wide assay that is amenable to high-throughput application . As qRT-PCR detects only preselected miRNAs, mostly the miRNAs that were shown to be differentially expressed in PDAC from normal tissue in other studies, it hinders the discovery of new miRNAs. Most importantly, the results of studies using qRT-PCR analysis [36–40] were consistent with those of microarray-based studies. In addition to the intra-platform deviations between microarray and qRT-PCR analyses , we excluded qRT-PCR-based studies and focused on studies using miRNA microarray platforms.
We identified a meta-signature of seven up- and three down-regulated miRNAs. To our knowledge, no meta-analysis of miRNA profiling studies has specifically investigated PDAC. Furthermore, this is the first study that used a combination of the two most commonly used methods in the meta-analysis of miRNA and gene profiling. To determine if the identified miRNAs could be used as diagnostic biomarkers, we experimentally validated the expression of these miRNAs in a set of PDAC samples.
There are several factors that must be considered when choosing miRNAs as candidate diagnostic biomarkers for PDAC. First, the fold-change of the biomarker should be significant enough to discriminate cancerous tissue from benign tissue. As is shown in Tables 2 and 3, the average fold changes of the 10 miRNAs identified in the microarray-based studies were all >2. In addition, the candidate miRNAs should be expressed in a majority of tissues. As was validated by qRT-PCR, the up-regulated miRNAs were all expressed in more than 85% of the samples tested (data not shown).
Second, the biological function of each individual miRNA should be thoroughly investigated. A single miRNA may have dozens of targets, and a specific mRNA may be regulated by multiple different miRNAs . A better understanding of the targets of the miRNAs would advance their use in clinical settings. As shown in Table 7, the ten most strongly enriched GO processes and pathways with respect to the meta-signature miRNA candidates were identified. In the GO processes list, regulation of transcription, DNA-dependent; positive regulation of transcription from RNA polymerase II promoter; and positive regulation of transcription, DNA-dependent were ranked as the top three, which is in accordance with the known primary functions of miRNAs [6, 7]. Pathways in cancer and Wnt signalling pathways were ranked first in the KEGG and Panther pathway lists, respectively, highlighting the essential roles of miRNAs in cancer development.
Third, there should be adequate information about the pattern of expression of the miRNAs in different types of specimens. It has been indicated that circulating miRNAs in plasma could be more tissue-specific than tumour-specific [41, 42]. In the context of the vast inconsistency between tissue-based and plasma-based results , we focused on studies that analysed miRNA expression between PDAC tissues and noncancerous pancreatic tissues in humans.
Last but not least, rigorous validation and demonstration of reproducibility in an independent cohort of patients are necessary to confirm the diagnostic value of miRNAs. With this in mind, we experimentally validated 10 candidate miRNAs in PDAC samples and confirmed that these 10 miRNAs were differentially expressed between PDAC tissues and noncancerous pancreatic tissues.
Considering that miRNA expression is able to successfully discriminate normal from cancerous pancreatic tissues, it is tempting to speculate that miRNAs could also predict cancer prognosis. However, our results do not exclude the possibility that other miRNAs are associated with prognosis, as we only studied a meta-signature of 10 miRNAs in a limited number of PDAC samples (n=78). The main reason for the possible association between miRNAs not within this meta-signature and prognosis may centre on the relatively small sample size in our study and others [25, 27]. It is quite unrealistic to include all the miRNAs in Kaplan-Meier survival analyses, as it would be very laborious and time-consuming. Thus, commonly, only the candidate miRNAs with the greatest fold changes are included. As mentioned above, although there were no strong disagreements between the individual miRNA profiling studies, the top lists varied considerably from study to study. To remedy this problem, it was critical to identify the most differentially expressed miRNAs. We used a meta-review approach, which combines the results of several individual studies to increase statistical power and to subsequently resolve the inconsistency among different profiling studies. A meta-signature of seven up- and three down-regulated miRNAs was identified. Then, in independent patient samples, miR-21, miR-31 and miR-375 were found to be associated with cancer prognosis.
From our point of view, great caution should be taken in future research in this field. To start, sample sizes should be increased to minimise random sampling error. Next, as it is impossible for every researcher to use the same platform, reliable microarray platforms should be employed in all experiments. Finally, it is advisable to obtain an integrated view of the candidate miRNAs from many studies to avoid one-sided opinions, as great discrepancies exist among the studies.
Our study presents a method to resolve the differences that exist among studies and might have some clinical significance for research on miRNAs in PDAC. The 10 identified miRNAs may be used as diagnostic biomarkers or even therapeutic targets. In addition to our meta-analysis, we performed further studies examining the expression of the candidate miRNAs in PDAC samples and confirmed miR-21, miR-31 and miR-375 as potential prognostic biomarkers for PDAC.
Hidalgo M: New insights into pancratic cancer biology. Ann Oncol. 2012, 23 (Suppl 10): 135-138. 10.1093/annonc/mds313.
Hidalgo M: Pancreatic cancer. N Engl J Med. 2010, 362: 1605-1617. 10.1056/NEJMra0901557.
Mardis ER: Applying next-generation sequencing to pancreatic cancer treatment. Nat Rev Gastroenterol Hepatol. 2012, 9: 477-486. 10.1038/nrgastro.2012.126.
Du Y, Liu M, Gao J, Li Z: Aberrant microRNAs expression patterns in pancreatic cancer and their clinical translation. Cancer Biother Radiopharm. 2013, 28: 361-369. 10.1089/cbr.2012.1389.
Costello E, Greenhalf W, Neoptolemos JP: New biomarkers and targets in pancreatic cancer and their application to treatment. Nat Rev Gastroenterol Heaptol. 2012, 9: 435-444. 10.1038/nrgastro.2012.119.
Gregory RI, Chendrimada TP, Cooch N, Shiekhattar R: Human RISC couples microRNA biogenesis and posttranscriptional gene silencing. Cell. 2005, 123: 631-640. 10.1016/j.cell.2005.10.022.
Yates LA, Norbury CJ, Gilbert RJ: The long and short of microRNAs. Cell. 2013, 153: 516-519. 10.1016/j.cell.2013.04.003.
Singh R, Mo YY: Role of microRNAs in breast cancer. Cancer Biol Ther. 2013, 14: 201-212. 10.4161/cbt.23296.
Baer C, Claus R, Plass C: Genome-wide epigenetic regulation of miRNAs in cancer. Cancer Res. 2013, 73: 473-477. 10.1158/0008-5472.CAN-12-3731.
Song S, Ajani JA: The role of microRNAs in cancers of upper gastrointestinal tract. Nat Rev Gastroenterol Hepatol. 2013, 10: 109-118.
Lee HK, Hsu AK, Sajdak J, Qin J, Pavlidis P: Coexpression analysis of human genes across many microarray data sets. Genome Res. 2004, 14: 1085-1094. 10.1101/gr.1910904.
Griffith OL, Melck A, Jones SJ, Wiseman SM: Meta-analysis and meta-review of thyroid cancer gene expression profiling studies identifies important diagnostic biomarkers. J Clin Oncol. 2006, 24: 5043-5051. 10.1200/JCO.2006.06.7330.
Chan SK, Griffith OL, Tai IT, Jones SJ: Meta-analysis of colorectal cancer gene expression profiling studies identifies consistently reported candidate biomarkers. Cancer Epidemiol Biomarkers Prev. 2008, 17: 543-552. 10.1158/1055-9965.EPI-07-2615.
Ma Y, Zhang P, Yang J, Liu Z, Yang Z, Qin H: Candidate microRNA biomarkers in human colorectal cancer: systematic review profiling studies and experimental validation. Int J Cancer. 2012, 130: 2077-2087. 10.1002/ijc.26232.
Guan P, Yin Z, Li X, Wu W, Zhou B: Meta-analysis of human lung cancer microRNA expression profiling studies comparing cancer tissues with normal tissues. J Exp Clin Cancer Res. 2012, 31: 54-10.1186/1756-9966-31-54.
Kolde R, Laur S, Adler P, Vilo J: Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012, 28: 573-580. 10.1093/bioinformatics/btr709.
Võsa U, Vooder T, Kolde R, Vilo J, Metspalu A, Annilo T: Meta-analysis of microRNA expression in lung cancer. Int J Cancer. 2013, 132: 2884-2893. 10.1002/ijc.27981.
Singh S, Chitkara D, Kumar V, Behrman SW: Mahato RI:miRNA profiling in pancreatic cancer and restoration of chemosensitivity. Cancer Lett. 2012, 12: 00596-4.
Munding JB, Adai AT, Maghnouj A, Urbanik A, Zöllner 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. Int J Cancer. 2012, 131: E86-E95. 10.1002/ijc.26466.
Ma Y, Yu S, Zhao W, Lu Z, Chen J: miR-27a regulates the growth, colony formation and migration of pancreatic cancer cells by targeting Sprouty2. Cancer Lett. 2010, 298: 150-158. 10.1016/j.canlet.2010.06.012.
Szafranska AE, Davison TS, John J, Cannon T, Sipos B, Maghnouj A, Labourier E, Hahn SA: MicroRNA expression alterations are linked to tumorigenesis and non-neoplastic processes in pancreatic ductal adenocarcinoma. Oncogene. 2007, 26: 4442-4452. 10.1038/sj.onc.1210228.
Piepoli A, Tavano F, Copetti M, Mazza T, Palumbo O, Panza A, di Mola FF, Pazienza V, Mazzoccoli G, Biscaglia G, Gentile A, Mastrodonato N, Carella M, Pellegrini F, di Sebastiano P, Andriulli A: Mirna expression profiles identify drivers in colorectal and pancreatic cancers. PLoS One. 2012, 7: e33663-10.1371/journal.pone.0033663.
Bauer AS, Keller A, Costello E, Greenhalf W, Bier M, Borries A, Beier M, Neoptolemos J, Büchler M, Werner J, Giese N, Hoheisel JD: Diagnosis of pancreatic ductal adenocarcinoma and chronic pancreatitis by measurement of microRNA abundance in blood and tissue. PLoS One. 2012, 7: e34151-10.1371/journal.pone.0034151.
Lee EJ, Gusev Y, Jiang J, Nuovo GJ, Lerner MR, Frankel WL, Morgan DL, Postier RG, Brackett DJ, Schmittgen TD: Expression profiling identifies microRNA signature in pancreatic cancer. Int J Cancer. 2007, 120: 1046-1054.
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. 10.1001/jama.297.17.1901.
Schultz NA, Werner J, Willenbrock H, Roslind A, Giese N, Horn T, Wøjdemann M, Johansen JS: MicroRNA expression profiles associated with pancreatic adenocarcinoma and ampullary adenocarcinoma. Mod Pathol. 2012, 25: 1609-1622. 10.1038/modpathol.2012.122.
Jamieson NB, Morran DC, Morton JP, Ali A, Dickson EJ, Carter CR, Sansom OJ, Evans TR, McKay CJ, Oien KA: MicroRNA molecular profiles associated with diagnosis, clinicopathologic criteria, and overall survival in patients with resectable pancreatic ductal adenocarcinoma. Clin Cancer Res. 2012, 18: 534-545. 10.1158/1078-0432.CCR-11-0679.
Panarelli NC, Chen YT, Zhou XK, Kitabayashi N, Yantiss RK: MicroRNA expression aids the preoperative diagnosis of pancreatic ductal adenocarcinoma. Pancreas. 2012, 41: 685-690.
Ali S, Saleh H, Sethi S, Sarkar FH, Philip PA: MicroRNA profiling of diagnostic needle aspirates from patients with pancreatic cancer. Br J Cancer. 2012, 107: 1354-1360. 10.1038/bjc.2012.383.
Zhang S, Hao J, Xie F, Hu X, Liu C, Tong J, Zhou J, Wu J, Shao C: Downregulation of miR-132 by promoter methylation contributes to pancreatic cancer development. Carcinogenesis. 2011, 32: 1183-1189. 10.1093/carcin/bgr105.
Nagao Y, Hisaoka M, Matsuyama A, Kanemitsu S, Hamada T, Fukuyama T, Nakano R, Uchiyama A, Kawamoto M, Yamaguchi K, Hashimoto H: Association of microRNA-21 expression with its targets, PDCD4 and TIMP3, in pancreatic ductal adenocarcinoma. Mod Pathol. 2012, 25: 112-121. 10.1038/modpathol.2011.142.
Sato F, Tsuchiya S, Terasawa L, Tsujimoto G: Intra-platform repeatability and inter-platform comparability of microRNA microarray technology. Plos One. 2009, 4: e5540-10.1371/journal.pone.0005540.
Wang B, Howel P, Bruheim S, Ju J, Owen LB, Fodstad O, Xi Y: Systematic evaluation of three microRNA profiling platforms: microarray, breads array, and quantative real-time PCR array. Plos One. 2011, 6: e17167-10.1371/journal.pone.0017167.
Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C: Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010, 16: 991-1006. 10.1261/rna.1947110.
Etienne W, Meyer MH, Peppers J, Meyer RA: Comparison of mRNA gene expression by RT-PCR and DNA microarray. Biotechniques. 2004, 36: 618-620.
Zhang XJ, Ye H, Zeng CW, He B, Zhang H, Chen YQ: Dysregulation of miR-15a and miR-214 in human pancreatic cancer. J Hematol Oncol. 2010, 3: 46-10.1186/1756-8722-3-46.
Tavano F, di Mola FF, Piepoli A, Panza A, Copetti M, Burbaci FP, Latiano T, Pellegrini F, Maiello E, Andriulli A, di Sebastiano P: Changes in miR-143 and miR-21 expression and clinicopathological correlations in pancreatic cancers. Pancreas. 2012, 41: 1280-1284. 10.1097/MPA.0b013e31824c11f4.
Greither T, Grochola LF, Udelnow A, Lautenschlager C, Wurl P, Taubert H: Elevated expression of microRNAs 155, 203, 210 and 222 in pancreatic tumors is associated with poorer survival. Int J Cancer. 2010, 126: 73-80. 10.1002/ijc.24687.
Papaconstantinou IG, Manta A, Gazouli M, Lyberopoulou A, Lykou dis PM, Polymeneas G, Voros D: Expression of microRNAs in patients with pancreatic cancer and its prognostic significance. Pancreas. 2013, 42: 67-71. 10.1097/MPA.0b013e3182592ba7.
Ikenaga N, Ohuchida K, Mizumoto K, Yu J, Kayashima T, Sakai H, Fujita H, Nakata K, Tanaka M: MicroRNA-203 expression as a new prognostic marker of pancreatic adenocarcinoma. Ann Surg Oncol. 2010, 17: 3210-3218.
Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, Ferracin M, Dumitru CD, Shimizu M, Zupo S, Dono M, Alder H, Bullrich F, Negrini M, Croce CM: An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci USA. 2004, 101: 9740-9744. 10.1073/pnas.0403293101.
Babak T, Zhang W, Morris Q, Blencowe BJ, Hughes TR: Probing microRNAs with microarrays: tissue specificity and functional inference. RNA. 2004, 10: 1813-1819. 10.1261/rna.7119904.
This work was supported by National Natural Science Foundation of China (grant no. 81272747). The funding sources had no role in the study design, data collection, analysis or interpretation, or the writing of this manuscript. The authors thank the Department of General Surgery of Ruijin Hospital for providing the PDAC tissue samples and Dr. Fei Yuan for the pathology assessments.
The authors declare that they have no competing interests.
MZM, XK and MZW conceived the study and participated in the data collection and analysis. MZM, XK and MZW performed the experiments. MZM and KX analysed the data. MZM, XK, ZWQ, WG and CHP wrote the paper. All authors read and approved the final manuscript.
Ming-Zhe Ma, Xiang Kong, Ming-Zhe Weng contributed equally to this work.
About this article
Cite this article
Ma, MZ., Kong, X., Weng, MZ. et al. Candidate microRNA biomarkers of pancreatic ductal adenocarcinoma: meta-analysis, experimental validation and clinical significance. J Exp Clin Cancer Res 32, 71 (2013). https://doi.org/10.1186/1756-9966-32-71