Correlative Analysis of miRNA Expression and Oncotype Dx Recurrence Score in Estrogen Receptor Positive Breast Carcinomas

Abstract

Altered expression of miRNAs has been observed in many types of cancer, including breast cancer, and shown to contribute to cancer growth, aggressiveness, and response to therapies. In this pilot study, we investigated the possible correlation of miRNAs with risk of recurrence of estrogen receptor positive, lymph node-negative mammary carcinomas as determined by the Oncotype DX^® Breast Cancer assay. To accomplish this, we extracted RNA from a collection of breast carcinomas that had previously been analyzed by Oncotype DX^®. Multiple Let-7 family members were negatively correlated with the recurrence score (RS), which is consistent with their tumor suppressor properties. Additional miRNAs were found to positively correlate with RS, including miR-377-5p, miR-633b, miR-548t and miR-3648. Pathway analysis of putative and validated targets suggests that these miRNAs may have a diverse range of functions that may contribute to tumor recurrence. Taken together, these findings provide evidence that a miRNA expression signature can be developed to aid existing methods to determine the risk of recurrence for women with estrogen receptor positive breast cancers treated with endocrine therapy.

About the Authors:

Rajyasree Emmadi

* E-mail: [email protected]

Affiliation: Department of Pathology, University of Illinois at Chicago, Chicago, Illinois, United States of America

Emanuele Canestrari

Affiliation: Center for Pharmaceutical Biotechnology, University of Illinois at Chicago, Chicago, Illinois, United States of America

Zarema H. Arbieva

Affiliation: Core Genomics Facility, University of Illinois at Chicago, Chicago, Illinois, United States of America

Wenbo Mu

Affiliation: Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois, United States of America

Yang Dai

Affiliation: Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois, United States of America

Jonna Frasor

Affiliation: Department of Physiology and Biophysics, University of Illinois at Chicago, Chicago, Illinois, United States of America

Elizabeth Wiley

Affiliation: Department of Pathology, University of Illinois at Chicago, Chicago, Illinois, United States of America

Introduction

Invasive breast tumors are biologically diverse with different rates of progression, treatment response, and outcomes. The current standard of care requires all newly diagnosed cases of invasive breast cancer to be routinely tested by immunohistochemistry (IHC) for the expression of estrogen receptor alpha (ER), progesterone receptor (PR) and the growth factor receptor HER2/neu (ERBB2). ER-positive tumors are treated with endocrine therapy, such as tamoxifen or aromatase inhibitors. However, about one third to one half of patients treated with tamoxifen or other endocrine therapy will develop recurrence, in part due to development of acquired endocrine resistance. The differential responsiveness of ER-positive breast cancers to endocrine therapy is explained by the heterogeneity of these tumors, which has been clearly demonstrated by genome-wide expression profiling studies. The two sub-classifications of ER-positive tumors, representing biologically distinct tumor types, luminal A and luminal B, have been identified, with the luminal A tumors having a good response to endocrine therapy and luminal B tumors responding poorly [1, 2].

