- Research article
- Open Access
MicroRNA expression as risk biomarker of breast cancer metastasis: a pilot retrospective case-cohort study
BMC Cancer volume 14, Article number: 739 (2014)
MicroRNAs (miRNAs) are small, non-coding RNA molecules involved in post-transcriptional gene regulation and have recently been shown to play a role in cancer metastasis. In solid tumors, especially breast cancer, alterations in miRNA expression contribute to cancer pathogenesis, including metastasis. Considering the emerging role of miRNAs in metastasis, the identification of predictive markers is necessary to further the understanding of stage-specific breast cancer development. This is a retrospective analysis that aimed to identify molecular biomarkers related to distant breast cancer metastasis development.
A retrospective case cohort study was performed in 64 breast cancer patients treated during the period from 1998–2001. The case group (n = 29) consisted of patients with a poor prognosis who presented with breast cancer recurrence or metastasis during follow up. The control group (n = 35) consisted of patients with a good prognosis who did not develop breast cancer recurrence or metastasis. These patient groups were stratified according to TNM clinical stage (CS) I, II and III, and the main clinical features of the patients were homogeneous. MicroRNA profiling was performed and biomarkers related to metastatic were identified independent of clinical stage. Finally, a hazard risk analysis of these biomarkers was performed to evaluate their relation to metastatic potential.
MiRNA expression profiling identified several miRNAs that were both specific and shared across all clinical stages (p ≤ 0.05). Among these, we identified miRNAs previously associated with cell motility (let-7 family) and distant metastasis (hsa-miR-21). In addition, hsa-miR-494 and hsa-miR-21 were deregulated in metastatic cases of CSI and CSII. Furthermore, metastatic miRNAs shared across all clinical stages did not present high sensitivity and specificity when compared to specific-CS miRNAs. Between them, hsa-miR-183 was the most significative of CSII, which miRNAs combination for CSII (hsa-miR-494, hsa-miR-183 and hsa-miR-21) was significant and were a more effective risk marker compared to the single miRNAs.
Women with metastatic breast cancer, especially CSII, presented up-regulated levels of miR-183, miR-494 and miR-21, which were associated with a poor prognosis. These miRNAs therefore represent new risk biomarkers of breast cancer metastasis and may be useful for future targeted therapies.
Breast cancer is the most frequent tumor after skin malignancies, representing the second most common cancer-related mortality in women . Although TNM staging provides important clinical prognostic information, mammary tumors are known to be biologically heterogeneous with regard to therapeutic responses as well as molecular profiling [2, 3]. For example, these tumors can be characterized as luminal A, luminal B, basal-like, human epidermal growth factor receptor 2 (HER-2)-overexpressing or claudin-low , and this profiling provides additional molecular prognostic markers. In addition, a recent study demonstrated that the microRNA expression signature appears to provide a better characterization of cancer subtypes than transcriptional profiles and may therefore represent a new classification system for breast cancer .
MicroRNAs (miRNAs) are small, non-coding RNAs of 19–25 nt that control a wide array of physiological and pathological processes by modulating the expression of their cognate target genes through cleaving mRNA molecules or inhibiting their translation . Most cancer tissues are archived as formalin-fixed, paraffin-embedded (FFPE) samples, and microRNAs are promising biomarkers because they are a stable form of RNA. Moreover, recent studies have shown a good correlation between microRNA samples from frozen and FFPE sections [7, 8].
MiRNAs regulate key biological processes such as development, differentiation, stress response, apoptosis and proliferation [9–11] and are consequently implicated in several diseases including cancer [12–14].
Iorio et al.  demonstrated the influence of miRNA deregulation in the development of breast cancer in several tissues and lineages, and other studies have correlated miRNA profiles with mRNA subtypes, particularly with regard to estrogen receptor (ER), progesterone receptor (PR) and HER2 status [4, 14]. In addition, specific miRNAs have been associated with steps of the metastasis cascade, such as micrometastasis, local invasion, intravasation and metastatic colonization . However, in practice, few miRNA expression signatures have been shown to correlate with breast cancer metastasis, and although the epithelial to mesenchymal transition (EMT) is a critical event for metastasis from carcinomas, few studies evaluating miRNAs during EMT have been performed in breast cancer [16, 17]. In breast cancer, one example is miR-183, which have been associated with migration and invasion [18, 19]. However, the results regard to miR-183 is controversial. In a recent study of meta-analysis, comparing cancer tissues with controls of several tumor types, this miRNA presented inconsistent regulation in breast cancer . However, miR-183 expression has not been evaluated in a metastatic breast cancer context.
The main goal of this study was to identify miRNA biomarkers of breast cancer metastasis. Using a collection of FFPE samples, we selected clinically homogeneous samples, and we paired metastatic and non-metastatic patients according to tumor grade. Using this strategy, we detected, with improved precision, miRNA biomarkers that could characterize metastasis irrespective of clinical staging as well as stage-specific biomarkers.
A retrospective case cohort study  was performed in 782 patients with invasive breast cancer (ductal or lobular), without metastasis at diagnosis (clinical stage I, II and III) who had previously received treatment at Barretos Cancer Hospital between 1998 and 2001. In this study, the case group consisted of patients with a poor prognosis who developed breast cancer recurrence and/or metastasis during a follow up of ten years. The rate of cases with metastasis according each clinical stage was: 12.8% to CSI (total of 117), 25% to CSII (total of 352), 51% to CSIII (total of 313).
The control group consisted of a random sample of patients who did not develop breast cancer recurrence and/or metastasis and had a good prognosis. The control group presented the same apparent risk and length of follow up period as the case group. The rate of cases with non-metastasis according each clinical stage was: 87.2% to CSI (117 total cases), 75% to CSII (total cases 352), 49% (313 total cases).
The groups also had a similar distribution of clinical staging (CS) according to TNM classification (TNM 7th edition) . All patients were treated at the same institution with same treatment protocol and received regular follow-up assessments at the Department of Mastology and Reconstructive Surgery at Barretos Cancer Hospital, Barretos, Sao Paulo, Brazil. Patients were excluded if they had a second primary tumor, an insufficient blocked tumor or the absence of high-quality miRNA for extraction. This study was approved by Barretos Cancer Hospital ethical committee, (protocol n°362/2010).
The same pathologist reviewed all of the medical records. Immunohistochemistry evaluation was performed in all cases, including the assessment of ER and PR status and the expression of Ki-67, Her2 and cytokeratin 5/6. ER status was evaluated using the Pathway anti-Her-2 790–2991 monoclonal antibody (Ventana Medical Systems, Roche Diagnosis, Tucson, Arizona 85755, USA) at a dilution of 1:200.
PR status was evaluated using the Rabbit monoclonal antibody clone clone SP1 at a dilution of 1:600. Ki-67 expression was evaluated using the monoclonal antibody MIB-1 (Dako, Sao Paulo, Brazil) at a dilution of 1:200. Her2 status was evaluated using the mouse antibody clone 4B5 (Ventana Medical Systems, Roche Diagnosis, Tucson, Arizona 85755, USA) at a 1:2,000 dilution.
The cytokeratin 5/6 status was evaluated using the mouse monoclonal D5/16B4 (Dako, Sao Paulo, Brazil) at a 1:100 dilution. ER and PR expression was considered positive when 1% of tumor cells showed positive staining. The Ki-67 cutoff value was 14%. For Her2 semi-quantitative immunohistochemistry (2+ and 3+), the DISH test was performed using Her2 Dako K5331 (Dako, Sao Paulo, Brazil). Basal-like tumors were characterized according to triple-negative receptor status and positive C5/6 status. We applied immunohistochemistry molecular characteristic subgroups based on previous reports [22, 23].
Patients and case selection
During the paring we considered a case-cohort ratio 1:1. After using the inclusion and exclusion criteria, 64 patients were selected including 29 in the case group and 35 in the control group. Ductal invasive carcinoma represented 76.6% of the group histology. The average patient age was 53.1 years (29–95), and the average tumor size was 3.1 cm (1.0-8.5 cm).
The groups were homogeneous, and there were no differences between the variables selected prior to pairing. Table 1 summarizes the main patient features, with the associated p-values from the Fisher test using SPSS software. Using a t-test for independent variables, there were no differences between groups related to tumor size (mean 3.24 × 3.04; p = 0.592) or age (mean 49.1 × 56.4; p = 0.07). After pairing and pathologic classification, we observed that only the PR status differed between the groups (Table 1).
The mean group follow-up duration was 82.0 months (5.1-162 months); excluding the case group, this period was 120.6 months. At the end of the study, 25/29 patients in the case group had died of cancer, while the others remained alive with tumor recurrence.
Total RNA isolation from FFPE sections
FFPE samples were submitted to a total RNA isolation protocol using the Recover All TM Total Nucleic Acid Isolation kit (Life Technologies). The samples were initially treated with xylene, followed by double washing with absolute ethanol and proteinase K treatment at 50°C for 3 hours. Quantification was performed using a nanodrop ND-1000 spectrophotometer (NanoDrop Products, Wilmington, DE), and RNA quality was assessed using an Agilent Small RNA chip with a Bioanalyzer device (Agilent Technologies).
The Agilent Human miRNA Microarray (8 × 15K - G4471A, Agilent Technologies) was used in all samples from FFPE sections. Additional microarrays of frozen samples were performed in comparison to those obtained from FFPE, as quality control (Additional file 1: Figure S2). A total of 100 ng of total RNA was hybridized using miRNA complete labeling and the Hyb Kit (Agilent Technologies), according to the manufacturer’s instructions. The reactions followed a 2-step preparation, represented by dephosphorylation and denaturation of the total RNA incorporated with Cy3 fluorochrome by the T4 ligase. The next steps included standard washing procedures and hybridization with microarrays slides. The images were scanned using an Agilent DNA microarray scanner with SureScan technology (Agilent Technologies).
miRNA microarray data analysis
The raw data were obtained using Feature Extraction software v.11.0 (Agilent Technologies) and submitted to R environment v. 2.15.0  for further analysis. The median signals (gMedianSignal and gBGMedianSignal) were used. Following background subtraction and log2 scale transformation, normalization was performed using the quantile method with the aroma light package . Differentially expressed microRNAs were obtained by rank product analysis using the RankProd package , considering p-values and positive false predictions (pfp) ≤ 0.05. The rank product analyses were separated according to clinical staging between groups (CSI-CSIM, CSII-CSIIM, CSIII-CSIIIM). Differentially expressed miRNAs were further ranked according to sensitivity and specificity to determine the best candidates between non-metastatic patients and metastatic patients in a stage-specific manner. Sensitivity and specificity are defined as the number of true positive decisions/the number of actually positive cases and the number of true negative decisions/the number of actually negative cases. The area under the ROC curve (AUC) is a measure for overall performance, which can be interpreted as the average value of sensitivity for all possible values of specificity . The criteria for biomarker selection included sensitivity and specificity values ≥ 80%, as determined using the ROCR package . Differentially expressed microRNAs were clustered by Euclidian distance and average linkage using the heatmap.2 function of the gplots package .
Real-time PCR and disease-free survival analysis
The criteria for miRNA selection for further RT-qPCR confirmation, and further analysis, was the high significance (pvalue and pfp) of the miRNA, biological relevance after literature search and which is present in (1) both groups and (2) at least 1 miRNA specific of each metastatic group compared to their specific primary patient of the same CS. We performed a Venn diagrams’ to represent this selection, using gplots package .
Taqman microRNA assays (Life Technologies, Foster City, CA, USA) were used to confirm the microarray data. In brief, these reactions consist of reverse transcription with miRNA-specific primers in a real-time PCR reaction with Taqman probes. The reverse transcriptase reactions used in this study contained 10 ng of total RNA and utilized the High Capacity cDNA Reverse Transcription Kit (Life Technologies), according to the manufacturer’s instructions, in a thermocycler (Eppendorf). All real-time PCR reactions were performed in triplicate in a 7900 HT Fast Real-time PCR System (Applied Biosystems USA). All analysis procedures were performed in R environment. The normalization step was performed according to the 2-ΔΔCt method . Cycle threshold (Ct) values from selected miRNA targets were subtracted from the Ct values of the endogenous small noncoding RNA control RNU48 (Control miRNA Assay, Applied Biosystems, Foster City, CA, USA). A subsequent ΔΔCt value was calculated by subtracting metastatic ΔCt values from non-metastatic ones. The data cutoff for modulation (up/down) in each patient was estimated according to the threshold obtained after receiver operating characteristic (ROC) curve analysis using the ROCR package. To assess the combination of biomarkers, a general logistic model (glm) was performed prior to ROC analysis. MiRNA modulation (up/down) was used, and data from the first clinical evaluation until recurrence were considered for the non-parametric estimation of disease-free survival using the Kaplan-Meier method. Risk curves were used, and comparisons were performed using the log-rank test and the Cox hazard model, considering p-values ≤0.05. In multiple Cox analysis, each clinical variable, such as molecular subtype, TNM and histology (ductal and lobular), was compared separately with the miRNA expression. The survival R package was used in this analysis .
Tissue microarray analysis
Considering that the three selected miRNAs shared PTEN as target we performed a tissue microarray analysis of these 64 cases of breast cancer, in duplicate, to evaluate the expression of this protein in both group of patients (metastatic and metastatic breast cancer). The TMAs were assembled using a manual tissue microarrayer (Beecher Instruments, Silver Springs, MD - USA). The TMA was done, with control tissues (testis and placenta tissue) and a series of 64 duplicate samples of breast ductal carcinoma and lobular histology were related to clinical stage I, II and III the same samples used in the screening by microarray, totaling 128 points in the receptor block.
Sections of 5 μm were cut from the TMA and further processed to immunostaning with PTEN antibody. For the immunohistochemistry of PTEN a Cell Signaling Monoclonal (reference number 9551P) antibody was used, diluted 1:250. Deparaffinised formalin-fixed, paraffin-embedded tissue sections were heated in the pretreatment module of the autostainer in Tris–HCl pH 8.5 buffer (for 20 minutes at 98°C). To detect the immune reaction we used the sequential system with peroxidase-conjugated secondary antibody (Abcam, USA) or the amplification system SS polymer Polymer Volume Mega-HRP Detection Kit (BioGenex, USA) according to the manufacturer's specifications. For the development of the tags the chromogen diaminobenzidine, DAB (Sigma, USA) and counter-stained with hematoxylin was used. Images of histological sections were digitized using Eclipse 50i microscope (Nikon, Japan) coupled to a digital video camera Sight DS system - Fi1 (Nikon, Japan) with an image analyzer Image-Pro Express version 6.0 (Media Cybernetics, USA).
Target prediction was performed using the mirDIP interface . In the present search, we selected at least 3 of 12 algorithms available for prediction. The targets were analyzed by functional analysis using the Database for Annotation, Visualization and Integrated Discovery (DAVID) version 6.7 [33, 34]. This approach was used to identify significant biological processes and pathways that could be shared between targets of miRNAs of interest. A biological process or pathway was considered significant if it contained a minimum of 3 genes per category, featuring score values less than 0.05 including the Benjamini-Hochberg correction.
The present study included 64 breast cancer patients, including non-metastatic patients (n = 35) and those with metastatic recurrence or metastasis at follow up (n = 29). The groups were stratified according to clinical staging (CS) as CSI, CSII and CSIII with the intention of identifying specific and shared miRNAs associated with metastasis. There were no differences between the groups before pairing with regard to CS; after pairing, the only difference between the groups was related to PR status (p = 0.02; Table 1).
miRNAs differentially expressed between non-metastatic and metastatic patients
The differentially expressed miRNAs (non-metastatic vs. metastatic groups) stratified according to clinical stage (CSI, CSII and CSIII) are shown in Additional file 2: Figure S1. The optimal biomarkers were ranked according to sensitivity and specificity values ≥ 80%. The clusters resulting from this analysis are represented in Figure 1, revealing miRNAs highly specific for the metastatic process. The miRNA hsa-miR-183 (CSII) was also selected for further analysis.Seven microRNAs were identified in the metastatic group irrespective of clinical staging (hsa-let-7a, hsa-let-7b, hsa-let-7c, hsa-miR-1308, hsa-miR-21, hsa-miR-494 and hsa-miR-923_v12.0, Figure 2). Of these, hsa-miR-494 and hsa-miR-21 were selected for further analysis using real-time PCR.
Relative risk analysis
The miRNA expression levels of the 64 patients evaluated in the microarray were assessed using real-time PCR for relative risk analysis. The primers used are shown in Additional file 3: Table S1.
The 3 microRNAs selected for further real-time PCR analysis shared important biological features, especially with regard to phosphoprotein- and kinase protein-associated functions, according to functional analysis.
The ROC curve, considering the non-metastatic group as the reference, was used to select the threshold value for miRNA modulation of quantitative values obtained by real-time PCR. The values for the area under the curve (AUC), sensitivity and specificity are presented in Table 2.
The relative risk curves using recurrence data as events are shown in Figure 3. Because of the small number of CSI patients, it was not possible to perform this analysis for this group. For CSII patients, hsa-miR-183 was the only miRNA that was significant when analyzed independently using the log-rank test (p = 0.03). In contrast, hsa-miR-21 and hsa-miR-494 showed no significance (p = 0.88 and 0.86, respectively). After combining the up-regulated levels of these 3 miRNAs our analysis, showed that were associated with metastatic events (p = 0.002). The combination of hsa-miR-21 with hsa-miR-494 was not significant after the log-rank test (p = 0.123), whereas other combinations such as hsa-miR-21 with hsa-miR-183 (p = 0.004) and hsa-miR-183 with hsa-miR-494 (p = 0.001) were significant. Among CSIII cases, we did not find significant results (data not shown). The miRNA expression and clinical variable data were included in the multiple Cox models, and the best results are shown in Table 3.
The main characteristics that provided some degree of risk included lobular histology, T2, N1 and the Her2 molecular subtype. However, none of these factors were significant in CSII patients and miRNA expression was a more effective independent prognostic factor in all cases.Finally, our results showed that PTEN protein could not be detected in almost all cases of breast cancer samples when the three miRNAs were induced by TMA methodology (Figure 4).
Despite recent findings regarding the role of microRNAs in metastasis, the molecular mechanisms of breast cancer progression remain incompletely understood. For example, the prognostic significance of tumor grade in this type of cancer remains unknown, as well as the molecular mechanisms for why small tumors from CSI (as well as those from CSII and CSIII) can lead to metastasis, which presented 10-years of follow-up. The CSI metastatic samples are extremely rare to obtain, and, in the present study it was possible to obtain four cases in hundreds of patients analyzed. Although there were difficulties associated with obtaining sufficient numbers of CSI metastatic samples, the similarities between tumors of different clinical stages were considered. We performed a retrospective study and identified new biomarkers of breast cancer metastasis that are both shared and specific to clinical stages I-III using FFPE samples. Despite low quality of FFPE samples, several studies have showed the stability of miRNAs [7, 8], because may be less affected than mRNA by formalin fixation and paraffin embedding perhaps due to their smaller size and lack of poly A tails . In the present study, an additional quality control was performed, showing high correlation of frozen samples and FFPE samples used (Additional file 1: Figure S2).
The first application of microarray-based gene expression profiling analysis to the study of breast cancer consisted of disease assessment at the molecular level. In addition, class-prediction studies aim to identify miRNA predictors that could be applicable to all patients with breast cancer, with the goal of separating patients according to prognosis and selecting candidate genes for metastasis during follow up .
The Mamaprint, Oncotype DX and Breast Cancer Index selected patients who were ER positive, and the Veridez 76-gene study evaluated patients without lymph node metastases . The present case cohort study evaluated patients with invasive breast cancer based on metastasis development and selected microRNAs related to metastasis development independent of lymph node or hormonal status.
The Oncotype DX study initially evaluated a retrospective cohort, and the genes of interest were selected using a univariate Cox analysis with a median of 15.1 years of follow up . The Mamaprint study evaluated the odds of developing distant metastases after a 5-year follow-up period . We performed a case cohort study with a control group that had a median follow-up duration of 10 years, which decreased the bias related to patient class migration.
The main prognostic factors related to breast cancer are summarized according to the TNM classification; therefore, at diagnosis, the tumor size, lymph node status and distant metastases represent the main prognostic factors. Another independent prognostic factor is the gene signature; however, this is not easy to evaluate in clinical practice. Although a semi-quantitative assessment of ER, PR, HER2 and Ki67 status-using immunohistochemistry is frequently used in clinical practice, this method does not provide the true gene signature .
In our study, we standardized the cases and controls according to TNM classification that is considered the main breast cancer prognostic factor and we did not observe any differences between the groups (Table 1). Other variables included in the analysis were related to histology and immunohistochemistry markers. Although the PR status was different between the groups, this did not serve as a prognostic factor, which may be due to the limited number of patients evaluated or a possible bias associated with the group selection.
Considering recent findings regarding breast cancer metastasis, Valastyan  reviewed the role of microRNAs according to the steps of the metastatic cascade. In our study, miRNAs shared by clinical stages I-III were found in concordance with processes such as cell motility (let-7 family) and distant metastases (hsa-miR-21). Of these miRNAs, we chose to confirm hsa-miR-21 using real-time PCR because it is the most frequently reported miRNA in several types of solid cancers, representing a potential oncomir [40, 41].
These miRNAs have been described as regulating oncogenes or tumor suppressor genes. The let-7 family appears to regulate the expression of RAS and HMGA2 in breast cancer cells and is associated with several mechanisms of carcinogenesis, including EMT [42, 43]. The other microRNA differentially expressed in both clinical staging, hsa-miR-494, targets several molecules relevant to cancer, including PTEN. Transfection studies have shown this miRNA to act as a master cell cycle regulator at the G1/S checkpoint by targeting CDK6 and at G2/M arrest by targeting PLK1, PTTG1, CCNB1, CDC2, CDC20 and TOP2A, and this miRNA also affects cell proliferation in A549 lung cells by regulating IGFBP1 and IGF2 and in gastrointestinal tumor cells by direct targeting KIT. In addition, miR-494 also appears to have a role in TRAIL-induced apoptosis  as well as the immune system via its regulation of key transcription factors, such as interferon ɣ and TGFβ1 . The evidence of its role in metastasis was demonstrated as an exossomal miRNA in pre-metastatic sites targeting cadherin-17 in mice . Together, these data as well as the findings of the present study highlight this miRNA as a new interesting candidate for verification in metastatic breast cancer.
Considering the evidences of the importance of miR-21 and miR-494, both confirmed by RT-qPCR in all CS. However, it was not clear the sensitivity and specificity of these biomarkers. After ROC curve analysis, despite the miR-21 was found as differentially expressed in all CS, only presented AUC ≥ 80% in CSI-CSIM and miR-494 were not significative. For these reasons, we also selected other miRNAs of CSII and CSIII to improve the sensibility and specificity.
Supervised analysis, based on clinical staging stratification, identified hsa-miR-183 as the best microRNA with regard to sensitivity and specificity in CSII (Figure 2). The hsa-miR-140-3p, the miRNA tested for CSIII, did not confirm microarray results (data not shown). Of these stage-specific microRNAs, miR-183, which was specific to CSII, is considered an oncogene because it targets DKK3, SMAD4, EGR1, PTEN and the PI3K pathway, and it is frequently described in breast cancer as well as other reproductive system related-cancers such as prostate, ovarian and urothelial carcinomas [52, 53]. Moreover, miR-183 has been considered a metastatic inhibitor by targeting ezrin  and lymph node metastasis in medullary thyroid carcinoma , and this miRNA was recently described as being involved in breast cancer progression . Despite these findings, our study is the first to report the association between miR-183 and CSII patients.
Despite the fact that microRNAs hsa-miR-494 and hsa-miR-21 share important targets such as PTEN, this combination was not significant after the log-rank test for CSII or all clinical stages pooled together (data not shown). These miRNAs presented an increased risk only when they were analyzed in combination with miR-183. Moreover, in regards to the molecular mechanisms shared by targets of these miRNAs, we identified phosphoproteins, specifically kinases, among both miRNA targets. For example, KIT and BCL6 were among the hsa-miR-494 proto-oncogenic targets, which represented a significant category with 34 genes, whereas tumor suppressors and apoptotic genes were more evident among hsa-miR-183 target genes, including PTEN, PDCD4 and BCL10, as well as tyrosine kinase signaling pathways. Protein expression of PTEN by TMA analysis showed no differences between metastatic and non metastatic breast tumors (Figure 4). The results showed that expression of PTEN protein was repressed in all cases of breast cancer while the three microRNAs are induced. Considering that some spots of TMA was lost, specially to CSI non metastatic tumor samples, probably this could be a limitation to see differences between metastatic group versus non metastatic group about the this protein expression. In our concern further studies using functional assays can be design to better explain the role of these microRNAs in PTEN regulation.
The 3 miRNAs selected (hsa-miR-21, hsa-miR-494 and hsa-miR-183) share PTEN as a target, and this combination demonstrated an increased risk for metastasis (Figure 3), which suggests some potentially shared mechanism of action. Moreover, the Cox regression analysis showed that the risk of breast cancer metastasis was more likely to be related to miRNA expression and appeared to be independent of clinic pathological variables. In this study, a homogeneous population was intentionally selected to evaluate the effect of miRNA deregulation with increased efficacy. However, it will be necessary to perform further studies in larger populations to validate these findings. Together, our findings indicate that miRNAs can be independently associated with patient prognosis in breast cancer and may represent risk biomarkers for the development of breast cancer metastasis.
Further studies are necessary to understand the role of these new candidate risk biomarkers and the effects of the combination of these miRNAs in breast cancer metastasis, especially in CSII.
Taken together, we demonstrate that miR-183, miR-494 and miR-21 were up-regulated in metastatic breast cancer tissues that was associated a poor prognosis. The TMA analysis showed that the expression of PTEN protein was repressed in all cases of breast cancer while the three miRNAs are induced. These data can indicate that these miRNAs represent new risk biomarkers of metastatic breast cancer and may be useful for future targeted studies.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA Cancer J Clin. 2011, 61: 69-90. 10.3322/caac.20107.
Sobin LH, Gospodarowicz M, Wittekind C: TNM Classification of Malignant Tumours. 2010, Oxford: Wiley-Blackwell, 7
Sandhu R, Parker JS, Jones WD, Livasy CA, Coleman WB: Microarray-based gene expression profiling for molecular classification of breast cancer and identification of new targets for therapy. Lab Med. 2010, 41: 364-372. 10.1309/LMLIK0VIE3CJK0WD.
Perou CM, Sørlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lønning PE, Børresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.
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.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Xi Y, Nakajima G, Gavin E, Morris CG, Kudo K, Hayashi K, Ju J: Systematic analysis of microRNA expression of RNA extracted from fresh frozen and formalin-fixed paraffin-embedded samples. RNA. 2007, 13: 1668-1674. 10.1261/rna.642907.
Li J, Smyth P, Flavin R, Cahill S, Denning K, Aherne S, Guenther SM, O’Leary JJ, Sheils O: Comparison of miRNA expression patterns using total RNA extracted from matched samples of formalin-fixed paraffin-embedded (FFPE) cells and snap frozen cells. BMC Biotechnol. 2007, 7: 36-10.1186/1472-6750-7-36.
Carleton M, Cleary MA, Linsley PS: MicroRNAs and cell cycle regulation. Cell Cycle. 2007, 6: 2127-2132. 10.4161/cc.6.17.4641.
Boehm M, Slack FJ: MicroRNA control of lifespan and metabolism. Cell Cycle. 2006, 5: 837-840. 10.4161/cc.5.8.2688.
Baltimore D, Boldin MP, O’Connell RM, Rao DS, Taganov KD: MicroRNAs: new regulators of immune cell development and function. Nat Immunol. 2008, 9: 839-845. 10.1038/ni.f.209.
Garzon R, Calin GA, Croce CM: MicroRNAs in Cancer. Annu Rev Med. 2009, 60: 167-179. 10.1146/annurev.med.59.053006.104707.
Iorio MV, Ferracin M, Liu C-G, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Ménard 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.
Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin S-F, Dunning MJ, Barbosa-Morais NL, Teschendorff AE, Green AR, Ellis IO, Tavaré 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.
Valastyan S: Roles of microRNAs and other non-coding RNAs in breast cancer metastasis. J Mammary Gland Biol Neoplasia. 2012, 17: 23-32. 10.1007/s10911-012-9241-9.
Yang J, Weinberg RA: Epithelial-mesenchymal transition: at the crossroads of development and tumor metastasis. Dev Cell. 2008, 14: 818-829. 10.1016/j.devcel.2008.05.009.
Sreekumar R, Sayan BS, Mirnezami AH, Sayan AE: MicroRNA control of invasion and metastasis pathways. Front Genet. 2011, 2: 58-
Lowery AJ, Miller N, Dwyer RM, Kerin MJ: Dysregulated miR-183 inhibits migration in breast cancer cells. BMC Cancer. 2010, 10: 502-10.1186/1471-2407-10-502.
Luo D, Wilson JM, Harvel N, Liu J, Pei L, Huang S, Hawthorn L, Shi H: A systematic evaluation of miRNA:mRNA interactions involved in the migration and invasion of breast cancer cells. J Transl Med. 2013, 11: 57-10.1186/1479-5876-11-57.
Zhang Q-H, Sun H-M, Zheng R-Z, Li Y-C, Zhang Q, Cheng P, Tang Z-H, Huang F: Meta-analysis of microRNA-183 family expression in human cancer studies comparing cancer tissues with noncancerous tissues. Gene. 2013, 527: 26-32. 10.1016/j.gene.2013.06.006.
Rundle AG, Vineis P, Ahsan H: Design options for molecular epidemiology research within cohort studies. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1899-1907. 10.1158/1055-9965.EPI-04-0860.
Carey LA, Perou CM, Livasy CA, Dressler LG, Cowan D, Conway K, Karaca G, Troester MA, Tse CK, Edmiston S, Deming SL, Geradts J, Cheang MCU, Nielsen TO, Moorman PG, Earp HS, Millikan RC: Race, breast cancer subtypes, and survival in the Carolina Breast Cancer Study. JAMA. 2006, 295: 2492-2502. 10.1001/jama.295.21.2492.
Cheang MCU, Chia SK, Voduc D, Gao D, Leung S, Snider J, Watson M, Davies S, Bernard PS, Parker JS, Perou CM, Ellis MJ, Nielsen TO: Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst. 2009, 101: 736-750. 10.1093/jnci/djp082.
The R Project for Statistical Computing. [http://www.r-project.org]
Bengtsson H, Hössjer O: Methodological study of affine transformations of gene expression data with proposed robust non-parametric multi-dimensional normalization method. BMC Bioinformatics. 2006, 7: 100-10.1186/1471-2105-7-100.
Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics. 2006, 22: 2825-2827. 10.1093/bioinformatics/btl476.
Park SH, Goo JM, Jo C-H: Receiver operating characteristic (ROC) curve: practical review for radiologists. Korean J Radiol. 2004, 5: 11-18. 10.3348/kjr.2004.5.1.11.
Sing T, Sander O, Beerenwinkel N, Lengauer T: ROCR: visualizing classifier performance in R. Bioinformatics. 2005, 21: 3940-3941. 10.1093/bioinformatics/bti623.
gplots: Various R programming tools for plotting data. [http://cran.r-project.org/web/packages/gplots/index.html]
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Therneau T, Lumley T: R: Survival analysis, including penalised likelihood. [http://stat.ethz.ch/R-manual/R-patched/library/survival/html/00Index.html]
mirDIP: microRNA Data Integration Portal. [http://ophid.utoronto.ca/mirDIP/]
DAVID Functional Annotation Bioinformatics Microarray Analysis. [http://david.abcc.ncifcrf.gov]
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008, 4: 44-57. 10.1038/nprot.2008.211.
Hall JS, Taylor J, Valentine HR, Irlam JJ, Eustace A, Hoskin PJ, Miller CJ, West CML: Enhanced stability of microRNA expression facilitates classification of FFPE tumour samples exhibiting near total mRNA degradation. Br J Cancer. 2012, 107: 684-694. 10.1038/bjc.2012.294.
Reis-Filho JS, Pusztai L: Gene expression profiling in breast cancer: classification, prognostication, and prediction. Lancet. 2011, 378: 1812-1823. 10.1016/S0140-6736(11)61539-0.
Cobleigh MA, Tabesh B, Bitterman P, Baker J, Cronin M, Liu M-L, Borchik R, Mosquera J-M, Walker MG, Shak S: Tumor gene expression and prognosis in breast cancer patients with 10 or more positive lymph nodes. Clin Cancer Res. 2005, 11 (24 Pt 1): 8623-8631.
Veer LJ V ’t, Dai H, van de Vijver MJ, He YD, Hart AAM, Mao M, Peterse HL, van der Kooy 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.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thürlimann B, Senn H-J: Panel members: Strategies for subtypes–dealing with the diversity of breast cancer: highlights of the St. Gallen international expert consensus on the primary therapy of early breast cancer 2011. Ann Oncol. 2011, 22: 1736-1747. 10.1093/annonc/mdr304.
Asangani IA, Rasheed SAK, Nikolova DA, Leupold JH, Colburn NH, Post S, Allgayer H: MicroRNA-21 (miR-21) post-transcriptionally downregulates tumor suppressor Pdcd4 and stimulates invasion, intravasation and metastasis in colorectal cancer. Oncogene. 2008, 27: 2128-2136. 10.1038/sj.onc.1210856.
Lu Z, Liu M, Stribinskis V, Klinge CM, Ramos KS, Colburn NH, Li Y: MicroRNA-21 promotes cell transformation by targeting the programmed cell death 4 gene. Oncogene. 2008, 27: 4373-4379. 10.1038/onc.2008.72.
Watanabe S, Ueda Y, Akaboshi S, Hino Y, Sekita Y, Nakao M: HMGA2 maintains oncogenic RAS-induced epithelial-mesenchymal transition in human pancreatic cancer cells. Am J Pathol. 2009, 174: 854-868. 10.2353/ajpath.2009.080523.
Mayr C, Hemann MT, Bartel DP: Disrupting the pairing between let-7 and Hmga2 enhances oncogenic transformation. Science. 2007, 315: 1576-1579. 10.1126/science.1137999.
Liu Y, Lai L, Chen Q, Song Y, Xu S, Ma F, Wang X, Wang J, Yu H, Cao X, Wang Q: MicroRNA-494 Is required for the accumulation and functions of tumor-expanded myeloid-derived suppressor cells via targeting of PTEN. J Immunol. 2012, 188 (11): 5500-5510. 10.4049/jimmunol.1103505.
Olaru AV, Ghiaur G, Yamanaka S, Luvsanjav D, An F, Popescu I, Alexandrescu S, Allen S, Pawlik TM, Torbenson M, Georgiades C, Roberts LR, Gores GJ, Ferguson-Smith A, Almeida MI, Calin GA, Mezey E, Selaru FM: MicroRNA down-regulated in human cholangiocarcinoma control cell cycle through multiple targets involved in the G1/S checkpoint. Hepatology. 2011, 54: 2089-2098. 10.1002/hep.24591.
Yamanaka S, Campbell NR, An F, Kuo SC, Potter JJ, Mezey E, Maitra A, Selaru FM: Coordinated effects of microRNA-494 induce G2/M arrest in human cholangiocarcinoma. Cell Cycle. 2012, 11: 2729-2738. 10.4161/cc.21105.
Ohdaira H, Sekiguchi M, Miyata K, Yoshida K: MicroRNA-494 suppresses cell proliferation and induces senescence in A549 lung cancer cells. Cell Prolif. 2012, 45: 32-38. 10.1111/j.1365-2184.2011.00798.x.
Kim WK, Park M, Kim Y-K, Tae YK, Yang H-K, Lee JM, Kim H: MicroRNA-494 Downregulates KIT and Inhibits Gastrointestinal Stromal Tumor Cell Proliferation. Clin Cancer Res. 2011, 17: 7584-7594. 10.1158/1078-0432.CCR-11-0166.
Romano G, Acunzo M, Garofalo M, Di Leva G, Cascione L, Zanca C, Bolon B, Condorelli G, Croce CM: MiR-494 is regulated by ERK1/2 and modulates TRAIL-induced apoptosis in non-small-cell lung cancer through BIM down-regulation. Proc Natl Acad Sci U S A. 2012, 109: 16570-16575. 10.1073/pnas.1207917109.
Ohdaira H, Sasaki T, Yoshida K: A subset of microRNAs potentially acts as a convergent hub for upstream transcription factors in cancer cells. Oncol Rep. 2010, 24: 1371-1381.
Rana S, Malinowska K, Zoller M: Exosomal tumor microRNA modulates premetastatic organ cells. Neoplasia. 2013, 15: 281-295.
Yamada Y, Enokida H, Kojima S, Kawakami K, Chiyomaru T, Tatarano S, Yoshino H, Kawahara K, Nishiyama K, Seki N, Nakagawa M: MiR-96 and miR-183 detection in urine serve as potential tumor markers of urothelial carcinoma: correlation with stage and grade, and comparison with urinary cytology. Cancer Sci. 2011, 102: 522-529. 10.1111/j.1349-7006.2010.01816.x.
Wang G, Mao W, Zheng S: MicroRNA-183 regulates Ezrin expression in lung cancer cells. FEBS Lett. 2008, 582: 3663-3668. 10.1016/j.febslet.2008.09.051.
Abraham D, Jackson N, Gundara JS, Zhao J, Gill AJ, Delbridge L, Robinson BG, Sidhu SB: MicroRNA profiling of sporadic and hereditary medullary thyroid cancer identifies predictors of nodal metastasis, prognosis, and potential therapeutic targets. Clin Cancer Res. 2011, 17: 4772-4781. 10.1158/1078-0432.CCR-11-0242.
Giricz O, Reynolds PA, Ramnauth A, Liu C, Wang T, Stead L, Childs G, Rohan T, Shapiro N, Fineberg S, Kenny PA, Loudig O: Hsa-miR-375 is differentially expressed during breast lobular neoplasia and promotes loss of mammary acinar polarity. J Pathol. 2012, 226: 108-119. 10.1002/path.2978.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/739/prepub
We thank the Researcher Support Center of Barretos Cancer Hospital, especially the statistician Zanardo C. for assisting in the statistical analysis. This study received financial support from Fundação de Amparo à Pesquisa do Estado de São Paulo (Fapesp, Proc: 10/ 16796–0, São Paulo, Brazil).
The authors declare that they have no competing interests.
ALFM conceived, developed and led the overall study. AFE participated in the data collection, analysis and helped in the manuscript preparation. TM helped with Real Time PCR experiments. RACV and LML provided clinical and pathology support in this study. LFAM helped in TMA analysis. ALFM and HCSS critically reviewed the manuscript. MMCM conceived, provided advice during the study development and prepared the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Figure S1: Heatmaps of all differentially expressed miRNAs stratified according to clinical stage. Figure S1A shows non-metastatic vs. metastatic patients in CSI; Figure S1B shows patients in CSII, and Figure S1C shows patients in CSIII. (TIFF 1 MB)
Additional file 1: Figure S2: Scatterplot comparing microRNA expression profile of FFPE sections (y-axis) against frozen tissue (x-axis). The R square is 0.781. (PDF 135 KB)
About this article
Cite this article
Marino, A.L., Evangelista, A.F., Vieira, R.A. et al. MicroRNA expression as risk biomarker of breast cancer metastasis: a pilot retrospective case-cohort study. BMC Cancer 14, 739 (2014). https://doi.org/10.1186/1471-2407-14-739
- Breast cancer
- Molecular profile
- Retrospective study