Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparison of MicroRNA Deep Sequencing of Matched Formalin-Fixed Paraffin-Embedded and Fresh Frozen Cancer Tissues

  • Wei Meng,

    Affiliation Department of Radiation Oncology, Wexner Medical Center, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

  • Joseph P. McElroy,

    Affiliation Center for Biostatistics, The Ohio State University, Columbus, Ohio, United States of America

  • Stefano Volinia,

    Affiliation Department of Molecular Virology, Immunology and Medical Genetics, The Ohio State University, Columbus, Ohio, United States of America

  • Jeff Palatini,

    Affiliation Microarray Shared Resource, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

  • Sarah Warner,

    Affiliation Microarray Shared Resource, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

  • Leona W. Ayers,

    Affiliation Department of Pathology, College of Medicine, The Ohio State University, Columbus, Ohio, United States of America

  • Kamalakannan Palanichamy,

    Affiliation Department of Radiation Oncology, Wexner Medical Center, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

  • Arnab Chakravarti ,

    Arnab.Chakravarti@osumc.edu (AC); TL@osumc.edu (TL)

    Affiliation Department of Radiation Oncology, Wexner Medical Center, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

  • Tim Lautenschlaeger

    Arnab.Chakravarti@osumc.edu (AC); TL@osumc.edu (TL)

    Affiliation Department of Radiation Oncology, Wexner Medical Center, Arthur G. James Comprehensive Cancer Center and Richard L. Solove Research Institute, The Ohio State University, Columbus, Ohio, United States of America

Abstract

MicroRNAs regulate several aspects of tumorigenesis and cancer progression. Most cancer tissues are archived formalin-fixed and paraffin-embedded (FFPE). While microRNAs are a more stable form of RNA thought to withstand FFPE-processing and degradation there is only limited evidence for the latter assumption. We examined whether microRNA profiling can be successfully conducted on FFPE cancer tissues using SOLiD ligation based sequencing. Tissue storage times (2–9 years) appeared to not affect the number of detected microRNAs in FFPE samples compared to matched frozen samples (paired t-test p>0.7). Correlations of microRNA expression values were very high across microRNAs in a given sample (Pearson’s r = 0.71–0.95). Higher variance of expression values among samples was associated with higher correlation coefficients between FFPE and frozen tissues. One of the FFPE samples in this study was degraded for unknown reasons with a peak read length of 17 nucleotides compared to 21 in all other samples. The number of detected microRNAs in this sample was within the range of microRNAs detected in all other samples. Ligation-based microRNA deep sequencing on FFPE cancer tissues is feasible and RNA degradation to the degree observed in our study appears to not affect the number of microRNAs that can be quantified.

Introduction

MicroRNAs are noncoding small RNAs with a length of 22 to 26 nucleotides (nt). MicroRNAs result in degradation of complementary messenger RNA in many species. MicroRNAs play an important role in cell development, cell death, cell proliferation [1], apoptosis [2] and angiogenesis [3]. The pioneer bead-based flow cytometric microRNA expression profiling method revealed that microRNA profiles reflect the developmental lineage and differentiation of tumors [4]. Normal tissues and cancer tissues were shown to have distinctive expression profiles and microRNA profiles can describe microRNA-driven pathways in solid tumors [5].

The current methods for microRNA quantification are: real-time reverse transcription-PCR [6], [7], microarray [8], Nanostring, and next generation sequencing [9][11]. Especially, the newly developed next generation sequencing technologies have the potential to discover novel microRNAs and other small RNAs. Most tissues processed in the health care setting will undergo formalin-fixation and paraffin-embedding (FFPE). Most clinical assays performed in pathology laboratories are optimized for FFPE tissues. FFPE tissues are typically stored at room temperature. Formalin preserves tissue samples by creating cross-linking between proteins, DNA and RNA. DNA extracted from FFPE tissue has been shown to be useful for copy number analysis and mutation analysis on at least one platform in some settings [12]. However, the chemical modification between macromolecules and RNAs caused by formalin fixation can accelerate the degradation of RNA [13] [14]. The stability of microRNAs is generally much more robust than messenger RNA [15]. Next generation sequencing based microRNA expression analysis has been developed on several platforms, including Roche/454 platform, Illumina’s Genome Analyzer and ABI's SOLiD platform [16] [17], and different commercial protocols for miRNA library preparation have been developed by the three companies [18]. There is currently only limited evidence available that microRNA sequencing results are reliable and valid. Ma et al. demonstrated the feasibility of miRNA profiling based on Sanger sequencing in 10-year-old archived tissue [19]. To our knowledge there are only two peer-reviewed studies published that compared microRNA sequencing results of matched frozen and FFPE specimens using the Illumina platform [9], [20]. The objective of this study was to determine if FFPE samples can be successfully characterized by deep miRNA sequencing. To this end we analyzed matched frozen and FFPE tissues of different histologies for which detailed processing and storage information was available.

