- Research
- Open access
- Published:
Hsa-miR-375 is a predictor of local control in early stage breast cancer
Clinical Epigenetics volume 8, Article number: 28 (2016)
Abstract
Background
A long-term analysis by the Early Breast Cancer Trialist Group (EBCTG) revealed a strong correlation between local control and cancer-specific mortality. MicroRNAs (miRs), short (20–25 nucleotides) non-coding RNAs, have been described as prognosticators and predictors for breast cancer in recent years. The aim of the current study was to identify miRs that can predict local control after breast conserving therapy (BCT) in early stage breast cancer.
Results
Clinical data of 46 early stage breast cancer patients with local relapse after BCT were selected from the institutional database. These patients were matched to 101 control patients showing identical clinical features but without local relapse. The study was conducted in two steps. (1) In the pilot study, 32 patients (16 relapses versus 16 controls) were screened for the most de-regulated microRNAs (= candidate microRNAs) in a panel of 1250 miRs by microarray technology. Eight miRs were found to be significantly de-regulated. (2) In the validation study, the candidate microRNAs were analyzed in an independent cohort of 115 patients (30 relapses versus 85 controls) with reverse transcription quantitative polymerase chain reaction (RT-qPCR). From these eight candidates, hsa-miR-375 could be validated. Its median fold change was 2.28 (Mann-Whitney U test, corrected p value = 0.008). In the log-rank analysis, high expression levels of hsa-miR-375 correlated with a significantly higher risk of local relapse (p = 0.003). In a multivariate analysis (forward stepwise regression) including established predictors and prognosticators, hsa-miR-375 was the only variable that was able to distinguish the statistical significance between relapse and control groups (raw p value = 0.000195 HR = 0.76, 95 % CI 0.66–0.88; corrected p value = 0.005).
Conclusions
Hsa-miR-375 predicts local control in patient with early stage breast cancer, especially in estrogen receptor α (ER-α)-positive patients. It can therefore serve as an additional molecular marker for treatment choice independently from known predictors and prognosticators. Validation in larger prospective studies is warranted.
Background
In 2014, the estimated number of deaths from breast cancer in Europe is 89,300 [1] and 40,000 in the US [2]. Every uncontrolled primary tumor comprises the risk of subsequent metastatic disease, and the younger the patient, the more likely she will die from it. A long-term analysis by the Early Breast Cancer Trialist Group (EBCTG) showed a 4:1 relation: if four local relapses can be avoided at year 10, one breast cancer-specific death can be avoided at year 15 [3, 4]. Due to advances in all disciplines related to breast conserving therapy (BCT), the local relapse rate is generally low, with age still appearing as the most important prognosticator and predictor for local relapse. Nowadays, the estimates for age adjusted annual in-breast recurrence rates are 0.4–0.7 % for patients >50 years, 0.72–1.2 % in the age group 41–50 years, and 0.72–2 % for patients younger than 40 years, respectively [5–7].
Additionally, many other patient- and therapy-related factors as well as molecular markers have been described to influence clinical outcome [8–10]. Seventy percent of invasive breast cancers show increased expression of the estrogen receptor α (ER-α). Up-regulation of ER-α is an early event in tumorigenesis and cancer progression [11, 12]. In recent years, microRNAs (miRs) were described as potential prognosticators and predictors for breast cancer. MiRs are short non-coding RNAs (21–25 nucleotides), which play a role in human cancer development both as oncogenes and tumor suppressors. More than 50 % of the known miR genes are located in genetically altered regions of the genome [13]. In 2005, Iorio et al. suggested a miR signature for breast cancer. This panel of de-regulated miRs correlated with known clinical and biological features such as hormonal receptor status, tumor size, lymph node status, vascular invasion, proliferation index, and p53 [14]. Recently, scientific interest has partially focused on the connection between miRs and hormone receptors [15–20]. About 15 miRs have been identified as regulators of ER-α [21].
Thus far, only one clinical study investigated the role of miRs with respect to local control (LC) in early stage breast cancer [22]. Zhou et al. showed that hsa-miR-9 could discriminate between patients with and without local relapse. Interestingly, the effect was most pronounced in ER-α positive patients. In order to further elucidate the importance of miRs for LC in early stage breast cancer, the aim of the current study was to identify miRs that can differentiate between patients with and without local relapse after BCT. This is especially important in the context of individualized local treatment selection.
Results
Patients
Potentially prognostic characteristics were distributed without significant differences between relapse and control groups (Table 1). The median follow-up in the pilot cohort was 126 months (range 26–200 months), and in the validation cohort, it was 136 months (range 17–207 months). The median time to local relapse was 37 months (range 15–123 months) in the pilot cohort and 74 months (range 35–185 months) in the validation cohort (Table 2).
Pilot phase
In the pilot study, five miRs showed a significant fold change between relapse patients and their controls (linear models for microarray analysis (LIMMA), raw p value <0.05): hsa-miR-362-3p, hsa-miR-660, hsa-miR-375, hsa-miR-223, and hsa-miR-125a-3p. Three miRs had a fold change (relapse versus control) that was close to significance (LIMMA, raw p value = 0.05): hsa-miR-532-3p, hsa-miR-487b, and hsa-miR-210. Hierarchical clustering of these eight miRs revealed two main groups with different miR-expression patterns (Fig. 1). One group only consisted of controls (pattern B), while relapse patients formed the majority of the other group (pattern A). The log-rank comparison revealed a significant p value of 0.002 (Fig. 2). The levels of hsa-miR-375 were significantly higher in the relapse group, with a raw p value of 0.009 (LIMMA) and a fold change of 2.97 between relapse and controls, which is illustrated in (Additional file 1: Figure S2). The microarray data were deposited in NCBI Gene Expression Omnibus under the accession number GSE69951; they can be accessed via the following link: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE69951.
Validation phase
For further validation, we selected the eight miRs that were able to differentiate between relapse and controls in the hierarchical cluster analysis. Of these eight candidate miRs, hsa-miR-375 could be validated by reverse transcription quantitative polymerase chain reaction (RT-qPCR), in the sense that it was able to differentiate between relapse and control group in an independent patient cohort (Table 3). The levels of hsa-miR-375 were—similar to the pilot phase—significantly higher in the relapse group compared to controls (Mann-Whitney U test, raw p value = 0.001; p value corrected according to Bonferroni = 0.008) with a median fold change of 2.28 (Table 3). The box plot in Fig. 3 shows the elevated levels of hsa-miR-375 in relapse patients compared to controls (Fig. 3). In the log-rank analysis (event = local relapse), high expression levels of hsa-miR-375 correlated to a significantly increased risk of local relapse (p value = 0.003; Fig. 4a). When tested in the multivariate analysis (forward stepwise regression) together with the patient- and treatment-related parameters listed in Table 1, hsa-miR-375 was the only variable that was able to distinguish the statistical significance between the relapse and control groups (raw p value = 0.000195 HR = 0.76, 95 % CI 0.66–0.88; p value corrected according to Bonferroni = 0.005).
In the subgroup of ER-α-positive patients, the difference in hsa-miR-375 levels between the relapse and control patients was also significant (p value = 0.001, Mann-Whitney U test), whereas in ER-α-negative patients, it was not (p = 0.495, Mann-Whitney U test). The log-rank analysis (event = local relapse) showed that high expression levels of hsa-miR-375 correlated to a higher probability of local relapse, which was more significant than in the whole patient cohort (p = 0.0005, see Fig. 4b).
The comparison of the hsa-miR-375 levels in ER-α-positive to triple negative patients showed a fold change of 4.17 (p value = 0.24, Mann-Whitney U test). Further subgroup comparisons showed that the levels of hsa-miR-375 were higher in the ER-α-positive patients compared to those in the ER-α-negative patients but without significance (p value 0.53). The same applies to her2-negative patients compared with her2-positive patients 1.52 (p value 0.66). For further details, please see Additional file 2: Table S2.
External validation
Hsa-miR-375 was cross-validated in randomly selected samples of five relapses and their 14 corresponding controls from the validation cohort. The result of the first validation step concerning hsa-miR-375 could be verified with a median fold change of 2.84 and a p value of 0.019 (Mann-Whitney U test).
Discussion
Although the clinical outcome in early stage breast cancer patients is generally good, some patients experience local relapse after BCT. Apart from clinical features, biomarkers can add to the optimal treatment choice. In the current study, we were able to identify and validate hsa-miR-375 as a potential predictor for LC after BCT. Our results allow us to assume that patients with high levels of hsa-miR-375 in the primary tumor have an increased risk of local relapse. Therefore, hsa-miR-375 is a new potential predictor for LC in early stage breast cancer independent of the known clinical and molecular markers. It could also become an additional selection tool for the choice of primary operative treatment. Patients with a high level of hsa-miR-375 show a higher risk of local relapse; however, threshold values have yet to be defined before clinical translation. These patients might be possible candidates for more aggressive local treatment, i.e., mastectomy or dose-intensified adjuvant radiotherapy following breast conserving surgery.
The first study to show that miRs in general play a role in breast cancer was published by Iorio and colleagues. The authors compared miR-expression levels by microarray technology in 76 primary tumors versus 34 normal tissues and identified a set of 15 miRs that could differentiate between tumor and normal tissue with 100 % accuracy [14].
Recent studies have focused on the relationship between ER-α and miR-expression levels. This is of importance for the current study since 80 % of the 147 patients we analyzed were ER-α positive. Profiling studies identified numerous miRs that are associated with either ER-α-positive or ER-α-negative tumors [14, 23–26]. The web of miRs involved in the regulation of the estrogen receptor activity is complex. One of these is hsa-miR-375, which plays a role in various types of cancer such as prostate [27], stomach [28], lung [29], and pancreas [30]. A connection between hsa-miR-375 and ER-α positivity was described in vitro [11, 31], in an interventional study in rodents [12] as well as in the clinical context of Tamoxifen resistance [31–34] and progression from normal tissue through lobular carcinoma in situ (LCIS) to invasive lobular carcinoma (ILC) [35].
Souza et al. compared the profiles of miRs and mRNA in cell culture. Inhibition of hsa-miR-375 reduced ER-α activation and cell proliferation. Dexamethasone-induced RAS-related protein 1 (RASD1) was found to be a direct target of hsa-miR-375, which negatively regulates ER-α expression. Hsa-miR-375 was up-regulated in ER-α-positive cells compared to ER-α-negative cells and normal cells. This led to the assumption that up-regulation of hsa-miR-375 is both a key driver of proliferation and an early event in tumorigenesis in ER-α positive tissues [34].
In an interventional study in rodents, Munagala et al. could verify the correlation between hsa-miR-375 and ER-α. ACI (August-Copenhagen-Irish) rats were treated with estradiol alone, its inhibitor ellagic acid alone, and estradiol + ellagic acid. ER-α was found to be regulated—among other miRs—by hsa-miR-375. In estradiol-treated cancers, up-regulation of hsa-miR-375 is common and can be reversed by ellagic acid [12, 36].
Ward et al. combined in vitro data with retrospective clinical results. Hsa-miR-375 was found to be down-regulated in Tamoxifen resistant cells. Its re-expression reverses both Tamoxifen resistance and accompanying epithelial-mesenchymal transition (EMT)-like properties in breast cancer. The levels of metadherin (MTDH), a direct target of hsa-miR-375, were measured in patients and correlated inversely with disease-free survival [31]. Hoppe et al. described the connection between miR profiles, ER-α status, and Tamoxifen resistance in postmenopausal patients. Low levels of hsa-miR-375 predicted worse relapse free time, but as opposed to hsa-miR-10a and hsa-miR-126, this effect was not statistically significant [32]. Lyng et al. conducted a similar study in 152 patients with ER-α-positive tumors. Distant recurrence was the statistical endpoint. Of note, patients in the pilot set were matched according to their clinical features, whereas patients in the two validation sets were randomly selected. None of the miRs in the pilot set overlapped with the miRs in the validation sets, and the only statistically highly significant difference between the group of recurrent and non-recurrent patients was the number of tumor-infiltrated lymph nodes [33]. Jonsdottir correlated miR profiles with known clinicopathological factors in 204 lymph node-negative patients with distant metastasis-free survival. Hsa-miR-375 levels were high in ER-α-positive patients. As opposed to hsa-miR-106, hsa-miR-375 was unable to differentiate between patients with or without distant metastases [34]. Giricz et al. investigated hsa-miR-375 in normal tissue, lobular carcinoma in situ and invasive lobular carcinoma (ILC) in formalin-fixed paraffin-embedded (FFPE) specimens of 31 patients and in vitro. The highest levels of hsa-miR-375 were detected in ILC and LCIS adjacent to ILC. Their gene transfer experiments revealed that induced elevation of this microRNA to levels found in human lobular breast cancer specimens changes the ability to maintain appropriate tissue organization. The authors conclude that high levels of hsa-miR-375 contribute to progression from normal tissue to ILC [35]. Since over 90 % [35] of ILC are ER-α-positive, there seems to be a functional connection between ER-α overexpression and elevated levels of hsa-miR-375 as described by Souza et al. [11].
Thus far, the only study with LC as statistical endpoint was performed by Zhou in a cohort of 68 patients. Only FFPE samples with a tumor burden of at least 70 % were included in the analysis. Hsa-miR-9 was able to differentiate between patients with and without local relapse. Hsa-miR-375 was one of the candidate molecules selected in the pilot phase but was finally unable to predict LC in the patients of the validation set. Nevertheless, the authors found a strong correlation between hsa-miR-375 and ER-α status [22].
Taken together, these studies show that to date, it remains a matter of debate whether hsa-miR-375 is of predictive and/or prognostic relevance for early stage breast cancer. In accordance with Jonsdottir et al., our results also show that high levels of hsa-miR-375 lead to a worse outcome [34]. The effect is particularly pronounced in ER-α-positive patients, which coincides with the mechanistic evidence provided by Souza et al. [11]. Discrepancies between our analysis and literature may be explained by our focus on LC in early stage breast cancer. Different clinical settings result in different statistical endpoints, such as disease-free survival and may therefore lead to diverging results. Additionally, sample-related differences may be important. (1) The current study compares tumor tissues with each other, which means that the difference in the miR levels of relapse versus controls must be large to be detected. Iorio et al., e.g., compare normal and tumor tissues [14]. (2) Some studies were performed after microdissection of the specimen [22, 34], whereas in the present analysis, tumor tissue contains stroma. This strategy is based on daily clinical procedures and therefore feasible in diagnostic routine without additional manipulation of the specimen, which adds to the robustness of our results. (3) Whole tissue sections were used in the present analysis since we believe that a given miR—if clinically relevant as a predictor and/or prognosticator—should be able to differentiate between patients even if the tumor cell content is small.
However, a limitation of our study concept is that the cellular signal is either from stromal or tumor tissue. Moreover, there is no validated endogenous reference gene available for breast cancer [37–39]. According to literature, the standard approach is to use small nuclear or small nucleolar RNAs as controls. The reference gene we used was SNORD-61, a small nucleolar RNA, which was validated by Spervelage for neuro-endocrine tumors of the ileum [38].
Beyond the scope of the current analysis, it is important to further elucidate the oncogenic interference of hsa-miR-375 in cancer pathways. This entails in particular the necessity to define mRNA molecules, which contain potential binding sites for hsa-miR-375. In theory, the result would be a down-regulation of a tumor suppressive protein. Our search in five different databases (TargetScan, miRDB, PITA, DIANA, DIANA Cancer) performed in May 2015 resulted in a large variety of potential target molecules. We chose the results given by TargetScan and PITA for two reasons. First, it is generally recommended to select overlapping results from databases that implement different algorithms [40]. Second, these two databases provide a reasonable number of potential targets. TargetScan (Release 6.2) identified 229 putative mRNA targets based on their conserved complementarities to the seed region of the miR [41]. PITA additionally includes secondary structures and binding energies for prediction [42], which consequently resulted in a smaller number of 116 potential targets (Additional file 3). The comparison of both databases revealed 15 overlaps. Among these, we identified one gene that has been validated experimentally by reporter gene assay in breast cancer cells: RASD1 [11]. It belongs to the Ras superfamily of small G proteins, but it does not share key residues in reactive domains with other Ras family members, such as V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS), transforming protein p21 (HRAS), and neuroblastoma ras viral oncogene homolog (NRAS). As opposed to these structurally related molecules, RASD1 inhibits clonal cell expansion suggesting a tumor suppressive role [43]. Additionally, we used the GeneCoDis software [44] to correlate this predicted target with cancer pathways (Panther classification system). According to this analysis, RASD1 is related to the PI3K-pathway (Additional file 4). The involvement in known cancer pathways and its experimental validation as a direct target of hsa-miR-375 makes RASD1 a promising candidate for further evaluation and underlines the potential biological significance of hsa-miR-375 in early stage breast cancer.
Conclusions
In summary, we could show that high levels of hsa-miR-375 are associated with a significantly higher probability of local relapse, especially in ER-α-positive patients. Hence, hsa-miR-375 is an independent predictor for local treatment outcome after breast conserving therapy. Although our results require validation in larger prospective studies, we suggest hsa-miR-375 as an additional molecular marker for treatment choice in early stage breast cancer.
Methods
Patients
Between October 1998 and October 2012, 5093 patients were treated for breast cancer at the Department of Radiation Oncology at the Paracelsus Medical University, Salzburg. For the current study, 147 breast cancer patients were selected. According to the 7th edition of the TNM/AJCC classification system, early stage breast cancer includes I, IIA, and IIB (T2 N1). 145/147 (99 %) patients in the current study are in one of this disease stage. 2/147 (1 %) patients (one relapse patient and her matched control) in the pilot cohort had N2 disease. All patients gave their informed consent for surgery, radiotherapy, and systemic treatment. Patients with local relapse were matched to controls without local relapse. Local relapse was defined as re-appearance of the invasive cancer in the same breast, regardless of whether the invasive tumor re-appeared within the former index quadrant, defined by the tumor bed with a 3-cm margin, or not [45–47]. Patients who showed no radiological and histological evidence of local relapse were defined as locally controlled. The study cohort comprised 46 in-breast relapses (in and out quadrant) and 101 matched controls. The patients in our study were selected based on the following criteria: (1) invasive carcinoma; (2) minimum follow-up of five years; (3) the tumor specimen had to be processed by the Department of Pathology at the Paracelsus Medical University; and (4) sufficient FFPE material. Patients were matched according to the year of diagnosis, the type of surgery (mastectomy or lumpectomy), the type of radiotherapy (whole-breast irradiation with percutaneous or intraoperative boost), age, tumor size, lymph node involvement, grading, histology, hormonal receptor status, her2 status, menopausal status, and Ki67 proliferation index. The study was approved by the local ethics committee in Salzburg (Ethikkommission für das Bundesland Salzburg 415-EP/73/85-2012).
Experimental design
The study was conducted in a pilot and a validation phase. In the pilot phase, primary tumors of 32 patients (16 patients who experienced a local relapse versus 16 matched controls who had no relapse) were screened for the most de-regulated miRs by means of a commercially available microarray (Agilent™ Sure Print). This part of the study was carried out by the Comprehensive Biomarker Center™ in Heidelberg. Candidate miRs were selected according to their fold changes (i.e., the ratio of a miR levels in patients with local relapse to those without), and their statistical significance was estimated with the moderated t test (LIMMA).
The validation study in an independent cohort was performed in primary tumor tissue of 115 patients and consisted of two steps. In the first step, the candidate miRs were analyzed in 30 patients with local relapse compared to 85 matched controls without relapse. This step was performed by the Laboratory of Molecular Pathology at the Department of Pathology, Paracelsus Medical University, Salzburg. The expression levels of these candidate miRs were assessed with RT-qPCR. In the second step, these results were cross-validated in 19 patients (5 relapses versus 14 matched controls) randomly selected from the validation cohort of 115 patients. This part of the analysis was done at the Department of Pathology of the Ludwig Maximilian University in Munich on a technically different RT-qPCR platform. The experimental design is summarized in the study flow-chart (Additional file 5).
Samples
Surgical samples were FFPE and archived in the tissue bank of the Department of Pathology at the Paracelsus Medical University, Salzburg. Apart from routine procedures applied to clinical samples, no additional processing was performed. For the current analysis, the whole tissue sections with a maximum proportion of tumor tissue were selected. The tumor burden in the biopsies ranged from 10 to 90 % (see Table 1).
Pilot phase: microarray
By means of microarray technology, a panel of 1250 miRs was screened. For this purpose, the FFPE samples of the primary tumors were retrieved from the tissue bank, and seven consecutive sections per patient were prepared (depending on tissue size, the thickness of the sections was 2 to 4 μm). Isolation of total miR and chip-based microarrays (Agilent’s Sure PrintG3 Human miRNA microarrays™) were performed according to standard procedures by the Comprehensive Biomarker CenterTM, Heidelberg.
Validation phase: RT-qPCR
RNA extraction and reverse transcription
Routinely, we used seven sections with a 10-μm thickness for RNA extraction with the miRNeasy FFPE kit (QiagenTM) according to the manufacturer’s instructions. Each sample was treated with DNase I (included in the kit) to avoid DNA contamination. The RNA was eluted in 20 μl water and stored at −80 °C. For RNA quantification, the Qubit™ RNA Assay Kit (InvitrogenTM) was used. One hundred thirty nanogram of RNA were reversely transcribed using miScript II RT Kit (QiagenTM) in a total volume of 20 μl. Each reaction was incubated at 37 °C for 60 min and heat-inactivated at 95 °C for 5 min. The complementary DNA (cDNA) was stored at −20 °C until further use.
Real-time PCR
Expression levels of candidate miRs were measured by miScript SYBR Green PCR Kit (QiagenTM). Specific PCR target information is provided in Additional file 6: Table S1. The primer assays for each target miR and SNORD61 (catalog number MS00033705) were purchased from QiagenTM. The reaction set up was 0.1 μl cDNA, 1 μl 10× miScript Universal Primer, 5 μl QuantiTect SYBR Green PCR master mix, and 1 μl Primer miScript Primer assay in a total reaction volume of 10 μl. Following the initial activation step at 95 °C for 15 min, the cycling conditions were 40 times 15 s 94 °C, 30 s at 55 °C, and 30 s at 70 °C. All PCR reactions were performed in technical duplicates.
qPCR validation
The PCR data were analyzed with the Rotor-Gene 2.0.2.4 software (QiagenTM) according to the suggestions of MIQE guidelines (minimum information for publication of quantitative real-time PCR experiments) [48]. To test the PCR specificity, every reaction was subjected to melting curve analyses. The amplification efficiency of each primer pair was determined by a qPCR assay using a series of 1:5 dilutions of cDNA from the control samples. Each reaction was performed in triplicates. The slope of the resulting standard curve was used to calculate the PCR efficiency according to the following formula: efficiency = 10(−1/slope).
Data analysis
The Ct value of the reference gene was subtracted from the Ct value of a given miR to receive the ∆Ct value. Subsequently, the median of all ∆Ct values was calculated. The median of the control group was subtracted from the median of the relapse group to receive a ∆∆Ct value. The fold change was calculated with the ∆∆Ct method [49]. The non-parametric Mann-Whitney U test was used to compare the ∆Ct values of relapse and control groups. Raw p values were corrected for multiple testing according to Bonferroni.
External validation
This second validation step on a technically different platform was performed with the sole purpose of verifying the results of the first validation step; therefore, only hsa-miR-375 was analyzed in this cohort of 19 randomly selected patients. SNORD61 was used as the reference gene.
SNORD61 expression analysis
Reverse transcription of 100 ng total RNA was carried out using the RevertAid First Strand cDNA Synthesis Kit (Thermo ScientificTM), according to the instructions of the manufacturer. RT-qPCR analyses were performed on a LightCycler 480 instrument using the TaqMan Gene Expression Assay (Life TechnologiesTM) following the instructions of the manual. Briefly, the final reaction mix contained 2 ng cDNA, 1× TaqMan Gene Expression Assay, and 1× TaqMan Gene Expression Mastermix II, no UNG. Single qPCRs analyses were run in duplicates.
Hsa-miR-375 expression analysis
Reverse transcription reactions were performed using the TaqMan™ Reverse Transcription Kit (Life TechnologiesTM) according to the manufacturer’s instructions. For single RT reactions, miR-specific RT primers were used and miRs were converted into cDNA in one reaction using 10 ng total RNA. cDNA templates were analyzed applying the TaqMan Small RNA Assay for hsa-miR-375, following the instructions of the manufacturer. Reaction Mixes contained 1× TaqMan Gene Expression Mastermix II, no UNG, 1× TaqMan miRNA assay, and a 1:15 dilution of the RT product. The temperature profile was as follows: 10 min at 95 °C, followed by 40 cycles of 95 °C for 15 s, and 60 °C for 1 min. All samples were analyzed in duplicate reactions. The relative expression of the miRs was calculated by the ΔΔCt comparative threshold method [49].
Statistics
In the pilot phase, differentially expressed miRs were selected as candidates for further validation according to their fold change and its statistical significance estimated with the moderated t test (LIMMA) [50]. Raw p values were adjusted for multiple testing according to Benjamini-Hochberg [51]. In order to avoid losing any biomarker with possibly discriminative potential, this rather weak correction method was applied on purpose in the pilot phase. We tried to detect possible clusters in rows (miRs) and columns (samples) of the normalized expression matrix by hierarchical clustering (bottom-up complete linkage clustering using the Euclidean distance as a measure). The dendrogram on top of the expression matrix demonstrates samples which cluster together. Thereby, a miR pattern could be defined that differentiates patients with and without local relapse.
In the validation phase, the Mann-Whitney U test was used to compare ∆Ct values, patient-, tumor-, and treatment-related characteristics between groups for statistically significant differences. LC was estimated with the Kaplan-Meier method. The comparison of time-to-event analyses in dependence of differentially expressed miR levels was performed with the log-rank test. Local relapse was defined as an event (= in-breast recurrence). Multivariate analyses were performed with the forward stepwise regression (Cox regression). A raw p value <0.05 was considered significant. To adjust for multiple testing, the rather strict Bonferroni correction was applied in this phase of the study in order to test potential markers for statistical robustness [52].
Abbreviations
- ACI-rats:
-
August-Copenhagen-Irish rats
- BCT:
-
breast conserving therapy
- EBCTG:
-
Early Breast Cancer Trialist Group
- EMT:
-
epithelial-mesenchymal transition
- ER-α:
-
estrogen receptor α
- FFPE:
-
formalin fixed paraffin embedded
- HRAS:
-
transforming protein p21
- KRAS:
-
V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog
- LC:
-
local control
- LIMMA:
-
linear models for microarray analysis (moderated t test)
- MIQE:
-
minimum information for publication of quantitative real-time PCR experiments
- miR (miRNA):
-
microRNA
- MTDH:
-
metadherin
- NRAS:
-
neuroblastoma ras viral oncogene homolog
- RASD1:
-
dexamethasone-induced ras-related protein 1
- RT-qPCR:
-
reverse transcription quantitative polymerase chain reaction
- TGF-β:
-
transforming growth factor
References
Malvezzi M, Bertuccio P, Levi F, La Vecchia C, Negri E. European cancer mortality predictions for the year 2014. Ann Oncol. 2014;25(8):1650–6. doi:10.1093/annonc/mdu138.
American Cancer Society. Cancer facts and figures 2014. Atlanta: American Cancer Society; 2014.
Sautter-Bihl ML, Sedlmayer F, Budach W, Dunst J, Feyer P, Fietkau R, et al. One life saved by four prevented recurrences? Update of the Early Breast Cancer Trialists confirms: postoperative radiotherapy improves survival after breast conserving surgery. Strahlenther Onkol. 2012;188(6):461–3. doi:10.1007/s00066-012-0092-5.
Darby S, McGale P, Correa C, Taylor C, Arriagada R, Clarke M, et al. Effect of radiotherapy after breast-conserving surgery on 10-year recurrence and 15-year breast cancer death: meta-analysis of individual patient data for 10,801 women in 17 randomised trials. Lancet. 2011;378(9804):1707–16. doi:10.1016/S0140-6736(11)61629-2.
Bartelink H, Horiot JC, Poortmans P, Struikmans H, Van den Bogaert W, Barillot I, et al. Recurrence rates after treatment of breast cancer with standard radiotherapy with or without additional radiation. N Engl J Med. 2001;345(19):1378–87. doi:10.1056/NEJMoa010874.
Whelan T, MacKenzie R, Julian J, Levine M, Shelley W, Grimard L, et al. Randomized trial of breast irradiation schedules after lumpectomy for women with lymph node-negative breast cancer. J Natl Cancer Inst. 2002;94(15):1143–50.
Bentzen SM, Agrawal RK, Aird EG, Barrett JM, Barrett-Lee PJ, Bliss JM, et al. The UK Standardisation of Breast Radiotherapy (START) Trial A of radiotherapy hypofractionation for treatment of early breast cancer: a randomised trial. Lancet Oncol. 2008;9(4):331–41. doi:10.1016/S1470-2045(08)70077-9.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thurlimann B, Senn HJ. 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(8):1736–47. doi:10.1093/annonc/mdr304.
Voduc KD, Cheang MC, Tyldesley S, Gelmon K, Nielsen TO, Kennecke H. Breast cancer subtypes and the risk of local and regional relapse. J Clin Oncol. 2010;28(10):1684–91. doi:10.1200/JCO.2009.24.9284.
Lazaridis G, Lambaki S, Karayannopoulou G, Eleftheraki AG, Papaspirou I, Bobos M, et al. Prognostic and predictive value of p-Akt, EGFR, and p-mTOR in early breast cancer. Strahlenther Onkol. 2014;190(7):636–8. doi:10.1007/s00066-014-0620-6. 40–5.
de Souza Rocha Simonini P, Breiling A, Gupta N, Malekpour M, Youns M, Omranipour R, et al. Epigenetically deregulated microRNA-375 is involved in a positive feedback loop with estrogen receptor alpha in breast cancer cells. Cancer Res. 2010;70(22):9175–84. doi:10.1158/0008-5472.CAN-10-1318.
Munagala R, Aqil F, Vadhanam MV, Gupta RC. MicroRNA ‘signature’ during estrogen-mediated mammary carcinogenesis and its reversal by ellagic acid intervention. Cancer Lett. 2013;339(2):175–84. doi:10.1016/j.canlet.2013.06.012.
Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 2004;101(9):2999–3004. doi:10.1073/pnas.0307323101.
Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, et al. MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005;65(16):7065–70. doi:10.1158/0008-5472.CAN-05-1783.
Adams BD, Furneaux H, White BA. The micro-ribonucleic acid (miRNA) miR-206 targets the human estrogen receptor-alpha (ERalpha) and represses ERalpha messenger RNA and protein expression in breast cancer cell lines. Mol Endocrinol. 2007;21(5):1132–47. doi:10.1210/me.2007-0022.
Castellano L, Giamas G, Jacob J, Coombes RC, Lucchesi W, Thiruchelvam P, et al. The estrogen receptor-alpha-induced microRNA signature regulates itself and its transcriptional response. Proc Natl Acad Sci U S A. 2009;106(37):15732–7. doi:10.1073/pnas.0906947106.
Bhat-Nakshatri P, Wang G, Collins NR, Thomson MJ, Geistlinger TR, Carroll JS, et al. Estradiol-regulated microRNAs control estradiol response in breast cancer cells. Nucleic Acids Res. 2009;37(14):4850–61. doi:10.1093/nar/gkp500.
Di Leva G, Piovan C, Gasparini P, Ngankeu A, Taccioli C, Briskin D, et al. Estrogen mediated-activation of miR-191/425 cluster modulates tumorigenicity of breast cancer cells depending on estrogen receptor status. PLoS Genet. 2013;9(3):e1003311. doi:10.1371/journal.pgen.1003311.
McCafferty MP, McNeill RE, Miller N, Kerin MJ. Interactions between the estrogen receptor, its cofactors and microRNAs in breast cancer. Breast Cancer Res Treat. 2009;116(3):425–32. doi:10.1007/s10549-009-0429-7.
Zhao JJ, Lin J, Yang H, Kong W, He L, Ma X, et al. MicroRNA-221/222 negatively regulates estrogen receptor alpha and is associated with tamoxifen resistance in breast cancer. J Biol Chem. 2008;283(45):31079–86. doi:10.1074/jbc.M806041200.
Tessel MA, Krett NL, Rosen ST. Steroid receptor and microRNA regulation in cancer. Curr Opin Oncol. 2010;22(6):592–7. doi:10.1097/CCO.0b013e32833ea80c.
Zhou X, Marian C, Makambi KH, Kosti O, Kallakury BV, Loffredo CA, et al. MicroRNA-9 as potential biomarker for breast cancer local recurrence and tumor estrogen receptor status. PLoS One. 2012;7(6):e39011. doi:10.1371/journal.pone.0039011.
Lowery AJ, Miller N, Devaney A, McNeill RE, Davoren PA, Lemetre C, et al. MicroRNA signatures predict oestrogen receptor, progesterone receptor and HER2/neu receptor status in breast cancer. Breast Cancer Res. 2009;11(3):R27. doi:10.1186/bcr2257.
Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, et al. Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer. 2006;5:24. doi:10.1186/1476-4598-5-24.
Blenkiron C, Goldstein LD, Thorne NP, Spiteri I, Chin SF, Dunning MJ, et al. MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype. Genome Biol. 2007;8(10):R214. doi:10.1186/gb-2007-8-10-r214.
Le Quesne J, Caldas C. Micro-RNAs and breast cancer. Mol Oncol. 2010;4(3):230–41. doi:10.1016/j.molonc.2010.04.009.
Selth LA, Townley S, Gillis JL, Ochnik AM, Murti K, Macfarlane RJ, et al. Discovery of circulating microRNAs associated with human prostate cancer using a mouse model of disease. Int J Cancer. 2012;131(3):652–61. doi:10.1002/ijc.26405.
Zhang X, Yan Z, Zhang J, Gong L, Li W, Cui J, et al. Combination of hsa-miR-375 and hsa-miR-142-5p as a predictor for recurrence risk in gastric cancer patients following surgical resection. Ann Oncol. 2011;22(10):2257–66. doi:10.1093/annonc/mdq758.
Hamamoto J, Soejima K, Yoda S, Naoki K, Nakayama S, Satomi R, et al. Identification of microRNAs differentially expressed between lung squamous cell carcinoma and lung adenocarcinoma. Mol Med Rep. 2013;8(2):456–62. doi:10.3892/mmr.2013.1517.
Lee KH, Lee JK, Choi DW, Do IG, Sohn I, Jang KT, et al. Postoperative prognosis prediction of pancreatic cancer with seven microRNAs. Pancreas. 2015;44(5):764–8. doi:10.1097/MPA.0000000000000346.
Ward A, Balwierz A, Zhang JD, Kublbeck M, Pawitan Y, Hielscher T, et al. Re-expression of microRNA-375 reverses both tamoxifen resistance and accompanying EMT-like properties in breast cancer. Oncogene. 2013;32(9):1173–82. doi:10.1038/onc.2012.128.
Hoppe R, Achinger-Kawecka J, Winter S, Fritz P, Lo WY, Schroth W, et al. Increased expression of miR-126 and miR-10a predict prolonged relapse-free time of primary oestrogen receptor-positive breast cancer following tamoxifen treatment. Eur J Cancer. 2013;49(17):3598–608. doi:10.1016/j.ejca.2013.07.145.
Lyng MB, Laenkholm AV, Sokilde R, Gravgaard KH, Litman T, Ditzel HJ. Global microRNA expression profiling of high-risk ER+ breast cancers from patients receiving adjuvant tamoxifen mono-therapy: a DBCG study. PLoS One. 2012;7(5):e36170. doi:10.1371/journal.pone.0036170.
Jonsdottir K, Janssen SR, Da Rosa FC, Gudlaugsson E, Skaland I, Baak JP, et al. Validation of expression patterns for nine miRNAs in 204 lymph-node negative breast cancers. PLoS One. 2012;7(11):e48692. doi:10.1371/journal.pone.0048692.
Giricz O, Reynolds PA, Ramnauth A, Liu C, Wang T, Stead L, et al. Hsa-miR-375 is differentially expressed during breast lobular neoplasia and promotes loss of mammary acinar polarity. J Pathol. 2012;226(1):108–19. doi:10.1002/path.2978.
Kovalchuk O, Tryndyak VP, Montgomery B, Boyko A, Kutanzi K, Zemp F, et al. Estrogen-induced rat breast carcinogenesis is characterized by alterations in DNA methylation, histone modifications and aberrant microRNA expression. Cell Cycle. 2007;6(16):2010–8.
Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, et al. A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009;10(6):R64. doi:10.1186/gb-2009-10-6-r64.
Sperveslage J, Hoffmeister M, Henopp T, Kloppel G, Sipos B. Establishment of robust controls for the normalization of miRNA expression in neuroendocrine tumors of the ileum and pancreas. Endocrine. 2014;46(2):226–30. doi:10.1007/s12020-014-0202-5.
Davoren PA, McNeill RE, Lowery AJ, Kerin MJ, Miller N. Identification of suitable endogenous control genes for microRNA gene expression analysis in human breast cancer. BMC Mol Biol. 2008;9:76. doi:10.1186/1471-2199-9-76.
Witkos TM, Koscianska E, Krzyzosiak WJ. Practical aspects of microRNA target prediction. Curr Mol Med. 2011;11(2):93–109.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120(1):15–20. doi:10.1016/j.cell.2004.12.035.
Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007;39(10):1278–84. doi:10.1038/ng2135.
Vaidyanathan G, Cismowski MJ, Wang G, Vincent TS, Brown KD, Lanier SM. The Ras-related protein AGS1/RASD1 suppresses cell growth. Oncogene. 2004;23(34):5858–63. doi:10.1038/sj.onc.1207774.
Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40(Web Server issue):W478–83. doi:10.1093/nar/gks402.
Smith TE, Lee D, Turner BC, Carter D, Haffty BG. True recurrence vs. new primary ipsilateral breast tumor relapse: an analysis of clinical and pathologic differences and their implications in natural history, prognoses, and therapeutic management. Int J Radiat Oncol Biol Phys. 2000;48(5):1281–9.
Haffty BG, Carter D, Flynn SD, Fischer DB, Brash DE, Simons J, et al. Local recurrence versus new primary: clinical analysis of 82 breast relapses and potential applications for genetic fingerprinting. Int J Radiat Oncol Biol Phys. 1993;27(3):575–83.
Huang E, Buchholz TA, Meric F, Krishnamurthy S, Mirza NQ, Ames FC, et al. Classifying local disease recurrences after breast conservation therapy based on location and histology: new primary tumors have more favorable outcomes than true local disease recurrences. Cancer. 2002;95(10):2059–67. doi:10.1002/cncr.10952.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22. doi:10.1373/clinchem.2008.112797.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. doi:10.1006/meth.2001.1262.
Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3. doi:10.2202/1544-6115.1027.
Benjamini Y, Hochberg Y. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57(1):289–300.
Shaffer JP. Multiple hypothesis-testing. Annu Rev Psychol. 1995;46:561–84. doi:10.1146/annurev.psych.46.1.561.
Acknowledgements
We thank Lisa Kronberger and Pia Luger for their participation in part of the molecular analysis. This study was funded by the Paracelsus Medical University Salzburg, grant no.: PMU-FFF E12/16/085-SED.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
FZ designed the study concept, helped carry out the molecular and bioinformatics analysis, and wrote the manuscript. CHK and BZ carried out the molecular analysis, wrote the methodological part of the manuscript, and critically reviewed it. FH and CS carried out the Munich part of the molecular studies, helped design the concept, and critically reviewed the manuscript. UB carried out the bioinformatics analysis and critically reviewed the manuscript. GF, RR, and TF provided part of the clinical data and critically reviewed the manuscript. HD and PS contributed to the study concept and the bioinformatics analysis and critically reviewed the manuscript. FS helped design the study, supervised it, and wrote the manuscript. All authors read and approved the final manuscript.
Additional files
Additional file 1: Figure S1.
Comparison of hsa-miR-375 levels in the pilot study. The most prominent single miR that could differentiate relapse from control patients was hsa-miR-375 (LIMMA, raw p value 0.009). The box plot shows relapse versus control group, and the expression values of hsa-miR-375 are shown on the y-axis. (PDF 24.3 kb)
Additional file 2: Table S2.
Comparison of subgroups according to receptor status. This table shows that in ER-α+ patients, the levels of hsa-miR-375 are higher than in TNBC and ERα− patients, with fold changes of 1.87 and 4.17. The p value was estimated with the non-parametric Mann-Whitney U test. None of the fold changes is significant. The group of her2-negative patients has slightly higher levels of hsa-miR-375 than the her2-positive patient group, again the difference is not significant. (PDF 124 kb)
Additional file 3: Figure S3.
Target prediction for hsa-miR-375. The Venn diagram shows the number of predicted targets from TargetScan (blue) http://www.targetscan.org/and PITA (red) http://genie.weizmann.ac.il/pubs/mir07/mir07_data.html; 15 genes were suggested by both algorithms. (PDF 44.5 kb)
Additional file 4: Figure S4.
Predicted targets of hsa-miR-375 and cancer pathways. Gene enrichment analysis of the predicted targets by GeneCoDis3 software (http://genecodis.cnb.csic.es/). RASD1 was assigned to the PI3K-pathway (http://www.pantherdb.org/pathway/). The numbers on the x-axis indicate how many of the 15 “overlapping” genes in Additional file 3: Figure S3 are found in a specific pathway. (PDF 20.2 kb)
Additional file 5: Figure S1.
Study flow-chart. This additional file shows the various steps of the study. (PDF 42.8 kb)
Additional file 6: Table S1.
RT-qPCR target information. This table provides information on the exact sequences of the primers that were used for RT-qPCR including the Qiagen™ catalog number. (PDF 29.9 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zehentmayr, F., Hauser-Kronberger, C., Zellinger, B. et al. Hsa-miR-375 is a predictor of local control in early stage breast cancer. Clin Epigenet 8, 28 (2016). https://doi.org/10.1186/s13148-016-0198-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13148-016-0198-1