- Research article
- Open Access
- Open Peer Review
Evaluation of a new high-dimensional miRNA profiling platform
- Julie M Cunningham†1Email author,
- Ann L Oberg†2,
- Pedro M Borralho3,
- Betsy T Kren4,
- Amy J French1,
- Liang Wang1,
- Brian M Bot2,
- Bruce W Morlan2,
- Kevin AT Silverstein5,
- Rod Staggs5,
- Yan Zeng6,
- Anne-Francoise Lamblin5,
- Christopher A Hilker1,
- Jian-Bing Fan7,
- Clifford J Steer4 and
- Stephen N Thibodeau1
© Cunningham et al; licensee BioMed Central Ltd. 2009
- Received: 12 January 2009
- Accepted: 27 August 2009
- Published: 27 August 2009
MicroRNAs (miRNAs) are a class of approximately 22 nucleotide long, widely expressed RNA molecules that play important regulatory roles in eukaryotes. To investigate miRNA function, it is essential that methods to quantify their expression levels be available.
We evaluated a new miRNA profiling platform that utilizes Illumina's existing robust DASL chemistry as the basis for the assay. Using total RNA from five colon cancer patients and four cell lines, we evaluated the reproducibility of miRNA expression levels across replicates and with varying amounts of input RNA. The beta test version was comprised of 735 miRNA targets of Illumina's miRNA profiling application.
Reproducibility between sample replicates within a plate was good (Spearman's correlation 0.91 to 0.98) as was the plate-to-plate reproducibility replicates run on different days (Spearman's correlation 0.84 to 0.98). To determine whether quality data could be obtained from a broad range of input RNA, data obtained from amounts ranging from 25 ng to 800 ng were compared to those obtained at 200 ng. No effect across the range of RNA input was observed.
These results indicate that very small amounts of starting material are sufficient to allow sensitive miRNA profiling using the Illumina miRNA high-dimensional platform. Nonlinear biases were observed between replicates, indicating the need for abundance-dependent normalization. Overall, the performance characteristics of the Illumina miRNA profiling system were excellent.
- miRNA Expression
- miRNA Target
- Average Bias
- Simple Linear Model
- miRNA Array
First identified nearly 15 years ago , microRNAs (miRNAs) are a family of short RNA molecules that predominantly inhibit gene expression at the post-transcriptional level in eukaryotes [2, 3]. In the nucleus, genes encoding primary miRNAs (pri-miRNA) are much longer than the mature form. These primary transcripts are processed by a nuclease (Drosha) and the double-stranded RNA binding protein, DGCR8, to short 60–70 nucleotide stem-loop structures (pre-miRNA). After export to the cytoplasm, pre-miRNAs are processed by interaction with the endonuclease Dicer. Mature miRNA are ~22 nucleotides in length and guide the RNA-induced silencing complex (RISC, the core components of which contain Argonaute proteins) to the target sites, usually located in the 3' untranslated region of gene transcripts. Binding of the RISC leads to suppression of translation and possibly degradation of target mRNAs.
The manner in which suppression of translation is mediated is poorly understood . However, it is becoming evident that this class of gene regulator controls the expression of 20–30% percent of all human genes . As miRNAs regulate the expression of a large number of protein-encoding genes [7–10], a wide range of biological processes are affected, such as metabolism, organogenesis, development, cell growth, cell death, and cell fate determination. Altered expression of miRNAs has also been associated with human disease, including cancer [11–18]. Indeed, some miRNAs have been coined "oncomirs", acting in the manner of oncogenes and tumor suppressors [11, 19], and oncogenes such as Myc have incorporated miRNA regulation into their tumorigenic potential .
There are now more than 541 human miRNAs identified in version 10.1 of the miRBase http://microrna.sanger.ac.uk/. For the Illumina's human miRNA BeadArray, 735 miRNA were incorporated into the beta version based on miRBase version 9.1 (February 2007 release), plus an additional 265 miRNAs derived from the literature . Other than the miRNA sample preparation, the chemistry utilized by the miRNA BeadArray is similar to that used in the DASL process (DNA Annealing, Selection and Ligation). The assay is highly multiplexed, using the universal Sentrix Array Matrices (SAM, Illumina).
To better understand the important physiological functions of miRNAs, high-throughput, miRNA microarray techniques have been employed to determine and compare global miRNA expression in different tissues and cell types and under different conditions [22–29]. These methods typically require micrograms of input RNA and often have a limited dynamic range. The Illumina miRNA profiling system has been described previously, showing the accuracy of the platform by rt-PCR and digital sequencing  using RNA from four cell lines (the same cell lines were used in this report) and purchased human tissues from Ambion. In this study, we evaluated the performance characteristics of this platform utilizing total RNA from five colon cancer cases and from four well-characterized cell lines. There is no other publication to date addressing reproducibility and performance of the Illumina miRNA profiling platform. We believe it to be highly relevant to evaluate the platform's reproducibility, and this forms the primary focus of this report. More specifically, the goals of this study were to understand variability due to plate, extraction, dilution and technical replication using cell line and clinical samples.
This study was approved by the Institutional Review Board of the Mayo Clinic (IRB 07-004516), and all patients provided informed consent. Tumor tissue was obtained from five patients with colon cancer. The samples were initially collected and stored in RLT buffer; total RNA was extracted using Trizol® LS (Invitrogen, Corp., Carlsbad, CA). For each case, frozen tumor tissue was cut on a cryostat to generate 10-micron-thick sections. The equivalent of ~2.5-square cm of tissue at 10 microns was placed into 400 μL of RLT buffer (QIAGEN, Chatsworth, CA) including 4 μL of β-mercaptoethanol. The tube was immediately flash frozen in liquid nitrogen and then stored at -80°C until utilized for RNA extraction. 1.2 ml of Trizol® LS was added to each of the vials containing the tissue sections in 400 μl of RLT, and the vials were shaken vigorously to thaw the tissue. After thawing, the 1.6 ml RLT-Trizol® LS solution was divided between three 1.5-ml tubes and homogenized for 10 sec/tube using a motor-driven, disposable RNase/DNase free pestle (Thermo Fisher Scientific, Inc., Chicago, IL). The three tubes were then combined into two, resulting in 0.8 ml RLT-Trizol tissue homogenate aliquots, and incubated at room temperature for 5 min. Following the addition of 160 μl of chloroform, the tubes were shaken vigorously for 15 sec, incubated for 5 min at room temperature, followed by centrifugation at 12,000 × g for 15 min at 4°C. From each tube, 400 μl of the aqueous phase was transferred to a fresh 1.5-ml tube, 400 μl of 100% isopropanol added, mixed by vortex and incubated at room temperature for 10 min. The tubes were then centrifuged at 7,500 × g for 10 min at 4°C and the supernatant removed by pipette and 1 ml of room-temperature 75% ethanol added. Following vortex mixing, the tubes were centrifuged at 7,500 × g for 10 min at 4°C. The supernatant was removed by pipette, and the RNA pellet air dried for 5 min prior to addition of 25 μl of RNase free water. Following resuspension, RNA was quantitated by UV spectrophotometry at 260 nm, aliquoted and flash frozen in liquid nitrogen and stored at -80°C until used. One μg was examined by agarose gel electrophoresis followed by ethidium bromide staining and visualized by UV light.
Total RNA from the four cell lines used in this study was provided by Illumina. These cell lines were the following: PC-3 (prostate adenocarcinoma), MCF-7 (breast adenocarcinoma), 293 (embryonic kidney), and HeLa (cervical adenocarcinoma).
Illumina RNA processing
The chemistry utilized by the miRNA BeadArray is similar to that used in the DASL process , and the workflow is available at http://www.illumina.com/downloads/MicroRNAAssayWorkFlow.pdf. Varying amounts of total RNA were polyadenylated and then converted to cDNA using a biotinylated oligo-dT primer with a universal PCR sequenced at its 5' end. This was followed by annealing of a miRNA-specific oligonucleotide pool (MSO) which consists of three parts: a universal PCR priming site at the 5' end, an address sequence complementary to a capture sequence on the BeadArray and a microRNA-specific sequence at the 3' end. Extension of MSO was facilitated by addition of a polymerase, but only if their 3' bases were complementary to the cognate sequence in the cDNA template. Common primers were used to amplify the cDNA templates; the primer complimentary to the BeadArray was fluorescently labelled. The single-stranded PCR product was hybridized to the Sentrix Array Matrix (SAM), where the labelled strand binds to the bead on the array containing the complementary address sequence. The SAMs were imaged using an Illumina BeadArray Reader, which measures the fluorescence intensity at each addressed bead location. Intensity files were analyzed using BeadStudio version 3.1.1. Expression levels were expressed as an average signal value.
Sample allocation on SAM 1 was designed to assess reproducibility between two separate total RNA extractions, between duplicate samples (technical replicates) and between two input RNA amounts. Sample allocation on SAM 2 was designed to assess technical reproducibility and the effect of six input levels of total RNA. Sample allocation on SAM 3 was designed to assess technical replication. The four cell lines, but none of the 5 patient RNA, were run on a fourth SAM. Reproducibility across all four SAMs was also assessed.
The goals of this study were to understand variability due to SAM, extraction, dilution and technical replication. Note that normalization by definition removes systematic variation from the data. Thus, we present an initial set of analyses using un-normalized data in order to fully understand the systematic effects of interest in this study. However, since most analyses will be performed only after normalization, we also present results on normalized data. The raw miRNA expression values were exported from the BeadStudio software without background subtraction or normalization for all 735 probes. All data are presented on the log base 2 scale, where a value of 1.0 corresponds to a 2-fold change.
A common method of evaluating agreement for microarrays is via scatter plot of all expression values for sample 1 versus all expression values for sample 2. A y = x identity line is frequently overlaid on this plot along with a loess smoother and the formula for the linear regression of sample 1 on sample 2 with corresponding correlation value. Another common method of evaluating agreement is the Bland-Altman plot , commonly called MVA (for Minus Versus Average) or MA  plot in the microarray literature. For each probe, the difference between the intensities on two arrays is plotted on the vertical axis and the average of the two intensities is plotted on the horizontal axis; there is a data point for each probe on the array. While scatter plots and correlations are useful, it is generally easier to visually detect and assess magnitude of disagreement and whether the variability in disagreement is a function of abundance in an MVA plot. If two replicates yield identical results, all points will lie on the y = 0 horizontal line (indicated on the plots for reference). A loess smoother, representing a moving average, is indicated on the plots in addition to demonstrating the average bias curve as a function of average intensity. A wholesale vertical shift away from zero indicates a global mean shift between replicates. Nonlinearities in the average bias curve indicate that the bias is a function of abundance. Axis scale is kept constant throughout all MVA plots to enable comparisons of biases between plots.
It is useful to evaluate variation about the average bias curve. As noted above, the MVA plots clearly demonstrate that variability in disagreement is a function of mean abundance. Thus, the standard deviations (SD) of the differences (the vertical axis on the MVA plots) are reported separately for the probes in the bottom 50% of the abundance distribution and the top 50% of the abundance distribution. The 50th percentile is generally around 9.5 to 10 on the horizontal axis. For comparison with other manuscripts, Spearman's correlations between pairs of arrays are displayed. The standard deviations and Spearman's correlations showed essentially no change from pre- to post-normalization; so only pre-normalization results are shown. Hierarchical clustering with a simple Euclidean distance metric was used as another measure of overall agreement for patient samples.
Normalizations were performed using fastlo , a model-based nonlinear normalization similar to cyclic loess . Fastlo can be conceptualized as a loess smooth together with a simple linear model. Envision a residual MVA plot for each sample where, for a given sample, the horizontal axis corresponds to the mean over all samples (the simplest linear model); the vertical axis corresponds to the difference of this sample from the mean of all samples (i.e., residuals from this simple linear model). A loess smooth fit to this residual MVA plot corresponds to the mean bias which is subtracted out in the algorithm for each probe, allowing the bias adjustment to be a function of mean expression as recommended in Bolstad et al  and justified by the data presented here. Like most normalization algorithms, this algorithm assumes that most probes are not differentially expressed and that there is symmetry in the number of differentially expressed probes. Ballman et al  show that this implementation is equivalent to the cyclic loess implementation in  and that the computational time required is far less. The linear model is easily tailored to allow normalization to be performed within groups if most probes are expected to be differentially expressed and/or to accommodate labeled technologies [33, 35] or other statistical experimental designs. Thus, the strengths of fastlo include an easily tailored, model-based estimation of intensity dependent biases. However, it does require knowledge of linear models in order to tailor the algorithm to specific cases. As applied here, the patient samples were normalized together since nearly all probes were expected to be similarly expressed in the five patients and between technical replicates, and no global effect due to total input RNA was observed; the cell line samples were normalized together in a separate normalization.
The detection calls generated by Illumina BeadStudio software were evaluated as part of the assessment of the effect of dilution. These are intended to be a measure of whether the signal for a probe is significantly greater than the average signal seen in the negative control probes. A cut-off of p = 0.01 was used where p < 0.01 indicates the signal was detected, as is commonly done. Venn diagrams and detection coverage plots were used to demonstrate detection overlap between dilutions and to understand differences in detection as a function of average probe abundance.
Validation of miRNA expression results
Quantitative PCR analysis of a select number of miRNA targets was performed on the four aforementioned cell lines by Illumina . Twelve probes were selected, including those with large and small differences in expression between the cell lines and those that had similar expression levels (miR-100, 125a, 125b, 135a, 146a, 150, 17-3p, 221, 26a, 31, 93 and 328). These data were produced independently from the miRNA array data and were used for validation purposes. Association of raw miRNA array expression (log2) versus the negative CT value was explored via Spearman's correlation to assess the absolute quantification capabilities of the array. All pairwise differences (log2) or fold changes (raw) between the four cell lines were also explored both within the array and qPCR data. These differences were compared in order to assess the relative quantification of the array.
Four miRNA Bead Array Plates (SAMs) were used to assess variability in patient and cell line samples due to separate total RNA extractions, technical replicates of each extraction and varying total RNA input (Figure 1).
Reproducibility between technical replicates
Reproducibility between SAMs
Reproducibility between extractions
Two separate total RNA extractions were carried out on five CRC samples. In order to compare between extractions, technical replicate 1 within an extraction, was chosen without loss of generality for presentation here. Bias curves comparing the intensities for extraction 1 versus the intensities of extraction 2 for each patient sample are displayed in panel E of Figure 3 (pre-normalization, left set; post-normalization, right set). These curves indicate relatively small average biases between extractions for most patients. Standard deviations for the differences between extractions (Figure 4, left box, panel E) in the low-abundant probes are similar to those within an extraction (panel B); values for the high-abundant data are remarkably similar to those within an extraction with similar median values. Likewise, the Spearman's correlations (Figure 4, right box, panel E) are high, ranging from 0.95 to 0.98. Post-normalization, the biases are again nearly completely removed. Clusters on pre-normalization data are not grouped by patient, extraction or dilution (100 or 200 ng) (see Additional file 7). However, post-normalization data are grouped by patient and not dilution or extraction, indicating similar expression profiles for each patient regardless of extraction or dilution (see Additional file 8). Between-extraction reproducibility is comparable to within-SAM reproducibility in these limited data.
Effect of varying input of RNA
As it is likely that many clinical samples will have limited amounts of tissue, and hence RNA, available for assessing miRNA expression profiling, we examined whether reliable results could be attained when input was as low as 25 ng. Two of the CRC samples (IDs 565 and 45) were run with one replicate of 25, 50 and 100 ng and two technical replicates of 200, 400 and 800 ng of total RNA on SAM 2. There is little difference between dilutions as evidenced by box and whisker plots; there is no global shift in intensity across input amount as would be expected if there was a dilution effect (Figure 5, right panel). Before normalization, the average bias curves of each dilution versus the 200 ng dilution in panel F of Figure 3 (left set) showed no evidence of a shift in abundance, since the curves are mostly centered about zero rather than the corresponding expected fold changes. Most curves show little nonlinearity. The distribution of standard deviations for both the low- and high-abundant data (Figure 4, left box, panel F) are shifted higher than for other systematic sources of variation indicating greater variability in disagreement. Likewise, the Spearman's correlations are shifted lower (Figure 4, right box, panel F) with a median value of 0.95 suggesting weaker agreement. After normalization, the nonlinear biases are mostly removed (Figure 3, panel F of right set). Clusters on both pre- and post-normalization data are not grouped by either patient or dilution (100 or 200 ng) (see Additional files 9 and 10). These data indicate that microRNA profiling using this chemistry can be achieved with as little as 25 ng of total RNA.
It is of interest to compare the detection rates between dilutions. The Illumina software defines a detection p-value for a probe by comparing the probe intensity to the distribution of the background probes. We set a cut-off of p = 0.01 to determine whether a probe was "significantly detected". For the two patients, 53% (386/735) and 51% (377/735) of the probes were detected in both 200 ng technical replicates, while 12% (88/735) and 5% (39/735) were detected in only one of the replicates. Venn diagrams show the overlap in detection with replicate 1 of 200 ng with all other dilutions (see Additional file A11). Overall, no systematic differences were noted in the number detected or not detected across the dilutions. Additionally, of the 461 probes detected in replicate 1 of the 200 ng dilution for one patient, 75%, 10%, 4%, 7%, 3%, and 2% were also detected in 5, 4, 3, 2, 1 or 0 of the other 5 dilutions for the one patient. Thus, 89% of the probes detected at 200 ng were detected in three or more of the other dilutions. Of the 274 probes not detected at 200 ng, 0%, 1%, 3%, 6%, 14% and 75% were detected in 5, 4, 3, 2, 1 or 0 other dilutions. Thus, of the probes not detected at 200 ng, only 4% were detected in three or more dilutions. These percentages are similar for the second patient.
Validation of miRNA data using qPCR
This report is aimed at understanding the performance characteristics of a newly developed, highly sensitive, high-dimensional miRNA expression platform, the Illumina human miRNA BeadArray. In particular, we wished to rigorously assess reproducibility of the assay and present raw, un-normalized data to better understand the performance of the assay. Overall, the Illumina miRNA array system performed very well. With respect to reproducibility (both within a plate and between plates), this assay was found to be similar to other high-dimensional gene expression platforms. The recommended input of RNA was 200 ng. Using this as the reference, however, we found no evidence for a dilution effect, with comparable results obtained from 25 ng to 800 ng RNA. Therefore, miRNA from small tissue samples can be reliably assayed. Additionally, this platform can test up to 96 samples at a time, providing high throughput capability. A description of the platform has been previously published  and included validation of 12 individual miRNA by qPCR, and on a more global scale, miRNA expression was validated using digital gene expression.
The experimental approach was to look at variability between extractions, between technical replicates (both within and between plates) and between varying inputs of RNA using RNA from five colon cancer specimens and four cell lines. The agreement between extractions appears to be similar to that between technical replicates within an extraction and was similar for all the samples tested. As one would expect, the between-plate agreement was slightly less than within-plate agreement, although most correlation values were >0.90. There were several outliers noted in one SAM. Processed at the same time as another SAM, these may be related to sample handling or perhaps to a few underperforming probes.
We present both scatter plots and MVA plots to evaluate the data; the latter are a useful tool for assessing biases in high-dimensional data. Just as with high-dimensional gene expression microarray platforms [32, 36, 34, 38, 37], the bias curves observed here demonstrate that bias is not constant across all abundance levels for high-dimensional microRNA data indicating a need for nonlinear normalization. There is debate in current literature as to whether normalization algorithms used for high density mRNA microarrays are applicable to data from microRNA microarrays [39–43]. Most commonly used mRNA normalization algorithms assume that only a small portion of probes are differentially expressed, that the distribution of differentially expressed probes is approximately symmetric about identity. In addition, there must be sufficient probes for estimation of bias without over-fitting. These assumptions must be evaluated specific to the experiment at hand. The data we present here suggest that there are sufficient numbers of probes expressed on this microRNA platform for estimation of biases without over-fitting using standard mRNA normalization algorithms in this experiment. In addition, there is no differential expression expected in the replicates analyzed herein. Thus, the assumptions hold for these data and nonlinear normalizations such as quantile [34, 44] and cyclic loess [32, 33] should be directly applicable to this experiment. Indeed, the fastlo normalization method utilized in this report removed the nonlinearity seen in the bias curves.
Finally, we selected 12 miRNA targets to further assess the accuracy of the Illumina array utilizing a second method to quantitate miRNA levels (qPCR). Using RNA from the four cells lines used in this study, the correlation between miRNA expression levels derived from the Illumina platform compared to that from the qPCR analysis was excellent for the 12 targets. These results were essentially the same as those previously reported .
Overall, the Illumina array appears to be quite specific for the mature form of the miRNA. Chen et al. hypothesize that the cDNA synthesis may be more complete for the short mature miRNAs than for the pre-miRNA templates. Also, the stem-loop structure of pre-miRNAs could hinder the cDNA synthesis and annealing of the oligonucleotides, resulting in relatively selective detection of expression of the mature miRNAs .
In summary, we found this high-dimensional miRNA profiling platform to be highly sensitive, providing reproducible data over a wide range of RNA input amounts. The variability between extraction, replicates, and SAMs was found to be acceptable. After comparison with independently produced qPCR data, validation of the absolute and relative quantification of the technology seemed adequate as well.
Illlumina's miRNA profiling application provides an excellent tool for determining miRNA expression in clinical and research samples. Small amounts of RNA may be used to generate highly reproducible data. The Illumina miRNA panel therefore, presents a robust tool for a variety of research applications, providing advantages over existing tools
We thank Malinda Woodward for excellent secretarial support and Stacy A. Johnson for technical support. The data have been deposited in the Gene Expression Omnibus (GEO ID GSE17346; http://www.ncbi.nlm.nih.gov/projects/geo/).
This work was supported by grants from the Minnesota Partnership for Biotechnology and Medical Genomics [2006-100412]; and by Mayo Clinic Comprehensive Cancer Center Shared Resource Grant [CA15083].
- 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: 843-854. 10.1016/0092-8674(93)90529-Y.View ArticlePubMedGoogle Scholar
- Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R: Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006, 20: 515-524. 10.1101/gad.1399806.View ArticlePubMedGoogle Scholar
- Zamore PD, Haley B: Ribo-gnome: the big world of small RNAs. Science. 2005, 309: 1519-1524. 10.1126/science.1111444.View ArticlePubMedGoogle Scholar
- Mattick JS, Makunin IV: Small regulatory RNAs in mammals. Hum Mol Genet. 2005, 14 (Spec No 1): R121-132. 10.1093/hmg/ddi101.View ArticlePubMedGoogle Scholar
- Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9: 102-114. 10.1038/nrg2290.View ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.View ArticlePubMedGoogle Scholar
- Berezikov E, Guryev V, Belt van de J, Wienholds E, Plasterk RH, Cuppen E: Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005, 120: 21-24. 10.1016/j.cell.2004.12.031.View ArticlePubMedGoogle Scholar
- Chan CS, Elemento O, Tavazoie S: Revealing posttranscriptional regulatory elements through network-level conservation. PLoS Comput Biol. 2005, 1: e69-10.1371/journal.pcbi.0010069.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu PW, Huang HD, Hsu SD, Lin LZ, Tsou AP, Tseng CP, Stadler PF, Washietl S, Hofacker IL: miRNAMap: genomic maps of microRNA genes and their target genes in mammalian genomes. Nucleic Acids Res. 2006, 34: D135-139. 10.1093/nar/gkj135.View ArticlePubMedGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121: 207-221. 10.1016/j.cell.2005.04.004.View ArticlePubMedGoogle Scholar
- Esquela-Kerscher A, Slack FJ: Oncomirs – microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6: 259-269. 10.1038/nrc1840.View ArticlePubMedGoogle Scholar
- Abbott AL, Alvarez-Saavedra E, Miska EA, Lau NC, Bartel DP, Horvitz HR, Ambros V: The let-7 MicroRNA family members mir-48, mir-84, and mir-241 function together to regulate developmental timing in Caenorhabditis elegans. Dev Cell. 2005, 9: 403-414. 10.1016/j.devcel.2005.07.009.View ArticlePubMedPubMed CentralGoogle Scholar
- Cui XS, Shen XH, Kim NH: Dicer1 expression in preimplantation mouse embryos: Involvement of Oct3/4 transcription at the blastocyst stage. Biochem Biophys Res Commun. 2007, 352: 231-236. 10.1016/j.bbrc.2006.11.009.View ArticlePubMedGoogle Scholar
- Abelson JF, Kwan KY, O'Roak BJ, Baek DY, Stillman AA, Morgan TM, Mathews CA, Pauls DL, Rasin MR, Gunel M, et al: Sequence variants in SLITRK1 are associated with Tourette's syndrome. Science. 2005, 310: 317-320. 10.1126/science.1116502.View ArticlePubMedGoogle Scholar
- Nakamoto M, Jin P, O'Donnell WT, Warren ST: Physiological identification of human transcripts translationally regulated by a specific microRNA. Hum Mol Genet. 2005, 14: 3813-3821. 10.1093/hmg/ddi397.View ArticlePubMedGoogle Scholar
- Schetter AJ, Leung SY, Sohn JJ, Zanetti KA, Bowman ED, Yanaihara N, Yuen ST, Chan TL, Kwong DLW, Au GKH, et al: MicroRNA Expression Profiles Associated With Prognosis and Therapeutic Outcome in Colon Adenocarcinoma. JAMA. 2008, 299: 425-436. 10.1001/jama.299.4.425.PubMedPubMed CentralGoogle Scholar
- Cheng Y, Ji R, Yue J, Yang J, Liu X, Chen H, Dean DB, Zhang C: MicroRNAs are aberrantly expressed in hypertrophic heart: do they play a role in cardiac hypertrophy?. Am J Pathol. 2007, 170: 1831-1840. 10.2353/ajpath.2007.061170.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang TH, Fan B, Rothschild MF, Hu ZL, Li K, Zhao SH: MiRFinder: an improved approach and software implementation for genome-wide fast microRNA precursor scans. BMC Bioinformatics. 2007, 8: 341-10.1186/1471-2105-8-341.View ArticlePubMedPubMed CentralGoogle Scholar
- Akao Y, Nakagawa Y, Kitade Y, Kinoshita T, Naoe T: Downregulation of microRNAs-143 and -145 in B-cell malignancies. Cancer Sci. 2007, 98: 1914-1920. 10.1111/j.1349-7006.2007.00618.x.View ArticlePubMedGoogle Scholar
- Blower PE, Chung JH, Verducci JS, Lin S, Park JK, Dai Z, Liu CG, Schmittgen TD, Reinhold WC, Croce CM, et al: MicroRNAs modulate the chemosensitivity of tumor cells. Mol Cancer Ther. 2008, 7: 1-9. 10.1158/1535-7163.MCT-07-0573.View ArticlePubMedGoogle Scholar
- Chen J, Lozach J, Garcia EW, Barnes B, Luo S, Mikoulitch I, Zhou L, Schroth G, Fan JB: Highly sensitive and specific microRNA expression profiling using BeadArray technology. Nucleic Acids Res. 2008, 36: e87-10.1093/nar/gkn387.View ArticlePubMedPubMed CentralGoogle Scholar
- Babak T, Zhang W, Morris Q, Blencowe BJ, Hughes TR: Probing microRNAs with microarrays: tissue specificity and functional inference. RNA. 2004, 10: 1813-1819. 10.1261/rna.7119904.View ArticlePubMedPubMed CentralGoogle Scholar
- Barad O, Meiri E, Avniel A, Aharonov R, Barzilai A, Bentwich I, Einav U, Gilad S, Hurban P, Karov Y, et al: MicroRNA expression detected by oligonucleotide microarrays: system establishment and expression profiling in human tissues. Genome Res. 2004, 14: 2486-2494. 10.1101/gr.2845604.View ArticlePubMedPubMed CentralGoogle Scholar
- Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA. 2005, 11: 241-247. 10.1261/rna.7240905.View ArticlePubMedPubMed CentralGoogle Scholar
- Esau C, Kang X, Peralta E, Hanson E, Marcusson EG, Ravichandran LV, Sun Y, Koo S, Perera RJ, Jain R, et al: MicroRNA-143 regulates adipocyte differentiation. J Biol Chem. 2004, 279: 52361-52365. 10.1074/jbc.C400438200.View ArticlePubMedGoogle Scholar
- Goff LA, Yang M, Bowers J, Getts RC, Padgett RW, Hart RP: Rational probe optimization and enhanced detection strategy for microRNAs using microarrays. RNA Biol. 2005, 2: 93-100.View ArticlePubMedGoogle Scholar
- Calin GA, Liu CG, Sevignani C, Ferracin M, Felli N, Dumitru CD, Shimizu M, Cimmino A, Zupo S, Dono M, et al: MicroRNA profiling reveals distinct signatures in B cell chronic lymphocytic leukemias. Proc Natl Acad Sci USA. 2004, 101: 11755-11760. 10.1073/pnas.0404432101.View ArticlePubMedPubMed CentralGoogle Scholar
- Miska EA, Alvarez-Saavedra E, Townsend M, Yoshii A, Sestan N, Rakic P, Constantine-Paton M, Horvitz HR: Microarray analysis of microRNA expression in the developing mammalian brain. Genome Biol. 2004, 5: R68-10.1186/gb-2004-5-9-r68.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu J, Carmell MA, Rivas FV, Marsden CG, Thomson JM, Song JJ, Hammond SM, Joshua-Tor L, Hannon GJ: Argonaute2 is the catalytic engine of mammalian RNAi. Science. 2004, 305: 1437-1441. 10.1126/science.1102513.View ArticlePubMedGoogle Scholar
- Steemers FJ, Gunderson KL: Illumina, Inc. Pharmacogenomics. 2005, 6: 777-782. 10.2217/146224184.108.40.2067.View ArticlePubMedGoogle Scholar
- Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1: 307-310.View ArticlePubMedGoogle Scholar
- Dudoit S, Yang Y, Callow M, Speed T: Statistical Methods for Identifying Differentially Expressed Genes in Replicated Cdna Microarray Experiments. Statistica Sinica. 2002, 12: 111-139.Google Scholar
- Ballman KV, Grill DE, Oberg AL, Therneau TM: Faster cyclic loess: normalizing RNA arrays via linear models. Bioinformatics. 2004, 20: 2778-2786. 10.1093/bioinformatics/bth327.View ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.View ArticlePubMedGoogle Scholar
- Eckel JE, Gennings C, Therneau TM, Burgoon LD, Boverhof DR, Zacharewski TR: Normalization of two-channel microarray experiments: a semiparametric approach. Bioinformatics. 2005, 21: 1078-1083. 10.1093/bioinformatics/bti105.View ArticlePubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.View ArticlePubMedPubMed CentralGoogle Scholar
- Kepler TB, Crosby L, Morgan KT: Normalization and analysis of DNA microarray data by self-consistency and local regression. Genome Biol. 2002, 3: RESEARCH0037-10.1186/gb-2002-3-7-research0037.View ArticlePubMedPubMed CentralGoogle Scholar
- Naef F, Socci ND, Magnasco M: A study of accuracy and precision in oligonucleotide arrays: extracting more signal at large concentrations. Bioinformatics. 2003, 19: 178-184. 10.1093/bioinformatics/19.2.178.View ArticlePubMedGoogle Scholar
- Hua YJ, Tu K, Tang ZY, Li YX, Xiao HS: Comparison of normalization methods with microRNA microarray. Genomics. 2008, 92: 122-128. 10.1016/j.ygeno.2008.04.002.View ArticlePubMedGoogle Scholar
- Pradervand S, Weber J, Thomas J, Bueno M, Wirapati P, Lefort K, Dotto GP, Harshman K: Impact of normalization on miRNA microarray expression profiling. RNA. 2009, 15: 493-501. 10.1261/rna.1295509.View ArticlePubMedPubMed CentralGoogle Scholar
- Rao Y, Lee Y, Jarjoura D, Ruppert AS, Liu CG, Hsu JC, Hagan JP: A comparison of normalization techniques for microRNA microarray data. Stat Appl Genet Mol Biol. 2008, 7: Article22.PubMedGoogle Scholar
- Risso D, Massa MS, Chiogna M, Romualdi C: A modified LOESS normalization applied to microRNA arrays: a comparative evaluation. Bioinformatics. 2009Google Scholar
- Sarkar D, Parkin R, Wyman S, Bendoraite A, Sather C, Delrow J, Godwin AK, Drescher C, Huber W, Gentleman R, Tewari M: Quality assessment and data analysis for microRNA expression arrays. Nucleic Acids Res. 2009, 37: e17-10.1093/nar/gkn932.View ArticlePubMedGoogle Scholar
- Astrand M: Contrast normalization of oligonucleotide arrays. J Comput Biol. 2003, 10: 95-102. 10.1089/106652703763255697.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/57/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.