|Home | About | Journals | Submit | Contact Us | Français|
microRNAs (miRNAs) are ~22nt RNAs that regulate target gene expression. Altered expression of miRNAs has been demonstrated in many different human cancers. Many studies using microarray technologies to characterize miRNA expression profiles have relied on fresh tissue to determine the miRNA signatures. In this study, we prepared total RNA from paired samples of formalin-fixed paraffin-embedded (FFPE) and fresh frozen malignant melanoma, and used that in microarray experiments to compare miRNA expression profiles between FFPE and fresh tissue with corresponding mRNA expression profiles from the same tissue sources. We demonstrate that miRNA expression profile from FFPE tissues closely resembles that from fresh tissues, and the correlation is significantly better than that for mRNA profiles from FFPE and fresh tissues. These results underscore the suitability of FFPE tissues as appropriate resources for molecular expression analyses and support the notion that miRNAs are more vigorous analytes for this purpose than mRNAs.
MicroRNAs (miRNAs) are non-coding RNAs of approximately 20-22 nucleotides in length in their active form. They function in a diverse set of biological processes, including differentiation, proliferation and apoptosis. They are transcribed as longer precursor RNAs of ~60-110 nucleotides in length and form a stem loop structure (pri-miRNAs and pre-miRNAs). Their maturation occurs through a series of nuclear and cytoplasmic cleavage events mediated by the Drosha and Dicer endonucleases, which result in the formation of a ~22 nucleotide RNA:RNA duplex. One strand of this duplex is incorporated into the RNA Induced Silencing Complex (RISC). RISC is a multi-protein complex that first facilitates the annealing between the active single stranded miRNA to its target mRNA, and then mediates the down-regulation of target mRNA gene expression. The latter occurs either through target mRNA cleavage or translational repression [1-3].
Dysregulation of miRNA expression has been linked to cancer. Alterations in the expression profiles of many miRNAs have been described in numerous human tumors [2, 4]. Larger scale miRNA analyses underscore the robustness of miRNA expression signatures as informative molecular surrogates, facilitating the distinction of benign tissues from their malignant counterparts, as well as the utility of miRNA expression signatures in clustering histologically ambiguous tumors according to their cell of origin [5-7].
Many different platforms have been devised to systematically interrogate miRNAs, including northern blots , qRT-PCR , and microarrays [10-12]. The gold standard of source tissue for miRNA expression profiling studies has traditionally been fresh-frozen tissue. However, given the enormous amount of physiologic information stored in archived formalin-fixed paraffin-embedded tissue (FFPE) and the abundance of clinical data retrievable in association with it, it will prove invaluable if miRNA expression levels could be routinely and systematically analyzed in FFPE tissues, particularly for retrospective studies and for the characterization of rare or small tumors. Indeed, numerous studies have documented that it is feasible and informative to characterize miRNA from frozen tissue, and an increasing number of studies have demonstrated this with FFPE using qRT-PCR [13-16], as well as microarrays [6, 12, 16-20]. Consistent among these studies is a strong correlation of miRNA expression levels between paired FFPE and frozen tissues that is largely independent of time of fixation as well as storage time in paraffin (up to approximately 10 years without significant compromise of miRNA expression detection) [15-17]. Although relatively few studies have directly compared miRNA and mRNA expression levels from paired FFPE and frozen samples, these similarly demonstrate a consistent greater correlation of miRNA expression levels compared to mRNA expression levels between paired FFPE and fresh tissues [13, 16].
In this study, we confirm the reproducibility of extracting adequate amounts of intact mRNA and miRNA from paired frozen and FFPE tissue samples of malignant melanoma to characterize the relative miRNA and mRNA expression profiles in these matched tissue samples. However, we demonstrate that miRNA is a far more vigorous analyte than mRNA for the interrogation of FFPE tissue. These results underscore the robustness of miRNA expression profiling and more importantly, support the applicability of FFPE tissues as important resources to characterize important molecular pathways and targets in human malignancy.
All tissues were collected according to the guidelines and policies of the Hospital of the University of Pennsylvania Institutional Board. Fresh frozen melanoma tissues were divided and half kept frozen and half fixed in formalin and processed for paraffin embedded. Over 90% of the tissue is composed of tumor cells. For frozen tissue, three (3) 10 μm sections per biological sample were obtained for the extraction of total RNA using mirVana (Ambion #1560) isolation kit for both mRNA and miRNA profiling, according to the manufacturer's instructions. Four (4) 20 μm sections per sample were obtained from the selected block for the extraction of total RNA using the RecoverAll™ Total Nucleic Acid Isolation system (Ambion, Cat# 1975) according to the manufacturer's instructions. The total RNA yield was determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Tech, Wilmington, DE). RNA integrity was assessed by Agilent 2100 Bioanalyzer electrophoresis compared to standard reference RNA as previously described for total RNA prepared from FFPE tissue [12, 19, 21].
The NuGEN WT-ovation FFPE system (Cat# 3400) was utilized for all amplifications, according to the manufacturer's instructions. Amplifications were performed starting with 100 ng of total RNA for both fresh frozen samples and FFPE. The amplified cDNA product was quantitated using a NanoDrop ND-1000. 5 μg of amplified cDNAwas labeled with the NuGEN FL-Ovation cDNA Biotin Module V2 (Cat# 4200) following the manufacturer's instructions. Fragmented and biotin-labeled target was hybridized to Affymetrix HG-U133A Plus2.0 GeneChip arrays. Hybridization, washing, staining and scanning were performed according to manufacturer's instructions.
Targets were prepared using the miRCURY LNA™ microRNA Power Labeling Kit (Exiqon #208032-A) for Exiqon miRNA microarrays according to the manufacturer's instructions. Labeling reactions were performed essentially as previously described  using equivalent 1μg aliquots of paired (frozen and FFPE) melanoma samples. The labeled miRNA was applied to miRCURY LNA microRNA array (Exqion # 208002-A) according to the manufacturer's instructions for hybridization and washing. A corresponding dye-swap sample pair was hybridized to a second microarray. For mRNA profiling, biotin-labeled cDNAs (5 μg per sample) were applied to Human Affymetrix (U133+V2.0) chips and hybridized according to the manufacturer's instructions.
GenePix results files (GPR) were imported into Partek Genomics Suite (v6.4, Partek, Inc., St. Louis, MO). Normalization of miRNA array data was performed as follows. The on-chip replicate signal intensities (quadruplicate spots of a particular miRNA sequence) were first averaged for each identifier (miRNA sequence). Next, for each identifier, we subtracted the background median signal from the averaged signal intensity. Any value less than 1 was assigned a value of 1 for ease of future manipulations (see below). LoESS normalization was applied to correct for intensity dependent dye effects, yielding log2 ratios for each probe. For each dye swapped pair of frozen and FFPE tissues, the average log2 ratio of FFPE:frozen signal was calculated, imported into an Excel spread sheet (normalized and averaged ratio of FFPE tissue signal/frozen tissue signal).
Affymetrix mRNA expression data was generated using a single color array platform. 34,321 discrete probesets were selected following a raw intensity filter across all arrays. In order to facilitate a comparison with the miRNA expression data which was generated from a two-color array platform, this expression data was also converted into a ratio. For each patient sample, the signal ratio (FFPE tissue/frozen tissue) was calculated for each of the 34,321 individual spots, imported into an Excel spread sheet and log2 transformed (normalized and averaged ratio of FFPE tissue signal/frozen tissue signal).
These data were then plotted. The standard deviation from a mean of zero (0) was calculated for both miRNA signals and mRNA signals.
Raw Affymetrix data files were processed as above with the exception of the ratio calculation. This yielded log2 transformed intensities for 34,321 probesets for each of the 6 samples. For the miRNA data, intensities were extracted as above, yielding 2 intensities (one for each dye) for each of the 6 biological samples. For each sample the dye-specific intensities were averaged, yielding a single average intensity for each sample. Any value less than one was assigned a value of 1 and the data were log2 transformed. Pairwise Pearson correlations were calculated for the miRNA data. These are reported in a table with the color of each cell corresponding to the r-value.
Both miRNA and Affymetrix intensities were clustered hierarchically using Euclidean distance and average linkage. The dendograms displaying the relative similarity relationship among the samples from each dataset are presented in Figure 4.
Total RNA was prepared from matched frozen and FFPE samples of malignant melanoma. All sample sources (matched frozen and FFPE) yielded adequate amounts of total RNA for labeling and microarray hybridization. Inspection of RNA by Agilent Bioanalyzer electrophoresis demonstrated a discrete miRNA peak from the frozen tissue and a similar, albeit slightly obscured peak from the FFPE tissue (Figure 1). These observations are consistent with previous studies, which described a similar profile and attributed this obscuring peak to degraded mRNA, tRNA and rRNA .
We first examined the similarity between miRNA expression profiles between paired fresh-frozen and FFPE samples from the same patient using intensity data from the microarrays (see Methods). As illustrated in Figure 2, the miRNA expression profiles from FFPE tissue correlate most highly with those generated from the matched fresh frozen tissue from the same original tissue source (r value ranges from 0.956 to 0.976).
We compared the relative miRNA and mRNA expression levels between paired fresh and FFPE samples of malignant melanomas. In order to facilitate a comparison between the miRNA expression data - generated from a two-color array platform - and the mRNA expression data - generated from a single color Affymetrix array platform, we chose to represent all data as a log2-ratio of FFPE:frozen signal for each spot. In this way, we could compare the relative distribution of these log2-transformed ratios as a reflection of the relative variance of signals generated from FFPE compared to matched frozen tissue for the complete set of miRNAs versus the complete set of mRNAs. Minimal variance between the FFPE and frozen sample would be predicted to generate a signal ratio of one (1) and a log2-(ratio) of zero (0). In contrast, discrepancies between these values would produce a signal ratio unequal to one and therefore, a log2-(ratio) unequal to zero (0). As shown in Figure 3, the log2-ratios for mRNA signals exhibit significantly greater dispersion about 0 than the log2-ratios for miRNA signals. Consistent with this, the standard deviation from the mean of zero (0) for the log2-miRNA ratios was 0.357, compared to a standard deviation of 1.001 for the log2-mRNA ratios from the same samples.
Using the intensity data, miRNA and Affymetrix signal intensities were subjected to hierarchical clustering (see Methods). The dendograms displaying the relative similarity relationship among the samples from each dataset are presented in Figure 4. This clustering analysis demonstrates that miRNA expression signatures cluster most closely according to the origin of the sample: the miRNA expression profiles are most similar between paired fresh-frozen and FFPE (i.e. from an individual patient). In contrast, mRNA expression signatures cluster according to the method by which the tissue was prepared (fresh-frozen versus FFPE, see Figure 4) and are independent of the individual patient from which a given tissue source was procured.
In this study, we extracted mRNA and miRNA from paired frozen and FFPE tissue samples of malignant melanoma for expression profiling experiments. To our knowledge, this is the first study directly comparing both miRNA as well as mRNA relative expression levels in paired fresh and FFPE samples of primary malignant melanoma. Using microarrays and total RNA derived from benign melanocytic nevi in paired fresh and FFPE tissue samples, Glud et al demonstrated a similarly strong correlation for miRNA signals (overall Spearman correlation coefficient of 0.80) . We demonstrate that although there is degradation observed in the total RNA prepared from FFPE tissue compared to concomitantly prepared matched fresh-frozen tissue, miRNA expression profiles from FFPE malignant melanoma samples faithfully recapitulate those from matched fresh frozen tissue; this correlation is significantly better comparing miRNA versus mRNA employing distinct statistical approaches. Indeed, the smaller standard deviation for the log2-miRNA ratios of fresh:FFPE (0.357) compared to that for log2-mRNA ratios (1.001) confirms a greater degree of correlation for miRNA - compared to mRNA - expression profiles in paired fresh-frozen and FFPE tissue samples. These results are consistent with previous studies demonstrating that miRNA is a more vigorous analyte than mRNA for the interrogation of FFPE tissue.
The standard of source tissue for most miRNA expression profiling studies has traditionally been fresh-frozen tissue. However, there is an increasing body of evidence demonstrating that FFPE is an equally reliable resource for this purpose [5, 6, 12-21]. A subset of these studies employed RT-PCR to directly compare miRNA expression levels between paired FFPE and frozen tissues to mRNA expression levels from the same tissue pairs [13, 16]. Using RT-PCR, Xi et al demonstrated significantly greater degrees of correlation for miRNA (r2=0.86-0.89) than for mRNA (r2=0.28-0.56) from paired FFPE and fresh samples of rat liver . Our results are in agreement with and expand upon these studies. Using a microarray platform, we demonstrated a high degree of correlation between miRNA expression signatures in paired fresh versus FFPE samples from the same patient sample (r value from 0.956 to 0.976, Figure 2). Using a different statistical approach to compare the miRNA data to the mRNA data, we demonstrated a significantly greater correlation of miRNA expression levels between paired FFPE and fresh tissue samples of malignant melanoma compared to mRNA expression levels from the same tissues (Figure 3). Most importantly, when the intensity data is subjected to hierarchical clustering analysis, miRNA expression profiles from fresh frozen and FFPE tissue sources cluster most closely according to the patient from which they were prepared. In contrast, mRNA expression profiles are most similar among the fresh frozen samples and FFPE, respectively independent of the source from which they were derived. These results argue that miRNA expression profiles derived from FFPE samples faithfully recapitulate the distinctive biology of each individual patient. In contrast, mRNA expression profiles cluster together according to the method by which the tissue was processed (fresh frozen versus FFPE) rather than according to the patient's identity. These data argue that the processing steps in FFPE potentially dampen biologically informative differences between individual patient samples. Instead, processing would appear to impart a more generalized, less unique “fixation mRNA signature”.
Many explanations have been proposed to explain these observations. Systematic analyses describing the chemical modifications RNA undergoes as a result of formalin fixation show that the principle change includes the addition of a hydroxymethyl group (termed “methylol”) to the nucleic acid backbones. This methylol group in turn facilitates the formation of a methylene bridge between two amino groups within the nucleic acid chain. Of note, adenosine is the preferred site for these additions. In addition, formaldehyde facilitates depurination of nucleic acids as well as hydrolysis of phosphodiester bonds which leads to shorter nucleic acid chains [23, 24]. Additional factors that could further exacerbate the impact of these modifications include the variability of FFPE specimens derived from routine clinical practice: a variable time before the tissue is subjected to fixation, a variable amount of time in fixation and a variable time embedded in paraffin wax. Indeed, delays of up to 2-3 hours prior to processing compromise mRNA abundance by ~50-60%, and mRNA abundance is significantly reduced by ~85-99% depending upon amplicon size after formalin fixation of 1 day and longer . Furthermore, capillary electrophoresis of RNA from FFPE breast tissue demonstrated that the RNA survives processing and extraction as ~300 bp fragments . Together these factors impede the ability to extract intact, biologically informative RNA from FFPE tissues. The chemical modifications would further interfere with necessary processing steps (reverse transcription into cDNA, for example). These studies explain in part why smaller nucleic acid targets, like miRNAs, are emerging as ideal substrates for the molecular characterization of FFPE tissue: their small size and close association with large protein aggregates more likely shields them to some degree from these chemical modification and degradation pathways. Furthermore, whereas conventional miRNA microarray labeling procedures employ a direct end-labeling procedure of the miRNAs , mRNA labeling protocols rely on reverse transcription. The latter would introduce further bias to the extent that these latter protocols depend upon oligo-dT labeling on a highly modified poly-A mRNA tail, since adenosine is the preferred site for metholyl addition.
There are already many well documented examples of altered miRNA expression patterns in malignant melanoma. miRNA-221 and miRNA-222 are up-regulated in melanoma cell lines and exhibit increasing expression levels along the continuum of dysplastic nevi and melanoma samples in situ. Furthermore, overexpression of miRNA-221 and miRNA-222 in melanoma cells correlates with a more rapid cell cycle progression (G1-S transition), increased invasion and migration in cell culture and enhanced tumorigenicity when these cells are injected into nude mice. At least in part, these miRNAs contribute to these phenotypes via the down-regulation of p27Kip1 and c-KIT mRNAs in the melanoma cells . Similarly, the miRNA let-7a is highly expressed in melanocytes, but exhibits reduced expression in malignant melanoma. One critical target of let-7a in melanoma appears to be β3–integrin, whose expression increases in melanoma in an inverse relationship with let-7a. Increased expression of β3–integrin, in turn, correlates with increasing metastatic potential in melanoma . Finally, array comparative genomic hybridization (aCGH) studies have shown a significant correlation (85.9%) between genomic regions demonstrating copy number alterations and miRNA loci in malignant melanoma as well as copy number alterations of the miRNA processing genes, including Dicer1 and Argonaute 2 . Together, these studies underscore the importance of a further systematic characterization of altered miRNA expression levels in the biology of pigmented lesions.
This and other studies have demonstrated that miRNAs are more vigorous analytes than mRNAs in the interrogation of FFPE tissues. Taken together with studies revealing the increased utility of miRNA expression profiles in the molecular characterization of tumors and ascribing the nature of tumors of unknown origin [2, 5-7], it is clear that miRNAs encode exceedingly important information in the archived tissues in which they are stored. Given the difficulty of procuring fresh primary melanoma specimens, FFPE is often the only resource for the study of primary melanomas. Therefore, miRNAs will likely emerge as critical resources in the further characterization of these lesions. As shown here, the ability to use FFPE primary melanoma specimens for this purpose will prove of significant value.
Support for this study was provided by a grant from the Specialized Program of Research Excellence on Skin Cancer (CA-093372).