Based on gene expression studies, the Oncotype DX^® assay was developed through the NSABP clinical trials [3]. Since 2007, the Association for Clinical Oncology (ASCO) Guidelines have used Oncotype DX^® testing for treatment stratification of ER-positive, lymph node negative breast carcinomas [4]. Oncotype DX^® is a reverse transcription polymerase chain reaction (RT-PCR) based assay that is performed on RNA isolated from formalin fixed paraffin embedded (FFPE) tumor tissue blocks. Based on the expression of 21 genes in the tumor, a Recurrence Score (RS) is issued. The 16 cancer genes included in the assay are associated with proliferation (Ki67, STK15, Survivin, Cyclin B1, MYBL2), estrogen action (ER, PR, Bcl-2, SCUBE2), HER2/cERBB2 action (GRB7, HER2), invasion (Stromelysin 3, Cathepsin L2), or other functions (GSTM1, CD68 and BAG1). Five reference genes (beta-actin, GADPH, RPLPO, GUS and TFRC) are also included in the assay. The recurrence score is represented by a number between 0 and 100, which correlates to a specific likelihood of breast cancer recurrence within 10 years of initial diagnosis. The scores are grouped into 3 risk categories (low: RS ≤17, intermediate: RS = 18–30, and high: RS>31). Oncologists use the RS to tailor the use of chemotherapies based on the finding that women with high RS benefit from chemotherapy whereas women with low RS do not [5]. Hence, women at low risk for recurrence may avoid the toxicity of chemotherapy and be successfully managed with hormonal therapy alone. Patients with high RS, on the other hand, are at an increased risk of recurrence and may receive one or both types of therapy. Little is known about what distinguishes the intermediate from the high RS groups so they generally receive similar therapies. In fact, in 2006, the National Cancer Institute (NCI) sponsored a prospective randomized clinical trial, the TAILORx trial, involving hormone-receptor-positive, HER2-negative, axillary node-negative breast cancers less than 5.0 cm in size, assigned to one of three treatment arms based on their OncotypeDX RS results. The results of this trial have recently been released, showing that patients with tumors that had favorable gene-expression profiles (RS of 0–10) had very low rates of recurrence at 5 years with endocrine therapy alone. Of note, however, is that new cut-offs were used to designate the three subgroups of RS compared to the original OncotypeDX assay. The new RS ranges were 0–10 (low risk), 11–25 (intermediate risk) and ≥26 (high risk), in contrast to prior ranges of <18, 18–30 and ≥31, respectively. The overall sample size of 10,253 women was driven by the need to include a sufficient number of patients in the ‘new’ intermediate risk range of RS of 11–25.

In addition to gene expression signatures, recent studies have revealed associations between miRNAs and tumor expression of ER, PR and HER2, as well as breast carcinoma subtypes and tumor grades [6, 7]. In addition to these correlative studies, specific miRNAs have been shown to play key roles in hormone-dependent breast cancer biology. For example, down-regulation of miR-451 can promote breast cancer cell survival and endocrine resistance [8], a miRNA196a2-TP63 circuit controls breast cancer proliferation and invasiveness properties [9], loss of miR-200c promotes breast cancer cell migration and invasion [10], and down-regulation of miR-29 contributes to progestin-induced stem cell expansion [11]. Let-7 family members have been found to inhibit breast cancer stem cell renewal [12] and have been implicated in tamoxifen response through modulation of ER levels [13]. Furthermore, miRNAs that are differentially expressed in tamoxifen-sensitive vs. resistant cell lines and tumors have been identified [14, 15]. These findings suggest that miRNAs may function as biomarkers of tumor response to therapy and/or risk of recurrence, as well as important drivers of these phenotypically different tumors. In this study we conducted miRNA profiling in a collection of breast tumors with previous testing for recurrence risk by the Oncotype DX^® assay in an attempt to identify additional miRNAs biomarkers of breast tumor recurrence and poor outcome in patients with ER-positive breast cancers. This is a pilot study with a small sample size. However, cases were selected to span the range of the Oncotype DX^® Recurrence Scores. We acknowledge the limitations of the study including retrospective analysis and a small sample size.

Materials and Methods

Patient Cohort and Tumor Characteristics

This research was a retrospective, non-interventional analysis, conducted following approval from the Institutional Review Boards (IRBs) of the University of Illinois Cancer Center and Provena Saint Joseph Medical Center in Illinois. The research was judged to qualify for waiver of informed consent based on the provisions under HHS regulations at 45 CFR 46.116(d). The research did not involve any risk to patients, data was anonymized and no patient identifiers were included. Formalin-fixed, paraffin-embedded (FFPE) tumor tissues from twenty-three cases of early stage breast carcinomas representing low, intermediate and high Oncotype DX Recurrence Score were obtained. Patient samples and clinical data were collected and processed in compliance with protocols approved by the University of Illinois Cancer Center and Provena Saint Joseph Medical Center Institutional Review Boards. Tumor tissue was macro-dissected by comparison to an adjacent H&E stained section to ensure that tissue used for miRNA analysis contained >70% tumor, in accordance with samples sent for Oncotype DX^® testing. Clinical data collected on each patient included information on ER, PR, p53, Ki67 and Her2 status for correlation with study results (Table 1).

