Skip to main content

mRNA expression profiles show differential regulatory effects of microRNAs between estrogen receptor-positive and estrogen receptor-negative breast cancer

Abstract

Background

Recent studies have shown that the regulatory effect of microRNAs can be investigated by examining expression changes of their target genes. Given this, it is useful to define an overall metric of regulatory effect for a specific microRNA and see how this changes across different conditions.

Results

Here, we define a regulatory effect score (RE-score) to measure the inhibitory effect of a microRNA in a sample, essentially the average difference in expression of its targets versus non-targets. Then we compare the RE-scores of various microRNAs between two breast cancer subtypes: estrogen receptor positive (ER+) and negative (ER-). We applied this approach to five microarray breast cancer datasets and found that the expression of target genes of most microRNAs was more repressed in ER- than ER+; that is, microRNAs appear to have higher RE-scores in ER- breast cancer. These results are robust to the microRNA target prediction method. To interpret these findings, we analyzed the level of microRNA expression in previous studies and found that higher microRNA expression was not always accompanied by higher inhibitory effects. However, several key microRNA processing genes, especially Ago2 and Dicer, were differentially expressed between ER- and ER+ breast cancer, which may explain the different regulatory effects of microRNAs in these two breast cancer subtypes.

Conclusions

The RE-score is a promising indicator to measure microRNAs' inhibitory effects. Most microRNAs exhibit higher RE-scores in ER- than in ER+ samples, suggesting that they have stronger inhibitory effects in ER- breast cancers.

Background

MicroRNAs (miRNAs) are a class of small noncoding (19- to 24-nucleotide) RNAs that regulate the expression of target mRNAs at the post-transcriptional level [1, 2]. In higher eukaryotic organisms, it is estimated that miRNAs account for about 1% of genes and regulate the expression of more than 30% of mRNAs [3].

It has been shown that miRNAs play critical roles in a variety of biological processes such as cell proliferation [4], apoptosis [5], development [6], and differentiation [7]. In humans, strong links between cancer and miRNA deregulation have been suggested by recent studies [8, 9]. A lot of known miRNAs are found to be located in the fragile sites (regions with high frequencies of copy number alterations in cancers) of human chromosomes, indicating that many miRNAs may be linked to carcinogenesis [10]. Furthermore, it has been shown that aberrant expression of miRNAs contributes to carcinogenesis by promoting the expression of proto-oncogenes or by inhibiting the expression of tumor suppressor genes. For instance, the down-regulation of let-7, which represses expression of the proto-oncogene RAS, has been found in a large proportion of lung cancer specimens [11]. Other examples are miR-15 and miR-16, which repress the anti-apoptotic factor gene BCL2 in chronic lymphocytic leukemia [12]. In addition, some recent studies suggest that expression profiles of miRNAs are informative for the classification of human cancers. Based on miRNA-expression profiles, Lu et al. [13] reported the classification of 334 leukemia and solid cancers that agrees well with the developmental lineage and differentiation state of the tumors. Rosenfield et al. [14] demonstrated that by using miRNA as biomarkers, tumors can be classified into subclasses according to their primary origins. Nowadays, miRNAs are thought of as promising biomarkers for cancer diagnosis and prognosis.

It has been proposed that animal miRNAs regulate gene expression mainly by inhibiting translation of their target mRNAs [15, 16]. More recent studies, however, have demonstrated that expression regulation at the mRNA level (via mRNA degradation or deadenylation) also serves as a critical mechanism for miRNA function in animals [17–23]. Over-expression of miRNA in cell lines cause moderate down-regulation of a large number of transcripts, many of which contain the complementary sequences of the over-expressed miRNA in their 3' untranslated regions (UTRs) [23]. Conversely, gene expression analysis from miRNA knockdown animals reveals that miRNA recognition motifs are strongly enriched in the 3' UTRs of up-regulated genes, but depleted in the 3' UTRs of down-regulated genes[20]. Motivated by these findings, several studies have demonstrated the effectiveness of investigating miRNA regulation by examining their target mRNA expression levels [24–27]. For example, Yu et al. [27] show that miRNA targets have lower expression levels in mature mouse and Drosophila tissues than in embryos via global analysis of miRNA target gene expression.

In this study, we investigate differential miRNA regulation between estrogen receptor (ER) positive (ER+) and negative (ER-) breast cancers by examining changes in the expression of the miRNAs' target genes. Breast cancer is a common disease, ranking first in terms of annual mortality in women worldwide [28]. According to the ER status and responsiveness to estrogen, breast cancer can be divided into two subtypes: ER+ and ER-. The links between miRNA expression and breast cancer have been shown using miRNA microarray techniques [13, 29]. Specifically, the differential expression of miRNAs between ER+ and ER- breast cancers has been investigated in [30–32]. In comparison with the large number of mRNA expression datasets [33–41], miRNA expression datasets for ER+ and ER- breast cancer are still limited. Moreover, results and conclusions from these studies are generally not consistent and sometimes even conflicting [30–32]. In this study, we take advantage of those mRNA expression datasets to investigate differential miRNA regulation between ER+ and ER- breast cancers.

For each miRNA, we calculate a regulatory effect (RE)-score, which measures the expression difference between the targets and non-targets of the miRNA in an expression profile. Then, we compare the RE-scores of miRNAs in ER+ tumor samples with their RE-scores in ER- samples to identify microRNAs with changing RE-scores (which we term RE-changing microRNAs). We applied our method to five independent microarray datasets that include gene expression profiles for both ER+ and ER- samples. In all of them, our results indicate that the majority of RE-changing miRNAs showed higher RE-scores in ER- than in ER+ samples, suggesting stronger inhibitory effects of miRNAs on their targets in ER- breast cancer. To check the robustness, we performed the same analyses using different miRNA target prediction methods, RE-score calculation methods, and RE-changing miRNA identification thresholds and obtained consistent results. Moreover, we examined the expression levels of genes in the miRNA biogenesis pathway and found that Ago1 and Ago2 (which encode argonautes, the key proteins forming the RNA-induced silencing complex (RISC)) had significantly higher expression levels in ER- than in ER+ breast cancer. This may suggest higher RISC activities and, therefore, that miRNAs down-regulate target gene expression in ER- breast cancer with higher efficiency.

Results and discussion

Identification of RE-changing miRNAs between ER+ and ER-breast cancers

To measure the inhibitory effect of a miRNA, we calculate the RE-score, denoted as the difference of average ranks between the miRNA's non-target and target genes. It should be noted that the RE-scores for different miRNAs may not be directly comparable because the miRNAs regulate different sets of target genes. However, we can compare the RE-scores for the same miRNA in different conditions (that is, using different expression profiles). A higher RE-score indicates lower expression levels of target genes and, thereby, a stronger inhibitory effect of the corresponding miRNA. Given a breast cancer microarray dataset, we calculate the RE-scores for each miRNA in all samples. Then, we compare the RE-scores in ER+ and ER- samples to identify miRNAs that show different regulatory effects between these two breast cancer subtypes. We refer to these miRNAs as RE-changing miRNAs. Using ER+ as the reference, some RE-changing miRNAs show stronger inhibitory effects, while others show weaker inhibitory effects in ER- breast cancer. The false discovery rate (FDR) was estimated using a similar method to the significance analysis of microarrays (SAM) method [42]. A flow diagram of our analysis is shown in Figure 1.

Figure 1
figure 1

