Ganepola GA, Rutledge JR, Suman P, Yiengpruksawan A, Chang DH. Novel blood-based microRNA biomarker panel for early diagnosis of pancreatic cancer. World J Gastrointest Oncol 2014; 6(1): 22-33
Corresponding Author of This Article
David H Chang, PhD, Research Scientist, Center for Cancer Research and Genomic Medicine, The Daniel and Gloria Blumenthal Cancer Center, The Valley Hospital, 1 Valley Health Plaza, Paramus, NJ 07652, United States. email@example.com
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
World J Gastrointest Oncol. Jan 15, 2014; 6(1): 22-33 Published online Jan 15, 2014. doi: 10.4251/wjgo.v6.i1.22
Novel blood-based microRNA biomarker panel for early diagnosis of pancreatic cancer
Ganepola AP Ganepola, John R Rutledge, Paritosh Suman, Anusak Yiengpruksawan, David H Chang
Ganepola AP Ganepola, John R Rutledge, Paritosh Suman, David H Chang, Center for Cancer Research and Genomic Medicine, The Valley Hospital, Paramus, NJ 07652, United States
Ganepola AP Ganepola, Anusak Yiengpruksawan, Department of Surgery, The Valley Hospital, Ridgewood, NJ 07450, United States
ORCID number: $[AuthorORCIDs]
Supported by The Valley Hospital Foundation Research Fund and private donations
Author contributions: Ganepola GAP, Suman P, Yiengprukswan A and Chang DH designed the research; Ganepola GAP and Chang DH performed the experiments; Ganepola GAP, Chang DH and Rutledge JR analyzed the data and wrote the manuscript.
Correspondence to: David H Chang, PhD, Research Scientist, Center for Cancer Research and Genomic Medicine, The Daniel and Gloria Blumenthal Cancer Center, The Valley Hospital, 1 Valley Health Plaza, Paramus, NJ 07652, United States. firstname.lastname@example.org
Telephone: +1-201-6345542 Fax: +1-201-6345383
Received: August 22, 2013 Revised: December 5, 2013 Accepted: December 12, 2013 Published online: January 15, 2014
AIM: To develop a panel of blood-based diagnostic biomarkers consisting of circulating microRNAs for the detection of pancreatic cancer at an early stage.
METHODS: Blood-based circulating microRNAs were profiled by high throughput screening using microarray analysis, comparing differential expression between early stage pancreatic cancer patients (n = 8) and healthy controls (n = 11). A panel of candidate microRNAs was generated based on the microarray signature profiling, including unsupervised clustering and statistical analysis of differential expression levels, and findings from the published literature. The selected candidate microRNAs were then confirmed using TaqMan real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) to further narrow down to a three-microRNA diagnostic panel. The three-microRNA diagnostic panel was validated with independent experimental procedures and instrumentation of RT-qPCR at an independent venue with a new cohort of cancer patients (n = 11), healthy controls (n = 11), and a group of high risk controls (n = 11). Receiver operating characteristic curve analysis was performed to assess the diagnostic capability of the three-microRNA panel.
RESULTS: In the initial high throughput screening, 1220 known human microRNAs were screened for differential expression in pancreatic cancer patients versus controls. A subset of 42 microRNAs was then generated based on this data analysis and current published literature. Eight microRNAs were selected from the list of 42 targets for confirmation study, and three-microRNAs, miR-642b, miR-885-5p, and miR-22, were confirmed to show consistent expression between microarray and RT-qPCR. These three microRNAs were then validated and evaluated as a diagnostic panel with a new cohort of patients and controls and found to yield high sensitivity (91%) and specificity (91%) with an area under the curve of 0.97 (P < 0.001). Compared to the CA19-9 marker at 73%, the three-microRNA panel has higher sensitivity although CA19-9 has higher specificity of 100%.
CONCLUSION: The identified panel of three microRNA biomarkers can potentially be used as a diagnostic tool for early stage pancreatic cancer.
Core tip: This study employed high throughput screening as a screening tool to identify blood-based circulating microRNA markers for detection of early stage pancreatic cancer. Two levels of confirmation were performed to ensure the validity of the identified microRNA targets. First, a panel of potential microRNA markers was generated and confirmed using a more specific and sensitive secondary assay, real-time quantitative reverse transcription polymerase chain reaction. Second, the confirmed panel of microRNA markers was independently validated with different experimental procedures and instruments, by independent personnel, and at a different institution, to diagnose a new cohort of patients and controls.
Citation: Ganepola GA, Rutledge JR, Suman P, Yiengpruksawan A, Chang DH. Novel blood-based microRNA biomarker panel for early diagnosis of pancreatic cancer. World J Gastrointest Oncol 2014; 6(1): 22-33
Pancreatic cancer is one of the most lethal human cancers and continues to be a major unsolved health problem. It has a five-year survival rate of only 6% and is estimated to have 43920 new cases and cause 37390 deaths in 2012 in the United States, a number that has been steadily increasing for more than a decade[2,3]. Conventional treatment approaches such as surgery, radiation, chemotherapy, or a combination thereof have had little impact on the course of this aggressive cancer. Collective studies from Japan indicate that those patients who were incidentally found through other imaging modalities to have early stage carcinoma have an improved five-year survival rate of 30% for those with a 2 cm carcinoma, 57% for those with a 1 cm or less “minute” carcinoma, and 100% for patients with a ductal epithelium tumor measuring less than 1 cm[4-6]. Therefore, these studies emphasize that the hope for better control of this disease is through diagnosis at its earlier stages when surgical resection may be curative.
The current most widely used biomarkers for pancreatic cancer are CA19-9, and, to a lesser degree, carcinoembryonic antigen (CEA) and genetic markers such as K-RAS and p53. Whereas clinicians may rely on CA19-9 levels as a prognostic tool when managing patients with late stage disease, or in determining operability or monitoring patients for recurrence, these markers have generally inadequate specificity and unreliable sensitivity to pancreatic cancer and are not recommended for screening and diagnosis of early disease[9,10]. Consequently, there is an urgent, unmet need for development of valid, reliable biomarkers for early detection and monitoring of pancreatic cancer.
MicroRNA (miRNA) are small non-coding RNA about 18-24 nucleotides in size. A large body of evidence indicates that miRNAs regulate gene expression at the post-translational level in almost every biological event and play important roles in tumorigenesis, cancer development, migration and metastasis. The differential expression of miRNAs has been related to various cancers, and efforts have been made to profile the global miRNA expression patterns associated with these cancers. Numerous investigations have evaluated the miRNA signature of pancreatic cancer, utilizing pancreatic tumor tissue and cell lines, searching for biomarkers and their association with tumorigenesis and patient survival[14-20]. However, tissue-based miRNA signature profiling is limited to the availability of tissue specimens. Therefore, it would be technologically challenging to detect cancer at its earlier stages when the tumor size is still small and proper tissue procurement may be difficult. Thus, a simple, noninvasive procedure, such as blood-based signature profiling, would be ideal for detecting pancreatic cancer at earlier stages.
Recent studies have shown that miRNAs are relatively stable and can be readily extracted and detected in bodily fluids such as blood plasma. Therefore, the presence of circulating extracellular miRNAs can potentially be used as markers for cancer detection in a noninvasive way. Studies investigating the potential role of circulating miRNAs as pancreatic cancer markers have shown some promise although the findings have been limited to a small number of predefined target miRNAs[23-25]. High throughput screening studies of known circulating miRNAs as biomarkers of pancreatic cancer have only recently begun to emerge, and by combining the efforts from all researchers in this area, a final panel of ideal markers can be developed to combat this dreadful disease[26-28].
This is a pilot study to employ high throughput screening microarray technology to screen for 1220 known human miRNAs (based on the miRBase version 16.0 database released in 2010), comparing the differential expression signature among eight early stage pancreatic cancer patients and eleven healthy controls. The identified panel of miRNAs were subsequently confirmed, and also validated for their diagnostic potential in an independent cohort of eleven early stage pancreatic cancer patients.
MATERIALS AND METHODS
This study was reviewed and approved by the Valley Hospital Institutional Review Board. Written informed consent was obtained from all study participants. The pancreatic cancer patient group included stage IIA/IIB patients whose stage was confirmed post-operatively by pathologists from the Valley Hospital. Eight patients diagnosed with ductal adenocarcinoma were used in the microarray screening and real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) confirmation study. A second group of eleven ductal adenocarcinoma patients was used in the validation study. The control group for all studies was comprised of eleven healthy participants with no family history of pancreatic cancer. The high risk group included eleven healthy participants who had a strong family history of pancreatic cancer, including ten participants with at least two “first degree relatives” and one participant with two “second degree relatives” having a pancreatic cancer diagnosis. The participant demographics are shown in Table 1.
Patients were excluded from the study if they had prior pancreatic cancer surgery, had other concomitant cancers other than non-melanoma basal cell skin cancer, or had a history of HIV infection.
Blood specimen collection and processing
Patient blood was drawn by peripheral venipuncture into BD Vacutainer® CPT™ (Cell Preparation Tubes) with Sodium Citrate (Becton Dickinson, Franklin Lakes, NJ). Blood processing was typically done within two hours of collection and performed according to the manufacturer’s protocol. Harvested plasma was stored at -80 °C.
MiRNA microarrays and analysis (performed at Ocean Ridge Biosciences, FL)
Plasma samples were processed at Ocean Ridge Biosciences (ORB, Palm Beach Gardens, FL) for analysis using custom multi-species microarrays containing 1209 probes covering 1220 human mature miRNAs present in the miRBase version 16.0 database released in 2010. The sensitivity of the microarray is such that it could detect as low as 20 amoles of synthetic miRNA being hybridized along with each sample. The microarrays were produced by Microarrays Inc. (Huntsville, Alabama), and consisted of epoxide glass substrates that have been spotted in triplicate with each probe.
Sample processing: Samples were isolated from 0.7 to 1.0 mL of plasma using TRI Reagent® BD (Molecular Research Center, Cincinnati, OH) as per manufacturer instructions. For quality control, a mixture of 10 synthetic miRNAs were added (spike-in) at a mass of 12.5 fmoles/mL of plasma to each plasma sample during isolation and one miRNA was added at 200 amoles per sample prior to labeling and hybridization. Total RNA was 3’-end labeled with Oyster-550 fluorescent dye using the Flash Tag RNA labeling Kit (Genisphere, Hatfield, PA). Labeled RNA samples were hybridized to the miRNA microarrays according to conditions recommended in the Flash Tag RNA labeling Kit manual. The microarrays were scanned on an Axon Genepix 4000B scanner, and data was extracted from images using GenePix v4.1 software.
Data pre-processing: Spot intensities were obtained for the 8816 features on each microarray by subtracting the median local background from the median local foreground for each spot. The 95th percentile of the negative control spots was also calculated for each array. The spot intensities and 95th percentile of negative controls (TPT95) were transformed by taking the Logarithm base 2 (indicated as log2) of each value. The normalization factor (N) for each microarray was determined by obtaining the 20% trimmed mean of the human probe intensities that were detected one log2 unit above TPT95 (TPT95 + 1) in all samples and with standard deviation of probe intensities among all samples less than 1.25. The log2-transformed spot intensities for all 8816 features were normalized, by subtracting N from each spot intensity, and scaled by adding the grand mean of N across all microarrays. The mean probe intensities for each of the 1209 human probes on each of the 20 arrays were then determined by averaging the triplicate spot intensities. Spots flagged as poor quality during data extraction were omitted prior to averaging. The 1209 human non-control log2-transformed, normalized, and averaged probe intensities were filtered to obtain a panel of 290 human miRNA probes showing probe intensity greater than one log2 unit above TPT95 (> TPT95 + 1) in at least 10% of the samples.
Quality control: Sensitivity of the microarray hybridization was confirmed by detection of hybridization signal for all 11 spikes that were added during isolation and labeling well above TPT95. The array also contains a set of specificity control probes complementary to three different miRNAs. Each specificity control includes a perfect match, single mismatch, double mismatch, and shuffled version of the probe. Specificity of the hybridization was confirmed by detection of hybridization signal on the microarray for the perfect match probes and not the double mismatch and shuffled version of the probes. Reproducibility of the arrays was determined by monitoring the hybridization intensity for the triplicate human spots on each array.
Differential expression analysis: For statistical analysis, samples were binned in two groups: Healthy Controls and Pancreatic Cancer. The log2-transformed and normalized spot intensities for the 290 detectable human probes were examined for differences between the groups by 1-way ANOVA using National Institute of Ageing Array Analysis software. This ANOVA was conducted using the Bayesian Error Model and 50 degrees of freedom. A total of 116 probes exhibited significant differences between the healthy controls and the pancreatic cancer patients. Statistical significance was determined using the False Discovery Rate (FDR) method; an FDR < 0.15 was considered significant in this study.
Hierarchical clustering of miRNA array data: Data for the 290 detectable human probes were clustered using Cluster 3.0 software. Genes were median-centered prior to hierarchical clustering. Hierarchical clustering was conducted using Centered Correlation as the similarity metric and Average Linkage as the clustering method. Intensity scale shown is arbitrary.
Confirmation with quantitative real-time RT-PCR (RT-qPCR, performed at ORB)
Total RNA was extracted from 0.9 to 2.1 mL of plasma using TRI Reagent® BD (Molecular Research Center, Cincinnati, OH) according to the manufacturer’s instructions with minor modifications. For down-stream quality control monitoring, a mixture of 10 synthetic miRNAs were added at a mass of 12.5 femtomoles/mL of biofluid following homogenization of the samples in Trizol-BD.
RNA isolated from plasma was diluted and used in 10 μL reverse transcription (RT) reactions using ABI (Applied Biosystems by Life Technology) miRNA-specific RT primers (Life Tech, Carlsbad, CA). For each miRNA and sample, RNA equivalent to 50 μL of biofluid was used in each RT reaction except for the probe sets for miR-642b-3p and miR-7, in which case RNA equivalent to 100 μL was used. The cDNA product was diluted 1:10 in water and 4.5 μL of the diluted product was combined in triplicate PCR-plate wells with 1 × ABI Universal PCR amplification mix and ABI miRNA-specific Taqman PCR primers in a final volume of 10 μL/well. The PCR plate was subjected to thermal cycling in an ABI StepOne Plus real-time PCR instrument. The cycling conditions were an initial incubation for 10 min at 95 °C, followed by 40 cycles of: 15 s at 95 °C, 1 min at 60 °C. The Cq (Quantification Cycle according to Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines, also known as threshold cycle (CT), or number of cycles required for the well to reach a specified threshold of fluorescence intensity) was determined for each well using ABI software and a threshold setting of 0.05. In the initial pilot experiment, only miRNA with Cq value at 36 or above were used for confirmation study. Triplicate Cq values were averaged to obtain the “mean Cq value.”
Independent validation RT-qPCR (performed at the Valley Hospital, NJ)
MiRNA was extracted from 2.0 to 4.0 mL of plasma using QIAamp circulating Nucleic Acid Kit (Qiagen, Germantown, MD) with QIAvac Connecting System (Qiagen), following manufacturer instructions. The concentration of extracted miRNA was determined using Agilent Small RNA kit with 2100 Bioanalyzer (Agilent, Santa Clara, CA), and was in the range of 10 to 100 ng.
RT-qPCR experiments were performed using TaqMan MiRNA Assays (Life Tech, Carlsbad, CA), for each miRNA, following manufacturer instructions. The TaqMan miRNA assays used were: hsa-miR-885-5p (ID#002296), hsa-miR-22-3p (ID#000398), hsa-miR-642b-3p (ID#462949_mat), and hsa-miR-3196 (ID#241941_mat). For RT, 1 to 10 ng RNA isolated from plasma was diluted and used in 5 μL of RT reactions with each miRNA-specific RT primers and MiRNA Reverse Transcription kit (Life Tech), without pre-amplification. The cDNA product was diluted 1:3 in water, and 5 μL of diluted cDNA was used in triplicate PCR-plate wells with miRNA-specific TaqMan PCR primer and Universal Master mix II (Life Tech). The qPCR reaction was performed on 7500 real-time PCR system (Life Tech). The Cq value required for the well to reach a specific threshold of fluorescence intensity was determined for each well using 7500 Software (Life Tech) with threshold setting at 0.05. All Cq values were 36 or less and were included for subsequent calculation. Triplicate Cq values were averaged to obtain the “mean Cq value”.
RT-qPCR and microarray data analysis using the comparative C(T) Method
The normalization control used for the RT-qPCR experiments was hsa-miR-3196. The endogenous control has been selected based on both microarray and pilot RT-qPCR experiments. MiR-3196 demonstrated relatively consistent and stable expression in microarray analysis and also in a series of pilot RT-qPCR. The use of the same extracted plasma miRNA (patients and controls, n = 19), with the same miRNA input, yielded consistent Cq values at 29.33 ± 0.50.
Fold change calculation was performed following the Comparative C(T) Method. The “mean Cq value” of hsa-miR-3196 was subtracted from the mean Cq value for each miRNA to obtain a “normalized mean ΔCq value” for each miRNA. The “normalized mean ΔCq value” for each miRNA was then converted to 2-(ΔCq). Fold change of each miRNA was determined by calculating the ratio of the mean of 2-(ΔCq) of all pancreatic cancer patients to the mean of 2-(ΔCq) of all healthy controls; and the mean of 2-(ΔCq) of all high risk controls to the mean of 2-(ΔCq) of all healthy controls. The fold change values were transformed to a log2 scale for the purpose of plotting the values on a continuum.
In order to further evaluate the RT-qPCR experimental results against those observed from the microarray analysis, microarray data for eight miRNAs (miR-642b-3p, miR-762, miR-4253, miR-885-5p, miR-18a, miR-486-5p, miR-7, and miR-22) were re-analyzed following the Comparative C(T) Method. The spot intensities of all data points (in log2) were normalized by subtracting the values for each miRNA from miR-3196 (as normalization control) to generate “normalized log2 values”, which were used as the equivalent of “normalized mean ΔCq value” in RT-qPCR for further calculation. The fold change values are then calculated as stated above for RT-qPCR data analysis.
The CA19-9 serum marker test was performed by the Valley Hospital Histology Laboratory following the routine diagnostic laboratory testing protocol using Tosoh A1A-360 Immunoanalyzer (Tosoh Bioscience, Inc.). CA19-9 range above 47 U/mL is considered to be a positive test result.
The Wilcoxon-Mann-Whitney test was used to compare the fold changes observed among the patient groups for each miRNA. Receiver operating characteristic (ROC) curves were generated to evaluate the sensitivity and specificity of predicting cancer cases and control cases for each miRNA and for the combination of miRNAs. Area under the curve and their respective 95%CI were calculated for all ROC curves. P-values ≤ 0.05 were considered to be statistically significant. All analyses were done using IBM-SPSS software (Version 19).
The miRNA expression signature was profiled using custom miRNA microarray chips covering 1220 human miRNAs derived from the miRBase database, version 16, released in 2010. The unsupervised hierarchical clustering of the 290 miRNAs with acceptable detection intensities is shown in Figure 1A. Remarkably, as depicted in the hierarchical dendrogram, the clustering pattern clearly separates the pancreatic cancer from most healthy controls.
Figure 1 The microarray signature profile of circulating microRNA in pancreatic cancer patients and healthy controls.
A: The unsupervised hierarchical clustering of 290 miRNAs that are differentially expressed among pancreatic cancer patients (n = 8) and healthy controls (n = 11). The dendrogram on top indicates the hierarchical clustering relationship between pancreatic cancer patients (in red) and healthy controls (in blue); B: Heat map depicting the subset of 42 miRNAs chosen for confirmation RT-qPCR study, arranged as indicated. Heat map color scale represents fold increase (red) or decrease (green) from baseline. MiRNA: MicroRNA.
To compile a set of targeted miRNA markers for further investigation and confirmation, the list of 290 miRNAs was narrowed down to 31 miRNAs with significant differential expression and combined with an additional 22 miRNAs, which have been shown to be candidate biomarkers for pancreatic cancer by other investigators[14,17-20,25]. The latter criterion was added to ensure the inclusion of miRNAs which may not have exhibited significant differential expression in this experimental setting but have otherwise been shown by others to be good candidates. As shown in Table 2, among the final set of 42 miRNAs, 11 miRNAs exhibiting the most significantly different expression levels (FDR ≤ 0.15) in this study (miR-194, miR-18a, miR-7, miR-26b, miR-301a, miR-106b, miR-16, miR-93, miR-106a, miR-19b, and let-7i) as well as in the previously published literature. The heat map of the 42 miRNAs in Figure 1B demonstrates a differential expression signature between pancreatic cancer patients and controls similar to the results observed from the unsupervised clustering in Figure 1A.
The 42 final target list based on the 290 miRNAs identified by miRNA microarray. The fold change shown is calculated as the ratio of pancreatic cancers over healthy controls (under “Up-regulated in Pancreatic Cancer”) or healthy controls over pancreatic cancers (under “Up-regulated in Healthy Controls”). MiRNA: MicroRNA; FDR: False discovery rate.
RT-qPCR confirmed three potential miRNA diagnostic markers
The panel of 42 miRNAs were subjected to further investigation using real time quantitative RT-PCR (RT-qPCR). In the initial pilot experiment using available miRNA TaqMan probes and testing the relative abundance of expression, 8 of the 42 miRNAs gave acceptable and most consistent signals (data not shown) and therefore were chosen for the subsequent confirmation study.
As shown in Figure 2, when comparing microarray and TaqMan RT-qPCR results side-by-side after normalizing both data sets with the same control, miR-3196, the relative expression level of each miRNA varied due to the differing nature of the experimental technology. However, three of the eight miRNAs, miR-642b-3p, miR-885-5p, and miR-22-3p, demonstrate consistent outcomes across the two methodologies with higher expression seen in pancreatic cancer patients than in healthy controls. Therefore, they were chosen for further validation study as a potential diagnostic panel.
Figure 2 Confirmation study of the expression profile of eight microRNAs.
Microarray (solid bar) and RT-qPCR (open bar) data sets were analyzed based on the “Comparative C(T) Method”, normalized to miR-3196, and calculated as the ratio of mean fold change of pancreatic cancer patients (n = 8) to healthy controls (n = 11). MiRNA: MicroRNA.
The validation study was performed on a new independent cohort of 11 pancreatic cancer patients, 11 healthy controls and 11 high risk controls. As shown in Figure 3A, the three miRNAs, plotted using the raw data of “normalized mean ΔCq value”, differentiate pancreatic cancer patients from both the healthy and high risk controls. Likewise, as shown in Figure 3B, when determining the relative fold expression change by calculating the ratio of pancreatic cancer to healthy controls, expression of the three miRNAs is clearly up regulated in pancreatic cancer patients. For the high risk controls, on the other hand, the expression level is either comparable or lower when compared to the healthy controls. Furthermore, the relative expression levels of all three miRNAs are remarkably consistent in direction and magnitude between confirmation RT-qPCR (Figure 2) and validation RT-qPCR (Figure 3B), despite the use of different patient samples and independent experimental procedures. This lends support to the suggestion that they could play potential roles as diagnostic biomarkers for pancreatic cancer.
Figure 3 Validation of three-microRNA panel using quantitative real-time reverse transcription polymerase chain reaction.
A: Box plot of relative expression of three miRNAs, based on normalized mean ΔCq values [mean quantification cycle (Cq) normalized to miR-3196] of pancreatic cancer patients (PC, n = 11), healthy controls (HC, n = 11), and high risk controls (HR, n = 11). The whiskers extend to the observations which are no more than 1.5 times the length of the box (interquartile range) away from the box. More extreme observations are considered outliers and are indicated as dots. "a" indicates P value ≤ 0.017 (0.05/3, a Bonferroni-adjusted α-level based on the 3 multiple comparisons performed among patient groups); B: Fold change expression levels of three miRNAs calculated as the ratio of mean fold change of PC/HC and HR/HC. MiRNA: MicroRNA.
Assessing diagnostic potential of the three-miRNA panel
To assess the potential use of the three miRNA as a diagnostic panel for pancreatic cancer, ROC analysis was performed on the validation data set for each of the individual miRNAs and the combination of the three. As summarized in Table 3 and shown in Figure 4B and 4C, miR-885-5p and miR-22-3p each demonstrated high sensitivity of 82% for cancer case identification and relatively high specificity of 73% and 82%, respectively, for identifying healthy controls. For miR-642b-3p (Figure 4A), a high sensitivity of 82% for identifying cancer cases was demonstrated, but a lower specificity for identifying healthy controls of 55% was observed. When using the optimal cut-point, as shown in Figure 4D, the composite of all three miRNAs yielded both a sensitivity and specificity of 91%.
Table 3 Receiver operating characteristic curve analysis of the 3 microRNAs validated by quantitative real-time reverse transcription polymerase chain reaction for diagnosing pancreatic cancer.
Composite of 3 miRNAs
AUC: Area under the curve; MiRNA: MicroRNA.
Figure 4 Receiver operating characteristic curve analysis of three microRNAs.
Validation study with 11 stage IIA/IIB pancreatic cancer patients and 11 healthy controls. A: Receiver operating characteristic (ROC) curve for the miR-642b-3p data alone [area under the curve (AUC) = 0.79]; B: ROC curve for the miR-885-5p data alone (AUC = 0.84); C: ROC curve for the miR-22-3p data alone (AUC = 0.86); D: ROC curve for all three miRNAs (miR-642b-3p, -885-5p, 22-3p) as a composite panel (AUC = 0.97). MiRNA: MicroRNA.
Given that the CA19-9 serum marker is the only marker currently available for routine diagnostic testing, we sought to compare CA19-9 results with the three-miRNA panel for predicting pancreatic cancer. The CA19-9 test has been performed on our cancer patient sample (n = 11) and non-cancer patient samples (healthy and high-risk controls, n = 22), and the sensitivity was observed to be 73% (8 out of 11 patients), and the specificity was observed to be 100% (all controls were below the reference point).
Early cancer detection remains a major challenge in pancreatic cancer but holds promise of resulting in a more favorable disease outcome. In light of the fact that current research progress into early detection of pancreatic cancer has resulted in limited actual clinical applications using various biomarkers such as tissue nucleic acid, proteins, tumor cells, and plasma proteins, we elected to focus on a relatively new source of potential biomarker, blood-based circulating miRNA. Given this test is blood-based, it would be noninvasive and ideal for diagnosing asymptomatic cancer. We performed an array-based high throughput screening process for all known human miRNA species (released by miRBase in 2010). We employed two levels of confirmation with RT-qPCR, using two independent samples of pancreatic cancer patients studied under two different sets of experimental conditions. We have identified and validated a panel of three miRNAs (miR-642b-3p, miR-885-5p and miR-22-3p) with high combined sensitivity of 91% and specificity of 91%.
Three prior studies have performed screening of circulating blood miRNAs for pancreatic cancer. Ali et al profiled plasma miRNAs based on a pooled plasma specimen from 50 newly diagnosed pancreatic cancer patients (without specifying cancer staging). They identified miR-21 to be significantly higher and the expression of the let-7 family (especially let-7d) and miR-146a to be significantly lower in cancer. Liu et al used Illumina Sequencing by Synthesis (SBS) technology and identified seven miRNAs (miR-20a, miR-21, miR-24, miR-25, miR-99a, miR-185, miR-191) as a potential panel of biomarkers, with high sensitivity of 83.6%. However, surprisingly, the SBS technology employed did not identify any new, as yet unknown miRNA markers, considering that more mature human miRNA continue to be identified and updated yearly (a total of 2578 mature human miRNA are included in miRBase version 20 issued in 2013). Also, Liu et al screened for pancreatic cancer patients from all stages (13.2% stage I, 24.4% stage II, 22.8% stage III, 33.5% stage IV, and 6.1% stage unknown). Therefore, it is unclear whether the seven-miRNA panel should be used as pan-pancreatic cancer prognostic markers or as early stage cancer screening markers, considering the gene expression profile can change dramatically from early to late stages[35,36]. Carlsen et al used the same strategy as presented in our study and identified circulating miR-375 as the sole potential marker although with relatively low accuracy (70%), and it did not outperform CA19-9 as a diagnostic marker. However, Carlsen et al used chronic pancreatitis patients as controls instead of normal healthy donors. It is uncertain why chronic pancreatitis patients were used as controls considering only a small population (8%) of pancreatic cancer patients have concomitant chronic pancreatitis. There is a weak link between pancreatitis and pancreatic cancer, and only 4% of patients within 20 years of chronic pancreatitis diagnosis have developed into pancreatic cancer.
The research presented here focuses solely on stage II pancreatic ductal adenocarcinoma patients, of which the gene expression profile should resemble early stage pancreatic cancer more closely than advanced stage cancer. The validation experiment was done by comparing pancreatic cancer patients to healthy controls, and in particular to high risk controls who might have inherited genetic susceptibility. Our three-miRNA panel can differentiate pancreatic cancer patients from healthy controls with high sensitivity and specificity. Nevertheless, the 11 miRNAs, as presented by the three prior studies mentioned above, were also found in our 290-miRNA panel to have acceptable detection intensities (presented in Figure 1A). It would be of great interest to conduct a future larger scale clinical trial comparing the three-miRNA panel identified here to the eleven-miRNAs identified previously to see how they fare in predicting pancreatic cancer, or if the combination of all 14 miRNAs could be developed into a diagnostic test for cancer.
Of the three miRNAs identified in this study, two miRNAs (miR-642b-3p and miR-885-5p) were shown to be significantly up-regulated in cancer patients by our screening process while the third miRNA (miR-22-3p) was shown to be up-regulated in cancer patients in the literature[19,20]. It should be noted that miR-642b is a relatively novel miRNA marker with no prior publication about its potential functional role or utility as a marker. However, miR-885-5p has been shown to be a potential serum marker for liver pathologies, including hepatocellular carcinoma, liver cirrhosis, and chronic hepatitis B. Functionally, miR-885-5p is found to be located in the 3p25.3 genomic region and is known to have a tumor suppressive function by triggering cell cycle arrest and senescence and/or apoptosis. MiR-885-5p activates the p53 pathway, causes down-regulation of cyclin-dependent kinase and mini-chromosome maintenance protein, and suppresses matrix metallopeptidase 9 expression and Caspase genes[42-44]. MiR-22, on the other hand, is one of the most common miRNAs in the colorectal cancer transcriptome and has been studied for its critical role in breast cancer and bone metastasis[45,46]. It is known for directly targeting the estrogen receptor α mRNA and is proposed to be a putative tumor suppressor by repressing the EVI1 oncogene expression. It inhibits cell cycle progression by repressing Max and ErbB3 expression post-transcriptionally, mediates the effects of the tumor-suppressor p53, and suppresses interferon gene expression by blocking interferon regulatory factor-5[49-51]. MiR-22 has been proposed as a potential serum marker for non-small cell lung cancer, esophageal squamous cell carcinoma, and nasopharyngeal carcinoma.
Our experimental strategy in this study focused on using hybridization-based microarray technology as the means to screen thousands of genes simultaneously, but it is also severely limited due to its issues with reproducibility and its tendency to produce a high rate of false positive and false negative results. This restricts its potential use as a reliable diagnostic tool for cancer. RT-qPCR, on the other hand, is based on sequence-specific amplification, which is highly specific and sensitive for individual testing targets. It has been developed for use in diagnostic/prognostic tests such as the Oncotype DX test for breast cancer and the Cervista® HPV HR assay. Furthermore, our approach using two layers of confirmation with RT-qPCR, utilizing different sets of patients, independent experimental procedures and instrumentation, still showed a remarkable consistency that suggests the potential future application of RT-qPCR-based diagnostic tests using circulating miRNA markers.
Although CA19-9 is not considered an ideal biomarker for the early diagnosis of pancreatic cancer, data from this study demonstrated a relatively high sensitivity of 73% and specificity of 100%. It should be noted, however, that all patients recruited for this study were “confirmed cancer cases” by pre-surgery imaging [magnetic resonance imaging (MRI) and computed tomography (CT)] and post-surgical examination by pathologists. Therefore, even with 100% confirmed cancer cases, CA19-9 has only a 73% sensitivity to detect them. The three-miRNA panel, on the other hand, exhibited a 91% sensitivity. It would be of great interest to test the three-miRNA panel, alongside CA19-9, in a future large scale clinical trial of suspected cancer cases to see how the two-marker system fares when they are compared to one another or when they work together.
In summary, we have identified three blood-based circulating miRNA targets, miR-642b-3p, miR-885-5p and miR-22-3p, which, when combined, provided a high level of diagnostic accuracy for early stage pancreatic cancer. Our plan is to study an expanded sample of patients to further develop and refine the diagnostic miRNA panel based on RT-qPCR. This new panel may work alone or in conjunction with other known immunoassays, such as CA19-9 and CEA, as a diagnostic test for early stage pancreatic cancer. We envision that the future miRNA biomarker panel can immediately apply to the category of patients at high risk for pancreatic cancer before more expensive and invasive modalities like CT, MRI, endoscopic ultrasound and endoscopic retrograde cholangio-pancreatography are used. A similar strategy can also be utilized to identify miRNA panels for other cancer types where early detection is crucial for a favorable disease outcome.
The authors are grateful for the generous support from the community of the Valley Hospital, including the Department of Oncology Clinical Trials for patient follow-up, sample collection and processing. The authors also want to thank Dr. Lawrence Harrison of the Valley Hospital for enrolling patients, and Drs. David Willoughby and Joseph Benito and the staffs at ORB for performing microarray screening of microRNA and real time quantitative RT-PCR.
Pancreatic cancer is one of the most lethal human cancers with a mere 6% 5-year survival rate. Studies have shown that early detection is the best option available for controlling this disease. The goal of this study, therefore, is to explore and compile a diagnostic biomarker panel, based on microRNA (miRNA) in the circulating blood, for detection of pancreatic cancer at earlier stages.
A diagnostic panel of biomarkers for cancer is highly desirable because it would help to fight cancer at the earliest possible stage when the disease is still curable. However, for the past fifteen to twenty years, despite numerous studies and publications, this goal remains elusive with only scanty numbers of risk assessment and prognostic panels eventually developing into clinical tests. One of the major problems is study reproducibility. The results of one study cannot be readily reproduced by another study. In addition, gene expression profiles can vary greatly among the different stages for each type of cancer, as well as between blood-based markers and tissue-specific markers. Therefore, it is critical for each investigator to be clear on their general experimental strategy which can, in turn, address issues pertaining to the development of future clinical tests.
Innovations and breakthroughs
The current study focuses on identifying a panel of blood-based biomarkers for pancreatic cancer. If it is effective, pancreatic cancer detection can be performed by a simple non-invasive blood draw instead of an invasive procedure. The approach used in this study is innovative due to the fact that, in addition to the general strategy of microarray screening followed by polymerase chain reaction-based confirmation, the authors employed a second layer of validation experiments, using different experimental procedures, instrumentation, and lab personnel at an independent location. The results of the two experimental strategies are remarkably similar even with a new cohort of patient specimens, suggesting high validity of this diagnostic panel. Furthermore, the authors have included a group of high risk individuals as controls. High risk controls are subjects with a strong family history of pancreatic cancer (at least two first degree relatives with the disease). Therefore, they have inherited genetic susceptibility to developing pancreatic cancer and hence are genotypically closer to individuals with disease than to normal healthy controls. Remarkably, the three-miRNA panel identified in this study can differentiate pancreatic cancer patients from both normal and high risk controls, demonstrating its high sensitivity and specificity for pancreatic cancer.
The identified panel of three miRNAs can potentially be used as a diagnostic detection set for early stage pancreatic cancer.
MiRNA are small non-coding RNA approximately 18-24 nucleotides in size. The abnormal expression of miRNA found in patient blood is known to be associated with cancer progression. Hence, the utilization of miRNA biomarkers was proposed as a way to potentially detect pancreatic cancer.
This is a well written paper examining a potentially useful means of screening for pancreatic carcinoma. The authors would recommend publication of the manuscript.
P- Reviewers: Li S, Morris-Stiff G, Safe S, Specchia ML S- Editor: Wen LL L- Editor: A E- Editor: Liu SQ
Siegel R, Naishadham D, Jemal A. Cancer statistics, 2012.CA Cancer J Clin. 2012;62:10-29.
Smith BD, Smith GL, Hurria A, Hortobagyi GN, Buchholz TA. Future of cancer incidence in the United States: burdens upon an aging, changing nation.J Clin Oncol. 2009;27:2758-2765.
Tsuchiya R, Noda T, Harada N, Miyamoto T, Tomioka T, Yamamoto K, Yamaguchi T, Izawa K, Tsunoda T, Yoshino R. Collective review of small carcinomas of the pancreas.Ann Surg. 1986;203:77-81.
Ishikawa O, Ohigashi H, Imaoka S, Nakaizumi A, Uehara H, Kitamura T, Kuroda C. Minute carcinoma of the pancreas measuring 1 cm or less in diameter--collective review of Japanese case reports.Hepatogastroenterology. 1999;46:8-15.
Ariyama J, Suyama M, Satoh K, Sai J. Imaging of small pancreatic ductal adenocarcinoma.Pancreas. 1998;16:396-401.
Steinberg WM, Gelfand R, Anderson KK, Glenn J, Kurtzman SH, Sindelar WF, Toskes PP. Comparison of the sensitivity and specificity of the CA19-9 and carcinoembryonic antigen assays in detecting cancer of the pancreas.Gastroenterology. 1986;90:343-349.
Schmidt C. Early detection tools for pancreatic cancer.J Natl Cancer Inst. 2012;104:1117-1118.
Poruk KE, Gay DZ, Brown K, Mulvihill JD, Boucher KM, Scaife CL, Firpo MA, Mulvihill SJ. The clinical utility of CA 19-9 in pancreatic adenocarcinoma: diagnostic and prognostic updates.Curr Mol Med. 2013;13:340-351.
Locker GY, Hamilton S, Harris J, Jessup JM, Kemeny N, Macdonald JS, Somerfield MR, Hayes DF, Bast RC. ASCO 2006 update of recommendations for the use of tumor markers in gastrointestinal cancer.J Clin Oncol. 2006;24:5313-5327.
Sayed D, Abdellatif M. MicroRNAs in development and disease.Physiol Rev. 2011;91:827-887.
Lee YS, Dutta A. MicroRNAs in cancer.Annu Rev Pathol. 2009;4:199-227.
Navon R, Wang H, Steinfeld I, Tsalenko A, Ben-Dor A, Yakhini Z. Novel rank-based statistical methods reveal microRNAs with differential expression in multiple cancer types.PLoS One. 2009;4:e8003.
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.
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.
Giovannetti E, van der Velde A, Funel N, Vasile E, Perrone V, Leon LG, De Lio N, Avan A, Caponi S, Pollina LE. High-throughput microRNA (miRNAs) arrays unravel the prognostic role of MiR-211 in pancreatic cancer.PLoS One. 2012;7:e49145.
Greither T, Grochola LF, Udelnow A, Lautenschläger C, Würl 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.
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.
Zhang Y, Li M, Wang H, Fisher WE, Lin PH, Yao Q, Chen C. Profiling of 95 microRNAs in pancreatic cancer cell lines and surgical specimens by real-time PCR analysis.World J Surg. 2009;33:698-709.
Wang X, Zhao J, Huang J, Tang H, Yu S, Chen Y. The regulatory roles of miRNA and methylation on oncogene and tumor suppressor gene expression in pancreatic cancer cells.Biochem Biophys Res Commun. 2012;425:51-57.
Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, Galas DJ, Wang K. The microRNA spectrum in 12 body fluids.Clin Chem. 2010;56:1733-1741.
Shen J, Stass SA, Jiang F. MicroRNAs as potential biomarkers in human solid tumors.Cancer Lett. 2013;329:125-136.
Wang J, Chen J, Chang P, LeBlanc A, Li D, Abbruzzesse JL, Frazier ML, Killary AM, Sen S. MicroRNAs in plasma of pancreatic ductal adenocarcinoma patients as novel blood-based biomarkers of disease.Cancer Prev Res (Phila). 2009;2:807-813.
Ho AS, Huang X, Cao H, Christman-Skieller C, Bennewith K, Le QT, Koong AC. Circulating miR-210 as a Novel Hypoxia Marker in Pancreatic Cancer.Transl Oncol. 2010;3:109-113.
Morimura R, Komatsu S, Ichikawa D, Takeshita H, Tsujiura M, Nagata H, Konishi H, Shiozaki A, Ikoma H, Okamoto K. Novel diagnostic value of circulating miR-18a in plasma of patients with pancreatic cancer.Br J Cancer. 2011;105:1733-1740.
Liu R, Chen X, Du Y, Yao W, Shen L, Wang C, Hu Z, Zhuang R, Ning G, Zhang C. Serum microRNA expression profile as a biomarker in the diagnosis and prognosis of pancreatic cancer.Clin Chem. 2012;58:610-618.
Carlsen AL, Joergensen MT, Knudsen S, de Muckadell OB, Heegaard NH. Cell-free plasma microRNA in pancreatic ductal adenocarcinoma and disease controls.Pancreas. 2013;42:1107-1113.
Ali S, Almhanna K, Chen W, Philip PA, Sarkar FH. Differentially expressed miRNAs in the plasma may provide a molecular signature for aggressive pancreatic cancer.Am J Transl Res. 2010;3:28-47.
Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data.Nucleic Acids Res. 2011;39:D152-D157.
Sharov AA, Dudekula DB, Ko MS. A web-based tool for principal component and significance analysis of microarray data.Bioinformatics. 2005;21:2548-2549.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing.J Roy Stat Soc B. 1995;57:289-300.
de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software.Bioinformatics. 2004;20:1453-1454.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments.Clin Chem. 2009;55:611-622.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method.Nat Protoc. 2008;3:1101-1108.
Lubezky N, Loewenstein S, Ben-Haim M, Brazowski E, Marmor S, Pasmanik-Chor M, Oron-Karni V, Rechavi G, Klausner JM, Lahat G. MicroRNA expression signatures in intraductal papillary mucinous neoplasm of the pancreas.Surgery. 2013;153:663-672.
Su GH, Kern SE. Molecular genetics of ductal pancreatic neoplasia.Curr Opin Gastroenterol. 2000;16:419-425.
Dzeletovic I, Harrison ME, Crowell MD, Pannala R, Nguyen CC, Wu Q, Faigel DO. Pancreatitis Before Pancreatic Cancer: Clinical Features and Influence on Outcome.J Clin Gastroenterol. 2013;Oct 22; Epub ahead of print.
Yadav D, Lowenfels AB. The epidemiology of pancreatitis and pancreatic cancer.Gastroenterology. 2013;144:1252-1261.
Steer ML, Waxman I, Freedman S. Chronic pancreatitis.N Engl J Med. 1995;332:1482-1490.
Konner J, O’Reilly E. Pancreatic cancer: epidemiology, genetics, and approaches to screening.Oncology (Williston Park). 2002;16:1615-1622, 1631-1632, discussion 1632-1633, 1637-1638.
Gui J, Tian Y, Wen X, Zhang W, Zhang P, Gao J, Run W, Tian L, Jia X, Gao Y. Serum microRNA characterization identifies miR-885-5p as a potential marker for detecting liver pathologies.Clin Sci (Lond). 2011;120:183-193.
Afanasyeva EA, Mestdagh P, Kumps C, Vandesompele J, Ehemann V, Theissen J, Fischer M, Zapatka M, Brors B, Savelyeva L. MicroRNA miR-885-5p targets CDK2 and MCM5, activates p53 and inhibits proliferation and survival.Cell Death Differ. 2011;18:974-984.
Yan W, Zhang W, Sun L, Liu Y, You G, Wang Y, Kang C, You Y, Jiang T. Identification of MMP-9 specific microRNA expression profile as potential targets of anti-invasion therapy in glioblastoma multiforme.Brain Res. 2011;1411:108-115.
Guan X, Liu Z, Liu H, Yu H, Wang LE, Sturgis EM, Li G, Wei Q. A functional variant at the miR-885-5p binding site of CASP3 confers risk of both index and second primary malignancies in patients with head and neck cancer.FASEB J. 2013;27:1404-1412.
Schee K, Lorenz S, Worren MM, Günther CC, Holden M, Hovig E, Fodstad O, Meza-Zepeda LA, Flatmark K. Deep Sequencing the MicroRNA Transcriptome in Colorectal Cancer.PLoS One. 2013;8:e66165.
Vimalraj S, Miranda PJ, Ramyakrishna B, Selvamurugan N. Regulation of breast cancer and bone metastasis by microRNAs.Dis Markers. 2013;35:369-387.
Pandey DP, Picard D. miR-22 inhibits estrogen signaling by directly targeting the estrogen receptor alpha mRNA.Mol Cell Biol. 2009;29:3783-3790.
Patel JB, Appaiah HN, Burnett RM, Bhat-Nakshatri P, Wang G, Mehta R, Badve S, Thomson MJ, Hammond S, Steeg P. Control of EVI-1 oncogene expression in metastatic breast cancer cells through microRNA miR-22.Oncogene. 2011;30:1290-1301.
Ting Y, Medina DJ, Strair RK, Schaar DG. Differentiation-associated miR-22 represses Max expression and inhibits cell cycle progression.Biochem Biophys Res Commun. 2010;394:606-611.
Ling B, Wang GX, Long G, Qiu JH, Hu ZL. Tumor suppressor miR-22 suppresses lung cancer cell progression through post-transcriptional regulation of ErbB3.J Cancer Res Clin Oncol. 2012;138:1355-1361.
Polioudakis D, Bhinge AA, Killion PJ, Lee BK, Abell NS, Iyer VR. A Myc-microRNA network promotes exit from quiescence by suppressing the interferon response and cell-cycle arrest genes.Nucleic Acids Res. 2013;41:2239-2254.
Franchina T, Amodeo V, Bronte G, Savio G, Ricciardi GR, Picciotto M, Russo A, Giordano A, Adamo V. Circulating miR-22, miR-24 and miR-34a as novel predictive biomarkers to pemetrexed-based chemotherapy in advanced non-small cell lung cancer.J Cell Physiol. 2014;229:97-99.
Zhang C, Wang C, Chen X, Yang C, Li K, Wang J, Dai J, Hu Z, Zhou X, Chen L. Expression profile of microRNAs in serum: a fingerprint for esophageal squamous cell carcinoma.Clin Chem. 2010;56:1871-1879.
Liu N, Cui RX, Sun Y, Guo R, Mao YP, Tang LL, Jiang W, Liu X, Cheng YK, He QM. A four-miRNA signature identified from genome-wide serum miRNA profiling predicts survival in patients with nasopharyngeal carcinoma.Int J Cancer. 2013;Sep 2; Epub ahead of print.