[Figure omitted. See PDF.]

Table 1. Patient and tumor characteristics.

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

miRNA Profiling

RNA was extracted from each sample using the FFPE RNA Purification Kit, according to the manufacturer’s instructions (Norgen Biotek Corp., Thorold, ON, Canada). The quality of the total RNA was verified by Agilent 2100 Bioanalyzer profile. Total RNA from each sample was labeled with Hy3 and a reference sample, consisting of equal amounts of total RNA combined from every individual sample, was labeled with Hy5 fluorescent label, using the miRCURY LNAmicroRNA Hi-Power Labeling Kit (Exiqon, Denmark) following the procedure described by the manufacturer. The Hy3-labeled samples and a Hy5-labeled reference sample were mixed pair-wise and hybridized to the miRCURY LNA microRNA Array 6th gen (Exiqon, Denmark), which contains capture probes targeting all microRNAs for human, mouse or rat registered in the miRBASE 16.0. The hybridization was performed using a Tecan HS4800 hybridization station (Tecan, Austria). After hybridization the microarray slides were scanned and stored in an ozone free environment (ozone level below 2.0 ppb) in order to prevent potential bleaching of the fluorescent dyes. The miRCURY LNA microRNA Array slides were scanned using the Agilent G2565BA Microarray Scanner System (Agilent Technologies, Inc., USA) and the image analysis was carried out using ImaGene 9 software (Exiqon, Denmark). The quantified signals were background corrected (Normexp with offset value 10) [16] and normalized using the global Lowess regression algorithm.

Data Analysis

Principal Components Analysis (PCA) was performed using the top 50 microRNAs that had the largest variation across all samples (Fig 1). A heatmap diagram was generated based on the two-way clustering method with the complete-linkage and Euclidean distance measurement (Fig 2). The intermediate RS group was included in the high RS group since the number of samples in the former was small. Student’s t-test was performed for differential expression analysis between low RS and high RS groups. Pearson correlation coefficient was used for correlation tests between the expression value of individual microRNAs and RS. All p-values were adjusted based on the Benjamini-Hochberg procedure. The threshold for the adjusted p-value <0.05 was used to determine the significant correlation (Table 2). The analysis was carried out using packages in R.

[Figure omitted. See PDF.]

Fig 1. Principal components analysis.

Principal components analysis performed on all samples using the top 50 miRNAs showing the highest standard deviation across all samples.

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

[Figure omitted. See PDF.]

Fig 2. Heat map of unsupervised hierarchical clustering.

Hierarchical clustering performed on all samples using the top 50 miRNAs showing the highest standard deviation across all samples.

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

[Figure omitted. See PDF.]

Table 2. Differential expression of miRNAs.

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

Predicted and Validated miRNA Targets

The miRNA story is an evolving one and the currently available databases are still incomplete. In order to maximize the information for correlative analysis, we looked at multiple miRNA databases. Predicted miRNA targets were obtained using DIANA-microT-CDS [17, 18]. Validated miRNA targets were obtained using miRWalk with a default p-value of 0.05 [19]. The findings of the predicted miRNA targets and pathways are presented in Fig 3 and validated targets and pathways are presented in Table 3 and Fig 4.

[Figure omitted. See PDF.]

Fig 3. Pathway analysis for putative miRNA target genes.

Pathways showing significant enrichment in putative target genes for miRNAs positively correlated with recurrence score (left) or negatively correlated with recurrence score (right) were clustered.

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

[Figure omitted. See PDF.]

Fig 4. Pathway analysis for validated miRNA target genes.

Pathways showing enrichment in validated targets genes for miR-377-5p and miR-663b.

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

[Figure omitted. See PDF.]