Materials and Methods

Ethics Statement

Human malignant tissue samples procured during period 2003–2009 were obtained from the Cooperative Human Tissue Network (CHTN/NCI), Midwestern Division, which is funded by the National Cancer Institute, based on an internal review board (IRB) approved research protocol (Ohio State University Human Subjects Protocol - 2011E0377). Other Investigators may have received specimens from the same subjects.

Tissue specimens

Eight paired human malignant tissue samples were received. Remnant surgical tissues were taken after diagnostic samples were secured from patients (five males and four females; median age 58 years, range 39–78 years) with breast invasive ductal carcinoma (2), renal clear cell carcinoma (2), lung adenocarcinoma (1), prostate adenocarcinoma (1), metastatic melanoma (1), and sarcoma of thigh (1). Surgical specimens were transported from the operating rooms to Tissue Procurement, examined under a pathologist’s supervision and preserved within 5 to 57 recorded minutes. Paired tissue samples selected for research were either immediately snap frozen in liquid nitrogen and placed in a –80°C freezer for monitored storage or fixed in 10% buffered formalin for up to 24 hours and then processed into a paraffin embedded block and stored at room temperature. When the research samples were received from CHTN, the coded samples were accompanied by a coded final pathology report and a quality assessment report verifying collection of tumor tissue. The specimens have been reviewed by a pathologist and 4 cases had 100% tumor, 3 cases 90–95% and 1 case 50–60% tumor, and matched FFPE and frozen samples were generally comparable in tumor content and size. The clinicopathologic features of the cases used for microRNA sequencing and expression analysis are listed in table 1.

thumbnail
Table 1. The clinicopathologic features of the eight cases used for microRNA sequencing and expression analysis.

https://doi.org/10.1371/journal.pone.0064393.t001

RNA Isolation, quantitation, and quality Assessment

RNA from FFPE samples were extracted using Recover All Total Nucleic Acid Isolation Kit (Life Technologies Corporation, Carlsbad, CA). Briefly, 5–10 mg samples were sliced from paraffin blocks and deparaffinizated by xylene at 50°C, followed by 100% ethanol wash. The air-dry tissue samples were digested by proteinase K for 24 hrs in a microtube shaking incubator set at 50°C. The digested samples were mixed with appropriate volume of isolation additive and 100% ethanol. After passing the mixture through the filter cartridge, the DNA and RNA were retained on the filter. The DNA was removed by on-filter DNase digestion. The RNA was purified by washing buffer and eluted with nuclease-free water.

The microRNA samples from the fresh frozen samples were extracted using PureLink™ miRNA Isolation Kit (Life Technologies Corporation, Carlsbad, CA). Briefly, 5 mg fresh frozen tissue samples were mixed with 300 µl Binding Buffer (L3) and homogenized using tissue homogenizer. The homogenized samples were centrifuged, and supernatant were mixed with 300 µl 70% ethanol. The solution containing the RNA was purified by two round of spin columns cleanup. The RNA was eluted with 50–100 µl sterile RNase-free water.

MicroRNA sequencing

Small RNAs from FFPE and fresh frozen samples were prepared for SOLiD sequencing as follows: the total RNA samples were processed by the flashPAGE fractionator (Ambion) and flashPAGE Clean-Up Kit (Ambion). The enriched small RNA was then processed according to the SOLiD Small RNA Expression Kit protocol (Applied Biosystems). The purified small RNAs were ligated with 5' and 3' adapter mix using RNA ligase. The ligated products (40–60 bases in length) were reverse transcribed and purified on Novex 10% TBE-Urea gel. Subsequently, 15–18 cycles of PCR were performed by amplifying the purified cDNA with barcoded PCR primer sets provided in the kit, which differed by a unique 6-nucleotide sequence. The amplified products were loaded on Novex 6% TBE gel (Invitrogen) and the gel bands containing 110 to 130 bp fragments were excised. The amplified products were purfied from excised gel band, amplified by emulsion PCR, then loaded on Applied Biosystems SOLiD 4 next generation high throughput sequencing system for data acquisition. The quality of the samples and libraries were verified on the Agilent Bioanalyzer [21], [22]. The raw SOLiD sequencing data were uploaded to the NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE45740 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45740).

