Decreased expression of key tumour suppressor microRNAs is associated with lymph node metastases in triple negative breast cancer
© Avery-Kiejda et al.; licensee BioMed Central Ltd. 2014
Received: 5 August 2013
Accepted: 29 January 2014
Published: 31 January 2014
Breast cancer is the most common malignancy that develops in women, responsible for the highest cancer-associated death rates. Triple negative breast cancers represent an important subtype that have an aggressive clinical phenotype, are associated with a higher likelihood of metastasis and are not responsive to current targeted therapies. miRNAs have emerged as an attractive candidate for molecular biomarkers and treatment targets in breast cancer, but their role in the progression of triple negative breast cancer remains largely unexplored.
This study has investigated miRNA expression profiles in 31 primary triple negative breast cancer cases and in 13 matched lymph node metastases compared with 23 matched normal breast tissues to determine miRNAs associated with the initiation of this disease subtype and those associated with its metastasis.
71 miRNAs were differentially expressed in triple negative breast cancer, the majority of which have previously been associated with breast cancer, including members of the miR-200 family and the miR-17-92 oncogenic cluster, suggesting that the majority of miRNAs involved in the initiation of triple negative breast cancer are not subtype specific. However, the repertoire of miRNAs expressed in lymph node negative and lymph node positive triple negative breast cancers were largely distinct from one another. In particular, miRNA profiles associated with lymph node negative disease tended to be up-regulated, while those associated with lymph node positive disease were down-regulated and largely overlapped with the profiles of their matched lymph node metastases. From this, 27 miRNAs were identified that are associated with metastatic capability in the triple negative breast cancer subtype.
These results provide novel insight into the repertoire of miRNAs that contribute to the initiation of and progression to lymph node metastasis in triple negative breast cancer and have important implications for the treatment of this breast cancer subtype.
KeywordsmiRNA Breast cancer Triple negative Metastases Tumour suppressor
Breast cancer is the most common malignancy that develops in women worldwide, its incidence continues to rise and it is responsible for the highest cancer-associated death rates . It is an extremely heterogeneous disease, made up of a number of different subtypes. Classification of breast cancer into subtypes can be partly attributed to the presence or absence of receptors for the hormones estrogen (ER), progesterone (PR) and human epidermal growth factor 2 (HER2). Triple negative breast cancer (TNBC) represents a particularly important clinical subtype, characterised by an absence of ER, PR and HER2 and which therefore lacks common targets used for anti-hormone therapies [2, 3]. Although TNBCs comprise only a small percentage of all breast cancers diagnosed (10-24%), they have been recently the subject of intense investigation because of their aggressive clinical behaviour. Patients who are diagnosed with TNBC are of younger age, tend to develop tumours of larger size, and have an increased likelihood of distant metastasis and death within 5 years of diagnosis [2, 3]. Thus, TNBCs represent a major clinical problem for which targeted therapies are currently not available.
microRNAs (miRNAs) are a class of small (~22 nucleotides) non-coding RNAs that control gene expression by targeting mRNAs and triggering either translational repression or RNA degradation . miRNAs have emerged as an attractive candidate for molecular biomarkers and novel therapeutic targets in cancer because of their stability, ease of detection and their ability to act as endogenous antisense regulators of entire gene sets that regulate cancer growth .
Several studies have identified critical roles for miRNAs in breast cancer. Iorio and colleagues were the first to report significant deregulation in miRNA expression profiles in breast cancer when compared to normal breast tissue, where they showed that the expression of several miRNAs was associated with breast cancer subtypes and clinico-pathological features including hormone receptor status, clinical stage and proliferation index , suggesting that miRNA expression in breast cancer may have diagnostic and prognostic value. Since then, several miRNAs have been identified that have prognostic significance, including miR-210, miR-126, miR-21 and miR-205 [7–11]. However, the clinical cohorts used to generate these miRNA profiles have been highly heterogeneous, in regards to clinicopathological variables (e.g. grade, tumour size, subtype). Moreover, because these studies have concentrated on using miRNA profiles to better classify breast cancer subtypes, there are few studies that have profiled all miRNAs within a particular breast cancer subtype and even fewer that have used matched normal tissue as controls in these analyses; and so it is plausible that miRNAs important in disease initiation and progression that are subtype specific have been missed. TNBCs are distinct from ER-, PR- positive tumours at the molecular level and this classification has clinical implications . Therefore, the identification of miRNAs that control TNBC initiation and progression could identify individuals that have more aggressive disease and may also help to identify subgroups of patients that are more responsive to particular treatments within these subtypes.
The dissemination of primary cancer cells to the lymphatic system represents one of the first signs of metastatic spread. In breast cancer, the number of positive lymph nodes (LN) is known to have an inverse linear correlation with prognosis and survival . This is not the case with TNBC, where it has been shown that any LN involvement is associated with worse disease-free and overall survival . Therefore, identification of miRNAs that are differentially expressed within LN metastases in this highly aggressive breast cancer subtype may serve as a better indicator of prognosis than miRNA profiles derived from primary breast cancers.
In this study, miRNA expression profiles of 31 primary TNBCs were examined and compared to the profiles of 23 matched normal breast tissues to reveal miRNAs that were differentially expressed in TNBC and that are potentially associated with its initiation. In addition, by comparing the miRNA expression profiles of lymph node positive primary tumours and matched LN metastases, we have identified a panel of miRNAs that are associated with metastatic capability in the TNBC subtype.
Study cohort and tissue sampling
Demographic data of triple negative breast cancer cases
Lymph node negative cases (16)
Lymph node positive cases (19)
58 ± 14
54 ± 14
Tumour size (mm)
32 ± 12
32 ± 21
No. of positive lymph nodes
-Micrometastases (<2 mm) (No LN available)
Normal adjacent tissue
Extraction of RNA
Total RNA was extracted using the miRNeasy FFPE kit (Qiagen, Doncaster, VIC, Australia). RNA was quantified using the Quant-it RiboGreen RNA Assay kit (Invitrogen, Mulgrave, VIC, Australia) and purity assessed by A260/A280 and A260/230 ratios (>1.8) using the Nanodrop. The RNA integrity of selected samples was analysed using the 2100 Bioanalyser and the RNA 6000 Nano kit (Agilent Technologies, Mulgrave, VIC, Australia).
100 ng of total RNA was dephosphorylated and directly labelled with Cy3 using the miRNA Complete Labelling and Hyb Kit (Agilent Technologies). Labelled RNA was hybridised to Human miRNA microarrays (Sanger Release 14.0) according to the manufacturers’ instructions (Agilent Technologies) and scanned on an Agilent High-resolution C scanner.
miRNA array analysis
These miRNA array results have been deposited in Gene Expression Omnibus (GEO) with Accession No. GSE38167. Data from 15,000 probe features representing 904 unique miRNAs was extracted using Agilent Feature Extraction software (v10.7.3.1) and converted to background subtracted signal intensities. The extracted data was imported into Genespring GX (Agilent Technologies) where it was log2 transformed and median normalised. Seven samples did not meet quality control measurements recommended by the manufacturer and were removed from the microarray analysis (1 IDC from a LN negative patient, 1 NAT and 3 IDCs from LN positive patients; and 2 LN metastases). Of all the probes interrogated in this analysis, those corresponding to 570 miRNA transcripts were present at a signal intensity threshold above background in at least one of the tissue samples. Unpaired t-tests were used to identify miRNAs with significantly altered expression (>2-fold, p < 0.05). To correct for false positive results, a Benjamini and Hochberg False Discovery Rate (FDR) of 5.0% was used for multiple testing. Supervised hierarchical cluster analysis was performed on miRNAs that were found to be significantly different (>2 fold, p < 0.05, FDR < 0.05). Similarity in the expression patterns between miRNAs was measured by Pearson’s correlation coefficient.
Biological targets of differentially expressed miRNAs were identified by searching for the presence of conserved 8mer and 7mer sites within genes that match the seed region of each miRNA. Non-conserved sites were also included in this analysis. This analysis was performed using sRNA Target Base (starBase, http://starbase.sysu.edu.cn/index.php)  which integrates data from 21 Ago or TNRC6 CLIP-Seq sequence data sets with the target prediction programs Target Scan, Pictar and miRanda. The number of genes identified by each of these programs and those in common between each of the programs is shown in Additional file 1: Figure S2. Only those target genes that were predicted by all three target prediction programs (560 genes) were used for further pathway analysis. PANTHER  was used to annotate the biological pathways that predicted miRNA target genes were involved in as previously described . Pathways with a p-value <0.05 were considered to be significantly regulated by the miRNAs.
To determine the significance of differentially expressed miRNA families and clustered miRNAs; or to annotate the function of differentially expressed miRNAs, a freely available web-based resource, TAM (tool for annotations of miRNAs), was used .
Semi-quantitative real-time PCR
Total RNA (5 ng) was reverse transcribed to generate cDNA using the Taqman MicroRNA Reverse Transcription Kit and Megaplex RT Human Primer Pools Set v3.0 (Applied Biosystems, Mulgrave, VIC, Australia) according to the manufacturers’ instructions. cDNA was amplified using the TaqMan PreAmp Master Mix with Megaplex Human PreAmp Primer Pools Set v3.0 (Applied Biosystems) according to the manufacturers’ instructions. Real-time PCR analysis was performed in triplicate on all samples (Table 1) using TaqMan Universal PCR mix No AmpErase UNG and TaqMan MicroRNA Assays (Applied Biosystems) according to the manufacturers’ instructions, with results quantified on a 7500 real-time PCR system (Applied Biosystems). The expression of the following miRNAs was analysed for array validation: hsa-let-7a (Assay ID: 377), hsa-let-7b (Assay ID: 2619), hsa-let-7c (Assay ID: 379), hsa-miR-100 (Assay ID: 437), hsa-miR-101 (Assay ID: 2253), hsa-miR-126* (Assay ID:451), hsa-miR-26a (Assay ID: 405), hsa-miR-26b (Assay ID: 407), hsa-miR-130a (Assay ID: 454), hsa-miR-29c (Assay ID: 587), hsa-miR-205 (Assay ID: 509), miR-210 (Assay ID: 512), RNU44 (Assay ID: 1094) and RNU49 (Assay ID: 1005). The relative miRNA expression was calculated by normalising the miRNA of interest to RNU44 (2-ΔCt). Relative expression of miRNAs was also calculated using a second normaliser, RNU49 (data not shown). The same relative expression patterns for the miRs analysed was observed when normalised with RNU49 (examples shown in Additional file 1: Figure S3). The relative expression of RNU44 and RNU49 was not significantly different between the subgroups analysed and therefore served as an appropriate normaliser for this analysis (Additional file 1: Figure S3). All pre-amplified multiplex miRNA assays were validated against uniplex miRNA assays to verify that the multiplex reaction did not affect miRNA quantitation (Additional file 1: Supplementary Methods, Table S1 and Figure S4). hsa-miR-126* and hsa-miR-205 were outside our range of acceptable PCR efficiencies and were not used for further validations.
The normality of the data distribution was tested using a D’Agostino and Pearson Omnibus test. The values were found not to have been sampled from a Gaussian distribution and thus, non-parametric statistical tests were used to compare the real-time PCR data. A two-tailed Mann–Whitney U test was used to determine if there was a statistically significant difference in the expression of miRNAs between any two subgroups. The Kruskal-Wallis rank test followed by a Dunn’s Multiple Correction test was used to determine the statistical significance of miRNA expression between multiple (>2) subgroups. Analysis of the correlation between miRNA expression and clinical parameters was performed using Spearman’s correlation test. All analysis was performed using GraphPad Prism (version 5.04) software (GraphPad software Inc., La Jolla, CA, USA).
miRNAs differentially expressed in TNBC
Unique miRNAs identified as being differentially expressed in triple negative breast cancer
miRNAs differentially expressed in LN negative IDCs compared to LN positive IDCs
To determine miRNAs that were differentially expressed in patients with LN metastases, the primary breast cancers from 15 LN negative patients and 16 LN positive patients were compared. In this analysis, there were no miRNAs that were found to be significantly different. Additionally, the comparison of NAT from node negative and node positive women yielded no significant differences in miRNA profiles.
miRNAs associated with LN metastasis
Of note, the five unique miRNAs that were differentially regulated when all tumours were compared to all normal tissues (Table 2), were not commonly regulated when lymph node positive and negative cases were compared separately to matched normal tissue. Instead, the miRs that were positively regulated (miR-17*- 3.03, miR-590-5p- 2.15, miR-1280- 2.10) when all tumours were analysed became more highly up-regulated (miR-17*- 4.66, miR-590-5p- 3.13, miR-1280- 2.40) in the comparison of lymph node negative tumours versus matched normal tissue, but were not regulated in lymph node positive tumours. In contrast, while miR-130a was down-regulated in the comparison of all tumour tissues (−2.21), it became more strongly down-regulated in lymph node positive tumours (−3.32) when they were compared to matched normal tissues, and was not regulated in lymph node negative cases. This supports our data, that there are intrinsic differences between lymph node positive and negative TNBCs: miRs are mainly down-regulated in patients with lymph node metastases, and up-regulated in patients who do not have lymph node metastases and suggests that it is these patients (i.e. lymph node positive or negative tumours) that are driving the differential regulation (negative and positive respectively) when all patients are combined.
To determine miRNAs that were altered in LN metastases, the miRNA profile of LN metastases were compared to those of NAT (9 matched and 5 unmatched) from LN positive patients. From this analysis, 63 miRNAs were found to be differentially expressed (Additional file 1: Table S6). The majority of these miRNAs (58/63 miRNAs) were down-regulated and largely overlapped with the profile of their primary tumour (IDC, LN+; 32/63 miRNAs, Figure 2C). Many of the miRNAs overlap with those identified by Cascione et al., with 33% (21/63) of the differentially expressed miRNAs identified by our study also showing differential regulation in normal versus LN metastases comparisons in their study of TNBC . No miRNAs were differentially expressed when LN positive primary IDCs and LN metastases were compared. This suggests that the profile of miRNAs differentially expressed in primary breast cancers are maintained in LN metastases from these tumours. This is also concordant with the results of Cascione et al. who found very few differentially expressed miRNAs in this comparison . In contrast, only five of the 63 miRNAs that were differentially expressed in LN metastases were also differentially expressed in LN negative primary IDC’s (Figure 2C) when compared to NAT, suggesting that their miRNA profiles are largely distinct.
miRNAs differentially expressed in triple negative breast cancer patients with lymph node metastases
Lymph node positive cases
Lymph node negative cases
IDC v NAT
LN met v NAT
IDC v NAT
Biological functions and pathways regulated by the 27 miRNAs associated with LN metastasis
To annotate the function of the 27 miRNAs, we performed TAM analysis . Through this analysis, the majority of miRNAs were found to function as tumour suppressor miRs (13/27 miRNAs) or were involved in human embryonic stem cell regulation (11/27 miRNAs) and these functions were significantly over-represented. A number of miRNAs in this list are also known to have functions in cell death and/or the cell cycle and these functions were also significantly over-represented (Additional file 1: Table S7).
Predicted signalling pathways regulated by 27 differentially expressed miRNAs in TNBCs with LN metastases
Homo sapiens reference genome (No. of genes)
miRNA target genes (No. of genes)
Over/under represented (+/−)
Wnt signaling pathway
TGF-beta signaling pathway
p53 pathway by glucose deprivation
p53 pathway feedback loops
PI3 kinase pathway
Heterotrimeric G-protein signaling pathway-Gi alpha and Gs alpha mediated pathway
EGF receptor signaling pathway
De novo pyrmidine ribonucleotides biosythesis
Correlation between miRNA expression and clinical parameters
Correlation between miRNA expression levels and clinical variables in triple negative breast cancers
Age at diagnosis
Tumour size (mm)
% + ve lymph nodes
This study has investigated miRNA profiles in TNBC cases to determine miRNAs whose expression is associated with the initiation of this breast cancer subtype and those associated with its metastasis to the LN. We have identified a panel of 27 miRNAs that are associated with breast cancer progression in TNBC.
miRNAs differentially expressed in TNBC
By comparing the miRNA profiles of primary breast tumours to matched normal adjacent breast tissue, 71 miRNAs were identified as being differentially expressed. The majority of the miRNAs identified in this analysis have previously been implicated in breast cancer, suggesting that these miRNAs have a general function in the tumourigenic process in all breast cancers and that they are not specific to TNBC. This is perhaps not surprising, given the weak separation of the TNBC subtype by hierarchical clustering in miRNA expression profiles when compared to other breast tumour subtypes . The results of this study are highly concordant with that of Cascione et al., with 70% (50/71) of the differentially expressed miRNAs identified by our study also showing differential regulation in normal versus tumour comparisons in their study of TNBC . Differentially expressed miRNAs included members of the miR-17-92 oncogenic cluster, the let-7 family and the miR-200 family which are known to be dysregulated in several solid tumours, including breast cancer [29, 30]. Of particular note, the miR-17-92 cluster was recently shown to be up-regulated in TNBC and the basal subtype [10, 31] while miR-145 and miR-199b-5p were strongly repressed  consistent with our results. miR-210 and miR-21, well known oncomiRs [7, 11] were found to be over-expressed in our study. We also observed that miRNAs known to be bona fide regulators of ERα (miR-18a, miR-18b) were over-expressed ; while those that have been previously reported to be differentially expressed between ER-positive and -negative breast cancers (including let-7b, miR-200a, miR-21, miR-25, miR-106b) were significantly altered in this study .
The immunohistochemical and molecular profiles of TNBCs are similar to that of hereditary breast cancers that have mutations in the BRCA1 gene . In this regard, several miRNAs identified in this study are known to participate in the BRCA1 signalling axis including miR-146a, miR-155 and miR-335 [33–35]. Deregulation of these miRNAs in TNBC may contribute to altered BRCA1 signalling and could partly explain the similarities of these tumours with those in which BRCA1 function is lost.
miRNAs involved in LN positive TNBC
A significant and important finding from this study is that the miRNA profiles of LN positive primary breast cancers were strikingly distinct from that of LN negative patients. In particular, there was an overall up-regulation of miRNAs in LN negative patients and a dampening of miRNA expression in LN positive patients, suggesting that oncogenic miRNAs are associated with the development of LN negative TNBC and in contrast, that deregulated expression of tumour suppressor miRs is involved in LN positive disease. Two enzymes, Drosha and Dicer, are pivotal in the processing of pri-miRNA into mature double stranded miRNA fragments . Interestingly, in breast cancer, reduced expression of Dicer has been associated with shorter metastasis-free survival and with the TNBC subtype, where it is observed in 60-78% of patients [36–38]. The correlation of Dicer expression with LN metastases in TNBC patients was not examined in these studies [36, 38]. It could be hypothesised that the overall down-regulation of miRNAs observed in LN positive TNBCs in the current study is a result of reduced Dicer expression, but this remains to be determined.
miRNAs as markers for metastasis
In this study, we were not seeking to identify miRNAs that were differentially expressed in the transition of breast cancer progression i.e. primary breast cancer to metastasis or from ductal carcinoma in situ (DCIS) to IDC, as previous studies have done [10, 27, 28, 39]. In contrast to other studies, our analysis yielded no discriminatory miRNAs when matched primary breast cancer and LN metastases were compared suggesting that they are remarkably similar in their miRNA profiles and supporting the validity of our approach. The premise of the current study was to determine altered miRNA expression patterns in primary breast cancers that were also present in LN metastases. We reasoned that deregulated expression of key miRNAs that promote a highly invasive phenotype would be an early event in breast cancer progression and that these changes would be present in both the primary and metastatic lesion. Using this approach, we identified a panel of 27 miRNAs that are associated with metastasis to the LN. Given that the majority of these miRNAs are down-regulated, their value as prognostic markers remains to be determined and is the subject of ongoing investigations. However, 17/27 (63%) miRNAs are either known to be tumour suppressor miRs and/or have been shown to play a role in human embryonic stem cell regulation (Additional file 1: Table S7), supporting that deregulated expression of this panel of miRNAs is likely to result in altered differentiation states and increased invasiveness. In this regard, let-7 family members, miR-34a and miR-205 are markedly reduced in breast tumour initiating cells and are involved in epithelial to mesenchymal transition [21, 40, 41]. In addition, several of the miRNAs in this list have already been shown to have prognostic value in breast cancer including miR-210, miR-126, miR-26a, miR-125b, miR-205 and miR-214 [7–10, 42–44]. Furthermore, our results are in agreement with Cascione et al., who found that both miR-125b and miR-497 were down-regulated and related to survival in TNBC .
Within our panel of metastasis-related miRNAs are a number of miRNAs whose functional role in breast cancer progression has not previously been described. These include miR-320c, miR-29c, miR-130a and miR-195 among others. There are no reports on miR-320c function to date and thus its possible role in promoting LN metastasis in TNBC remains unknown. miR-29c has previously been reported to be up-regulated in a small cohort of breast cancer cases , however, we found miR-29c to be down-regulated and associated with LN metastasis. In agreement with our study, miR-29c is down-regulated in chronic lymphocytic leukaemia and a range of solid tumours including colorectal cancer where it has been reported to decrease with cancer progression and is a predictor of survival and early recurrence [45, 46]. miR-130a is down-regulated in a panel of ovarian cancer and lung cancer cell lines that are resistant to chemotherapies [23, 24]. Moreover, it has recently been shown to be decreased in prostate cancer and its over-expression in prostate cancer cell lines causes the repression of key oncogenic pathways such as the MAPK pathway . miR-195 was markedly down-regulated in our study. Its re-expression in breast cancer cell lines has recently been reported to reduce cellular proliferation and invasion, suggesting that it plays a key role in breast cancer progression .
A recent study by Buffa et al., described three miRNAs (miR-342, miR-27b, miR-150) that were prognostic for relapse-free survival in TNBC . In our hands, these results were not replicated for LN metastasis and moreover, these miRNAs were not differentially expressed when breast tumour was compared to normal tissue. In addition, miR-10b has been shown to be strongly expressed in metastatic breast cancer cells, where it regulates invasion and migration . However, we observed miR-10b to be strongly down-regulated in lymph-node positive and negative primary breast cancers and in LN metastases, suggesting it does not contribute to breast cancer invasiveness, and this is in strong agreement with several recent studies [6, 10, 20].
We used starBASE to predict signalling pathways regulated by the 27 miRNAs associated with LN metastasis. Interestingly, the pathways identified by this analysis have been widely reported to be deregulated in TNBC. In particular Wnt/β-Catenin, EGFR1 and the TGFβ receptor signalling axis are all known to be highly expressed in a large proportion of cases within the TNBC subtype and have been proposed as alternative treatment targets for TNBC [2, 3, 51]. A number of chemotherapeutic agents which target these pathways are currently under investigation in patients with TNBC, in both the neoadjuvant and metastatic setting . The findings from our study- decreased expression of key miRNAs which target these pathways- provides a biological mechanism for up-regulation of these pathways in this breast cancer subtype. Additionally, TP53, PIK3CA and EGFR are among the four most commonly mutated genes associated with aberrant expression of interacting genes in TNBC , exemplifying their importance in this subtype. The predicted regulation of these pathways by miRNAs, as found in our study, reveals another mechanism for disruption of these pathways in TNBC.
A disadvantage of our study design is that samples were acquired by punch biopsy, rather than isolating individual normal breast and tumour epithelial cells by laser capture microdissection. We recognise that tissue heterogeneity and contamination by other non-neoplastic cell types including adipocytes, stroma and lymphocytes was unavoidable and that it may have contributed to the miRNA profiles generated. Our approach of selecting miRNAs that were only differentially expressed in primary tumours and matched LN metastases and that were not altered in LN negative primary tumours would have culled a large proportion of those miRNAs that are specifically altered in inflammatory responses. Moreover, those miRNAs that were differentially expressed simply because of differences in cell type composition between normal tissue (higher content of adipocytes and stroma) when compared to tumour tissue (higher content of epithelium) would have shown overlapping miRNA profiles when both LN negative and LN positive primary tumours were compared to matched normal tissues (10 miRNAs, Figure 2C) and were not included in our list of miRNAs related to metastasis (Table 3). Therefore, because of this careful process of elimination, we consider the miRNAs identified in this study (Table 3) are strongly representative of the invasive capacity of the primary tumour.
Given that the sample size used in this study is relatively small, conclusions cannot be drawn about the relationship of these miRNAs to prognosis or their specificity to the TNBC subtype. However, the tumours are homogenous with regards to size, hormone receptor status and histological grade. Furthermore, the size of this cohort (n = 35) is in line with several recent studies investigating the prognostic value of miRNAs in TNBC (Buffa et al., n = 37; Farazi et al., n = 48; Radojocic et al., n = 49) [20, 28, 49] and represents the only study to date that has analysed the full repertoire of miRNAs in lymph node positive and negative TNBC samples separately compared to matched normal adjacent tissue. The uniqueness of our study design was highly advantageous for delineating critical differences in miRNA profiles between LN positive and negative primary tumours, which would have been missed if: 1) matched normal adjacent tissue was not used and 2) only primary tumours from LN positive versus LN negative patients were compared to one another. In a recent study, Farazi and colleagues showed that there was no difference in miRNA profiles of primary breast tumours between TNBC patients who relapsed when compared to those who did not . Indeed, we also found that no miRNAs were differentially expressed when LN negative versus positive primary tumours were directly compared using our statistical cut-offs. Further to this, in our study no miRNAs were found to be differentially expressed when normal tissue from LN positive patients was compared to that of LN negative patients. As such, the distinct miRNA profiles generated when LN negative and LN positive primary tumours were each separately compared to matched normal adjacent tissues are likely to reflect small, but critical changes in the regulation of miRNAs that are essential for the development of each of these respective phenotypes. Given that any LN involvement (regardless of number) is associated with worse disease-free and overall survival in TNBC , identification of molecular determinants such as miRNAs, that contribute to this highly invasive phenotype is urgently needed.
This study has provided novel insight into the repertoire of miRNAs that contribute to the initiation of and progression to LN metastasis in TNBC. These miRNAs may serve as markers for metastasis or treatment targets in the future.
Micro ribonucleic acid
Human epidermal growth factor receptor 2
Triple negative breast cancer
Messenger ribonucleic acid
Invasive ductal carcinoma
Normal adjacent tissue
Complementary deoxyribonucleic acid.
The authors would like to thank Dr. Ricardo Vilain for pathological review of all tumour and normal tissue specimens used in this analysis and Ms. Tina Hope for assistance with archival specimens. This work was supported by funding from the National Breast Cancer Foundation and the Hunter Medical Research Institute. Dr. Avery-Kiejda is supported by a Hunter Translational Cancer Research Unit Fellowship from the Cancer Institute NSW.
- Kamangar F, Dores GM, Anderson WF: Patterns of cancer incidence, mortality, and prevalence across five continents: defining priorities to reduce cancer disparities in different geographic regions of the world. J Clin Oncol. 2006, 24 (14): 2137-2150. 10.1200/JCO.2005.05.2308.View ArticlePubMedGoogle Scholar
- Podo F, Buydens LM, Degani H, Hilhorst R, Klipp E, Gribbestad IS, Van Huffel S, van Laarhoven HW, Luts J, Monleon D, Postma GJ, Schneiderhan-Marra N, Santoro F, Wouters H, Russnes HG, Sorlie T, Tagliabue E, Borresen-Dale AL: Triple-negative breast cancer: present challenges and new perspectives. Mol Oncol. 2010, 4 (3): 209-229. 10.1016/j.molonc.2010.04.006.View ArticlePubMedGoogle Scholar
- Carey L, Winer E, Viale G, Cameron D, Gianni L: Triple-negative breast cancer: disease entity or title of convenience?. Nat Rev Clin Oncol. 2010, 7 (12): 683-692. 10.1038/nrclinonc.2010.154.View ArticlePubMedGoogle Scholar
- Jackson RJ, Standart N: How do microRNAs regulate gene expression?. Sci STKE. 2007, 2007 (367): re1-View ArticlePubMedGoogle Scholar
- Andorfer CA, Necela BM, Thompson EA, Perez EA: MicroRNA signatures: clinical biomarkers for the diagnosis and treatment of breast cancer. Trends Mol Med. 2011, 17 (6): 313-319. 10.1016/j.molmed.2011.01.006.View ArticlePubMedGoogle Scholar
- 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 (16): 7065-7070. 10.1158/0008-5472.CAN-05-1783.View ArticlePubMedGoogle Scholar
- Foekens JA, Sieuwerts AM, Smid M, Look MP, de Weerd V, Boersma AW, Klijn JG, Wiemer EA, Martens JW: Four miRNAs associated with aggressiveness of lymph node-negative, estrogen receptor-positive human breast cancer. Proc Natl Acad Sci USA. 2008, 105 (35): 13021-13026. 10.1073/pnas.0803304105.View ArticlePubMedPubMed CentralGoogle Scholar
- Le Quesne JL, Jones J, Warren J, Dawson SJ, Ali R, Bardwell H, Blows F, Pharoah P, Caldas C: Biological and prognostic associations of miR-205 and let-7b in breast cancer revealed by in situ hybridisation analysis of micro-RNA expression in arrays of archival tumour tissue. J Pathol. 2012, 227 (3): 306-314. 10.1002/path.3983.View ArticleGoogle Scholar
- Rothe F, Ignatiadis M, Chaboteaux C, Haibe-Kains B, Kheddoumi N, Majjaj S, Badran B, Fayyad-Kazan H, Desmedt C, Harris AL, Piccart M, Sotiriou C: Global microRNA expression profiling identifies MiR-210 associated with tumor proliferation, invasion and poor clinical outcome in breast cancer. PLoS One. 2011, 6 (6): e20980-10.1371/journal.pone.0020980.View ArticlePubMedPubMed CentralGoogle Scholar
- Volinia S, Galasso M, Sana ME, Wise TF, Palatini J, Huebner K, Croce CM: Breast cancer signatures for invasiveness and prognosis defined by deep sequencing of microRNA. Proc Natl Acad Sci USA. 2012, 109 (8): 3024-3029. 10.1073/pnas.1200010109.View ArticlePubMedPubMed CentralGoogle Scholar
- Yan LX, Huang XF, Shao Q, Huang MY, Deng L, Wu QL, Zeng YX, Shao JY: MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA. 2008, 14 (11): 2348-2360. 10.1261/rna.1034808.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn 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 Nat Acad Sci USA. 2001, 98 (19): 10869-10874. 10.1073/pnas.191367098.View ArticlePubMedPubMed CentralGoogle Scholar
- Carter CL, Allen C, Henson DE: Relation of tumor size, lymph node status, and survival in 24,740 breast cancer cases. Cancer. 1989, 63 (1): 181-187. 10.1002/1097-0142(19890101)63:1<181::AID-CNCR2820630129>3.0.CO;2-H.View ArticlePubMedGoogle Scholar
- Hernandez-Aya LF, Chavez-Macgregor M, Lei X, Meric-Bernstam F, Buchholz TA, Hsu L, Sahin AA, Do KA, Valero V, Hortobagyi GN, Gonzalez-Angulo AM: Nodal status and clinical outcomes in a large cohort of patients with triple-negative breast cancer. J Clin Oncol. 2011, 29 (19): 2628-2634. 10.1200/JCO.2010.32.1877.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH: starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011, 39 (Database issue): D202-D209.View ArticlePubMedGoogle Scholar
- PANTHER Classification system. http://www.pantherdb.org,
- Avery-Kiejda KA, Bowden NA, Croft AJ, Scurr LL, Kairupan CF, Ashton KA, Talseth-Palmer BA, Rizos H, Zhang XD, Scott RJ, Hersey P: P53 in human melanoma fails to regulate target genes associated with apoptosis and the cell cycle and may contribute to proliferation. BMC Cancer. 2011, 11 (1): 203-10.1186/1471-2407-11-203.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu M, Shi B, Wang J, Cao Q, Cui Q: TAM: a method for enrichment and depletion analysis of a microRNA category in a list of microRNAs. BMC Bioinformatics. 2010, 11: 419-10.1186/1471-2105-11-419.View ArticlePubMedPubMed CentralGoogle Scholar
- Altuvia Y, Landgraf P, Lithwick G, Elefant N, Pfeffer S, Aravin A, Brownstein MJ, Tuschl T, Margalit H: Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33 (8): 2697-2706. 10.1093/nar/gki567.View ArticlePubMedPubMed CentralGoogle Scholar
- Radojicic J, Zaravinos A, Vrekoussis T, Kafousi M, Spandidos DA, Stathopoulos EN: MicroRNA expression analysis in triple-negative (ER, PR and Her2/neu) breast cancer. Cell Cycle. 2011, 10 (3): 507-517. 10.4161/cc.10.3.14754.View ArticlePubMedGoogle Scholar
- Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, Vadas MA, Khew-Goodall Y, Goodall GJ: The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol. 2008, 10 (5): 593-601. 10.1038/ncb1722.View ArticlePubMedGoogle Scholar
- Majid S, Dar AA, Saini S, Shahryari V, Arora S, Zaman MS, Chang I, Yamamura S, Chiyomaru T, Fukuhara S, Tanaka Y, Deng G, Tabatabai ZL, Dahiya R: MicroRNA-1280 inhibits invasion and metastasis by targeting ROCK1 in bladder cancer. PLoS One. 2012, 7 (10): e46743-10.1371/journal.pone.0046743.View ArticlePubMedPubMed CentralGoogle Scholar
- Acunzo M, Visone R, Romano G, Veronese A, Lovat F, Palmieri D, Bottoni A, Garofalo M, Gasparini P, Condorelli G, Chiariello M, Croce CM: miR-130a targets MET and induces TRAIL-sensitivity in NSCLC by downregulating miR-221 and 222. Oncogene. 2012, 31 (5): 634-642.PubMedGoogle Scholar
- Sorrentino A, Liu CG, Addario A, Peschle C, Scambia G, Ferlini C: Role of microRNAs in drug-resistant ovarian cancer cells. Gynecol Oncol. 2008, 111 (3): 478-486. 10.1016/j.ygyno.2008.08.017.View ArticlePubMedGoogle Scholar
- Shan X, Miao Y, Fan R, Qian H, Chen P, Liu H, Yan X, Li J, Zhou F: MiR-590-5P Inhibits Growth of HepG2 Cells via Decrease of S100A10 Expression and Inhibition of the Wnt Pathway. Int J Mol Sci. 2013, 14 (4): 8556-8569. 10.3390/ijms14048556.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiao X, Tang C, Xiao S, Fu C, Yu P: Enhancement of proliferation and invasion by MicroRNA-590-5p via targeting PBRM1 in clear cell renal carcinoma cells. Oncol Res. 2013, 20 (11): 537-544.View ArticlePubMedGoogle Scholar
- Cascione L, Gasparini P, Lovat F, Carasi S, Pulvirenti A, Ferro A, Alder H, He G, Vecchione A, Croce CM, Shapiro CL, Huebner K: Integrated microRNA and mRNA signatures associated with survival in triple negative breast cancer. PLoS One. 2013, 8 (2): e55910-10.1371/journal.pone.0055910.View ArticlePubMedPubMed CentralGoogle Scholar
- Farazi TA, Horlings HM, Ten Hoeve JJ, Mihailovic A, Halfwerk H, Morozov P, Brown M, Hafner M, Reyal F, van Kouwenhove M, Kreike B, Sie D, Hovestadt V, Wessels LF, van de Vijver MJ, Tuschl T: MicroRNA sequence and expression analysis in breast tumors by deep sequencing. Cancer Res. 2011, 71 (13): 4443-4453. 10.1158/0008-5472.CAN-11-0608.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 (10): R214-10.1186/gb-2007-8-10-r214.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 Nat Acad Sci USA. 2006, 103 (7): 2257-2261. 10.1073/pnas.0510565103.View ArticlePubMedPubMed CentralGoogle Scholar
- Enerly E, Steinfeld I, Kleivi K, Leivonen SK, Aure MR, Russnes HG, Ronneberg JA, Johnsen H, Navon R, Rodland E, Makela R, Naume B, Perala M, Kallioniemi O, Kristensen VN, Yakhini Z, Borresen-Dale AL: miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors. PLoS One. 2011, 6 (2): e16915-10.1371/journal.pone.0016915.View ArticlePubMedPubMed CentralGoogle Scholar
- Klinge CM: miRNAs and estrogen action. Trends Endocrinol Metab. 2012, 23 (5): 223-233. 10.1016/j.tem.2012.03.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Garcia AI, Buisson M, Bertrand P, Rimokh R, Rouleau E, Lopez BS, Lidereau R, Mikaelian I, Mazoyer S: Down-regulation of BRCA1 expression by miR-146a and miR-146b-5p in triple negative sporadic breast cancers. EMBO Mol Med. 2011, 3 (5): 279-290. 10.1002/emmm.201100136.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang S, Wang RH, Akagi K, Kim KA, Martin BK, Cavallone L, Haines DC, Basik M, Mai P, Poggi E, Isaacs C, Looi LM, Mun KS, Greene MH, Byers SW, Teo SH, Deng CX, Sharan SK: Tumor suppressor BRCA1 epigenetically controls oncogenic microRNA-155. Nat Med. 2011, 17 (10): 1275-1282. 10.1038/nm.2459.View ArticlePubMedPubMed CentralGoogle Scholar
- Heyn H, Engelmann M, Schreek S, Ahrens P, Lehmann U, Kreipe H, Schlegelberger B, Beger C: MicroRNA miR-335 is crucial for the BRCA1 regulatory cascade in breast cancer development. Int J Cancer. 2011, 129 (12): 2797-2806. 10.1002/ijc.25962.View ArticlePubMedGoogle Scholar
- Dedes KJ, Natrajan R, Lambros MB, Geyer FC, Lopez-Garcia MA, Savage K, Jones RL, Reis-Filho JS: Down-regulation of the miRNA master regulators Drosha and Dicer is associated with specific subgroups of breast cancer. Eur J Cancer. 2011, 47 (1): 138-150. 10.1016/j.ejca.2010.08.007.View ArticlePubMedGoogle Scholar
- Grelier G, Voirin N, Ay AS, Cox DG, Chabaud S, Treilleux I, Leon-Goddard S, Rimokh R, Mikaelian I, Venoux C, Puisieux A, Lasset C, Moyret-Lalle C: Prognostic value of Dicer expression in human breast cancers and association with the mesenchymal phenotype. Br J Cancer. 2009, 101 (4): 673-683. 10.1038/sj.bjc.6605193.View ArticlePubMedPubMed CentralGoogle Scholar
- Passon N, Gerometta A, Puppin C, Lavarone E, Puglisi F, Tell G, Di Loreto C, Damante G: Expression of Dicer and Drosha in triple-negative breast cancer. J Clin Pathol. 2012, 65 (4): 320-326. 10.1136/jclinpath-2011-200496.View ArticlePubMedGoogle Scholar
- Baffa R, Fassan M, Volinia S, O'Hara B, Liu CG, Palazzo JP, Gardiman M, Rugge M, Gomella LG, Croce CM, Rosenberg A: MicroRNA expression profiling of human metastatic cancers identifies cancer gene targets. J Pathol. 2009, 219 (2): 214-221. 10.1002/path.2586.View ArticlePubMedGoogle Scholar
- Kim NH, Kim HS, Li XY, Lee I, Choi HS, Kang SE, Cha SY, Ryu JK, Yoon D, Fearon ER, Rowe RG, Lee S, Maher CA, Weiss SJ, Yook JI: A p53/miRNA-34 axis regulates Snail1-dependent cancer cell epithelial-mesenchymal transition. J Cell Biol. 2011, 195 (3): 417-433. 10.1083/jcb.201103097.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C, Huang Y, Hu X, Su F, Lieberman J, Song E: let-7 regulates self renewal and tumorigenicity of breast cancer cells. Cell. 2007, 131 (6): 1109-1123. 10.1016/j.cell.2007.10.054.View ArticlePubMedGoogle Scholar
- Zhang Y, Yan LX, Wu QN, Du ZM, Chen J, Liao DZ, Huang MY, Hou JH, Wu QL, Zeng MS, Huang WL, Zeng YX, Shao JY: miR-125b is methylated and functions as a tumor suppressor by regulating the ETS1 proto-oncogene in human invasive breast cancer. Cancer Res. 2011, 71 (10): 3552-3562. 10.1158/0008-5472.CAN-10-2435.View ArticlePubMedGoogle Scholar
- Jansen MP, Reijm EA, Sieuwerts AM, Ruigrok-Ritstier K, Look MP, Rodriguez-Gonzalez FG, Heine AA, Martens JW, Sleijfer S, Foekens JA, Berns EM: High miR-26a and low CDC2 levels associate with decreased EZH2 expression and with favorable outcome on tamoxifen in metastatic breast cancer. Breast Cancer Res Treat. 2012, 133 (3): 937-947. 10.1007/s10549-011-1877-4.View ArticlePubMedGoogle Scholar
- Schwarzenbach H, Milde-Langosch K, Steinbach B, Muller V, Pantel K: Diagnostic potential of PTEN-targeting miR-214 in the blood of breast cancer patients. Breast Cancer Res Treat. 2012, 134 (3): 933-941. 10.1007/s10549-012-1988-6.View ArticlePubMedGoogle Scholar
- Kuo TY, Hsi E, Yang IP, Tsai PC, Wang JY, Juo SH: Computational analysis of mRNA expression profiles identifies microRNA-29a/c as predictor of colorectal cancer early recurrence. PLoS One. 2012, 7 (2): e31587-10.1371/journal.pone.0031587.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatopoulos B, Meuleman N, Haibe-Kains B, Saussoy P, Van Den Neste E, Michaux L, Heimann P, Martiat P, Bron D, Lagneaux L: microRNA-29c and microRNA-223 down-regulation has in vivo significance in chronic lymphocytic leukemia and improves disease risk stratification. Blood. 2009, 113 (21): 5237-5245. 10.1182/blood-2008-11-189407.View ArticlePubMedGoogle Scholar
- Boll K, Reiche K, Kasack K, Morbt N, Kretzschmar AK, Tomm JM, Verhaegh G, Schalken J, von Bergen M, Horn F, Hackermuller J: MiR-130a, miR-203 and miR-205 jointly repress key oncogenic pathways and are downregulated in prostate carcinoma. Oncogene. 2013, 32 (3): 277-285. 10.1038/onc.2012.55.View ArticlePubMedGoogle Scholar
- Li D, Zhao Y, Liu C, Chen X, Qi Y, Jiang Y, Zou C, Zhang X, Liu S, Wang X, Zhao D, Sun Q, Zeng Z, Dress A, Lin MC, Kung HF, Rui H, Liu LZ, Mao F, Jiang BH, Lai L: Analysis of MiR-195 and MiR-497 expression, regulation and role in breast cancer. Clin Cancer Res. 2011, 17 (7): 1722-1730. 10.1158/1078-0432.CCR-10-1800.View ArticlePubMedGoogle Scholar
- Buffa FM, Camps C, Winchester L, Snell CE, Gee HE, Sheldon H, Taylor M, Harris AL, Ragoussis J: microRNA-associated progression pathways and potential therapeutic targets identified by integrated mRNA and microRNA expression profiling in breast cancer. Cancer Res. 2011, 71 (17): 5635-5645. 10.1158/0008-5472.CAN-11-0489.View ArticlePubMedGoogle Scholar
- Ma L, Teruya-Feldstein J, Weinberg RA: Tumour invasion and metastasis initiated by microRNA-10b in breast cancer. Nature. 2007, 449 (7163): 682-688. 10.1038/nature06174.View ArticlePubMedGoogle Scholar
- King TD, Suto MJ, Li Y: The Wnt/beta-catenin signaling pathway: a potential therapeutic target in the treatment of triple negative breast cancer. J Cell Biochem. 2012, 113 (1): 13-18. 10.1002/jcb.23350.View ArticlePubMedGoogle Scholar
- Rakha EA, Chan S: Metastatic triple-negative breast cancer. Clin Oncol. 2011, 23 (9): 587-600. 10.1016/j.clon.2011.03.013.View ArticleGoogle Scholar
- Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G, Bashashati A, Prentice LM, Khattra J, Burleigh A, Yap D, Bernard V, McPherson A, Shumansky K, Crisan A, Giuliany R, Heravi-Moussavi A, Rosner J, Lai D, Birol I, Varhol R, Tam A, Dhalla N, Zeng T, Ma K, Chan SK, et al: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486 (7403): 395-399.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/51/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.