Table 3. Validated gene targets for miR-377-5p and miR-663b using miRWalk.

https://doi.org/10.1371/journal.pone.0145346.t003

Pathway Analysis

Predicted targets.

miRNA targets were predicted based on DIANA-microT-CDS. The combined effect of groups of miRNAs on pathways was determined based on the union of target genes for enrichment analysis using DIANA miRPath v.2.0 (Vlachos). Heatmaps were generated using the log-transformed enrichment p-values as features for each miRNA. [20].

Validated targets.

Validated miRNA targets were available only for miR-377-5p (miR-377*) and miR-663b using miRWalk. No gene targets were elucidated for miR-3648 and miR-748t. The gene targets for miR-377-5p and miR-663b were combined and pathway enrichment analysis was performed using the Panther classification system v9.0 and Panther pathway, enriched with a default p-value of 0.05 [21, 22].

Results

Genome wide profiling of miRNA expression in 23 human breast tumors was carried out using a locked nucleic acid (LNA) based microRNA array that contained ~2800 capture probes complementary to mature microRNAs based on miRBase v16 (Exiqon). All tumors were ER-positive, lymph node-negative, invasive mammary carcinomas. Clinical and laboratory data collected for each patient are shown in Table 1 and include the Recurrence Score (RS) based on the Oncotype DX^® assay and the estimated risk of recurrence (R-Rate). The tumors profiled included 12 low, 4 intermediate and 7 high RS scores.