Schematic diagram showing the method for identifying RE-changing miRNAs between ER- and ER+ breast cancer samples. For each miRNA in each sample, a RE-score is calculated by comparing average ranks of its target and non-target genes. RE-changing miR (ER- > ER+) and RE-changing miR (ER- < ER+) represent miRNAs that have significantly higher and lower RE-scores in ER- compared to ER+ samples, respectively. RE-invariant miR represents miRNAs that show no significant difference in RE-scores between these samples. Note that many miRNAs share the same target mRNA, while many mRNAs can also be targeted by the same miRNA, which constitutes a complex miRNA-mRNA network.

Most miRNAs show stronger inhibitory effects in ER- than in ER+breast cancer

We applied our analysis to 5 carefully selected large scale microarray datasets, each containing at least 30 expression profiles for both ER+ and ER- breast cancer samples. Among these datasets, four were measured by one-channel Affymetrix GeneChips and one was measured by two-channel cDNA arrays (see Materials and methods for details about these datasets). For each dataset, we calculated the RE-scores of each miRNA in all samples. To do this, we needed to determine the target and non-target gene sets for miRNAs. Several computational methods have been developed to identify microRNA targets and predictions using these can be considerably different (Additional data file 1, the distribution of miRNA target gene numbers for different prediction tools). In our analysis, the target genes for miRNAs were predicted using the PITA algorithm, which has been shown to have high prediction accuracy [43]. Subsequently, we computed t-scores (ER- versus ER+) to measure the difference between RE-scores for ER- and ER+ samples. A positive t-score for a miRNA suggests that this miRNA has higher overall RE-scores and, thereby, stronger inhibitory effects on its targets in ER- samples. Conversely, a negative ER-/ER+ t-score indicates a stronger inhibitory effect of a miRNA in ER+ samples. For example, to estimate the RE-score of miR-371 in a sample from the HE (Hess et al. [44]) dataset, we first grouped the total 14,327 genes in the HE dataset into two sets, one with 2,054 target genes and the other with 12,273 non-target genes. Second, we sorted the expression levels of the 14,327 genes and computed the average ranks of the 2,054 targets and 12,273 non-targets, respectively. The RE-score for miR-371 in each sample was calculated as the average rank of the non-targets minus the average rank of the targets. We performed the RE-score calculation for 82 ER+ samples and 51 ER- samples and found that the RE-scores for the ER- samples are significantly higher than those for ER+ samples (t-test, P = 3.74E-15). We also compared the RE-scores for ER+ samples with those for ER- samples in the other four datasets. As shown in Figure 2a, in all of the five datasets, the RE-scores for miR-371 are significantly higher in ER- samples. Namely, miR-371 represses the expression of its target mRNAs more efficiently in ER- breast cancers. In the next section we discuss the results based on other miRNA target prediction methods.

Figure 2
figure 2

Comparison of RE-scores between ER+ and ER- samples from five breast cancer expression datasets. (a) Box plots of RE-scores for miR-371. miR-371 shows significantly higher RE-scores in ER- than in ER+ samples for all five datasets. The statistical significance level of difference (FDR) for each dataset is also shown. (b) Distributions of the t-scores for RE-score comparison between ER- and ER+ samples. The t-scores for 470 miRNAs were calculated by comparing their RE-scores for ER- samples with those for ER+ samples. The t-score distributions for the five datasets are shown in different colors. Most t-scores are positive, indicating that most miRNAs exhibit higher RE-scores in ER- than in ER+ samples. (c) Proportion of RE-changing miRNAs with higher inhibitory effect in ER- samples (red) and RE-changing miRNAs with lower inhibitory effect in ER- samples (green) at three different significance levels (FDR ≤ 0.05, FDR ≤ 0.10, and FDR ≥ 0.20). The number on the top of a bar represents how many RE-changing miRNAs were identified from the corresponding mRNA microarray dataset. HE, MI, MN, VA and WA represent the microarray data published by Hess et al. [44], Miller et al. [38], Minn et al. [39], van't Veer et al. [34], and Wang et al. [40], respectively.