RT-qPCR for microRNA validation

The RNA from fresh frozen and FFPE tissues were validated by real-time (RT) quantitative PCR. microRNA expression was normalized to the small nuclear RNA U6 using TaqMan MicroRNA Assay kits (Applied Biosystems, Darmstadt, Germany). A 30 ng sample of RNA was processed by the TaqMan MicroRNA Reverse Transcription kit (Applied Biosystems, Darmstadt, Germany). In brief, 5 µl RNA was mixed with 1 mmol/l of each deoxyribonucleotide triphosphate, 50 units of Multiscribe Reverse Transcriptase, 5× reaction buffers, 4 units RNase inhibitor, and 5× gene-specific RT primers mix in a final reaction volume of 15 µl. Reactions were then incubated at 16°C for 30 min, 42°C for 30 min, 85°C for 15 min with a final hold at 4°C. Afterward, 15 µl of the cDNA solution was diluted by nuclease-free water to a final volume of 100 µl. Quantitive PCR was run on a 7900 HT PCR machine with denaturation step at 95°C for 10 minutes, followed by 40 cycles of a denaturation step at 95°C for 15 seconds, and an annealing/elongation step at 60°C for 60 seconds. Fold change for microRNA in tissue samples was then calculated using the equation 2-ΔCt test/2-ΔCt control.

Data analysis

