Small nucleolar RNA U91 is a new internal control for accurate microRNAs quantification in pancreatic cancer
BMC Cancer volume 15, Article number: 774 (2015)
RT-qPCR quantification of miRNAs expression may play an essential role in pancreatic ductal adenocarcinoma (PDAC) diagnostics. RT-qPCR-based experiments require endogenous controls for the result normalization and reliability. However, expression instability of reference genes in tumors may introduce bias when determining miRNA levels.
We investigated expression of 6 miRNAs, isolated from FFPE samples of pancreatic adenocarcinomas. Four internal controls were utilized for RT-qPCR result normalization: artificial miR-39 from C. elegans, U6 snRNA, miR-16 and snoRNA U91.
We found miR-21, miR-155 or miR-217 expression values in tumors may differ up to several times, depending on selected internal controls. Moreover, different internal controls can produce controversial results for miR-96, miR-148a or miR-196a quantification. Also, expression of our endogenous controls varied significantly in tumors. U6 demonstrated variation from −1.03 to 8.12-fold, miR-16 from −2.94 up to 7.38-fold and the U91 from −3.05 to 4.36-fold respectively. On the other hand, the most stable gene, determined by NormFinder algorithm, was U91. Each miRNA normalized relatively to the spike or U91, demonstrated similar expression values. Thus, statistically significant and insignificant differences between tumors and normal tissues for miRNAs were equal for the spike and the U91. Also, the differences between the spike and U91 were statistically insignificant for all of miRs except miR-217. Among three endogenous controls, U91 had the lowest average expression values and standard deviation in cancer tissues.
We recommend U91 as a new normalizer for miRNA quantification in PDACs.
Pancreatic ductal adenocarcinoma (PDAC) is the most common and the most aggressive primary pancreatic neoplasm. The majority of patients are diagnosed by the time the tumor had already invaded peripancreatic structures or has metastasized . Therefore, there is a need for biomarkers enabling early detection of asymptomatic PDACs. miRNAs are stable in tissues and blood plasma ; consequently they are ideal molecules to be utilized as biomarkers. miRNAs are involved in oncogenesis, apoptosis and cell growth; thereby functioning as tumor suppressors or oncogenes [3–6]. A large number of miRNAs are proven to be overexpressed in pancreatic cancer [7–11]. On the other hand, the expression of miRNA-coding genes, which act as tumor suppressors, could be inhibited in cancer cells [12–16]. Alterations in the miRNAs expression profile of cancer in comparison with normal tissues could be used in pancreatic cancer diagnostics. The high sensitivity of reverse transcription quantitative PCR (RT-qPCR) has made it a popular method in the measurement of tumor miRNA expression. RT-qPCR-based experiments require endogenous controls for result normalization, reliability and reproducibility. The endogenous control helps to correct differences between sample quality and variations during RNA extraction or reverse transcription procedures. Housekeeping genes, ribosomal, small nuclear or nucleolar RNAs can play the role of such internal controls. However, according to experimental data, expression levels of these genes may differ in neoplastic and normal tissues [17–19]. These variations may introduce bias to experiment results.
In this study we compared the expression of selected miRNAs in samples from pancreatic cancer and normal pancreatic parenchyma and evaluated the influence of different internal controls on the expression data alterations.
Patients and tissue specimens
FFPE blocks of pancreatic ductal adenocarcinomas were retrieved from the archive of the Department of Pathology of the 3rd Faculty of Medicine of the Charles University and University Hospital Kralovske Vinohrady in Prague. The samples were collected from 24 patients, who had undergone pancreatoduodenectomy, distal pancreatectomy or total pancreatectomy between 2007 and 2012. Participants signed a written informed consent before the study. The study was performed according to the Declaration of Helsinki and approved by the Ethics Committee of the Third Faculty of Medicine (Charles University in Prague, Czech Republic). The resolution 1006/2012 was signed by Dr. Marek Vacha, Ph.D, Head of the Ethics Committee.
In the selected FFPE blocks the tumor occupied the majority of the slide. As negative control, FFPE blocks containing normal pancreatic tissue of the respective patients were selected.
The age of patients with resected pancreatic adenocarcinoma ranged from 36 to 83 years, with a median of 65.5 years. In total, 11 patients were women and 13 patients were men. Genetic syndromes were described in none of the patients. Grossly, 18 tumors were located in the head of the pancreas, 1 in the body of the pancreas and 5 in the tail of the pancreas.
The tumors showed in all of the selected cases the features of conventional ductal pancreatic adenocarcinoma. According to the guidelines of the WHO Classification of Tumors of the Gastrointestinal Tract, 3rd and 4th edition, 1 tumor was well differentiated, 14 tumors were moderately differentiated and 9 tumors were poorly differentiated. In one patient, a synchronous mucinous cystic neoplasm (MCN) was identified in the cauda of the pancreas. In another patient the tumor originated from an MCN. In 3 patients the resected tumor was described as pT1, in 5 patients pT2, in 15 patients pT3 and in one patient pT4. Additionally, lymph node metastases were confirmed in the resected specimens of 18 patients.
RNA isolation and reverse transcription
One to three 6 μm thick unstained paraffin embedded tissue sections were procured for RNA extraction, using the miRNeasy FFPE kit (Qiagen), following the manufacturer’s instructions. Two microliters of isolated RNA were used for RNA quantity and purity analysis. Optical density at 260 and 280 nm was measured with a multi-detection microplate reader Synergy HT (BioTek), including Take3 micro-volume plate. RNA integrity was assessed with denaturing agarose gel electrophoresis and GeneTools 3.08 software (SynGene).
A mix of 10 stem-loop primers was used for miRNA reverse transcription. Stem-loop primers were selected for the analysis, because their structure reduces annealing of the primer to pre- and pri-miRNAs, therefore increasing the specificity of the assay. Primers were designed with miRNA primer designer software, kindly provided by Dr. Fuliang Xie, East Carolina University. The stem-loop primer sequences for the internal controls, including the alien spike (miR-39 from C. elegans), and the examined pancreatic miRNAs are listed in Tables 1 and 2. The spike RNA was added to the reaction mix directly before the reverse transcription. Alien spike can’t be used as a normalizer for differences between samples during the RNA isolation, because tissue sections may contain different amounts of tissue. Therefore, the addition of spike before RNA isolation may introduce bias, because a ratio between amount of the RNA and the alien spike concentration may vary from sample to sample.
Reverse transcription was carried out, using RevertAid Reverse Transcriptase (Thermo Scientific), in a 50 μl reaction mixture, containing the following reagents: 1 μg of DNA-free RNA, reaction buffer [50 mM Tris–HCl (pH 8.3 at 25 °C), 50 mM KCl, 4 mM MgCl2 and 50 mM DTT]; 1 mM of dATP, dTTP, dCTP, dGTP; 20 IU rRNasin ribonuclease inhibitor; 100 IU of moloney murine leukemia virus reverse transcriptase (M-MuLV RT) and the primer mix, including 20 pmol of each stem-loop primer. Artificial spike RNA (miR-39 from C. elegans, 5 × 108 copies) was also added to the reaction as internal control. After initial denaturation (5 min at 70 °C, then cooling samples on ice), the reactions were incubated at 25 °C (10 min), and then at 42 °C for 1 h. To stop the reaction, the mixture was heated at 70 °C for 10 min.
cDNA samples were amplified in duplicates, using the Applied Biosystems 7500 Fast real-time PCR system and Hot FirePol EvaGreen qPCR Mix Plus (Solis BioDyne). The reaction mix included 10 pmol of each primer (miRNA specific and the universal (Table 3)) and 2 μl of cDNA.
Amplification of the cDNAs was performed at the following thermal conditions: denaturation at 94 °C for 15 min, followed by 40 cycles consisting of denaturation at 94 °C for 15 s, annealing at 48 °C for 60 s and DNA synthesis at 72 °C for 40 s. Reaction product specificity was controlled with their respective melting curves.
The expression of miRNAs in neoplastic and normal tissues was compared utilizing a paired two-tailed Student’s t test as well as a one-way ANOVA analysis. P-values below 0.05 were regarded as statistically significant. RT-qPCR data (threshold cycles) were linearized, and the NormFinder algorithm was used to calculate the most stable gene among the internal controls.
Evaluation of miRNA expression levels in PDAC samples
We investigated the expression of 6 miRNAs isolated from FFPE samples of pancreatic adenocarcinomas from 24 patients. The following microRNAs were selected: miR-21, which promotes cell proliferation and may accelerate tumorigenesis [8, 9, 20]; miR-155, which interacts with TP53 INP1 and transforming growth factor β (TGF-β) [11, 21, 22]; miR-96 and miR-217, which may act as a tumor suppressors, inhibiting the KRAS-signaling pathway [13, 14] also miR-148a and miR-196a, which are frequently included in experimental panels for pancreatic carcinoma diagnosis [23–29].
Four internal controls were utilized for qRT-PCR result normalization: an alien spike (artificial miR-39 from C. elegans) and three endogenous controls – U6 snRNA, miR-16 and snoRNA U91. miRNA expression values were normalized relative to each of these controls, and significant variations for the same miRNAs were found (Fig. 1, Table 4). In comparison with normal pancreatic tissue, miR-21 was significantly overexpressed, up to 14.56-fold (p < 0.01) in the case of the alien spike. However, for other internal controls, fold change values were shifted to 5.44 for U6 (p < 0.01), 7.03 for miR-16 (p < 0.01) and 17.71 for U91 (p < 0.01), respectively (Table 4, Fig. 1). The miR-155 also demonstrated increased expression levels with great variations between internal controls: 15.1-fold for the spike (p < 0.01); 5.05-fold for U6 (p < 0.01); 6.39-fold for miR-16 (p < 0.01) and 13.36-fold for U91 (p < 0.01). miRNA miR-96 in pancreatic carcinoma did not show significant differences in comparison with normal tissues, when normalizing to the alien spike (−1.04-fold, p > 0.05), as well as to U91 (−1.17-fold, p > 0.05). But, this miRNA was significantly down-regulated, when the expression was measured relative to U6 (−3.22-fold, p < 0.01) or miR-16 (−2.32-fold, p < 0.01). Also, no significant differences were found for miR-148a, normalized to spike (1.25 fold, p > 0.05) and U91 (1.06, p > 0.05). But, this miRNA was significantly inhibited for U6 (−1.33 fold, p < 0.01) and miR-16 (−2.04 fold, p < 0.01). Expression of miR-196a was slightly up-regulated relatively the alien spike (1.09-fold) and U91 (1.13-fold), without the results being statistically significant (p > 0.05). On the other hand, miR-196a was significantly down-regulated for U6 (−2.22-fold, p < 0.01), as well as not statistically significant for the miR-16 (−1.35, p > 0.05). The expression of miR-217 was significantly lower in all PDACs than in normal pancreatic tissues for all the examined internal controls (p < 0.01) (Table 4, Fig. 1).
NormFinder algorithm was used to calculate the most stable pairing of internal controls. According to our results, the best combination was U6 and U91 (stability value = 0036; Table 5) for the miRNAs expression evaluation. Normalized to this most stable pair, miRNAs miR-21 and miR-155 demonstrated significant upregulation (9.67-fold and 8.79-fold, p < 0.01). Activity of miR-96, miR-196a and miR-217 was significantly inhibited (−1.85-fold, −1.34-fold and −7.19-fold, p < 0.01). The miR-148a expression also was down-regulated, but the decrease was not statistically significant (−1.27-fold, p > 0.05).
Determination of the best normalizers for the miRNAs expression measuring
To find the best normalizer, we compared miRNA expression levels, normalized relatively the artificial spike and other endogenous controls, including the combination of U6 + U91, with one-way ANOVA analysis for each individual miRNA. The null hypothesis (Ho) was, that average fold change values, calculated for the each individual miRNA, are the same for all the internal controls. However, the differences were significant in the case of miR-21, miR-96, miR-148a, miR-155 and miR-196a (p < 0.01). For miR-217, the difference was not statistically significant (p > 0.05). Consequently, we compared miRNA expression normalized to the spike, with their normalization to the other individual examined endogenous controls, using a paired two-tailed Student’s t-test. Differences of miRNA expression, normalized to U6 or the combination of U6 and U91, in comparison with the alien spike, were statistically significant for all miRNAs (p < 0.01; Table 6). In the case of miR-16, the difference was not significant for miR-217 only (p > 0.05; Table 6). On the other hand, the difference between spike and the U91 was statistically insignificant for all miRNAs (p > 0.05; Table 6), except for miR-217 (p < 0.05; Table 6). Thus, one endogenous control was found which demonstrated a behavior very similar to the alien spike.
We investigated the causes of miRNA expression variations and their dependence on certain normalizers, and thus attempted to find the most suitable normalizer. The 2-ΔΔCT method was used for the miRNAs expression quantification, where CT is cycle threshold and ΔΔCT = ((CTmiRNA)tumor-(CTcontrol)tumor)-((CTmiRNA)normal-(CTcontrol)normal) For accurate miRNA quantification, in theory, CT values for the internal control gene should be very close for tumors and normal tissues, ideally (CTcontrol)tumor = (CTcontrol)normal. However, this CTcontrol value may be shifted up to several cycles in tumors (for example, up to ± n cycles), if the expression of endogenous control differs in tumor and normal tissue. This difference may introduce a bias to the miRNA fold change calculations:
While analyzing the amplification curves of the different internal controls, almost in all tumor samples a cycle threshold (CT) shift of n = 5 or even 6 cycles upwards in comparison with the normal tissue was apparent (Table 7). For example, CT values of the spike were very similar for PDAC and the normal tissues, they differed less than n = 0.8 cycle (Table 7). Nevertheless, for other normalizers these values varied from n = −6.20 up to n = 5.8 cycles (Table 7). We measured expression levels of our endogenous controls, using the alien spike for normalization. As expected, U6 expression in tumors varied from −1.03 to 8.12-fold, miR-16 showed variations from −2.94 up to 7.38-fold in different tumors, and the U91 from −3.05 to 4.36-fold respectively. The difference in expression was statistically significant for all endogenous controls (p < 0.01; Table 8). Also, U6 was overexpressed in 22 tumors from 24, miR-16 in 18 tumors and U91 in 14 correspondently. miR-16 was downregulated in 6 tumors and U91 in 5 tumor samples (Table 9). Thus, all selected endogenous controls demonstrated expression instability in tumor samples.
To identify the most stable internal control, the NormFinder algorithm was used. RT-qPCR expression data for all internal controls were linearized and compared in two groups: tumors and normal pancreatic tissues (Table 5). The most stable gene was U91 (stability value = 0.056), but stability values for all internal controls were close (0.085; 0.056 and 0.078 for the spike, U6 and miR-16, respectively; Table 5). U6 had the same stability value as U91 (0.056), but it demonstrated higher levels of intergroup variation (Table 5). Surprisingly, the most unstable control was the artificial spike (0.085; Table 5). The NormFinder can calculate variations between two groups, including all normal or cancer samples, but it is unable to evaluate the differences between normal and cancer tissues among individual patients. This may be the reason for the “instability” of the alien spike. The most stable pair of internal controls was a combinations of two genes, U6 and U91 (stability value = 0.036; Table 5).
MicroRNA expression values depend on a selected internal control
Pancreatic ductal adenocarcinoma is one of the most frequently occurring solid cancers and it carries an extremely poor prognosis. As such, an extensive search for biomarkers of early disease is undergoing, miRNAs may have the ideal characteristics to fulfil this role. Due to their stability and resistance against RNase degradation, they are viable in a wide range of samples. Viable miRNAs for PDAC diagnosis may be isolated from frozen and paraffin embedded tissue samples, stool, blood plasma, or pancreatic juice, [24, 28, 30, 31].
For our analysis we have selected miRNAs, frequently described to be dysregulated in various types of PDACs samples. Studies mapping microRNA expression using microarrays have proven considerable heterogeneity in pancreatic carcinomas. Zhang et al. have demonstrated relative expression values miRNAs spanning 6-logs (from 0.01–10,000) among individual cases . Among tumor samples we determined up to 45-fold variability in both miR-21 and miR-155 levels. During our brief review of literature we have noticed that the mean values of miRNA-levels in tumors varied among studies. There are many factors including differences in reagents/materials for miRNAs quantification protocols and data-processing algorithms, which can contribute to the variation. One of these factors is a variety of controls, which were used for normalization. Thus, the differences in the mean expression of miRNAs may be at least partially explained by the choice of controls for normalization.
For example, when normalizing with snoRNA U6, Bloomston et al. measured a median 2.2-fold increase in tumor miR-21 levels . Zhang et al., using the same internal control, found that expression of miR-21 was up-regulated up to 6888-fold in several tumors . Hong et al. reported about up to 550-fold increase in PDACs, normalizing relative to U6 . When using both U6 and 5S as endogenous controls, du Rieu et al. detected a 20.1-fold tumor miR-21 up-regulation . In our study, when normalizing with U6, a mean 5.5-fold increase in miR-21 in tumors was present. However, when normalizing to miR-16 a 7.03-fold increase was present (p < 0.01, Table 4). Wang et al. detected in plasma samples with miR-16 only a mean 2.42-fold up-regulation . On the other hand, when normalizing to the artificial spike or with U91 we detected a mean 14.56-fold and 17.71-fold increase (p < 0.01, Table 4).
The data about miR-96 expression in PDAC are controversial. Several groups of authors reported about miR-96 expression fold increase in experiments with microarray [24, 32]. For example, Bloomston et al. measured an average 1.77-fold increase, when determining miR-96 levels in PDACs . Kent et al. also demonstrated 2.7-fold upregulation of miR-96 in pancreatic cancer cell lines . On the other hand, miR-96 has been shown to be frequently down-regulated in experiments, utilizing Northern blot or RT-qPCR [13, 15, 25, 31, 33]. Szafranszka et al. determined in their study a −4.35-fold decrease in tumor miR-96 expression, when normalizing to miR-24 . Hong et al. as well as Feng et al. showed that miR-96, normalized to U6, was downregulated in PDAC samples up to −8-fold [31, 33]. With U6, miR-16 and combination of U6 + U91 respectively, we demonstrated a statistically significant mean −3.22-fold, −2.32-fold and −1.85-fold decrease in tumor tissue (p < 0.01, Table 4). However, expression analysis with the artificial spike and U91 alone yielded a statistically insignificant alteration in miR-96 expression in tumors in comparison with normal tissues (p > 0.05, Table 4).
miRNA miR-148a expression is described to be down-regulated in PDAC due to promoter hypermethylation, which represents an early event in pancreatic carcinogenesis . Bloomston et al. as well as Jamieson et al. measured an average −5.5-fold and −7.14-fold decrease respectively, when determining miR-148 expression with a microarray [24, 29]. In experiments, utilizing RT-qPCR, Szafranszka et al. demonstrated −6.15-fold decrease of miR-148a levels, with miR-24 as normalizer . However, Ma et al. and Zhang et al., normalizing to U6, determined a respective −2.86-fold and −2.5-fold downregulation in PDAC samples [34, 35]. Hanoun et al. also reported about miR-148a down-regulation, using U6 like the endogenous control . In our study, tumor miR-148a levels were −2.04-fold and −1.33-fold decreased with U6 and miR-16 as a normalizers, respectively (p < 0.01 and p < 0.05, Table 4). On the other hand, analysis of miR-148a expression, normalized to the alien spike, U91 and a combination of U6 and U91 did not determine statically significant differences in expression between cancerous and non-cancerous tissues (p > 0.05, Table 4).
The miR-155 is an onco-miR, overexpressed in early pancreatic adenocarcinoma precursors and invasive PDAC . The miR-155 expression in PDACs and pancreatic cancer cell lines, measured by microarray, ranged from 1.8 to 2.9-fold in different studies [24, 28, 29, 36, 37]. On the other hand, Habbe et al. found a mean 11.6-fold increase in intraductal papillary mucinous neoplasms, which was measured by RT-qPCR relative to U6 . Zhang et al. also used U6 like a normalizer in their study. They reported about up to 52-fold increase in individual cases . In our pancreatic carcinomas a mean 5.05-fold increase was present, when normalizing to U6 (p < 0.01, Table 4). However, Ma et al. measured only a 2.11-fold increase with the same endogenous control . Wang et al. determined a 3.74-fold increase in serum miR-155 levels in cancer, when normalizing with miR-16 . Our PDAC samples showed, on the other hand, a mean 6.39-fold increase with miR-16 as internal control (p < 0.01, Table 4). However, the expression levels were several times higher, measured relative to the alien spike or U91 - 15.1 and 13.36-fold respectively (p < 0.01, Table 4).
The miR-196a is an onco-miR reported to be frequently dysregulated in PDAC [27, 30]. Zhang et al. measured, when normalizing to U6, 0.35-1557-fold variations in tumor miR-196 expression . In our tumors we determined a mean −2.2-fold decrease, when normalizing to U6 (p < 0.01, Table 7). Wang et al. demonstrated 16.05-fold increase in plasma samples with miR-16 as the endogenous control . On the other hand, for miR-16, as well as for the alien spike or U91 we did not find significant differences in miR196a expression between cancer and normal tissues (p > 0.05, Table 4).
The miR-217 inhibits in vitro tumor cell growth and it functions as a potential tumor-suppressor by influencing the Akt/KRAS signaling pathway, therefore, miR-217 is frequently down-regulated in PDAC. MicroRNA miR-217 was down-regulated 10-fold in the study by Szafranszka et al., normalized relative to miR-24 . However, Greither et al. determined only a mean −2-fold decrease with 18S as internal control . Ma et al. demonstrated −3.91-fold decrease, using U6 for normalization . On the other hand, Hong et al. found, that expression of miR-217 was down-regulated up to −62.5-fold in PDACs. They also used U6 like internal control . In our samples, miR-217 expression was significantly down-regulated across all internal controls, with a maximum −24.39-fold decrease with U6 and a minimum −7.19-fold decrease with U6 + U91 combination (p < 0.01, Table 4).
Thus, for miRNAs with high positive or negative expression levels, such as miR-21, mir-155 or miR-217, fold change values may differ up to several times, depending on selected internal controls. Moreover, different internal controls can produce controversial results for miRNAs quantification, as it was demonstrated for miR-96, miR-148a or miR-196a.
Comparing internal controls: U91 is a new endogenous control for microRNAs quantification in pancreatic cancer
RT-qPCR quantification of tumor miRNA expression may play an essential role in PDAC diagnostics, chemotherapy resistance and survival prediction. RT-qPCR-based experiments require endogenous controls for the result normalization, reliability and reproducibility. U6 small nuclear RNA [8–11, 14, 15, 27, 30, 37, 38], 18S  and 5S ribosomal RNAs [8, 15, 39, 40], small nucleolar RNAs RNU48, RNU43, RNU44 – commercial available Applied Biosystems assays , or miR-16 [30, 42, 43] were often used as the endogenous controls for miRNAs expression evaluation. However, there are data indicating, that expression levels of these reference genes may differ significantly in neoplastic and normal tissues [17–19]. This expression instability may introduce bias, when determining miRNA dysregulation in tumors. For example, U6 small nuclear RNA was the most common internal control [8–11, 14, 15, 27, 30, 38] for the quantification of miRNAs expression in PDAC. However, there are data, implying that U6 expression may be unstable in breast and cervical cancers [17, 19, 42]. Also, the amount of U6 may vary significantly in serum samples from patients with breast and colorectal cancers [18, 42]. According to our findings, U6 expression may show as high as an 8-fold difference in PDAC and normal pancreatic tissue (Table 8). On the other hand, U6 was determined as the second most stable gene by the NormFinder algorithm (Table 5). U6 also demonstrated greater expression stability in breast carcinoma tissue samples when compared with the snoRNAs RNU44, RNU48 and RNU43. Furthermore, changes in levels of these snoRNAs correlated with tumor morphology and patient prognosis . However, U6 alongside 5S and miR-16 showed remarkable expression variability in tissue samples from patients with breast carcinoma .
The data about miR-16 expression in serum samples from the breast cancer patients are controversial. On the first look, this miRNA demonstrated significant expression variations [18, 42]. On the other hand, analysis with the geNorm algorithm has identified miR-16 as well as miR-425 as the most stable normalizers . According to our measurements, expression of miR-16 varied significantly in pancreatic carcinomas (p < 0.01, Table 7). In addition, miR-16 was marked by the NormFinder algorithm as the least stable of the analyzed endogenous controls (Table 5).
Another possibility for RT-qPCR result normalization is the use of alien spike RNAs, such as miR-39 from C. elegans [18, 44], as internal controls. Also, these spike RNAs should be selected while taking into consideration that the same RNA sequences may already exist in the human genome. Surprisingly, according to the NormFinder analysis, the artificial spike was the least stable control (stability 0.085; Table 5). It must be taken into consideration, that the NormFinder algorithm can calculate variations between two groups including all normal and cancer samples, but it is unable to evaluate the differences between normal and cancer tissues from the individual patients. Accordingly, this may be the reason for the “instability” of the alien spike.
In this study, we compared the expression of 4 internal controls to determine the most stable of them. On the first look, the best internal control is the artificial spike, due to its amplification curves and threshold cycles, which have demonstrated to be very close for cancers and normal tissues (Table 7). On the other hand, according to the results, yielded by the NormFinder analysis, the best normalizer is the combination of U6 and U91. This combination has the best stability value, but as normalization results show, it differs significantly from the artificial spike (p < 0.01, Table 6). The most stable gene, determined by NormFinder, was U91 (Table 5). Each miRNA normalized relatively to the spike or U91, demonstrated similar expression values. Thus, statistically significant and insignificant differences between tumors and normal tissues for miRNAs were equal for the spike and the U91 (Table 4). Also, the differences between the spike and U91 were statistically insignificant for all of miRs except of miR-217 (Table 6). Among three endogenous controls, the U91 had the lowest average expression values and standard deviation in cancer tissues (Table 8).
Thus, we recommend U91 as a new normalizer of miRNA expression in pancreatic adenocarcinoma.
We found expression of traditional endogenous controls, such as U6 and miR-16 can be unstable in pancreatic tumors and may vary up to several times. This instability may introduce bias to the miRNAs quantification. On the other hand, U91, the new stable internal control for miRNAs expression evaluation in pancreatic cancers was found.
This study was carried out in compliance to the Minimum Information for Publication of Quantitative Real–Time PCR Experiments (MIQE; ).
Availability of data
Data files, including raw CT values or fold change tables are available on request, please, contact the correspondence author.
Pancreatic ductal adenocarcinoma
Reverse transcription quantitative polymerase chain reaction
Formalin-fixed, paraffin-embedded (tissues)
Hidalgo M. Pancreatic cancer. The New England journal of medicine. 2010;362(17):1605–17. doi:10.1056/NEJMra0901557.
Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. P Natl Acad Sci USA. 2008;105(30):10513–8. doi:10.1073/pnas.0804549105.
Sassen S, Miska EA, Caldas C. MicroRNA - implications for cancer. Virchows Arch. 2008;452(1):1–10. doi:10.1007/s00428-007-0532-2.
Shenouda SK, Alahari SK. MicroRNA function in cancer: oncogene or a tumor suppressor? Cancer Metast Rev. 2009;28(3–4):369–78. doi:10.1007/s10555-009-9188-5.
Zhang BH, Pan XP, Cobb GP, Anderson TA. microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007;302(1):1–12. doi:10.1016/j.ydbio.2006.08.028.
Lynam-Lennon N, Maher SG, Reynolds JV. The roles of microRNA in cancer and apoptosis. Biol Rev. 2009;84(1):55–71. doi:10.1111/j.1469-185X.2008.00061.x.
Lee EJ, Gusev Y, Jiang JM, Nuovo GJ, Lerner MR, Frankel WL, et al. Expression profiling identifies microRNA signature in pancreatic cancer. Int J Cancer. 2007;120(5):1046–54. doi:10.1002/Ijc.22394.
du Rieu MC, Torrisani J, Selves J, Al Saati T, Souque A, Dufresne M, et al. MicroRNA-21 Is Induced Early in Pancreatic Ductal Adenocarcinoma Precursor Lesions. Clin Chem. 2010;56(4):603–12. doi:10.1373/clinchem.2009.137364.
Moriyama T, Ohuchida K, Mizumoto K, Yu J, Sato N, Nabae T, et al. MicroRNA-21 modulates biological functions of pancreatic cancer cells including their proliferation, invasion, and chemoresistance. Mol Cancer Ther. 2009;8(5):1067–74. doi:10.1158/1535-7163.Mct-08-0592.
Ryu JK, Hong SM, Karikari CA, Hruban RH, Goggins MG, Maitra A. Aberrant MicroRNA-155 Expression Is an Early Event in the Multistep Progression of Pancreatic Adenocarcinoma. Pancreatology. 2010;10(1):66–73. doi:10.1159/000231984.
Habbe N, Koorstra JBM, Mendell JT, Offerhaus GJ, Ryu JK, Feldmann G, et al. MicroRNA miR-155 is a biomarker of early pancreatic neoplasia. Cancer Biol Ther. 2009;8(4):340–6. doi:10.4161/Cbt.8.4.7338.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, et al. RAS is regulated by the let-7 MicroRNA family. Cell. 2005;120(5):635–47. doi:10.1016/j.cell.2005.01.014.
Yu SN, Lu ZH, Liu CZ, Meng YX, Ma YH, Zhao WG, et al. miRNA-96 Suppresses KRAS and Functions as a Tumor Suppressor Gene in Pancreatic Cancer. Cancer Res. 2010;70(14):6015–25. doi:10.1158/0008-5472.Can-09-4531.
Zhao WG, Yu SN, Lu ZH, Ma YH, Gu YM, Chen J. The miR-217 microRNA functions as a potential tumor suppressor in pancreatic ductal adenocarcinoma by targeting KRAS. Carcinogenesis. 2010;31(10):1726–33. doi:10.1093/carcin/bgq160.
Hanoun N, Delpu Y, Suriawinata AA, Bournet B, Bureau C, Selves J, et al. The Silencing of MicroRNA 148a Production by DNA Hypermethylation Is an Early Event in Pancreatic Carcinogenesis. Clin Chem. 2010;56(7):1107–18. doi:10.1373/clinchem.2010.144709.
Mardin W, Mees ST. MicroRNAs: Novel Diagnostic and Therapeutic Tools for Pancreatic Ductal Adenocarcinoma? Ann Surg Oncol. 2009;16(11):3183–9. doi:10.1245/s10434-009-0623-1.
Peltier HJ, Latham GJ. Normalization of microRNA expression levels in quantitative RT-PCR assays: Identification of suitable reference RNA targets in normal and cancerous human solid tissues. Rna. 2008;14(5):844–52. doi:10.1261/Rna.939908.
Xiang MQ, Zeng Y, Yang RR, Xu HF, Chen Z, Zhong J, et al. U6 is not a suitable endogenous control for the quantification of circulating microRNAs. Biochem Bioph Res Co. 2014;454(1):210–4. doi:10.1016/j.bbrc.2014.10.064.
Hansen CN, Ketabi Z, Rosenstierne MW, Palle C, Boesen HC, Norrild B. Expression of CPEB, GAPDH and U6snRNA in cervical and ovarian tissue during cancer development. Apmis. 2009;117(1):53–9. doi:10.1111/j.1600-0463.2008.00015.x.
Qi LQ, Bart J, Tan LP, Platteel I, van der Sluis T, Huitema S et al. Expression of miR-21 and its targets (PTEN, PDCD4, TM1) in flat epithelial atypia of the breast in relation to ductal carcinoma in situ and invasive carcinoma. Bmc Cancer. 2009;9. doi: 10.1186/1471-2407-9-163.
Gironella M, Seux M, Xie MJ, Cano C, Tomasini R, Gommeaux J, et al. Tumor protein 53-induced nuclear protein 1 expression is repressed by miR-155, and its restoration inhibits pancreatic tumor development. Proc Natl Acad Sci U S A. 2007;104(41):16170–5. doi:10.1073/pnas.0703942104.
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(1):73–80. doi:10.1002/ijc.24687.
Szafranska AE, Davison TS, John J, Cannon T, Sipos B, Maghnouj A, et al. MicroRNA expression alterations are linked to tumorigenesis and non-neoplastic processes in pancreatic ductal adenocarcinoma. Oncogene. 2007;26(30):4442–52. doi:10.1038/sj.onc.1210228.
Bloomston M, Frankel WL, Petrocca F, Volinia S, Alder H, Hagan JP, et al. MicroRNA expression patterns to differentiate pancreatic adenocarcinoma from normal pancreas and chronic pancreatitis. Jama. 2007;297(17):1901–8. doi:10.1001/jama.297.17.1901.
Szafranska AE, Doleshal M, Edmunds HS, Gordon S, Luttges J, Munding JB, et al. Analysis of microRNAs in pancreatic fine-needle aspirates can classify benign and malignant tissues. Clin Chem. 2008;54(10):1716–24. doi:10.1373/clinchem.2008.109603.
Huang F, Tang J, Zhuang X, Zhuang Y, Cheng W, Chen W, et al. MiR-196a promotes pancreatic cancer progression by targeting nuclear factor kappa-B-inhibitor alpha. PloS one. 2014;9(2):e87897. doi:10.1371/journal.pone.0087897.
Zhang Y, Li M, Wang H, Fisher WE, Lin PH, Yao Q, et al. Profiling of 95 microRNAs in pancreatic cancer cell lines and surgical specimens by real-time PCR analysis. World journal of surgery. 2009;33(4):698–709. doi:10.1007/s00268-008-9833-0.
Link A, Becker V, Goel A, Wex T, Malfertheiner P. Feasibility of Fecal MicroRNAs as Novel Biomarkers for Pancreatic Cancer. PloS one. 2012;7(8):e42933. doi:10.1371/journal.pone.0042933.
Jamieson NB, Morran DC, Morton JP, Ali A, Dickson EJ, Carter CR, et al. MicroRNA molecular profiles associated with diagnosis, clinicopathologic criteria, and overall survival in patients with resectable pancreatic ductal adenocarcinoma. Clinical cancer research : an official journal of the American Association for Cancer Research. 2012;18(2):534–45. doi:10.1158/1078-0432.CCR-11-0679.
Wang J, Chen J, Chang P, LeBlanc A, Li D, Abbruzzesse JL, et al. MicroRNAs in plasma of pancreatic ductal adenocarcinoma patients as novel blood-based biomarkers of disease. Cancer prevention research. 2009;2(9):807–13. doi:10.1158/1940-6207.CAPR-09-0094.
Hong TH, Park IY. MicroRNA expression profiling of diagnostic needle aspirates from surgical pancreatic cancer specimens. Annals of surgical treatment and research. 2014;87(6):290–7. doi:10.4174/astr.2014.87.6.290.
Kent OA, Mullendore M, Wentzel EA, Lopez-Romero P, Tan AC, Alvarez H, et al. A resource for analysis of microRNA expression and function in pancreatic ductal adenocarcinoma cells. Cancer Biol Ther. 2009;8(21):2013–24.
Feng J, Yu JB, Pan XL, Li ZL, Chen Z, Zhang WJ, et al. HERG1 functions as an oncogene in pancreatic cancer and is downregulated by miR-96. Oncotarget. 2014;5(14):5832–44.
Ma MZ, Kong X, Weng MZ, Cheng K, Gong W, Quan ZW, et al. Candidate microRNA biomarkers of pancreatic ductal adenocarcinoma: meta-analysis, experimental validation and clinical significance. J Exp Clin Canc Res. 2013;32:71. doi:10.1186/1756-9966-32-71.
Zhang R, Li M, Zang W, Chen X, Wang Y, Li P, et al. MiR-148a regulates the growth and apoptosis in pancreatic cancer by targeting CCKBR and Bcl-2. Tumour biology : the journal of the International Society for Oncodevelopmental Biology and Medicine. 2014;35(1):837–44. doi:10.1007/s13277-013-1115-2.
Schultz NA, Werner J, Willenbrock H, Roslind A, Giese N, Horn T, et al. MicroRNA expression profiles associated with pancreatic adenocarcinoma and ampullary adenocarcinoma. Modern Pathol. 2012;25(12):1609–22. doi:10.1038/modpathol.2012.122.
Wang J, Raimondo M, Guha S, Chen JY, Diao LX, Dong XQ, et al. Circulating microRNAs in Pancreatic Juice as Candidate Biomarkers of Pancreatic Cancer. J Cancer. 2014;5(8):696–705. doi:10.7150/Jca.10094.
Nakata K, Ohuchida K, Mizumoto K, Kayashima T, Ikenaga N, Sakai H, et al. MicroRNA-10b is overexpressed in pancreatic cancer, promotes its invasiveness, and correlates with a poor prognosis. Surgery. 2011;150(5):916–22. doi:10.1016/j.surg.2011.06.017.
Sun M, Estrov Z, Ji Y, Coombes KR, Harris DH, Kurzrock R. Curcumin (diferuloylmethane) alters the expression profiles of microRNAs in human pancreatic cancer cells. Mol Cancer Ther. 2008;7(3):464–73. doi:10.1158/1535-7163.MCT-07-2272.
Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, et al. Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004;64(11):3753–6. doi:10.1158/0008-5472.CAN-04-0637.
Gee HE, Buffa FM, Camps C, Ramachandran A, Leek R, Taylor M, et al. The small-nucleolar RNAs commonly used for microRNA normalisation correlate with tumour pathology and prognosis. Brit J Cancer. 2011;104(7):1168–77. doi:10.1038/sj.bjc.6606076.
Appaiah HN, Goswami CP, Mina LA, Badve S, Sledge GW, Liu YL, et al. Persistent upregulation of U6:SNORD44 small RNA ratio in the serum of breast cancer patients. Breast Cancer Res. 2011;13(5):R86. doi:10.1186/Bcr2943.
McDermott AM, Kerin MJ, Miller N. Identification and Validation of miRNAs as Endogenous Controls for RQ-PCR in Blood Specimens for Breast Cancer Studies. PloS one. 2013;8(12):e83718. doi:10.1371/journal.pone.0083718.
Liu JQ, Gao J, Du YQ, Li ZS, Ren Y, Gu JJ, et al. Combination of plasma microRNAs with serum CA19-9 for early detection of pancreatic cancer. Int J Cancer. 2012;131(3):683–91. doi:10.1002/Ijc.26422.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22. doi:10.1373/clinchem.2008.112797.
This work was technically supported by the project OPPK No. CZ.2.16/3.1.00/24024, awarded by European Fund for Regional Development (Prague & EU – We invest for your future) and by the project PRVOUK Oncology P27, awarded by Charles University in Prague.
The authors declare that they have no competing interests.
AP carried out the molecular genetic studies, performed the statistical analysis and drafted the manuscript. AS provided technical and material support, helped in acquisition of data and drafted the manuscript. VM conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Popov, A., Szabo, A. & Mandys, V. Small nucleolar RNA U91 is a new internal control for accurate microRNAs quantification in pancreatic cancer. BMC Cancer 15, 774 (2015). https://doi.org/10.1186/s12885-015-1785-9