- Research article
- Open Access
- Open Peer Review
Differential long-term stability of microRNAs and RNU6B snRNA in 12–20 year old archived formalin-fixed paraffin-embedded specimens
BMC Cancervolume 17, Article number: 32 (2017)
The quantitative analysis of microRNA (miRNA) gene expression in archived formalin-fixed, paraffin embedded (FFPE) tissues has been instrumental to identifying their potential roles in cancer biology, diagnosis, and prognosis. However, it remains unclear whether miRNAs remain stable in FFPE tissues stored for long periods of time.
Here we report Taqman real-time RT-PCR quantification of miR-21, miR-141, miR-221, and RNU6B small nuclear RNA (snRNA) levels from 92 radical prostatectomy specimens stored for 12–20 years in FFPE blocks. The relative stability of each transcript over time was assessed using general linear models. The correlation between transcript quantities, sample age, and RNA integrity number (RIN) were determined utilizing Spearman rank correlation.
All transcript levels linearly decreased with sample age, demonstrating a clear loss of miRNA stability and RNU6B snRNA stability over time. The most rapid rates of degradation were observed for RNU6B and miR-21, while miR-141 and miR-221 were more stable. RNA quality was not correlated with sample age or with miR-21, miR-221, or RNU6B snRNA levels. Conversely, miR-141 levels increased with RNA quality.
MiRNA and snRNA levels gradually decreased over an eight year period in FFPE tissue blocks. Sample age was the most consistent feature associated with miRNA stability. The reference snRNA, RUN6B, was more rapidly degraded when compared to miR-141 and miR-221 miRNAs. Various miRNAs demonstrated differential rates of degradation. Quantitative miRNA studies from long-term archived FFPE tissues may therefore benefit from epidemiologic study design or statistical analysis methods that take into account differential storage-dependent transcript degradation.
MicroRNAs are small non-coding transcripts that post-transcriptionally inhibit the translation and stability of specific mRNAs . Deregulated miRNA expression has been reported in most every human cancer . Moreover, numerous miRNAs are aberrantly expressed during disease development or progression, making miRNAs attractive as diagnostic or prognostic biomarkers . Many of the miRNAs commonly deregulated in cancer have tumor suppressive or oncogenic properties when inhibited or over-expressed in cellular and animal models . This collective study of miRNA gene expression in clinical tissues and miRNA function in cancer models has been crucial to our current understanding of the roles miRNAs may play in human cancer.
Archival formalin-fixed, paraffin-embedded (FFPE) tissues have been an invaluable resource for gene expression studies in cancer biology. These samples provide a historical record of tissue histology, protein and nucleic acid expression that can be correlated with long-term follow up of disease outcomes. Properly fixed and embedded samples can be stored indefinitely at room temperature without loss of structural integrity. However, processing and long-term storage can result in the damage and degradation of proteins and RNA transcripts. The stability of coding mRNA in FFPE specimens is dependent upon many factors, including processing methods, storage time and transcript abundance [5, 6]. RNA quality measurements, such as RNA integrity number (RIN) score, are currently applied to identify more reliable specimens for mRNA gene expression analyses . It remains unclear, however, how well analyses of transcript size represent the quality or stability of small non-coding RNAs, such as miRNAs.
Comparative analyses of FFPE and rapidly frozen tissues suggests that miRNAs are unusually stable when compared to longer mRNA transcripts [8, 9]. In point of fact, the intentional degradation of RNA samples by exposure to high temperatures reduces RIN and mRNA stability without apparent influence on miRNA levels . However, the literature has reached variable conclusions regarding miRNA stability in FFPE specimens stored for long periods of time [9–17]. There is therefore a need to further study miRNA stability in aged FFPE samples. Here we present new data on miRNA and RNU6B snRNA stability from 92 FFPE radical prostatectomy specimens that were processed at a single institution and stored for 12–20 years. The results indicate steady degradation of miRNAs over time and indicate that different miRNAs have different stabilities. The commonly applied reference snRNA transcript, RNU6B, showed differential stability from some miRNAs over time. These results signify a need to consider sample age and potential differences between miRNA and other snRNA stabilities in miRNA expression analyses from aged FFPE tissue samples.
Archival samples and RNA extractions
Total RNA from FFPE radical prostatectomy specimens was previously obtained from a cohort of 92 patients who underwent surgery at Johns Hopkins between 1993 and 2001, and were followed for recurrence until 2004 . Briefly, malignant regions were identified by a pathologist and two 0.6 mm tumor cores were obtained from each sample. Samples were deparaffinized with xylene, extracted with ethanol, and protease digested. Total RNA was extracted with the RecoverAll Total Nucleic Acid Isolation Kit (Ambion, Grand Island, NY) and stored at −80 °C. RNA quantity was assessed by Nanodrop spectroscopy (Thermo Scientific, Wilmington, DE) and RNA integrity score were obtained with Agilent RNA Pico kit and the 2100 Bioanalyzer Instrument (Agilent Technologies, Wilmington, DE) as per manufacturer instructions.
Reverse transcription and quantitative PCR were previously performed using the Taqman MicroRNA Reverse Transcription Kit and Applied Biosystems qRT-PCR probe and primer sets (Applied Biosystems, Grand Island, NY) . Briefly, 2 ng of RNA from FFPE samples was applied to RT and QPCR was performed with the Bio-Rad MyiQ single color RT-PCR detection system (Bio-Rad Laboratories, Hercules, CA). Standard curves (2 ng/μl to 0.016 ng/μl) for RNU6B and for each miRNA transcript were generated from cell lines known to express each miRNA (PC3 for miR-221; LNCaP for miR-141; and MCF-7 for miR-21). The resulting standard curves and results were used to determine assay linearity, qRT-PCR efficiency, and SQmean.
General linear models were utilized to compare RNA transcript stability rates in models with and without interaction. If the interaction term was deemed insignificant (parallel slopes) comparisons were further assessed utilizing analysis of covariance (ANCOVA). Statistical comparisons with RNU6B were applied using levels determined from each cell line standard curve (PC3 for miR-221; LNCaP for miR-141; and MCF-7 for miR-21). The correlation between transcript quantity and RIN or FFPE block age (sample age) were determined by Spearman rank correlation, and RNU6B levels derived from all cell line standard curves reached the same conclusions. For ease of reporting only RNU6B derived from PC3 standard curves are reported for the correlation analysis. Two-sided tests were conducted using R (version 3.2.3), and in all analyses, p-values <0.05 were considered to be statistically significant.
Loss of miRNA and snRNA stability in archival FFPE specimens
Over the last several decades the Johns Hopkins retropubic radical prostatectomy program has treated thousands of men with localized prostate cancer [19–21]. Archival specimens from these patients provide a unique resource to study the effects of long term FFPE tissue block storage on miRNA stability. As a single institution study of a single type of cancer, these samples underwent standardized processing, fixation, and storage procedures. Therefore, they present a unique and well-controlled resource for miRNA analyses. Tissue blocks from a previously defined nested-case control cohort of 92 radical prostatectomy patients, treated at Johns Hopkins Hospital between 1993 and 2001, were identified and total RNA was extracted from pathologist-defined tumor tissue . The characteristics of patient and tumor samples are provided in Table 1. Transcript levels of miR-21, miR-141, miR-221, and the commonly applied reference gene, RNU6B snRNA, were quantified by Taqman quantitative RT-PCR assay. The results demonstrate a clear and statistically significant inverse association between transcript quantity and FFPE block age for all four small non-coding RNA transcripts by Spearman rank correlation analysis (Fig. 1, Additional file 1: Table S1). These results indicate a continuous loss of miRNA and RNU6B snRNA transcript stability over storage time.
Small non-coding RNA stability and RNA Quality
The quality of RNA from older FFPE samples is generally considered to be poor. To investigate whether RNA quality is indicative of miRNA and snRNA transcript levels, we calculated the RIN for each sample. RIN scores ranged from 1 to 3.9 (Average 2.1 ± 0.5). No association was found between sample RIN and FFPE block age (Additional file 1: Figure S1). Moreover, no statistically significant correlation was found between RIN and RNU6B snRNA, miR-21, or miR-221 transcript levels by Spearman rank correlation (Fig. 2a, b and d). However, a significant correlation was observed between RIN and miR-141 transcript levels (Fig. 2c), indicating that RNA quality can be associated with the stability and detectability of some miRNAs. Collectively, these results identify FFPE block age as the most consistent feature associated with miRNA and snRNA transcript detectability and stability. The results of these statistical analyses are summarized in Additional file 1: Table S2.
Differential stability between small non-coding RNA transcripts
Mature miRNAs are approximately 22 nucleotides in length while common endogenous reference controls, such as RNU6B snRNA, are slightly longer (Table 2). Given the differences in length and function between miRNAs and reference snRNAs, we investigated whether there may be corresponding differences in signal stability over time by general linear models. RNU6B and miR-21 were found to have comparable stabilities over time (Fig. 3a). However, the rate of RNU6B snRNA degradation was significantly more rapid than that of miR-141 and miR-221 (Fig. 3b–c). Accordingly, miR-21 was also found to have a more rapid degradation rate when compared to miR-141 and miR-221 (Fig. 3d). These data support that commonly studied microRNAs and reference controls can have significantly different stabilities in FFPE tissues stored over long periods of time. These results may reflect differences in small non-coding RNA sequence, function, or association with proteins. They further suggested that quantitative miRNA studies from long-term archived FFPE tissues may benefit from epidemiologic study design or statistical analysis methods that take into account differential storage-dependent transcript degradation.
The first reported miRNA gene, lin-4, was described in C. elegans in 1993 ; the same year that the first surgical samples in this study were isolated. In a perfect world these samples would remain unchanged over time; however, it is well known that processing and storage can lead to significant RNA degradation [5, 6]. The influence of such parameters on miRNA detection and stability remains debatable.
A recent study of six miRNAs in colorectal tissue blocks stored for up to 28 years found no significant effects of sample block age on miRNA detection . Similar observations were reported for miR-181b and 5S ribosomal RNA from blocks as old as ten years . Deep sequencing analyses of miRNAs from multiple different types of tissue, stored for 2 to 9 years, have also found no significant change in miRNA detection with sample age . In contrast, others have reported significant miRNA loss with extended FFPE block storage times. Comparison of miRNA expression from 1 to 11 year old FFPE samples reported a clear and gradual loss of miRNA signal with storage time . Gradual loss of miRNA signal was also reported from FFPE samples of tongue carcinoma , with detectable signal loss after only one year of storage. In yet another study, miRNA levels from human tissues processed and stored as FFPE blocks for more than 10–20 years showed a nearly 50% decrease in miRNA accessibility . These few examples of conflicting results underscore the need to better understand if and how FFPE block storage time may influence miRNA detection.
Here we report the effect of FFPE block storage time on the stability of three frequently studied miRNAs, and one of the most commonly utilized miRNA reference genes, in prostate cancer specimens stored for up to 20 years. Importantly, this study was performed on samples of a single tissue type that were isolated and processed by a single institution. Therefore tissue type, processing methodology and storage have been consistent. Our results support that miRNAs and snRNAs are not stable over time in FFPE samples stored for over 12 years (Fig. 1).
There are several factors that can contribute to the loss of RNA stability in FFPE blocks including fixation time, fixation method, or exposure to oxidation, extreme temperatures, or light [10, 13, 17]. Capillary electrophoresis analyses and RIN score are standard methods for evaluating RNA quality from such clinical specimens. These methods primarily focus on nucleotide fragment size distribution; therefore, they may not distinguish small RNA transcripts from mRNA degradation products. In our study, RIN score was not associated with small non-coding RNA stability in three of the four transcripts studied (Fig. 2). Moreover, two identically sized miRNAs, miR-21 and miR-141, had significantly different stabilities in this sample set (Fig. 3). On the other hand, miR-141 levels were found to be associated with RIN score. These results reflect a common need to identify better tools for assessing miRNA quality in clinical specimens. It is notable that our sample set had a limited range of RIN scores and did not include any samples younger than 12 years of age. Consequently, a more clear correlation between miRNA stability, sample age, and RIN score may have been observed with a broader range of sample ages and qualities.
One interesting observation from this study is the apparent differential stability between some miRNAs. MiR-221 and miR-141 were significantly more stable over time when compared to miR-21 and RNU6B (Fig. 3). This supports previous reports of differential stability between different miRNAs . Perhaps differences in miRNA function, localization, or extent and stability of binding to other macromolecules help account for these observations. Further studies are needed to characterize these features and to identify mechanisms that contribute to or inhibit miRNA degradation over time.
We report a linear loss of miRNA and RNU6B snRNA signal with sample age in FFPE blocks stored for 12–20 years. Some miRNAs were more stable than others. FFPE block age was the most consistent feature associated with miRNA and snRNA stability. It is important to emphasize the poor stability of RNU6B snRNA, when compared to some miRNAs, because it is one of the most widely utilized reference genes for miRNA expression normalization. Our results suggest that it would be beneficial to consider sample block age, rather than RNA quality, in miRNA expression analyses from older FFPE samples. For example, we have previously noted trends in miRNA and RNU6 snRNA stability and have adjusted analyses by calendar year of prostatectomy . Similar adjustments have been made in other studies . Future research is needed to identify alternative methods for assessing miRNA and snRNA quality in older FFPE samples, and for miRNA expression comparisons and normalization.
Analysis of covariance
Formalin-fixed, paraffin embedded
RNA integrity number
Hammond SM. An overview of microRNAs. Adv Drug Deliv Rev. 2015;87:3–14.
Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6(11):857–66.
Di Leva G, Croce CM. miRNA profiling of cancer. Curr Opin Genet Dev. 2013;23(1):3–11.
Kent OA, Mendell JT. A small piece in the cancer puzzle: microRNAs as tumor suppressors and oncogenes. Oncogene. 2006;25(46):6188–96.
Turashvili G, Yang W, McKinney S, Kalloger S, Gale N, Ng Y, Chow K, Bell L, Lorette J, Carrier M, et al. Nucleic acid quantity and quality from paraffin blocks: defining optimal fixation, processing and DNA/RNA extraction techniques. Exp Mol Pathol. 2012;92(1):33–43.
Cronin M, Pho M, Dutta D, Stephans JC, Shak S, Kiefer MC, Esteban JM, Baker JB. Measurement of gene expression in archival paraffin-embedded tissues: development and performance of a 92-gene reverse transcriptase-polymerase chain reaction assay. Am J Pathol. 2004;164(1):35–42.
Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7:3.
Xi Y, Nakajima G, Gavin E, Morris CG, Kudo K, Hayashi K, Ju J. Systematic analysis of microRNA expression of RNA extracted from fresh frozen and formalin-fixed paraffin-embedded samples. RNA. 2007;13(10):1668–74.
Meng W, McElroy JP, Volinia S, Palatini J, Warner S, Ayers LW, Palanichamy K, Chakravarti A, Lautenschlaeger T. Comparison of microRNA deep sequencing of matched formalin-fixed paraffin-embedded and fresh frozen cancer tissues. PLoS One. 2013;8(5):e64393.
Jung M, Schaefer A, Steiner I, Kempkensteffen C, Stephan C, Erbersdobler A, Jung K. Robust microRNA stability in degraded RNA preparations from human tissue and cell samples. Clin Chem. 2010;56(6):998–1006.
Hall JS, Taylor J, Valentine HR, Irlam JJ, Eustace A, Hoskin PJ, Miller CJ, West CM. Enhanced stability of microRNA expression facilitates classification of FFPE tumour samples exhibiting near total mRNA degradation. Br J Cancer. 2012;107(4):684–94.
Bovell L, Shanmugam C, Katkoori VR, Zhang B, Vogtmann E, Grizzle WE, Manne U. miRNAs are stable in colorectal cancer archival tissue blocks. Front Biosci. 2012;4:1937–40.
Kakimoto Y, Kamiguchi H, Ochiai E, Satoh F, Osawa M. MicroRNA Stability in Postmortem FFPE Tissues: Quantitative Analysis Using Autoptic Samples from Acute Myocardial Infarction Patients. PLoS One. 2015;10(6):e0129338.
Rentoft M, Fahlen J, Coates PJ, Laurell G, Sjostrom B, Ryden P, Nylander K. miRNA analysis of formalin-fixed squamous cell carcinomas of the tongue is affected by age of the samples. Int J Oncol. 2011;38(1):61–9.
Siebolts U, Varnholt H, Drebber U, Dienes HP, Wickenhauser C, Odenthal M. Tissues from routine pathology archives are suitable for microRNA analyses by quantitative PCR. J Clin Pathol. 2009;62(1):84–8.
Szafranska AE, Davison TS, Shingara J, Doleshal M, Riggenbach JA, Morrison CD, Jewell S, Labourier E. Accurate molecular characterization of formalin-fixed, paraffin-embedded tissues by microRNA expression profiling. J Mol Diagn. 2008;10(5):415–23.
Becker C, Hammerle-Fickinger A, Riedmaier I, Pfaffl MW. mRNA and microRNA quality control for RT-qPCR analysis. Methods. 2010;50(4):237–43.
Zheng Q, Peskoe SB, Ribas J, Rafiqi F, Kudrolli T, Meeker AK, De Marzo AM, Platz EA, Lupold SE. Investigation of miR-21, miR-141, and miR-221 expression levels in prostate adenocarcinoma for associated risk of recurrence after radical prostatectomy. Prostate. 2014;74(16):1655–62.
Partin AW, Pound CR, Clemens JQ, Epstein JI, Walsh PC. Serum PSA after anatomic radical prostatectomy. The Johns Hopkins experience after 10 years. Urol Clin North Am. 1993;20(4):713–25.
Pound CR, Partin AW, Epstein JI, Walsh PC. Prostate-specific antigen after anatomic radical retropubic prostatectomy. Patterns of recurrence and cancer control. Urol Clin North Am. 1997;24(2):395–406.
Han M, Partin AW, Pound CR, Epstein JI, Walsh PC. Long-term biochemical disease-free and cancer-specific survival following anatomic radical retropubic prostatectomy. The 15-year Johns Hopkins experience. Urol Clin North Am. 2001;28(3):555–65.
Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75(5):843–54.
This work was supported by shared resources of the Sidney Kimmel Comprehensive Cancer Center including the Oncology Tissue Services (TMA) Lab Core and Specimen Accessioning Core.
Funding for this project was supported by grants from the National Institutes of Health 5R01CA143299, 5P50CA058236 SPORE in Prostate Cancer Pathobiology Core and, the Patrick C. Walsh Prostate Cancer Research Fund, the Department of Defense Prostate Cancer Research Program Award numbers W81XWH-10-2-0056, W81XWH-10-2-0046 PCRP Prostate Cancer Biorepository Network (PCBN) and DAMD 17-03-273.
Availability of data and materials
Data is available upon request.
SBP observed trends in miRNA stability. SBP, JRB, and EAP performed statistical analyses. QZ performed RNA isolations and qRT-PCR analyses, AKM, AM, EAP, and SEL contributed to study design, data acquisition and data analysis. All authors contributed to the writing, editing and conclusions of the final manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was approved under exemption 4 of the U.S. Department of Health and Human Services Human Subject Regulations because it applied de-identified pathological specimens.
Spearman’s rank correlation of sample value and FFPE block age. Table S2. Spearman’s rank correlation of sample value and RNA Integrity Number. Figure S1. Association of RIN and Sample Block Age. No association was observed between sample age and RNA Integrity Number (RIN). (DOCX 258 kb)