Principal components analysis (PCA), which was conducted for the top 50 miRNAs with the largest variation in expression across all samples, revealed that the tumors with high or low RS largely clustered independently of each other but that tumors with intermediate RS were more widely distributed (Fig 1). Of the samples with low RS, one (#S7) was more closely related to the high RS samples than other tumors with low RS. Two-way hierarchical clustering by miRNA expression levels confirmed the relationship between samples and recurrence score groups (Fig 2).

To identify specific miRNAs that are associated with RS, two types of analyses were performed. First, differential expression of miRNAs between the low and high RS groups was carried out using Student’s T test followed by Benjamini Hochberg adjustment for multiple tests. Twenty-five miRNAs were found to be statistically different (P<0.05) between the two groups (Table 2). However, when all samples were analyzed as three groups (low, intermediate, and high) using ANOVA and Benjamini Hochberg adjustment, only three of miRNAs remained as significantly different among the groups. In the second approach, correlation analysis between miRNA expression levels and the actual RS was conducted. Expression levels for seven miRNAs were significantly correlated with RS (Fig 5), six of which overlapped with other analyses as indicated in Table 2.

[Figure omitted. See PDF.]

Fig 5. Correlation analysis between miRNA expression level and recurrence score.

The relationship between miRNA expression and recurrence score is shown for seven miRNAs with a statistically significant correlation.

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

miRNAs showing negative correlation with RS comprise several members of the Let-7 family, which have been found to inhibit breast cancer stem cell renewal [12] and have been implicated in Tamoxifen response through modulation of ER levels [13]. However, little is known about the other miRNAs (miR-377-5p, miR-663b, miR-548T and miR-3648) that we found to be higher in breast tumors with high RS.

To determine pathways that may be controlled by these miRNAs, DIANA miRPath was first used to identify putative targets (see S1 Data for complete results). The union of the putative targets from Let-7a, 7b and 7c (all negatively correlated with RS) and the union of putative targets from miR-663b, miR-377-5p and miR-3648 (all positively correlated with RS) were separately used for pathway enrichment analysis. The miRNAs were then clustered based on the log-transformed enrichment p-values and visualized in heatmaps. As shown in Fig 3, target genes for miRNAs positively correlated with RS are fairly divergent with chemical carcinogenesis and xenobiotic metabolism associated miR-3648, and glycan degradation and systemic lupus erythematosus with miR-633b. ErbB signaling, on the other hand, was enriched for both miR-377-5p and miR-633b. In contrast, pathways well characterized in breast cancer progression (PI3K/AKT, MAPK, cytokines, and WNT), are highly enriched for all three Let-7 family members showing positive correlation with RS score.

Validated gene targets obtained for miR-377-5p and miR-663b using miRWalk (Table 3) were then used for pathway enrichment analysis using PANTHER 9.0 and showed association with multiple pathways, in particular those involving angiogenesis, EGF receptor signaling, inflammation mediated by chemokines and cytokines, integrin signaling, PI3 kinase, TGF-beta signaling and Ras and p53 pathway feedback loops (Fig 4).

Discussion

Several publications have looked at the miRNA expression profiles in breast carcinoma cell lines and tumor tissues. Most of the studies have varied representation of the different subtypes of breast carcinoma (luminal, basal-like, Her2-enriched) and with or without lymph node involvement. As a result, mixed (non subtype-specific) miRNA signatures emerge and have been reported in ranges of 1–21 miRNAs [23–27]. Commonly reported differential expressions include upregulation of miR-21, miR-155, miR-206, miR-122a and miR-210 and downregulation of let-7, miR-10b, miR-125a, miR-125b and miR-145. Putative or validated target gene and pathway analyses demonstrate involvement of oncogenes and tumor suppressors, such as ErbB2, Akt, NF-κB, Myc, p53, Rb and PTEN with consequent proliferation and cell survival effects as well as inhibition of apoptosis. While cell cycle and proliferation signals would be expected to correlate with outcomes, other signatures described are associated with hypoxia, cytokines, stroma and chromosomal instability. In essence it appears that different reported signatures may represent different aspects of ultimate multifactorial mechanisms. Several of the miRNAs associated with RS in ER+ breast tumors have been described to play a role in cancer progression. For example, Let-7 family members, which in this study were found to be down-regulated in tumors with high recurrence score, have been found to inhibit breast cancer stem cell renewal [12] and have been implicated in tamoxifen response through modulation of ER levels [13]. On the other hand, little is known about the other miRNAs that we find to be higher in breast tumors with high RS. The inferences we have made by correlating the results of our analyses and literature are presented.

1. miRNA profiling captures basic biological relationship between high vs low risk of recurrence. More refined risk prediction is needed—both Oncotype DX^® and the MammaPrint^® assays have a discordance rate between the assay prediction and clinical-pathologic risk category of approximately 30% [28].

2. The intermediate group does not represent a unique group biologically but appears to be distributed along a continuum of high to low. This is similar to findings reported previously by Ivshina et al, who found 6 markers that could effectively separate grade 1 and grade 3 tumors and could also separate grade 2 tumors into two highly discriminate classes [29].

3. Suppression of the Let-7 family is associated with increased RS score. Of particular interest is the reduced expression of Let-7 family members in cases with high RS. Six family members were found to be significantly associated with low risk of recurrence using at least one statistical analysis. The hsa-let-7 family is well established as a tumor suppressor with interactions with cell cycle proliferative genes including Aurora Kinase A (AURKA), which it downregulates [30]. A recent meta-analysis shows a significant correlation between increased AURKA expression and decreased metastasis-free survival in ER-positive breast cancers [31]. The hsa-let-7 family involvement in breast carcinoma is well documented and serves to confirm the validity of our findings.

4. Novel miRNAs involved in breast carcinoma. The three miRNAs (hsa-miR-377-5p, hsa-miR-663b and hsa-miR-3648) that are positively correlated with high Oncotype DX^® RS scores have not been previously detailed in association with breast carcinoma outcomes. Very little specific data exists regarding the interactions and roles of these miRNAs with regard to ER-positive breast carcinoma.

Information available via the tools listed suggest that multiple pathways are involved by miR-377-5p and miR-663b, affecting tumor suppressor functions (LATS2, TXNIP, wtKRAS, CDKN1A, PPM1A), cytoskeleton regulation (PAK1, PKN1) and a proinflammatory and reactive oxygen species adaptive response (TNF, IFNG, IL1B, HMOX1, SOD1 and SOD2). We infer that multiple factors are involved in separating those tumors that have a high RS score including specific proliferation advantages as well as suppression of tumor response mechanisms.

While these genes and pathways may represent possible therapeutic targets for tumor suppression, these miRNAs also are of particular interest as they may provide novel approaches for specific, targeted, less toxic systemic treatment protocols in the future. Progress in microRNA-directed therapeutic approaches offers hope that such strategies might prove useful in reversing endocrine resistance and reducing breast cancer recurrence.

Supporting Information

[Figure omitted. See PDF.]

S1 Data.

https://doi.org/10.1371/journal.pone.0145346.s001

(XLSX)

Acknowledgments

This research was supported by the University of Illinois at Chicago through the Chancellor’s Discovery Fund (JF, EW, YD) and the National Institute of Health (R01 CA130932-05 to JF).

The research was conducted in collaboration with Provena St. Joseph Medical Center after obtaining IRBs from each institution. We thank Dr. J. Natarelli (Provena St. Joseph Medical center, Pathology) for reviewing the cases and selecting tumor tissue blocks for analysis.

Author Contributions

Conceived and designed the experiments: RE JF YD. Performed the experiments: RE EC ZA WM YD. Analyzed the data: RE JF ZA WM YD. Contributed reagents/materials/analysis tools: RE JF ZA YD. Wrote the paper: RE JF EW YD.

Citation: Emmadi R, Canestrari E, Arbieva ZH, Mu W, Dai Y, Frasor J, et al. (2015) Correlative Analysis of miRNA Expression and Oncotype Dx Recurrence Score in Estrogen Receptor Positive Breast Carcinomas. PLoS ONE 10(12): e0145346. https://doi.org/10.1371/journal.pone.0145346

References

1. Voduc KD, Cheang M, Tyldesley S, Gelmon K, Nielsen TO, Kennecke H. Breast Cancer Subtypes and the Risk of Local and Regional Relapse. J Clin Oncol. 2010; 28:1684–1691. pmid:20194857

2. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thürliman B, et al. Panel members. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Annals of Oncology. 2013; 24: 2206–2223 pmid:23917950

3. Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004; 351:2817–2826 pmid:15591335

4. Harris L, Fritsche H, Mennel R, Norton L, Ravdin P, Taube S, et al. American Society of Clinical Oncology 2007 update of recommendations for the use of tumor markers in breast cancer. J Clin Oncol. 2007; 25:5287–5312 pmid:17954709

5. Paik S, Tang G, Shak S, Kim C, Baker J, Kim W, et al. Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol. 2006; 24:3726–3734 pmid:16720680

6. Lowery AJ, Miller N, Devaney A, McNeill RE, Davoren PA, Lemetre C, et al. MicroRNA signatures predict oestrogen receptor, progesterone receptor and HER2/neu receptor status in breast cancer. Breast Cancer Res. 2009; 11(3):R27 pmid:19432961

7. Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin S- F, Dunning MJ, et al. MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biology. 2007; 8:R214 pmid:17922911

8. Bergamaschi A, Katzenellenbogen BS. Tamoxifen downregulation of miR-451 increases 14-3-3zeta and promotes breast cancer cell survival and endocrine resistance. Oncogene. 2012; 31:39–47 pmid:21666713

9. Kim K, Madak-Erdogan Z, Ventrella R, Katzenellenbogen BS. A MicroRNA196a2* and TP63 Circuit Regulated by Estrogen Receptor-alpha and ERK2 that Controls Breast Cancer Proliferation and Invasiveness Properties. Horm Cancer 2013 Apr;4(2):78–91 pmid:23250869

10. Cochrane DR, Spoelstra NS, Howe EN, Nordeen SK, Richer JK 2. MicroRNA-200c mitigates invasiveness and restores sensitivity to microtubule-targeting chemotherapeutic agents. Mol Cancer Ther. 2009; 8:1055–1066 pmid:19435871

11. Cittelly DM, Finlay-Schultz J, Howe EN, Spoelstra NS, Axlund SD, Hendricks P, et al. Progestin suppression of miR-29 potentiates dedifferentiation of breast cancer cells via KLF4. Oncogene. 2013 May 16;32(20):2555–64 pmid:22751119

12. Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C, et al. let-7 Regulates Self Renewal and Tumorigenicity of Breast Cancer Cells. Cell 2007; 131(6):1109–1123 pmid:18083101

13. Zhao Y, Deng C, Lu W, Xiao J, Ma D, Guo M, et al. let-7 MicroRNAs Induce Tamoxifen Sensitivity by Downregulation of Estrogen Receptor α Signaling in Breast Cancer. Mol Med. 2011; 17(11–12)1233–1241 pmid:21826373

14. Miller TE, Ghoshal K, Ramaswamy B, Roy S, Datta J, Shapiro CL, et al. MicroRNA-221/222 confers tamoxifen resistance in breast cancer by targeting p27Kip1. J Biol Chem. 2008 Oct 31; 283(44):29897–903. pmid:18708351

15. Manavalan TT, Teng Y, Appana SN, Datta S, Kalbfleisch TS, Li Y, et al. Differential expression of microRNA expression in tamoxifen-sensitive MCF-7 versus tamoxifen-resistant LY2 human breast cancer cells. Cancer Letters, 2011; 313(1):26–43. pmid:21955614

16. Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, et al. A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007; 23(20):2700–2707. pmid:17720982

17. Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, at al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013 Jul; 41(Web Server issue):W169–73. pmid:23680784

18. Reczko M, Maragkakis M, Alexiou P, Grosse I, Hatzigeorgiou AG. Functional microRNA targets in protein coding sequences. Bioinformatics. 2012; 28(6):771–776. pmid:22285563

19. Dweep H, Sticht C, Pandey P, Gretz N. miRWalk—database: prediction of possible miRNA binding sites by "walking" the genes of 3 genomes, Journal of Biomedical Informatics. 2011; 44: 839–7. pmid:21605702

20. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, et al. DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Research 2012 Jul; 40(Web server issue): W498–504. pmid:22649059

21. Mi H, Muruganujan A, Thomas PD. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucl. Acids Res. 2013; 41(D1):D377–D386.

22. Mi H, Thomas PD. Protein Networks and Pathway Analysis. Methods in Molecular biology, 2009; 563 Part 2: 123–140.

23. Wang V and Wu W. MicroRNA: A New Player in Breast Cancer Development. Journal of Cancer Molecules. 2007; 3(5): 133–138.

24. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, et al. MicroRNA expression profiles classify human cancers. Nature. 2005; 435:834–838. pmid:15944708

25. Iorio M, Ferracin M, Liu C-G, Veronese A, Spizzo R, Sabbioni S, et al. MicroRNA Gene Expression Deregulation in Human Breast Cancer. Cancer Res. 2005; 65(16):7065–70. pmid:16103053

26. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006; 103(7):2257–61. pmid:16461460

27. Si M-L, Zhu S, Wu H, Lu Z, Wu F and Mo Y-Y. miR-21-mediated tumor growth. Oncogene. 2007; 26:2799–2803. pmid:17072344

28. Albain KS, Paik S, van’t Veer L. Prediction of adjuvant chemotherapy benefit in endocrine responsive, early breast cancer using multigene assays. The Breast. 2009; S*, S141–S145. pmid:19914534

29. Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, et al. Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006 Nov 1; 66 (21):10292–301. pmid:17079448

30. Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, et al. The let-7 MicroRNA Represses Cell Proliferation Pathways in Human Cells. Cancer Res. 2007;67(16):7713–22 pmid:17699775

31. Weier H-U G, Mao J-H. Meta-analysis of Aurora Kinase A (AURKA) Expression Data Reveals a Significant Correlation between Increased AURKA Expression and Distant Metastases in Human ER-positive Breast Cancers. J Data Mining Genomics Proteomics. 2013; 4(1):127. pmid:24729942