We calculated the ER-/ER+ t-scores (measuring the difference between RE-scores for ER- versus ER+ samples) for 470 human miRNAs in all of the 5 datasets. Interestingly, we found that most miRNAs exhibit higher RE-scores in ER- than in ER+ samples, as suggested by the distributions of their t-scores in Figure 2b. We calculated the significance of the t-scores based on the permutation test using a similar method to SAM [42] (see Materials and methods for detail). At the 0.05 significance level (FDR ≤ 0.05), we identified 109, 188, 15 and 306 RE-changing miRNAs from a total of 475 miRNAs in the HE (Hess et al. [44]), MI (Miller et al. [38]), MN (Minn et al. [39]) and VA (van't Veer et al. [34]) datasets, respectively, and all of them show higher inhibitory effects in ER- breast cancer. In the WA (Wang et al. [40]) dataset, we identified 377 RE-changing miRNAs, of which 373 have higher inhibitory effects and only 4 lower inhibitory effects in ER- breast cancer. This suggests that most miRNAs exhibit stronger inhibitory effects on the expression of their targets in ER- compared to ER+ breast cancer. This conclusion could still be made when we relaxed the FDR threshold to 10% and 20%, as illustrated in Figure 2c. The t-score, P-value and FDR of each miRNA for all datasets are provided in Additional data file 2.

Use of other miRNA target prediction algorithms

Next, we investigated whether similar results can be obtained using other miRNA target prediction methods. It has been shown previously that distinct miRNA prediction methods may result in considerably different target gene sets (Additional data file 1, the distribution of miRNA target numbers for different prediction tools). To rule out the possible bias introduced by PITA, we repeated our analysis using three other miRNA target prediction methods: TargetScan [3], PicTar [45] and miRanda [46]. We chose these three out of a handful of miRNA target prediction methods not only because they have been prevalently used but also because they are, in some sense, complementary to the PITA method. Almost all miRNA target prediction methods first scan the 3' UTR of transcripts for potential miRNA binding sites that are complementary to the seed region of miRNAs. TargetScan and PicTar meet stringent seed pairing criteria, whereas the criteria are moderately stringent in PITA and miRanda. To further increase the prediction accuracy, PITA takes into account the local accessibility of the potential binding sites, whereas miRanda and PicTar apply a different strategy: they filter out those miRNA binding sites in non-conserved regions. TargetScan, the most widely used prediction method, considers both site conservation and context accessibility.

The results based on PicTar and miRanda are illustrated in Figure 3a, b. As shown, the t-scores for RE-score comparisons for ER- versus ER+ samples are more likely to be positive values for all five datasets, suggesting that miRNAs have stronger inhibitory effects on their targets in ER- breast cancer. Since both PITA and miRanda can require moderately stringent miRNA seed:target complementarity, in order to obtain more reliable target and non-target gene sets for miRNAs, we also tried another strategy: combining the prediction results of PITA and miRanda methods. For each miRNA, we define its target genes as those predicted by both methods and its non-target genes as those predicted by neither. This will presumably decrease both false positive and false negative prediction rates. Based on this target and non-target gene set definition, we again obtain similar results, as shown in Figure 3c.

Figure 3
figure 3

Distributions of the t-scores for comparison of RE-scores based on distinct miRNA target prediction algorithms. (a) PicTar algorithm. (b) miRanda algorithm. (c) Intersection of miRanda and PITA. (d-f) TargetScan algorithm where the site is conserved (d), the site context score is above -0.20 (e), and all potential targets are included (f). The t-score distributions for the five datasets are shown in different colors. The t-scores are more likely to be positive values in all five datasets, suggesting that miRNAs have stronger inhibitory effects on their targets in ER- breast cancer.

TargetScan is currently the most widely used microRNA target prediction tool, which relies on strict miRNA seed region complementarity [3, 47]. In addition, the conservation of binding site, the context of the miRNA-binding site, the proximal AU composition, and proximity to sites for co-clustered miRNAs can enhance the targeting efficacy of a binding site [48]. Choosing different parameters for target prediction results in quite different performance [49]. Among the parameters, site conservation and site accessibility (measured as context score) are the two most important [50, 51]. To evaluate the performance of different TargetScan cutoffs in the RE-score comparison, we chose three target sets - one in which the members have a conserved binding site, one in which the members have a context score greater than -0.20, and one that includes all potential targets - which we refer to as ConservedTS, ContextTS, and AllTS, respectively. These three TargetScan predictions are quite different. On average, 210, 765, and 2,026 targets per miRNA are predicted in ConservedTS, ContextTS and AllTS, respectively. After integrating mRNA expression data with all three target sets to compare the RE-scores in ER- and ER+ samples, we found again that the t-scores for ER- versus ER+ samples are more likely to be positive for all five datasets, as illustrated in Figure 3d-f. This demonstrates that the observation of higher RE-scores in ER- breast cancer, for most miRNAs, is not likely caused by a bias from the miRNA prediction method. Complete results based on three TargetScan predictions, miRanda, PicTar, and the intersection of PITA and miRanda can be found in Additional data file 2.

Use of alternative methods to compare miRNA inhibitory effects

To further substantiate our findings, we also used two alternative methods to investigate the inhibitory effects of miRNAs in ER+ and ER- breast cancers. The first method is similar to the one described above, but we use a different way to calculate the RE-scores for miRNAs in an expression profile. Instead of computing the average rank difference between the target and non-target gene sets for a miRNA, we calculate the RE-score as follows: first, calculate the relative expression levels of each gene across all of the samples by subtracting the mean and then dividing by the standard deviation; second, calculate the RE-score of a miRNA by comparing the relative expression levels of its target and non-target genes. For clarity, we will refer to these two RE-score calculation methods as rank comparison and expression comparison. Similar to what we found using the rank comparison method, RE-scores obtained using expression comparison tend to be higher in ER- samples as indicated by the t-score (ER- versus ER+) distribution (Figure 4a). These results are not dependent on the miRNA target prediction method because similar results are obtained using PITA and miRanda (complete results are given in Additional data file 3). As a matter of fact, the t-scores obtained by using expression comparison and rank comparison are highly correlated. For example, for the VA dataset, these two methods yield two sets of t-scores with a correlation coefficient of 0.928 (Figure 4b). As shown, 432 out of 466 miRNAs have positive t-scores from both methods, confirming stronger inhibitory effects of miRNAs in ER- breast cancer.

Figure 4
figure 4

Results obtained from an alternative RE-score calculation method based on expression comparison. (a) Distributions of the t-scores calculated by comparing the RE-scores from the expression comparison method. The employed target prediction algorithm was PITA. The t-score distributions for the five datasets are shown in different colors. (b) Correlation between the t-scores obtained from the two different RE-score calculation methods. The microarray dataset used was VA, published by van't Veer et al. [34]. The correlation coefficient (R) is 0.928, indicating that the t-scores obtained using expression comparison and rank comparison are highly correlated.

The other method, referred to as ARR (adapted ranked ratio), is similar to the ranked ratio (RR) method proposed by Yu et al. [27]. First, the expression levels of each gene in ER+ and ER- samples were compared and a t-score (ER+/ER-) was calculated to measure the expression differentiation of the gene in the two breast cancer subtypes. The t-scores for all genes were then ranked and genes were divided into two groups, those with high t-scores and those with low t-scores. For each miRNA, the ARR value was calculated by dividing the number of target genes in the 'low' ranked group by the 'high' ranked group. The ARR value is an indicator of the distribution of a miRNA's targets within all genes. A low ARR value (ARR < 1) indicates that a miRNA has more targets in genes with higher t-scores, that is, genes that are lowly expressed in ER- samples; the target genes of this miRNA tend to have lower expression levels in ER- breast cancer. We calculated the ARR values for all miRNAs in each of these five datasets. The numbers of miRNAs with ARR < 1 and ARR > 1 are listed in Table 1. As shown, more miRNAs have ARR < 1 in all datasets, indicating their stronger inhibitory effect in ER- breast cancer.

Table 1 Number of miRNAs with ARR < 1 and ARR > 1 in each dataset

Although the ARR method is similar to the RR method described by Yu et al. [27], they differ in some ways. The RR value for a miRNA in a tissue is calculated by dividing the number of targeted genes with 'low' expression by the number of target genes with 'high' expression after the expression levels of each gene across a series of tissues are ranked and split into 'high' and 'low' groups, with half the genes in each. In our ARR method, we first performed a t-test to compare the expression levels of each gene in ER+ and ER- samples. The t-scores were ranked and genes were divided into two groups corresponding to high ranked and low ranked genes, each containing half the genes. The ARR value of each miRNA was then calculated by dividing the number of targets with high rank by the number of targets with low rank. Compared with the RR method described by Yu et al. [27], our method is different in three aspects. First, for each gene, the expression levels were compared between ER+ and ER- samples. To reveal the expression difference between two groups, the t-score is more effective than the ranks across all samples. Second, the ARR value from our method is actually an indicator of difference between the expression distribution of a microRNA's target genes and that of all genes. Therefore, it directly reflects the regulatory effect of a microRNA on its target genes. Third, for a microRNA, only one ARR value is obtained based on the whole dataset with our method, and the ARR value facilitates a global inspection of the inhibitory activity differences of a microRNA between two sample groups.

Although the calculations of RE-score and ARR value are completely different, the results from each are highly consistent. We compared the RE-scores determined by expression comparison methods with the ARR results. First, we computed the Spearman correlation of the RE-scores and the ARR values for each microarray dataset. As illustrated in Table 2, the inhibitory activities calculated by these two different methods are highly correlated, with the correlation coefficients ranging from 0.578 to 0.861, which provides further confirmation that more microRNAs show higher inhibitory effects in ER- breast cancers. Second, we overlapped the microRNAs with higher or lower inhibitory activity in ER- cancers predicted by the RE-score and ARR values (Table 3). If a microRNA has a t-score (ER-/ER+) > 0 in the RE-score comparison and ARR < 1 in the ARR calculation, it is predicted to have higher inhibitory activity in ER- cancer by both methods, whereas a microRNA with a t-score < 0 and RR > 1 shows consistently higher activity in ER+ cancer. More than 80% of the miRNAs overlap, indicating that these two methods are in strong agreement. Furthermore, the number of miRNAs with consistently higher activity in ER- samples is much higher than the number with consistently lower activity in ER+ samples, again indicating that most miRNAs exhibit higher regulatory effects in ER- than in ER+ samples. Some significant miRNAs are identified by both methods. For example, it has been reported that miR-206, which regulates the estrogen receptor, has higher activity in ER- than ER+ cancers [52]. In our calculations, for all five microarray datasets, the ARR values of this microRNA are all <1, and the t-scores for RE-score comparison between ER- and ER+ cancers are all >0 (Table 3). These results are consistent with the activity difference between ER+ and ER- cancer reported by Adams et al. [52].

Table 2 Correlation between the results obtained using the ARR and RE-score calculation methods
Table 3 Regulatory activity of miR-206 predicted by the RE-score and ARR methods

Differential regulatory effects of miRNAs can not be explained by miRNA expression differences between ER+ and ER-cancer

To understand why miRNAs tend to have stronger inhibitory effects on their targets in ER- samples, we asked whether they are more highly expressed in ER- breast cancers. Using miRNA microarray technology, expression levels of miRNAs have been previously measured and compared in ER- and ER+ samples in three different studies [30–32]. Iorio et al. [31] identified 11 miRNAs that were differentially expressed between ER+ and ER- samples, of which 8 were down-regulated in the ER- samples. In contrast, many more miRNAs were reported to be differentially expressed by Blenkiron et al. [30] and Mattie et al. [32]. Specifically, Blenkiron et al. identified 35 differentially expressed miRNAs, of which 11 were up-regulated and 24 were down-regulated in the ER- samples. Mattie et al., however, reported that the majority of differentially expressed miRNAs were down-regulated in ER- samples (40 out of 43). These three miRNA expression studies do not support the idea that miRNAs tend to be more highly expressed in ER- than ER+ breast cancer. It should be noted that the three studies obtained substantially different results due to the technological issues of miRNA microarray experiments.

In addition, to measure the correlation between miRNAs' inhibitory effects and their expression levels, we calculated the Spearman correlations of the t-scores for the miRNA expression comparisons and those for the miRNA RE-score comparisons. As illustrated in Table 4, there is only a very weak positive correlation between them; particularly, the miRNA expression data published by Mattie et al. [32] shows almost no correlation with the miRNA regulatory effects predicted from all five mRNA expression datasets. This further indicates that the stronger inhibitory effect of miRNAs in ER- cancer cannot be explained by their expression levels.

Table 4 Correlation between microRNA RE-scores and their expression levels

Some microRNAs have large inconsistencies between their expression levels and RE-scores. For example, many studies have suggested that the expression levels of Dicer, the key gene in the generation of microRNAs, vary in different cancer subtypes [53–55]. In our study, Dicer is significantly down-regulated in ER- compared to ER+ cancers (see next section for details). A possible mechanism for this is that it is regulated epigenetically [56]. Six microRNAs, miR-103, miR-122a, miR-130a, miR-148a, miR-19a, and miR-29a, are commonly predicted to target Dicer by the prediction methods PITA, miRanda, PicTar and Targetscan. We investigated the expression levels of these microRNAs in two distinct datasets published by Blenkiron et al. [30] and Mattie et al. [32]. The expression levels of these microRNAs are mostly lower in ER- samples (Figure 5), which is opposite to our inference that they may be up-regulated to transcriptionally repress Dicer in ER- cancer. We then compared the RE-scores of these microRNAs in ER+ and ER- cancers. To our surprise, almost all microRNAs show stronger inhibitory effects in ER- cancers (Figure 5), which may explain why Dicer is expressed less in ER- cancer. Especially, miR-122a, which was reported to target Dicer and function in various cellular stresses [57, 58], is expressed at significantly lower levels but shows significantly higher inhibitory activity in ER- cancer, strongly indicating that the differential regulatory effects of miRNAs can not be explained by miRNA expression differences between ER+ and ER- cancer.

Figure 5
figure 5

The expressions and RE-scores of microRNAs predicted to target Dicer. BL and MA represent the microRNA microarray data published by Blenkiron et al. [30] and Mattie et al. [32]. HE, MI, MN, VA and WA represent the mRNA microarray data published by Hess et al. [44], Miller et al. [38], Minn et al. [39], van't Veer et al. [34], and Wang et al. [40], which were used to calculate the microRNA RE-scores. If the difference between ER+ and ER- samples is significant, the plot is flagged with three asterixes. The expression levels of these six microRNAs are mostly lower in ER- samples; however, almost all the RE-scores in ER- samples are higher, suggesting that the differential regulatory effects of miRNAs can not be explained by miRNA expression difference between ER+ and ER- cancers.

Several studies have reported that good classification of cancer subtypes can be achieved using the expression levels of miRNAs [13, 14]. Because striking differences in the RE-scores for a set of miRNAs between ER+ and ER- samples are observed, the RE-score of an miRNA could be a promising predictor for breast cancer subtype classification. We used the RE-scores of the top eight significantly RE-changing miRNAs in the MN dataset [39] to classify the ER+ and ER- subtypes. As expected, the accuracy was up to 89.29%. The RE-score profiles of these miRNAs are plotted in Figure 6. The classification accuracy was comparable or even better (85.76%) when estimated using the expression levels of the top 35 differentially expressed miRNAs in the dataset published by Blenkiron et al. [30], suggesting that the prediction of ER status of breast cancer based on miRNA regulatory effect or miRNA targeted mRNA expression is an alternative to that based on miRNA expression.

Figure 6
figure 6

RE-score profiles of microRNAs for the classification of ER+ and ER- breast tumors. The figure demonstrates unsupervised hierarchical clustering of 57 ER+ and 42 ER- samples in the MN dataset [39] using the top 8 RE-changing miRNAs. A dendrogram of the tumors is shown at the top, with ER+ samples in red and ER- samples in yellow. For hierarchical clustering, RE-scores of each miRNA were mean centered and normalized, and tumors were clustered using Pearson correlation (uncentered) and average linkage (CLUSTER and TREEVIEW software) [73].

Differential expression of miRNA processing genes between ER+ and ER-breast cancers

In addition to miRNA abundance, post-transcriptional regulation of miRNA expression may also be important for the inhibitory effect of miRNAs on their targets. Deregulation of genes required for miRNA biogenesis may be expected to lead to global changes in miRNA expression as well as the inhibitory effects of miRNAs. Therefore, we examined whether miRNA processing genes are differentially expressed in ER+ and ER- breast cancers.

We found that among the miRNA processing genes, Ago1 and Ago2 were significantly up-regulated in ER- compared to ER+ samples in all datasets, with combined P-values of 4.0E-8 and 2.0E-10, respectively, whereas Dicer and TRBP were significantly down-regulated, with combined P-values of 8.8E-6 and 2.9E-10, respectively (Figure 7b; Additional data file 4). Differential expression of Ago1, Ago2 and Dicer between ER+ and ER- breast cancer has been previously investigated and consistent results were reported by Blenkiron et al. [30]. As shown in Figure 7a, several proteins play a critical role in the miRNA processing pathway. DROSHA, a double-stranded RNA-specific ribonuclease, digests the pri-miRNA in the nuclease to release hairpin, precursor miRNA (pre-miRNA) [7]; then DICER, a member of the RNase III nuclease, cleaves the pre-miRNA into a single-stranded mature miRNA with the assistance of TRBP [59]; finally, the mature miRNA is incorporated into RISC consisting of DICER, TRBP, AGO and several other proteins [60–62]. Among the eight human AGO proteins, AGO1 and AGO2 are known to play the most important roles in transcriptional silencing mediated by miRNAs or small interfering RNAs. Assembly of human RISC minimally requires AGO2, DICER, and TRBP, among which AGO2 is the catalytic engine owing to its endonuclease activity and the DICER-TRBP complex acts simply as a platform [60, 63, 64]. The relatively lower abundance of AGO1 and AGO2 proteins in ER+ breast cancer may limit the activity of functional RISC, which would in turn lower the inhibitory effect of miRNAs on their targets. Moreover, since the expression levels of Dicer and Ago genes are anti-correlated in ER+ and ER- cancer, there is no necessary link between the mature miRNA expression levels and RISC activity. This may also explain the global up-regulation of miRNA expression levels in ER+ cancer observed by Blenkiron et al. since Dicer is significantly up-regulated [30].

Figure 7
figure 7

Differential expression of miRNA processing genes between ER+ and ER- breast cancer samples. (a) miRNA biogenesis and function pathway. Genes significantly up-regulated and down-regulated in ER- compared to ER+ samples are shown in yellow and cyan, respectively. A gene - for example, Drosha - is marked grey to denote that it shows no significant differential expression between ER- and ER+ samples. (b) Expression levels of Ago1, Ago2, Dicer, and TRBP in ER+ (red) and ER- (green) samples. The mean and the standard deviation of the expression levels for each gene are shown as a bar and vertical line, respectively. Data for a gene are not shown if it is missing from a dataset. The combined P-value for each gene is also shown.

It seems that the key genes in the microRNA biogenesis pathway are subjected to delicate regulation and their differential expression is likely to be associated with distinct tumor subtypes. More interestingly, genes in this pathway are not consistently regulated: Dicer and TRBP, which are involved in miRNA maturation and RISC assembly, are down-regulated whereas the catalytic engine of RISC is up-regulated in ER- relative to ER+ breast cancer. As a result, the capability of miRNAs (or more precisely RISC) to repress their targets may not be reflected by their expression levels. Using microRNA microarray experiments, Blenkiron et al. [30] found that the most differentially expressed miRNAs between ER+ and ER- cancers are down-regulated in the latter. They also examined the correlation between miRNA expression and changes in the mRNA levels of their direct targets but failed to detect enrichment for down- or up-regulation of predicted target miRNAs consistent with miRNA expression differentiation in most cases. This can be explained by the hypothesis that many miRNAs act at the level of translation rather than mRNA stability; nevertheless, this can also be explained by discordance in changes of expression between the key miRNA processing genes. Our results demonstrate that miRNAs tend to have stronger inhibitory effect on their mRNA targets in ER- breast cancer, suggesting that the AGO proteins (up-regulated in ER- cancer at the mRNA level), the catalytic engine of RISC, may eventually determine the efficiency of miRNAs to down-regulate their targets. In addition, deregulation of the key genes in the miRNA biogenesis pathway may be related to tumorigenesis of certain cancer types, as has been suggested by the fact that down-regulating DICER expression promoted tumorigenesis in vitro and in a mouse lung cancer model [65].

It has been reported that Ago2 is expressed more in ER- than in ER+ breast cancer cell lines, and that this is correlated with active ER signaling [66]. AGO2 enhances cell proliferation, reduces cell-cell adhesion, and increases cell migratory ability, which contribute to the tumor phenotype transformation from ER+ to ER- through overexpression of Ago2. Either gene amplification or activation of cell signaling cascades elevates Ago2 expression in ER- cancer cells. Up to now, no clear evidence or experimental data have shown that Ago2 is amplified in ER- cancer. The epidermal growth factor receptor (EGFR) and mitogen activated protein kinase (MAPK) signaling cascades are the major signal transduction pathways in ER- breast cancers [67, 68]. One of the frequent and remarkable features in ER- cancer is the up-regulated EGFR gene [69]. Adams et al. [66] proposed and confirmed that epidermal growth factor stimulated Ago2 expression in ER- cancers and that this was primarily regulated by the MAPK pathway. In addition, with the overexpression of Ago2, the inhibition activity of miR-206 was elevated, whereas without Ago2 the activity of miR-206 remained unchanged even with the overexpression of miR-206, suggesting that formation of Ago2-miRNA complexes is the main factor influencing miR-206 inhibitory activity [70]. This is consistent with our finding that the activity of a microRNA cannot be explained merely by its expression level. Based on this suggestion, a hypothesis can be provided that, with elevated Ago2 expression, an miRNA's inhibitory activity accordingly increases, which leads to low expression levels of genes involved in ER+ cell types and the predominant expression of genes involved in the oncogenic pathways leading to ER- cancer.

Despite growing evidence that Dicer mRNA levels vary between different tumor subtypes and that these variations are correlated with cancer progression [53–55], the regulation of Dicer remains unclear. Weisen et al. [56] reported that type I interferon represses Dicer. As already reported, MAPK signaling pathways comprise a major cascade in ER- cancers [68]. Type I interferon signals can be transduced by the MAPK pathway [71], and the activated MAPK pathway in ER- cancers may enhance the signal of type I interferon, which results in the inhibition of Dicer expression. Another possible explanation of the low expression of Dicer in ER- cancers may be the regulatory effect of miRNAs. DICER's epigenetic regulation could also occur via specific mechanisms involving the DICER 3' UTR and the binding of microRNAs [56]. In this study, we have shown that the activity of miRNAs is stronger in ER- than ER+ cancer and that Dicer is targeted and suppressed to a lower level in ER- compared to ER+ cancers.

Conclusions

In this study, we created the RE-score to measure the inhibitory effect of a miRNA on its targets. Based on RE-score calculations, we compared the inhibitory effects of miRNAs on their targets between two breast cancer subtypes, ER+ and ER-. miRNAs that showed significantly different inhibitory effects were identified for five independent datasets. We found that, for most miRNAs, the target genes were more repressed in ER- than ER+ breast cancer, suggesting that miRNAs have stronger inhibitory abilities in the former. The exact identity of the miRNA targets does not seem important since these findings are robust to several distinct methods of miRNA target prediction and are further consolidated by another two methods for comparing miRNA regulation. To seek the potential mechanisms contributing to the inhibitory effects of miRNAs, we explored miRNA abundance measured by miRNA microarrays and expression levels of genes involved in miRNA biogenesis and function. Our analysis indicates that a high inhibitory ability is not necessarily associated with high miRNA expression levels, because previous miRNA expression data do not suggest prevalent over-expression in ER- breast cancer. However, it is interesting to find that several key miRNA processing genes are significantly differentially expressed between ER+ and ER- breast cancer. Ago1 and Ago2 are significantly up-regulated in ER- cancer, while Dicer and TRBP are significantly down-regulated. These results imply that the miRNA processing pathway is subject to subtle regulation and that deregulation of key genes in this is involved in the cancer pathology. This method is easily applied and can be used to investigate the miRNA regulation underlying other microarray datasets.

Materials and methods

Breast cancer microarray datasets

All the microarray data used in this study were downloaded from public databases or from the websites provided by the original publications. Over ten breast cancer datasets have been generated in previous studies [33–41, 44]. From these datasets, we chose five according to the following criteria: contain at least 30 samples for both ER+ and ER- breast cancer; and expression of ER+ and ER- samples is measured using the same platform. The first criterion is to ensure a high power of statistical analysis, while the second criterion is to avoid bias introduced by platform effect. Among these five datasets, one used cDNA arrays and the other four used oligonucleotide arrays produced by Affymetrix. Numbers of ER+ and ER- samples in each dataset are listed in Table 5. The expression values are represented by normalized log ratios for cDNA microarrays or by log-transformed intensities after Robust Multichip Average normalization for Affymetrix oligonucleotide microarrays [72]. The probe or probeset IDs are mapped to NCBI Refseq IDs. When multiple probe sets are mapped to the same Refseq ID, their values are averaged to represent the expression level of this Refseq gene.

Table 5 Breast cancer gene expression datasets used in this study

miRNA target predictions

A number of miRNA target prediction approaches have been suggested in the past few years [43, 45–47]. In this paper, we utilize four sets of miRNA target prediction data derived using PITA [43], miRanda [46], PicTar [45] and TargetScan [3, 47], respectively. All miRNA target prediction datasets were downloaded from the most recently updated websites. To facilitate the analysis, the target gene IDs were also converted into NCBI Refseq IDs. For each miRNA, the target genes are defined as those presented in the microarray data and predicted to contain at least one binding site at their 3' UTR; the non-target genes are defined as those presented in the microarray data but not predicted to be regulated by the miRNA.

Measuring a miRNA's inhibitory effect with the average rank difference between its targets and non-targets

To measure the inhibitory effect for a miRNA, we defined the RE-score, which measures the difference in expression levels between its target and non-target genes. The RE-score can be calculated in two ways: one is based on rank comparison and the other is based on expression comparison.

The RE-score based on rank comparison is calculated as follows. We denote the number of a miRNA's targets and non-targets as N t and N n , respectively. After sorting the expression levels of all genes, the ranks of target genes and a non-target genes are denoted as R t and R n , respectively. The RE-score of a miRNA is defined as the difference of the average rank between its targets and non-targets:

where and represent the mean target and non-target ranks, respectively. The RE-score is essentially a transformation of the sum rank statistic (the sum of ranks for target genes) used in the Wilcoxon rank sum test. Since genes with high absolute expression values have high rank values, a positive RE-score indicates that the non-target genes of a miRNA tend to be expressed at higher levels than its target genes, presumably due to the inhibitory effect of the miRNA on its target genes. The higher the RE-score, the stronger the inhibitory effect of a miRNA on its targets.

Identifying microRNAs with significantly changed RE-scores between ER+ and ER-breast cancer

To investigate the difference of a miRNA regulatory effect between ER+ and ER- breast cancer, a two sample t-test was performed to compare RE-scores and determine whether the RE-scores of a miRNA are significantly different between ER+ and ER- cancer.

Since usually hundreds of miRNAs are examined simultaneously, multiple testing corrections needed to be considered. We calculated the FDR based on permutations similar to the method used in SAM [42]. If there were N1 ER+ and N2 ER- samples, the t-scores obtained from comparing RE-scores for each miRNA were calculated in the original data, denoted as TRES(r) for the rth miRNA. We then permutated the samples; at each permutation, N1 samples were randomly selected to form one permuted ER+ group, and the rest of the samples were used as the permutated ER- group. The permutated t-score, denoted as TRES(r, k), for the rth miRNA in the kth permutation, is recalculated. We then considered the histogram of all TRES(r, k) over all r and k, and used this null distribution to compute an FDR value for a given t-score TRES(r) = TRES*. If TRES* ≥ 0, the FDR is the ratio of the percentage of all (r, k) with TRES(r, k) ≥ 0, whose TRES(r, k) ≥ TRES*, divided by the percentage of miRNAs with TRES(r) ≥ 0, where TRES(r) ≥ TRES*, and similarly if TRES* < 0.

If the FDR for a miRNA is below a predefined threshold, we call this miRNA as a significantly RE-changing miRNA.

ER+ and ER-cancer subtype classification

The miRNA RE-score is a promising feature to classify tumor subtypes as well as microRNA expression. In this study, we constructed a multi-miRNA signature and used the algorithm of a linear support vector machine. We ranked the P-values that were derived from the comparison of RE-scores or expressions. The top N significant RE-score changing microRNAs or differentially expressed microRNAs were chosen to perform the classification analysis. To estimate the effect of the classifier, we adapted a leave one out cross validation strategy. Generally, the number of ER+ samples is different from the number of ER- samples. If there were N1 ER+ and N2 ER- samples, assuming that N1 > N2, in order to balance the sample effect, we randomly selected N2 ER+ samples. The total 2* N2 (N2 ER+ and N2 ER-) samples were used in the leave one out validation. The classification accuracy was determined by averaging the accuracies of the leave one out validations repeated 100 times.

Additional data files

The following additional data are available with the online version of this paper: a figure showing distributions of miRNA target numbers determined using different prediction tools (Additional data file 1); a table listing RE-score results for five breast cancer expression datasets (Additional data file 2); a table listing RE-score results calculated using the expression comparison method (Additional data file 3); a table listing expression levels of several miRNA processing genes in ER+ and ER- samples (Additional data file 4).

Abbreviations

ARR:

adapted ranked ratio

ER:

estrogen receptor

FDR:

false discovery rate

MAPK:

mitogen activated protein kinase

miRNA:

microRNA

RE:

regulatory effect

RISC:

RNA-induced silencing complex

RR:

ranked ratio

SAM:

significance analysis of microarrays

UTR:

untranslated region.

References

  1. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.

    Article  PubMed  CAS  Google Scholar 

  2. Ambros V: microRNAs: tiny regulators with great potential. Cell. 2001, 107: 823-826. 10.1016/S0092-8674(01)00616-X.

    Article  PubMed  CAS  Google Scholar 

  3. 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.

    Article  PubMed  CAS  Google Scholar 

  4. Cheng AM, Byrom MW, Shelton J, Ford LP: Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis. Nucleic Acids Res. 2005, 33: 1290-1297. 10.1093/nar/gki200.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Xu P, Guo M, Hay BA: MicroRNAs and the regulation of cell death. Trends Genet. 2004, 20: 617-624. 10.1016/j.tig.2004.09.010.

    Article  PubMed  CAS  Google Scholar 

  6. Karp X, Ambros V: Developmental biology. Encountering microRNAs in cell fate signaling. Science. 2005, 310: 1288-1289. 10.1126/science.1121566.

    Article  PubMed  CAS  Google Scholar 

  7. Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303: 83-86. 10.1126/science.1091903.

    Article  PubMed  CAS  Google Scholar 

  8. Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6: 857-866. 10.1038/nrc1997.

    Article  PubMed  CAS  Google Scholar 

  9. Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6: 259-269. 10.1038/nrc1840.

    Article  PubMed  CAS  Google Scholar 

  10. Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, Croce CM: Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci USA. 2004, 101: 2999-3004. 10.1073/pnas.0307323101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014.

    Article  PubMed  CAS  Google Scholar 

  12. Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci USA. 2005, 102: 13944-13949. 10.1073/pnas.0506654102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.

    Article  PubMed  CAS  Google Scholar 

  14. Rosenfeld N, Aharonov R, Meiri E, Rosenwald S, Spector Y, Zepeniuk M, Benjamin H, Shabes N, Tabak S, Levy A, Lebanony D, Goren Y, Silberschein E, Targan N, Ben-Ari A, Gilad S, Sion-Vardy N, Tobar A, Feinmesser M, Kharenko O, Nativ O, Nass D, Perelman M, Yosepovich A, Shalmon B, Polak-Charcon S, Fridman E, Avniel A, Bentwich I, Bentwich Z, et al: MicroRNAs accurately identify cancer tissue origin. Nat Biotechnol. 2008, 26: 462-469. 10.1038/nbt1392.

    Article  PubMed  CAS  Google Scholar 

  15. Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403: 901-906. 10.1038/35002607.

    Article  PubMed  CAS  Google Scholar 

  16. Pillai RS, Bhattacharyya SN, Filipowicz W: Repression of protein synthesis by miRNAs: how many mechanisms?. Trends Cell Biol. 2007, 17: 118-126. 10.1016/j.tcb.2006.12.007.

    Article  PubMed  CAS  Google Scholar 

  17. Farh KK, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP: The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 2005, 310: 1817-1821. 10.1126/science.1121158.

    Article  PubMed  CAS  Google Scholar 

  18. Stark A, Brennecke J, Bushati N, Russell RB, Cohen SM: Animal MicroRNAs confer robustness to gene expression and have a significant impact on 3'UTR evolution. Cell. 2005, 123: 1133-1146. 10.1016/j.cell.2005.11.023.

    Article  PubMed  CAS  Google Scholar 

  19. Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434: 338-345. 10.1038/nature03441.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Krutzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, Stoffel M: Silencing of microRNAs in vivo with 'antagomirs'. Nature. 2005, 438: 685-689. 10.1038/nature04303.

    Article  PubMed  Google Scholar 

  21. Bagga S, Bracht J, Hunter S, Massirer K, Holtz J, Eachus R, Pasquinelli AE: Regulation by let-7 and lin-4 miRNAs results in target mRNA degradation. Cell. 2005, 122: 553-563. 10.1016/j.cell.2005.07.031.

    Article  PubMed  CAS  Google Scholar 

  22. Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci USA. 2006, 103: 4034-4039. 10.1073/pnas.0510928103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433: 769-773. 10.1038/nature03315.

    Article  PubMed  CAS  Google Scholar 

  24. Wang X: Systematic identification of microRNA functions by combining target prediction and expression profiling. Nucleic Acids Res. 2006, 34: 1646-1652. 10.1093/nar/gkl068.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Arora A, Simpson DA: Individual mRNA expression profiles reveal the effects of specific microRNAs. Genome Biol. 2008, 9: R82-10.1186/gb-2008-9-5-r82.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Cheng C, Li LM: Inferring microRNA activities by combining gene expression with microRNA target prediction. PLoS ONE. 2008, 3: e1989-10.1371/journal.pone.0001989.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Yu Z, Jian Z, Shen SH, Purisima E, Wang E: Global analysis of microRNA target gene expression reveals that miRNA targets are lower expressed in mature mouse and Drosophila tissues than in the embryos. Nucleic Acids Res. 2007, 35: 152-164. 10.1093/nar/gkl1032.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Hartmann LC, Sellers TA, Schaid DJ, Nayfield S, Grant CS, Bjoraker JA, Woods J, Couch F: Clinical options for women at high risk for breast cancer. Surg Clin North Am. 1999, 79: 1189-1206. 10.1016/S0039-6109(05)70068-2.

    Article  PubMed  CAS  Google Scholar 

  29. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006, 103: 2257-2261. 10.1073/pnas.0510565103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin SF, Dunning MJ, Barbosa-Morais NL, Teschendorff AE, Green AR, Ellis IO, Tavare S, Caldas C, Miska EA: MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biol. 2007, 8: R214-10.1186/gb-2007-8-10-r214.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Menard S, Palazzo JP, Rosenberg A, Musiani P, Volinia S, Nenci I, Calin GA, Querzoli P, Negrini M, Croce CM: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070. 10.1158/0008-5472.CAN-05-1783.

    Article  PubMed  CAS  Google Scholar 

  32. Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, Fedele V, Ginzinger D, Getts R, Haqq C: Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer. 2006, 5: 24-10.1186/1476-4598-5-24.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Perou CM, Sorlie T, Eisen MB, Rijn van de M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.

    Article  PubMed  CAS  Google Scholar 

  34. van't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.

    Article  Google Scholar 

  35. Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, Wong JE, Liu ET, Bergh J, Kuznetsov VA, Miller LD: Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006, 66: 10292-10301. 10.1158/0008-5472.CAN-05-4414.

    Article  PubMed  CAS  Google Scholar 

  36. Pawitan Y, Bjohle J, Amler L, Borg AL, Egyhazi S, Hall P, Han X, Holmberg L, Huang F, Klaar S, Liu ET, Miller L, Nordgren H, Ploner A, Sandelin K, Shaw PM, Smeds J, Skoog L, Wedren S, Bergh J: Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohorts. Breast Cancer Res. 2005, 7: R953-964. 10.1186/bcr1325.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, Rijn van de M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lonning P, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  38. Miller LD, Smeds J, George J, Vega VB, Vergara L, Ploner A, Pawitan Y, Hall P, Klaar S, Liu ET, Bergh J: An expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. Proc Natl Acad Sci USA. 2005, 102: 13550-13555. 10.1073/pnas.0506230102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, Giri DD, Viale A, Olshen AB, Gerald WL, Massague J: Genes that mediate breast cancer metastasis to lung. Nature. 2005, 436: 518-524. 10.1038/nature03799.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.

    Article  PubMed  CAS  Google Scholar 

  41. Richardson AL, Wang ZC, De Nicolo A, Lu X, Brown M, Miron A, Liao X, Iglehart JD, Livingston DM, Ganesan S: X chromosomal abnormalities in basal-like human breast cancer. Cancer Cell. 2006, 9: 121-132. 10.1016/j.ccr.2006.01.013.

    Article  PubMed  CAS  Google Scholar 

  42. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39: 1278-1284. 10.1038/ng2135.

    Article  PubMed  CAS  Google Scholar 

  44. Hess KR, Anderson K, Symmans WF, Valero V, Ibrahim N, Mejia JA, Booser D, Theriault RL, Buzdar AU, Dempsey PJ, Rouzier R, Sneige N, Ross JS, Vidaurre T, Gomez HL, Hortobagyi GN, Pusztai L: Pharmacogenomic predictor of sensitivity to preoperative chemotherapy with paclitaxel and fluorouracil, doxorubicin, and cyclophosphamide in breast cancer. J Clin Oncol. 2006, 24: 4236-4244. 10.1200/JCO.2006.05.6861.

    Article  PubMed  CAS  Google Scholar 

  45. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet. 2005, 37: 495-500. 10.1038/ng1536.

    Article  PubMed  CAS  Google Scholar 

  46. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human microRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.

    Article  PubMed  CAS  Google Scholar 

  48. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455: 64-71. 10.1038/nature07242.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136: 215-233. 10.1016/j.cell.2009.01.002.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  51. Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105. 10.1101/gr.082701.108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Adams BD, Furneaux H, White BA: The micro-ribonucleic acid (miRNA) miR-206 targets the human estrogen receptor-alpha (ERalpha) and represses ERalpha messenger RNA and protein expression in breast cancer cell lines. Mol Endocrinol. 2007, 21: 1132-1147. 10.1210/me.2007-0022.

    Article  PubMed  CAS  Google Scholar 

  53. Merritt WM, Lin YG, Han LY, Kamat AA, Spannuth WA, Schmandt R, Urbauer D, Pennacchio LA, Cheng JF, Nick AM, Deavers MT, Mourad-Zeidan A, Wang H, Mueller P, Lenburg ME, Gray JW, Mok S, Birrer MJ, Lopez-Berestein G, Coleman RL, Bar-Eli M, Sood AK: Dicer, Drosha, and outcomes in patients with ovarian cancer. N Engl J Med. 2008, 359: 2641-2650. 10.1056/NEJMoa0803785.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  54. Chiosea S, Jelezcova E, Chandran U, Luo J, Mantha G, Sobol RW, Dacic S: Overexpression of Dicer in precursor lesions of lung adenocarcinoma. Cancer Res. 2007, 67: 2345-2350. 10.1158/0008-5472.CAN-06-3533.

    Article  PubMed  CAS  Google Scholar 

  55. Karube Y, Tanaka H, Osada H, Tomida S, Tatematsu Y, Yanagisawa K, Yatabe Y, Takamizawa J, Miyoshi S, Mitsudomi T, Takahashi T: Reduced expression of Dicer associated with poor prognosis in lung cancer patients. Cancer Sci. 2005, 96: 111-115. 10.1111/j.1349-7006.2005.00015.x.

    Article  PubMed  CAS  Google Scholar 

  56. Wiesen JL, Tomasi TB: Dicer is regulated by cellular stresses and interferons. Mol Immunol. 2009, 46: 1222-1228. 10.1016/j.molimm.2008.11.012.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Bhattacharyya SN, Habermacher R, Martine U, Closs EI, Filipowicz W: Relief of microRNA-mediated translational repression in human cells subjected to stress. Cell. 2006, 125: 1111-1124. 10.1016/j.cell.2006.04.031.

    Article  PubMed  CAS  Google Scholar 

  58. Marsit CJ, Eddy K, Kelsey KT: MicroRNA responses to cellular stress. Cancer Res. 2006, 66: 10843-10848. 10.1158/0008-5472.CAN-06-1894.

    Article  PubMed  CAS  Google Scholar 

  59. Ketting RF, Fischer SE, Bernstein E, Sijen T, Hannon GJ, Plasterk RH: Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans. Genes Dev. 2001, 15: 2654-2659. 10.1101/gad.927801.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  60. Chendrimada TP, Gregory RI, Kumaraswamy E, Norman J, Cooch N, Nishikura K, Shiekhattar R: TRBP recruits the Dicer complex to Ago2 for microRNA processing and gene silencing. Nature. 2005, 436: 740-744. 10.1038/nature03868.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  61. Diederichs S, Haber DA: Dual role for argonautes in microRNA processing and posttranscriptional regulation of microRNA expression. Cell. 2007, 131: 1097-1108. 10.1016/j.cell.2007.10.032.

    Article  PubMed  CAS  Google Scholar 

  62. Haase AD, Jaskiewicz L, Zhang H, Laine S, Sack R, Gatignol A, Filipowicz W: TRBP, a regulator of cellular PKR and HIV-1 virus expression, interacts with Dicer and functions in RNA silencing. EMBO Rep. 2005, 6: 961-967. 10.1038/sj.embor.7400509.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  63. MacRae IJ, Ma E, Zhou M, Robinson CV, Doudna JA: In vitro reconstitution of the human RISC-loading complex. Proc Natl Acad Sci USA. 2008, 105: 512-517. 10.1073/pnas.0710869105.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  64. Meister G, Landthaler M, Patkaniowska A, Dorsett Y, Teng G, Tuschl T: Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol Cell. 2004, 15: 185-197. 10.1016/j.molcel.2004.07.007.

    Article  PubMed  CAS  Google Scholar 

  65. Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T: Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nat Genet. 2007, 39: 673-677. 10.1038/ng2003.

    Article  PubMed  CAS  Google Scholar 

  66. Adams BD, Claffey KP, White BA: Argonaute-2 expression is regulated by epidermal growth factor receptor and mitogen-activated protein kinase signaling and correlates with a transformed phenotype in breast cancer cells. Endocrinology. 2009, 150: 14-23. 10.1210/en.2008-0984.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  67. Salomon DS, Bianco C, De Santis M: Cripto: a novel epidermal growth factor (EGF)-related peptide in mammary gland development and neoplasia. Bioessays. 1999, 21: 61-70. 10.1002/(SICI)1521-1878(199901)21:1<61::AID-BIES8>3.0.CO;2-H.

    Article  PubMed  CAS  Google Scholar 

  68. McKay MM, Morrison DK: Integrating signals from RTKs to ERK/MAPK. Oncogene. 2007, 26: 3113-3121. 10.1038/sj.onc.1210394.

    Article  PubMed  CAS  Google Scholar 

  69. Yarden Y, Sliwkowski MX: Untangling the ErbB signalling network. Nat Rev Mol Cell Biol. 2001, 2: 127-137. 10.1038/35052073.

    Article  PubMed  CAS  Google Scholar 

  70. Adams BD, Cowee DM, White BA: The role of miR-206 in the epidermal growth factor (EGF) induced repression of estrogen receptor-alpha (ER{alpha}) signaling and a luminal phenotype in MCF-7 breast cancer cells. Mol Endocrinol. 2009, 23: 1215-1230. 10.1210/me.2009-0062.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  71. Thyrell L, Hjortsberg L, Arulampalam V, Panaretakis T, Uhles S, Dagnell M, Zhivotovsky B, Leibiger I, Grander D, Pokrovskaja K: Interferon alpha-induced apoptosis in tumor cells is mediated through the phosphoinositide 3-kinase/mammalian target of rapamycin signaling pathway. J Biol Chem. 2004, 279: 24152-24162. 10.1074/jbc.M312219200.

    Article  PubMed  CAS  Google Scholar 

  72. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.

    Article  PubMed  Google Scholar 

  73. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge support from the NIH and from the AL Williams Professorship funds.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mark Gerstein.

Additional information

Authors' contributions

CC and MG conceived and designed the study. CC extracted the gene expression data. CC and XF preformed the full analysis. CC, XF, PA and MG wrote the manuscript.

Chao Cheng, Xuping Fu contributed equally to this work.

Electronic supplementary material

13059_2009_2240_MOESM1_ESM.pdf

Additional data file 1: Distribution of miRNA target numbers for four prediction tools, PITA, miRanda, PicTar, and TargetScan. In addition, three different parameters in TargetScan were chosen and are denoted as 'Conserved', 'ContextScore ≥ -0.20' and 'All', respectively. On average, 6,949, 2,026, 1,563, 765, 426, and 210 targets per miRNA were predicted by PITA, TargetScan(All), miRanda, TargetScan(ContextScore ≥ -0.20), PicTar, and TargetScan(Conserved), respectively. (PDF 15 KB)

13059_2009_2240_MOESM2_ESM.xls

Additional data file 2: Includes seven sheets containing the complete results based on RE-scores from rank comparison and seven miRNA target predictions from PITA, miRanda, PicTar, the intersection of PITA and miRanda, TargetScan(Conserved), TargetScan(ContextScore ≥ -0.20), and TargetScan(All), respectively. The t-score, P-value, and adjusted P-value (FDR) of all miRNAs in the five breast cancer datasets are provided. (XLS 960 KB)

13059_2009_2240_MOESM3_ESM.xls

Additional data file 3: Additional data file 3 includes two sheets containing the complete results based on RE-scores calculated from expression comparison. In the two sheets, miRNA target predictions determined by the miRanda and PITA tools are used, respectively. (XLS 298 KB)

Additional data file 4: Expression levels of several miRNA processing genes in ER+ and ER- samples. (XLS 21 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Cheng, C., Fu, X., Alves, P. et al. mRNA expression profiles show differential regulatory effects of microRNAs between estrogen receptor-positive and estrogen receptor-negative breast cancer. Genome Biol 10, R90 (2009). https://doi.org/10.1186/gb-2009-10-9-r90

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2009-10-9-r90

Keywords