The linker sequences of small RNA libraries were removed using cutadapt software (http://code.google.com/p/cutadapt/) under default parameters. The length distribution of microRNA was constructed by plotting reads counts against 0–30 nucleotide lengths. The small RNA Analysis Pipeline Tool from Applied Biosystems was used with the following parameters to include only sufficiently high quality reads in our study: 1) Minimal average read quality threshold was 20 (The QV value is calculated using a phred like score q = –10× log10 (p) where q is the quality value and (p) is the predicted probability that the color call is incorrect.), 2) to be stringent we allowed only 1 mismatch between reads and miRBase version 16 hits for each counted read 3) the minimal length of aligned reads had to be 17 bases. Data were normalized as reads per million (RPM). Sequences with expression value below 5 RPMs were considered not detected for the correlation between frozen and FFPE analysis part.

Results

Effects of FFPE sample degradation: microRNA read length distribution

Formalin fixation and long storage times are known to result in RNA degradation [23]. This degradation results in a shorter average RNA length [24]. To identify possible degradation effects of microRNAs, we analyzed the read length distribution of matched fresh frozen and FFPE cancer tissue samples. After the 3’ adaptor sequences of each read were trimmed, sequence reads length distribution was examined (figure 1). Read lengths were observed to be centered at 21 nt in all frozen samples and all but one FFPE samples. There was one FFPE sample whose peak microRNA length was found at 17 nt, indicating significant degradation. We then carefully examined the small RNA classification profiles and correlation of microRNA expression values of fresh frozen and FFPE samples. We found that there are negligible or no effects on those analyses. Storage times or preparation methods were not different between this degraded FFPE sample and the other samples analyzed in this study. Also, the proportion of small RNAs in the 20–22 nt range was lower for several cases in the FFPE compared to the matched frozen specimens (one renal cell carcinoma case, the lung adenocarcinoma case, and the prostate cancer case). This is also likely due to degradation in FFPE samples, although small variations during the library preparation could partially contribute to the observed differences.

thumbnail
Figure 1. Sequence reads length distribution among eight paired fresh frozen and FFPE samples.

https://doi.org/10.1371/journal.pone.0064393.g001

Effects of FFPE sample degradation: storage time

The storage time of the samples ranged from 2 to 9 years. 624 micoRNAs and their precursors were identified in the 8 pairs of tissue samples. Some microRNAs could not be detected in either fresh frozen or FFPE samples. Figure 2 shows the percentages of microRNAs that were detected in both fresh frozen and FFPE samples, only in fresh frozen samples and only in FFPE samples over the year of specimen collection. There were no storage-time-dependent changes in those percentages, indicating that there is no significant association between storage time and the number of microRNAs present during this 9 year time window.

thumbnail
Figure 2. Percentages of microRNA reads detected over the years of specimen collection.

https://doi.org/10.1371/journal.pone.0064393.g002

We next examined the percentage of microRNA reads among total reads (figure 3) to understand if the purity of microRNAs was reduced in the FFPE samples. Approximately 70% and 85% of reads in FFPE and frozen samples, respectively, were identified as microRNA sequences. A possible reason for this difference could be the presence of fragments of lncRNAs and mRNAs in the fraction of RNA smaller than 40 nucleotides (nt) in FFPE tissues (figure 4).

thumbnail
Figure 3. Percentage of microRNA reads among total reads.

https://doi.org/10.1371/journal.pone.0064393.g003

thumbnail
Figure 4. The classification of small RNA populations in fresh frozen samples and FFPE samples.

https://doi.org/10.1371/journal.pone.0064393.g004

Effects of FFPE sample degradation: classification of small RNAs

To better understand if the difference in the percentage of small RNA sequences among all reads between frozen and FFPE tissue data we examined the distribution of RNA classes in one of the breast cancer sample. There were a total of 11563414 and 10273837 sequencing reads, in the two examined matched samples and the reads were mapped to the human genome (hg19, National Center of biotechnology information building) using custom Python scripts. The reads were classified in different RNA groups according to the EnsEMBL database. The majority of small RNAs in both fresh frozen and FFPE libraries were microRNAs (figure 4). Comparing data from matched fresh frozen samples to data from FFPE samples there were lincRNA and sequences of unknown classes in data obtained from FFPE samples. Given that only RNAs smaller than 40 nt were extracted and processed for sequencing this suggests that there is fragmentation of lincRNA and other long RNA in FFPE tissues, which are longer than 40 nt.

A total of 469 to 548 (mean 519) different microRNAs and their precursors were identified in both fresh frozen and FFPE samples. We determined the percentages of microRNAs that have more than a 2 fold difference of expression between matched fresh frozen and FFPE samples (figure 5). Among the different cases, 69 to 145 microRNA transcripts were more than 2 fold overexpressed in fresh frozen samples. 73 to 186 microRNA transcripts were more than 2 fold overexpressed in FFPE tissues (figure 5).

thumbnail
Figure 5. The number of microRNAs that have more than a 2 fold expression change between matched fresh frozen and FFPE samples.

https://doi.org/10.1371/journal.pone.0064393.g005

Correlation of microRNA expression values of matched FFPE and frozen samples

We next determined the correlation of microRNA expression values of matched FFPE samples and fresh frozen samples. There were 250 microRNAs detected in all samples in this study. There was a strong correlation of microRNA expression values in matched FFPE and snap fresh frozen tissues across all microRNAs and all cases (r  = 0.85; figure 6). Correlations across detected microRNAs from the same tumor/patient were high, with r ranging from 0.71–0.95 (figure 7) [25].

thumbnail
Figure 6. Regression of Frozen miR values onto FFPE miR values across all individuals and miRs.

Only the 250 miRs with no missing data were used.

https://doi.org/10.1371/journal.pone.0064393.g006

thumbnail
Figure 7. Scatter plots of detected microRNAs (250 miRs with no missing data) profiles obtained from matched fresh frozen-FFPE samples by SOLiD 4 system.

Each dot represents the expression values of one microRNA in a fresh frozen-FFPE pair.

https://doi.org/10.1371/journal.pone.0064393.g007

The mean of Pearson correlation coefficients of individual microRNA expression values of FFPE and frozen tissues was 0.46 (figure 8), with only 41 mature microRNAs and precursors having correlation coefficients higher than 0.8 (table 2). We sought out to identify features that are associated with higher correlation coefficients indicating invariance to the tissue preservation method. To test if the level of expression of a given microRNA is associated with the correlation between fresh frozen and FFPE tissue, Fisher transformed (arctanh) correlations were regressed onto mean FFPE expression (figure 9). There was no association between expression values and correlations between data from frozen and FFPE tissue observed. We then sought out to understand if a possible lack of variation of individual microRNAs among the biological replicates is resulting in poor correlations. Fisher transformed correlations were regressed onto log transformed FFPE variances. The variance analysis showed that microRNAs with higher variance have significantly higher correlations of data from FFPE and frozen tissues (figure 10; arctanh(r)  = 0.47+0.59×log10 (FFPE variance)), which may be due to the presence of true biological variation present in these microRNAs.

thumbnail
Figure 8. MicroRNA Pearson correlation coefficient distribution of eight matched fresh frozen and FFPE samples among the 250 microRNAs detected in all samples.

https://doi.org/10.1371/journal.pone.0064393.g008

thumbnail
Figure 9. Scatter plot of microRNAs expression mean values versus Pearson correlation coefficients (Fisher transformed) from matched fresh frozen-FFPE samples.

https://doi.org/10.1371/journal.pone.0064393.g009

thumbnail
Figure 10. Scatter plot of microRNAs variance (log10 transformed) versus Pearson correlation coefficients (Fisher transformed) from matched fresh frozen-FFPE samples.

https://doi.org/10.1371/journal.pone.0064393.g010

thumbnail
Table 2. 41 microRNAs had correlation coefficients equal to or greater than 0.8 between FFPE and fresh frozen pairs.

https://doi.org/10.1371/journal.pone.0064393.t002

PCR validation of microRNA expression

We selected two microRNAs (miR-21 and miR-19a) that were expressed higher in frozen samples than in FFPE samples in our microRNA sequencing results for validation using PCR (Taqman). Both miR-21 and miR-19a were found to also be higher expressed in frozen tissue than in matched FFPE samples using a PCR assay (figure 11), consistent with our sequencing results.

thumbnail
Figure 11. Validation of differentiated microRNA expression profile by RT-qPCR method.

RNA samples from fresh frozen (n = 8) and FFPE (n = 8) tissues were analyzed. Small nuclear RNA U6 was used as reference. MiR-21 and miR-19a had higher expression in fresh frozen samples than in FFPE samples, consistent with sequencing results.

https://doi.org/10.1371/journal.pone.0064393.g011

Discussion

MicroRNAs as gene expression regulators appear to have an important role in cancer cell differentiation [4], proliferation [26], metastasis [27] and apoptosis [28]. There have been many microRNA profiling efforts using frozen cancer tissues [29] [30] [31]. Cancer tissue specimens are typically archived after formalin fixation and paraffin embedding (FFPE) and frozen tissue is only available for a minority of cases in pathology archives and tissue banks. A better understanding of the quality of microRNA sequencing results of FFPE tissues would allow for a more informed decision as to when use of FFPE tissue for small RNA sequencing experiments would be appropriate.

In this study, we analyzed small RNA sequencing profiles of eight pairs of matched fresh frozen and FFPE cancer tissue samples using the ABI SOLiD 4 platform. Storage time up to nine years appeared to not affect the number of detectable microRNAs in either frozen of FFPE tissue in this study. The presence of higher percentages of non-microRNAs in the fraction of nucleotides smaller than 40 nt suggested that FFPE tissues are more degraded, though. However, this apparent degradation did not result in fewer microRNAs being detected in FFPE tissues, suggesting that a sequencing approach can overcome sample degradation to some degree. While small RNAs in one of eight FFPE samples appeared degraded as suggested by the reduced read lengths centered at 17 nt compared to the other eight FFPE and all eight frozen samples centered at 21 nt, the number of aligned reads and the number of detected microRNAs in this degraded sample were similar to the data from other samples, which suggests that data obtained from samples with the above described degree of degradation would not necessarily need to be removed from all types of analysis.

Overall, there was excellent correlation between miRNA expression values between matched frozen and FFPE when combining all samples and all microRNAs (fig. 6), or when analyzing all microRNAs within individual patient sample pairs (fig. 7). Those findings are consisted with previous reports using different sequencing platforms suggesting that sequencing approaches can result in good correlations between FFPE and frozen samples when including several hundred microRNAs [32] [33] [9] [34].

Correlations of expression values of individual microRNAs (see figs. 8, 9, 10) were lower than correlations including all microRNAs (figs 6, 7). The underlying reason for this observation likely originates in the identified association of correlation coefficients of individual microRNAs and variance of microRNA expression values (fig 10). Very similar expression values of an individual microRNA among the eight cases (equal to low variance in fig. 10) typically result in low correlation coefficients because the variation between frozen and FFPE data is larger than the variation of expression between the eight cases included in this study. Therefore it is not necessarily surprising to have a wide distribution of correlations of individual microRNAs as analyzed in detail in figures 8, 9, and 10.

Limitations of this analysis include limited small sample size and potential PCR drifts and PCR selection biases during amplification in the multiplex PCR based library construction [35] [36]. This disadvantage should be overcome by the 3rd generation of sequencing technology, which is thought to not require amplification of the target sequence anymore but detect single nucleic acid molecules [37]. Furthermore, the use of adjacent tissue to create matched specimens for this analysis introduces some tissue heterogeneity effects and likely contributes to the observed differences between matched frozen and FFPE samples. In summary, microRNA sequencing of FFPE tissues is feasible using ligation-based sequencing and microRNA expression values of individual patients are highly correlated with expression data from matched frozen tissues. MicroRNAs with higher variance appear to have a higher correlation of their expression values in matched frozen and FFPE tissue data.

Acknowledgments

We thank Carlo M. Croce for helpful discussions.

Author Contributions

Conceived and designed the experiments: WM TL. Performed the experiments: WM JP SW. Analyzed the data: WM JPM SV TL. Contributed reagents/materials/analysis tools: JPM SV JP LWA KP AC TL. Wrote the paper: WM JPM LWA TL.

References

  1. 1. Ambros V (2004) The functions of animal microRNAs. Nature 431: 350–355.
  2. 2. Chang TC, Wentzel EA, Kent OA, Ramachandran K, Mullendore M, et al. (2007) Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis. Mol Cell 26: 745–752.
  3. 3. Dews M, Homayouni A, Yu D, Murphy D, Sevignani C, et al. (2006) Augmentation of tumor angiogenesis by a Myc-activated microRNA cluster. Nat Genet 38: 1060–1065.
  4. 4. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, et al. (2005) MicroRNA expression profiles classify human cancers. Nature 435: 834–838.
  5. 5. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, et al. (2006) A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A 103: 2257–2261.
  6. 6. Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, et al. (2005) Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 33: e179.
  7. 7. Shi R, Chiang VL (2005) Facile means for quantifying microRNA expression by real-time PCR. Biotechniques 39: 519–525.
  8. 8. Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, et al. (2004) An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci U S A 101: 9740–9744.
  9. 9. Weng L, Wu X, Gao H, Mu B, Li X, et al. (2010) MicroRNA profiling of clear cell renal cell carcinoma by whole-genome small RNA deep sequencing of paired frozen and formalin-fixed, paraffin-embedded tissue specimens. J Pathol 222: 41–51.
  10. 10. Morozova O, Marra MA (2008) Applications of next-generation sequencing technologies in functional genomics. Genomics 92: 255–264.
  11. 11. Mendes ND, Freitas AT, Sagot MF (2009) Current tools for the identification of miRNA genes and their targets. Nucleic Acids Res 37: 2419–2433.
  12. 12. Schweiger MR, Kerick M, Timmermann B, Albrecht MW, Borodina T, et al. (2009) Genome-wide massively parallel sequencing of formaldehyde fixed-paraffin embedded (FFPE) tumor tissues for copy-number- and mutation-analysis. PLoS One 4: e5548.
  13. 13. Scicchitano MS, Dalmas DA, Bertiaux MA, Anderson SM, Turner LR, et al. (2006) Preliminary comparison of quantity, quality, and microarray performance of RNA extracted from formalin-fixed, paraffin-embedded, and unfixed frozen tissue samples. J Histochem Cytochem 54: 1229–1237.
  14. 14. Xi Y, Nakajima G, Gavin E, Morris CG, Kudo K, et al. (2007) Systematic analysis of microRNA expression of RNA extracted from fresh frozen and formalin-fixed paraffin-embedded samples. RNA 13: 1668–1674.
  15. 15. Jung M, Schaefer A, Steiner I, Kempkensteffen C, Stephan C, et al. (2010) Robust microRNA stability in degraded RNA preparations from human tissue and cell samples. Clin Chem 56: 998–1006.
  16. 16. Soares AR, Pereira PM, Santos MA (2012) Next-generation sequencing of miRNAs with Roche 454 GS-FLX technology: steps for a successful application. Methods Mol Biol 822: 189–204.
  17. 17. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, et al. (2008) Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res 18: 610–621.
  18. 18. Tian G, Yin X, Luo H, Xu X, Bolund L, et al. (2010) Sequencing bias: comparison of different protocols of microRNA library construction. BMC Biotechnol 10: 64.
  19. 19. Ma Z, Lui WO, Fire A, Dadras SS (2009) Profiling and discovery of novel miRNAs from formalin-fixed, paraffin-embedded melanoma and nodal specimens. J Mol Diagn 11: 420–429.
  20. 20. Li S, Yang C, Zhai L, Zhang W, Yu J, et al. (2012) Deep sequencing reveals small RNA characterization of invasive micropapillary carcinomas of the breast. Breast Cancer Res Treat 136: 77–87.
  21. 21. Norden-Krichmar TM, Allen AE, Gaasterland T, Hildebrand M (2011) Characterization of the small RNA transcriptome of the diatom, Thalassiosira pseudonana. PLoS One 6: e22870.
  22. 22. Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, et al. (2008) Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol 26: 407–415.
  23. 23. Ferruelo A, El-Assar M, Lorente JA, Nin N, Penuelas O, et al. (2011) Transcriptional profiling and genotyping of degraded nucleic acids from autopsy tissue samples after prolonged formalin fixation times. Int J Clin Exp Pathol 4: 156–161.
  24. 24. Zhang X, Chen J, Radcliffe T, Lebrun DP, Tron VA, et al. (2008) An array-based analysis of microRNA expression comparing matched frozen and formalin-fixed paraffin-embedded human tissue samples. J Mol Diagn 10: 513–519.
  25. 25. Lake SL, Kalirai H, Dopierala J, Damato BE, Coupland SE (2012) Comparison of formalin-fixed and snap-frozen samples analysed by multiplex ligation-dependent probe amplification for prognostic testing in uveal melanoma. Invest Ophthalmol Vis Sci.
  26. 26. Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T (2007) Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nature Genetics 39: 673–677.
  27. 27. Tavazoie SF, Alarcon C, Oskarsson T, Padua D, Wang QQ, et al. (2008) Endogenous human microRNAs that suppress breast cancer metastasis. Nature 451: 147–U143.
  28. 28. Welch C, Chen Y, Stallings RL (2007) MicroRNA-34a functions as a potential tumor suppressor by inducing apoptosis in neuroblastoma cells. Oncogene 26: 5017–5022.
  29. 29. Lin J, Teo S, Lam DH, Jeyaseelan K, Wang S (2012) MicroRNA-10b pleiotropically regulates invasion, angiogenicity and apoptosis of tumor cells resembling mesenchymal subtype of glioblastoma multiforme. Cell Death Dis 3: e398.
  30. 30. Zhang W, Zhang J, Yan W, You G, Bao Z, et al.. (2012) Whole-genome microRNA expression profiling identifies a 5-microRNA signature as a prognostic biomarker in Chinese patients with primary glioblastoma multiforme. Cancer.
  31. 31. Zinn PO, Sathyan P, Mahajan B, Bruyere J, Hegi M, et al. (2012) A novel volume-age-KPS (VAK) glioblastoma classification identifies a prognostic cognate microRNA-gene signature. PLoS One 7: e41522.
  32. 32. Li J, Smyth P, Flavin R, Cahill S, Denning K, et al. (2007) Comparison of miRNA expression patterns using total RNA extracted from matched samples of formalin-fixed paraffin-embedded (FFPE) cells and snap frozen cells. BMC Biotechnol 7: 36.
  33. 33. Glud M, Klausen M, Gniadecki R, Rossing M, Hastrup N, et al. (2009) MicroRNA expression in melanocytic nevi: the usefulness of formalin-fixed, paraffin-embedded material for miRNA microarray profiling. J Invest Dermatol 129: 1219–1224.
  34. 34. de Biase D, Visani M, Morandi L, Marucci G, Taccioli C, et al. (2012) miRNAs expression analysis in paired fresh/frozen and dissected formalin fixed and paraffin embedded glioblastoma using real-time pCR. PLoS One 7: e35596.
  35. 35. Moltzahn F, Olshen AB, Baehner L, Peek A, Fong L, et al. (2011) Microfluidic-based multiplex qRT-PCR identifies diagnostic and prognostic microRNA signatures in the sera of prostate cancer patients. Cancer Res 71: 550–560.
  36. 36. Elnifro EM, Ashshi AM, Cooper RJ, Klapper PE (2000) Multiplex PCR: optimization and application in diagnostic virology. Clin Microbiol Rev 13: 559–570.
  37. 37. Kumar S, Tao C, Chien M, Hellner B, Balijepalli A, et al. (2012) PEG-Labeled Nucleotides and Nanopore Detection for Single Molecule DNA Sequencing by Synthesis. Sci Rep 2: 684.