- Research
- Open Access
- Published:
Bioinformatics-based identification of miR-542-5p as a predictive biomarker in breast cancer therapy
Hereditas volume 155, Article number: 17 (2018)
Abstract
Background
Tamoxifen is the first-line hormone therapy for estrogen receptor alpha positive (ERα+) breast cancer. However, about 40% of patients with ERα + breast cancer who receive tamoxifen therapy eventually develop resistance resulting in a poor prognosis. The aim of this study was to mine available data sets in the Gene Expression Omnibus (GEO) database, including in vitro (cell lines) and in vivo (tissue samples), and to identify all miRNAs associated with tamoxifen resistance (TamR) in breast cancer. Secondly, this study aimed to predict the key gene regulatory networks of newly found TamR-related miRNAs and evaluate the potential role of the miRNAs and targets as potential prognosis biomarkers for breast cancer patients.
Result
Microarray data sets from two different studies were used from the GEO database: 1. GSE66607: miRNA of MCF-7 TamR cells; 2. GSE37405: TamR tissues. Differentially expressed microRNAs (miRNAs) were identified in both data sets and 5 differentially expressed miRNAs were found to overlap between the two data sets. Profiles of GSE37405 and data from the Kaplan-Meier Plotter Database (KMPD) along with Gene Expression Profiling Interactive Analysis (GEPIA) were used to reveal the relationship between these 5 miRNAs and overall survival. The results showed that has-miR-542-5p was the only miRNA associated with overall survival of ERα + breast cancer patients who received adjuvant tamoxifen. Targets of has-miR-542-5p were predicted by miRanda and TargetScan, and the mRNA expression of the three 3 target gene, Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Beta (YWHAB), Lymphocyte Antigen 9 (LY9), and Secreted Frizzled Related Protein 1 (SFRP1) were associated with overall survival in 2 different databases. Copy-number alterations (CNAs) of SFRP1 confer survival disadvantage to breast cancer patients and alter the mRNA expression of SFRP1 in cBioPortal database.
Conclusion
This study indicates that miRNA has-miR-542-5p is associated with TamR and can predict prognosis of breast cancer patients. Furthermore, has-miR-542-5p may be acting through a mechanism involving the target genes YWHAB, LY9, and SFRP1. Overall, has-miR-542-5p is a predictive biomarker and potential target for therapy of breast cancer patients.
Background
Breast cancer is the most common malignancy among women and second leading cause of cancer death in the USA [1,2,3]. Tamoxifen is a selective estrogen receptor modulator and one of the most effective adjuvant treatments for estrogen receptor alpha positive (ERα+) breast cancer patients in clinical practice. Unfortunately, about 40% of breast cancer patients with ERα + who receive tamoxifen therapy develop tamoxifen resistance (TamR) [4, 5].
The mechanisms of TamR include aberrant ERα expression, altered signal transduction pathways, imbalance of co-regulatory proteins, genetic polymorphisms in tamoxifen metabolism [6, 7], and epigenetic modifications including expression changes of microRNAs (miRNAs) [8]. miRNAs are a class of small non-coding RNAs that post-transcriptionally control the translation and stability of mRNAs [9] by forming specific base-pairing interactions between the 5′ end of the miRNA (seed region) and the miRNA response elements within the coding region or untranslated regions (UTRs) of mRNAs, leading to mRNA destabilization and/or translational inhibition [10]. Fundamentally, miRNAs provide an extra level of regulation to gene expression.
The role of miRNAs in the progression of endocrine-resistant breast cancer is of intense interest as they are promising biomarkers and therapeutic targets to counter metastatic disease [11]. Consequently, studies are rapidly emerging on the role of miRNAs in endocrine resistant breast cancer [12, 13]. Several miRNAs have been associated with TamR in MCF-7 cells exposed to tamoxifen (1 μM) for >12 months [14]. For example, expressions of miR-190b and miR-516a-5p were altered in TamR cells and were predictive of treatment outcome in a cohort of ERα + breast cancer patients receiving adjuvant tamoxifen mono-therapy [13]; miR-519a confers TamR by targeting a network of tumor-suppressor genes in ERα + breast cancer [15]; miR-27b is epigenetically downregulated in TamR breast cancer cells due to promoter methylation and also regulates tamoxifen sensitivity by targeting High Mobility Group Box 3 [16]; miR-320a sensitizes TamR breast cancer cells to tamoxifen by targeting CAMP Regulated Phosphoprotein 19 and ERRγ in breast cancer cells and tissues [17]; and miR-378a-3p regulates tamoxifen sensitivity in MCF-7 cells through targeting Golgi Transport 1A [18].
In this study, we used a meta-analysis approach to discover which miRNAs are associated with TamR across multiple studies in order to reveal strong miRNA candidates and down-stream target genes that can be used as predictive biomarkers and novel targets for therapy. To accomplish this, we mined the Gene Expression Omnibus (GEO) database and analyzed all the available in vitro (cell lines) and in vivo (tissue samples) data sets for miRNAs associated with TamR in breast cancer, identified all the miRNAs related to TamR, predicted the key gene regulatory network of a novel miRNA associated with TamR, and evaluated the potential role of this miRNA and target genes as prognosis biomarkers for breast cancer patients. This is the first study to indicate and validate the potential association of miR-542-5p and its downstream genes Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Beta (YWHAB), Lymphocyte Antigen 9 (LY9), and Secreted Frizzled Related Protein 1 (SFRP1) with the prognosis of breast cancer.
Methods
Microarray data
From the GEO database we searched all available data on tamoxifen treated samples of breast cancer tissue (with resistance status) and all TamR breast cancer cell samples (Additional file 1), and found two sets of miRNA microarrays: GSE66607 and GSE37405. GSE66607 employed an Affymetrix GeneChip® miRNA 3.0 microarray to identify differentially expressed miRNAs between tamoxifen sensitive MCF7 parent cells and tamoxifen-induced resistant cells. The dataset of GSE37405 employed an Exiqon mercury LNA microRNA array on 152 ERα + primary tumors from high-risk breast cancer patients, all patients had received adjuvant tamoxifen as mono-therapy (median clinical follow-up: 4.6 years). Patients were divided into two groups, patients without recurrence (denoted N n = 76), and recurrence (denoted R n = 76) [12].
Quality control
To ensure accuracy and reliability of this experimental design, the median value of NUSE (Normalized Unscaled Standard Errors) and the median value of RLE (Relative Log Expression) were used as evaluation criteria for each microarray. This method has been previously published [19].
Identification of differentially expressed genes
In GSE66607, differentially expressed genes were analyzed by Gene Cloud of Biotechnology Information (GCBI) working platform (GMINIX Informatics Ltd. Co, Shanghai, P. R. China), the rank is sorted by the absolute value of the fold change. Differentially expressed genes in GSE37405 were analyzed using GEO2R, an R-based web application that helps users analyze GEO data [20]. Data points with a p-value <0.05 were retained (https://www.ncbi.nlm.nih.gov/geo/geo2r/?acc=GSE37405). All data are displayed as log2 transformed.
Strategy
As shown in Fig. 1 we used a multi-step strategy to identify miRNAs that become differentially expressed in TamR breast cancer cells and tissues. First, significantly changed miRNAs in GSE66607 (cells) and GSE37405 (tissues) were identified. Second, GSE37405 and Kaplan-Meier Plotter Database (KMPD) - an online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1809 patients who accept endocrine therapy [21] with miRpower of breast cancer as the validated set [21] were used to screen overlapped miRNAs in step 1 that related to overall survival. Third, Targets of miRNAs were predicted by miRanda and TargetScan, to reveal the relationship between miRNA target genes and overall survival, the following tools and databases were used: KMPD, GEO (Affymetrix HGU133A and HGU133 + 2 microarrays) [22], and Gene Analysis Differential genes analysis (GEPIA). GEPIA is a newly developed interactive web server using a standard processing pipeline for analyzing the RNA sequencing expression data of 9736 tumors and 8587 normal samples in the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression projects [23]. In addition, Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and TCGA Provisional data in cBioPortal [24] were used to investigate the relationship between mutation/ Copy-number alterations (CNAs) and patients overall survival as well as the mRNA expression of target genes in breast cancer.
The multi-step strategy used in this study. First, significantly changed miRNAs in GSE66607 (cells for miRNA in MCF-7/TamR vs MCF-7) and GSE37405 (tissues for breast cancer tissues) were identified. Second, GSE37405 and Kaplan-Meier Plotter Database (KMPD) for miRpower of breast cancer as the validated set were used to screen overlapped miRNAs in step 1 that related to overall survival. Third, for the target mRNA genes, KMPD, and Gene Expression Profiling Interactive Analysis (GEPIA), were used to reveal the relationship between target genes and overall survival
Statistical analysis
For each triplicate of microarray data, analysis of variance for comparisons between more than two groups was used in analyzing the CNAs and mRNA expression. A p-value of <0.05 was considered statistically significant. Differences in dichotomous variables between groups were compared using Pearson’s chi-square statistic or Fisher’s exact test. Differences in continuous variables between groups were compared using the median score test. Kaplan-Meier overall survival curves were computed and compared using the log-rank statistic. Adjusted p-values were computed using Sidak’s procedure for pairwise comparisons after a significant log-rank test in overall survival value.
Result
Identification of differentially expressed genes
There were 91 differentially expressed genes (48 up-regulated and 43 down-regulated miRNA genes) between normal and TamR MCF-7 cells in GSE66607 (Additional file 2). In GSE37405, 325 miRNAs were substituted as differentially expressed genes (Additional file 3). The two differentially expressed gene sets were overlapped and 5 miRNAs (Fig. 2a, hsa-miR-324-3p, hsa-miR-486-5p, hsa-miR-940, hsa-miR-542-5p, hsa-miR-421) were identified(Additional file 4).
Survival of tamoxifen resistance specific markers. a. The two Gene Expression Omnibus (GEO) sets (GSE66607 and GSE37405) were overlapped and 5 miRNAs (Fig. 2a, hsa-miR-324-3p, hsa-miR-486-5p, hsa-miR-940, hsamiR-542-5p, hsa-miR-421) were obtained. b. Association of has-miR-542-5p expression with tamoxifen resistance outcome in GSE37405. Survival curves were generated using the Kaplan-Meier method and the log rank test was used to assess the statistical significance of differences. (Log Rank =0.019), Red: high expression; black: low expression. c. In the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) data set, the high expression of has-miR-542-5p was also associated with overall survival in patients with ER- positive endocrine therapy in 695 patients (logrankP = 0.00041, HR = 0.6(0.45–0.8). Red: high expression; black: low expression
Survival analysis of TamR specific markers
Kaplan-Meier analysis of GSE37405 showed that high expression of has-miR-542-5p in TamR breast cancer tissues conferred a survival advantage to patients with tamoxifen systemic treatment [logRank (Mantel-Cox) p = 0.019, HR = 0.31(0.11–0.875), Fig. 2b]. In the METABRIC data set of KMPD, the high expression of has-miR-542-5p was also associated with overall survival in 695 patients with ERα + breast cancer receiving endocrine therapy [logrankP = 0.00041, HR = 0.6 (0.45–0.8), Fig. 2c]. Expression of the other 4 TamR-related miRNAs was not associated with overall survival in the GSE37405 dataset, but 2 of them [hsa-miR-486-5p: logrankP = 0.04900, HR = 0.74 (0.55–1); hsa-miR-421: logrankP = 0.05800, HR = 1.33 (0.99–1.8)] seemed to be conferred a survival relationship to patients with tamoxifen systemic treatment in the KMPD dataset (Additional file 4).
miRNA target gene prediction and overall survival
1 miRNA target prediction
Since miRNAs negatively regulate gene expression, upregulated miRNAs may result in downregulated target mRNAs. Targets of miRNAs were predicted by miRanda and TargetScan(Additional file 5). A total of three has-miR-542-5p target mRNAs were closely related to overall survival (Fig. 3) of breast cancer patients in both KMPD and GEPIA, and these are: YWHAB, LY9, and SFRP1.
The mRNA expression of has-miR-542-5p target genes and overall survival value. The mRNA expression of has-miR-542-5p target genes including Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Beta (YWHAB), Lymphocyte Antigen 9 (LY9), and Secreted Frizzled Related Protein 1 (SFRP1) associated with overall survival in breast cancer. The left panel (YWHAB a, LY9 c, SFRP1 e) were verified in Gene Expression Profiling Interactive Analysis (GEPIA) database Red: high expression; blue: low expression; and the right panel (YWHAB b, LY9 d, SFRP1 f) results comes from KMPD. Red: high expression; black: low expression
2 survival analysis of target genes
High expression of YWHAB confers a survival disadvantage to breast cancer patients in both GEPIA with TCGA samples [logRank p = 0.0057, HR (high) =1.6, Fig. 3a] and KMPD with GEO samples [217718_s_at, logRank p = 8e-04, HR = 1.44(1.16–1.79), Fig. 3b]. While high expression of the other two target genes LY9 [GEPIA: logrankp = 0.038, HR (high) = 0.71, Fig. 3c; KMPD: 210370_s_at, logRankp = 0.00014, HR = 0.66(0.53–0.82), Fig. 3d] and SFRP1 [GEPIA: logrankP = 0.028, HR(high) = 0.7, Fig. 3e; KMPD: 202037_s_at, logRank p = 0.00031, HR = 0.72(0.58–0.9), Fig. 3f] conferred a survival advantage to breast cancer patients.
The cBioPortal was used to identify CNAs associated with survival in YWHAB, LY9 and SFRP1 and the relationship between CNAs and mRNA expression. There are 10% cases had genetic alterations in SFRP1 (Fig. 4a), and these alterations confer survival disadvantage to breast cancer patients in TCGA and METABRIC cohort [Fig. 4d logrank P = 0.011, HR = 1.245 (1.051–1.474)]. The mRNA expression of SFRP1 was lower in cases with deep deletion (p = 0.502), shallow deletion (p = 0.000), gain (p = 0.007) and amplification (p = 0.046) than that without change in METABRIC samples (Fig. 5e), which was similar to the results in TCGA (Fig. 5f). Alterations of the other 2 genes were not associated with overall survival in the cBioPortal (and Fig. 4b and c). And the difference between the mRNA expression and CNAs of these 2 genes were not obvious (Fig. 5a and b for YWHAB, c and d for LY9). In addition, the mutations of the predicted genes were also investigated in the cBioPortal, and they are all rare ones and little data is available for further analysis.
Copy-number alterations of has-miR-542-5p target genes and overall survival value. Has-miR-542-5p target genes expression associated with copy-number alterations in breast cancer. a: The overall alteration situation of the 3 genes. b-d: The relationship between Copy-number alterations (CNAs) and overall survival value verified in cBioPortal (TCGA and Metabric) database (b for YWHAB, c for LY9, d for SFRP1) black: case without alterations; green: case with alterations
CNAs and mRNA expression of has-miR-542-5p target genes. Has-miR-542-5p target genes expression associated with CNAs in breast cancer. The left panel (a for YWHAB, c for LY9, e for SFRP1) were verified in Metabric database; and the right panel (b for YWHAB, d for LY9, f for SFRP1) results comes from KMPD. These samples were separated in to five categories, 1-Deep deletion (homozygous deletion), 2-Shallow deletion (heterozygous deletion), 3-diploid, 4-Gain (low level amplification), and 5-amplification (high level amplification). p < 0.05 compared to diploid group
Discussion
In the present study, we mined all miRNA profiling data of breast cancer available in the GEO database to examine miRNA changes associated with TamR, and identified a novel target for TamR diagnosis and prognosis of breast cancer patients. It is interesting to note that we found highly expressed has-miR-542-5p is associated with better overall survival in patients with ERα + breast cancer receiving endocrine therapy. In addition, three targets of has-miR-542-5p (YWHAB, LY9 and SFRP1) were also associated with prognosis in these patients.
Understanding the molecular mechanism of TamR is of critical importance for diagnosis and treatment of breast cancer patients [25, 26]. The results for 4 of the 5 TamR-related miRNAs found in this study corroborate with past findings. For instance, down regulation of miR-27b-3p enhances TamR in breast cancer by increasing NR5A2 and CREB1 expression [27]; deep sequencing analysis of the breast cancer cell line MCF7 over-expressing miR-335-5p and miR-335-3p revealed that this miRNA duplex represses gene expressions of factors involved in the ERα signaling pathway; thus, enhancing TamR in MCF7 cells [28]; miR-29b-1 has tumor suppressor activity in TamR cells [29]; and miR-221/222 is involved in cell cycle deregulation in breast cancer drug resistance [30]. Therefore, the novel strategy used in our study provides a reliable and robust meta-analysis approach to find strong miRNA candidates for prognosis biomarkers or novel therapeutic targets. Moreover, in the current study all the available miRNA data of breast cancer in GEO were included in the analysis and investigated with a comprehensive examination and validation. With this approach we identified a novel and very promising TamR-associated miRNA.
Examining the relationship between TamR-related miRNAs and overall survival using data from GSE37405 and METABRIC, we identified a new candidate factor for prognosis, has-miR-542-5p. Two of the other 4 miRNAs, has-miR-486-5p and has-miR-421, might also be associated with overall survival in ERα + breast cancer patients receiving endocrine therapy. A significant change in their expressions was detected in the larger sample set (Additional file 4); however, there are discrepancies in the expressions of these 2 miRNAs between the GEO and KMPD databases, which might be caused by different sample sizes of the datasets. The role that has-miR-486-5p plays in the prognosis of breast cancer patients has also been indicated in other investigations [31].
Studies have described the function of has-mir-542-5p in cancers. Kaplan-Meier survival analysis revealed significantly reduced event-free survival in patients with low miR-542-5p expression (p < 0.001, log-rank test) in neuroblastoma, which is consistent with the results of our study [32, 33]. The difference of transcript/translation expression for has-miR-542-5p between local and metastatic disease is very obvious in neuroblastoma [32], breast cancers [34], and endometrial carcinosarcoma [35]. Our research shows that high expression of has-miR-542-5p confers a survival advantage to patients receiving tamoxifen treatment (Fig. 2b), though has-miR-542-5p is overexpressed in breast cancer with lymph node metastasis [34]. Together, these studies suggest that has-miR-542-5p might be closely associated with epithelial-mesenchymal transition (EMT), a process by which epithelial cells lose their cell polarity and cell-cell adhesion ability, thus gaining migratory and invasive properties. Has-miR-421, another differentially expressed miRNA found in our study to be related to overall survival in the KMPD dataset, has previously been found to inhibit breast cancer metastasis by target E-cadherin [36] or metastasis associated 1 [37]. Recurring tumors after treatment with chemotherapy have an increased proportion of cells expressing EMT-associated genes, thus highlighting the importance of EMT in therapeutic resistance in breast cancer [38]. Other processes might also be involved in TamR. For instance, inhibition of miR-542-5p promoted efficient DNA repair and activated expression of Notch reporters [39]. In another report, miR-542-5p induced double-strand DNA breaks and reactive oxygen species accumulation in transfected cells [40], which might be related to the significantly increased sensitivity of therapy-resistant derivatives of breast cancer cells to a complex of DNA ligase III inhibitors that increased the number of DNA double-strand breaks [41].
MiRNAs affect the development of various diseases by inhibiting the expression of their target mRNAs [42]. Understanding the relevance of miRNA and mRNA expression patterns in TamR breast cancer is important to elucidate the relationship between the pathophysiological process and gene transcription. We analyzed the relationship of the downstream target genes of miR-542-5p with survival of the patients in the KMPD and GEPIA databases, and found 3 target genes (YWHAB, LY9 and SFRP1) associated with overall survival. YWHAB plays a key role in cellular proliferation and oncogenic transformation and has been reported to be involved in the development of breast cancer [43] and the activation of tumor/metastasis pathways and inhibition of apoptosis [44]. Knockdown of YWHAB upregulated connectivity hubs, successfully inhibited in vitro proliferation, colony formation, anchorage independence, migration and invasion in MDA-MB-231 cells [45]. This is consistent with our results that high YWHAB expression is associated with poor overall survival in breast cancer patients (Fig. 3a, b), which could also be explained by low expression of has-mir-542-5p for these patients shown in Fig. 2b, c.
The hazard ratio of the other two target genes (LY9 and SFRP1) could not be explained by negative regulation of miR-542-5p. This may be due to the different mechanism or other regulating factors in breast cancer. LY9 belongs to the signaling lymphocyte activation molecules family of immunomodulatory receptors and interacts with the adaptor molecule signaling lymphocyte activation molecules-associated protein [46], which is involved in hematopoietic stem cell differentiation pathways and lineage-specific markers [47]. Diseases associated with LY9 include Systemic Lupus Erythematosus. To our knowledge, there is no previously published report associating LY9 with breast cancer; thus, further studies are needed to evaluate the role of LY9 in the development of TamR breast cancer. SFRPs act as soluble modulators of Wnt-pathway signaling, which is associated with cancer proliferation, migration and invasion of breast cancer cells. Studies have reported that miR-27a affects in vivo and in vitro proliferation, migration, and invasion of breast cancer cells through targeting the SFRP1 gene via Wnt/β-catenin signaling pathway [48]. Expression of SFRP1 leads to antitumor synergy of combined HDAC and methyltransferase inhibitors in chemo-resistant cancers [49]. Furthermore, SFRP1 DNA methylation was reported to accumulate with age in normal-appearing kidney tissues and be associated with increased renal cancer risk [50], which indicates the importance of SFRP1 in cancer. In the present study, CNAs of SFRP1 gene is also closely related to prognosis (Fig. 4d), which might be related to the decreased mRNA expression in patients with these genetic alterations (Fig. 5 e and f). The mechanism of how has-mir-542-5p regulates TamR in breast cancer and how it influences its downstream target genes warrants further study.
Conclusion
In conclusion, our study provides a comprehensive bioinformatics analysis of TamR-related miRNAs and reveals the association of certain miRNAs with the overall survival of breast cancer patients. The current study indicates that has-miR-542-5p is related to TamR and affects the prognosis of breast cancer patients. Regulating the target genes YWHAB, LY9 and SFRP1 may be a possible mechanism of has-miR-542-5p in TamR of breast cancer.
Abbreviations
- CNAs:
-
Copy-number alterations
- EMT:
-
Epithelial-mesenchymal transition
- ERα + :
-
Estrogen receptor alpha positive
- GCBI:
-
Gene Cloud of Biotechnology Information
- GEO:
-
Gene Expression Omnibus
- GEPIA:
-
Gene Expression Profiling Interactive Analysis
- KMPD:
-
Kaplan-Meier Plotter Database
- LY9:
-
Lymphocyte Antigen 9
- METABRIC:
-
Molecular Taxonomy of Breast Cancer International Consortium
- MIRNA:
-
MicroRNA
- SFRP1:
-
Secreted Frizzled Related Protein 1
- TamR:
-
Tamoxifen resistance
- TCGA:
-
The Cancer Genome Atlas
- UTRs:
-
Untranslated regions
- YWHAB:
-
Tyrosine 3-Monooxygenase/Tryptophan 5-Monooxygenase Activation Protein Beta
References
Gradishar WJ, Anderson BO, Balassanian R, Blair SL, Burstein HJ, Cyr A, Elias AD, Farrar WB, Forero A, Giordano SH, et al. Invasive breast cancer version 1.2016, NCCN clinical practice guidelines in oncology. Journal of the National Comprehensive Cancer Network : JNCCN. 2016;14:324–54.
Sivasubramaniam PG, Zhang B-L, Zhang Q, Smith JS, Zhang B, Tang Z-H, Chen G-J, Xie X-M, X-Z X, Yang H-J, et al. Breast cancer disparities: a multicenter comparison of tumor diagnosis, characteristics, and surgical treatment in China and the US. Oncologist. 2015;20:1044–50.
Zheng S, Bai JQ, Li J, Fan JH, Pang Y, Song QK, Huang R, Yang HJ, Xu F, Lu N, Qiao YL. The pathologic characteristics of breast cancer in China and its shift during 1999-2008: a national-wide multicenter cross-sectional image over 10 years. Int J Cancer. 2012;131:2622–31.
Ring A, Dowsett M. Mechanisms of tamoxifen resistance. Endocr Relat Cancer. 2004;11:643–58.
Song QK, Li J, Huang R, Fan JH, Zheng RS, Zhang BN, Zhang B, Tang ZH, Xie XM, Yang HJ, et al. Age of diagnosis of breast cancer in china: almost 10 years earlier than in the United States and the European union. Asian Pacific journal of cancer prevention : APJCP. 2014;15:10021–5.
Viedma-Rodriguez R, Baiza-Gutman L, Salamanca-Gomez F, Diaz-Zaragoza M, Martinez-Hernandez G, Ruiz Esparza-Garrido R, Velazquez-Flores MA, Arenas-Aranda D. Mechanisms associated with resistance to tamoxifen in estrogen receptor-positive breast cancer (review). Oncol Rep. 2014;32:3–15.
Liu R, Guo CX, Zhou HH. Network-based approach to identify prognostic biomarkers for estrogen receptor-positive breast cancer treatment with tamoxifen. Cancer biology & therapy. 2015;16:317–24.
Badia E, Oliva J, Balaguer P, Cavailles V. Tamoxifen resistance and epigenetic modifications in breast cancer cell lines. Curr Med Chem. 2007;14:3035–45.
Kasinski AL, Slack FJ. Epigenetics and genetics. MicroRNAs en route to the clinic: progress in validating and targeting microRNAs for cancer therapy. Nat Rev Cancer. 2011;11:849–64.
Brummer A, Hausser J. MicroRNA binding sites in the coding region of mRNAs: extending the repertoire of post-transcriptional gene regulation. BioEssays : news and reviews in molecular, cellular and developmental biology. 2014;36:617–26.
Muluhngwi P, Klinge CM. Roles for miRNAs in endocrine resistance in breast cancer. Endocr Relat Cancer. 2015;22:R279–300.
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:e36170.
Joshi T, Elias D, Stenvang J, Alves CL, Teng F, Lyng MB, Lykkesfeldt AE, Brunner N, Wang J, Gupta R, et al. Integrative analysis of miRNA and gene expression reveals regulatory networks in tamoxifen-resistant breast cancer. Oncotarget. 2016;7:57239–53.
Coser KR, Wittner BS, Rosenthal NF, Collins SC, Melas A, Smith SL, Mahoney CJ, Shioda K, Isselbacher KJ, Ramaswamy S, Shioda T. Antiestrogen-resistant subclones of MCF-7 human breast cancer cells are derived from a common monoclonal drug-resistant progenitor. Proc Natl Acad Sci U S A. 2009;106:14536–41.
Ward A, Shukla K, Balwierz A, Soons Z, Konig R, Sahin O, Wiemann S. MicroRNA-519a is a novel oncomir conferring tamoxifen resistance by targeting a network of tumour-suppressor genes in ER+ breast cancer. J Pathol. 2014;233:368–79.
Li X, Wu Y, Liu A, Tang X. MiR-27b is epigenetically downregulated in tamoxifen resistant breast cancer cells due to promoter methylation and regulates tamoxifen sensitivity by targeting HMGB3. Biochem Biophys Res Commun. 2016;477:768–73.
Lu M, Ding K, Zhang G, Yin M, Yao G, Tian H, Lian J, Liu L, Liang M, Zhu T, Sun F. MicroRNA-320a sensitizes tamoxifen-resistant breast cancer cells to tamoxifen by targeting ARPP-19 and ERRgamma. Sci Rep. 2015;5:8735.
Ikeda K, Horie-Inoue K, Ueno T, Suzuki T, Sato W, Shigekawa T, Osaki A, Saeki T, Berezikov E, Mano H, Inoue S. miR-378a-3p modulates tamoxifen sensitivity in breast cancer MCF-7 cells through targeting GOLT1A. Sci Rep. 2015;5:13170.
Collisson EA, Sadanandam A, Olson P, Gibb WJ, Truitt M, Gu S, Cooc J, Weinkle J, Kim GE, Jakkula L, et al. Subtypes of pancreatic ductal adenocarcinoma and their differing responses to therapy. Nat Med. 2011;17:500–3.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41:D991–5.
Lanczky A, Nagy A, Bottai G, Munkacsy G, Szabo A, Santarpia L, Gyorffy B. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160:439–46.
Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, Szallasi Z. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123:725–31.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al: Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Science signaling 2013, 6:pl1.
Egeland NG, Lunde S, Jonsdottir K, Lende TH, Cronin-Fenton D, Gilje B, Janssen EA, Soiland H. The role of MicroRNAs as predictors of response to tamoxifen treatment in breast cancer patients. Int J Mol Sci. 2015;16:24243–75.
Sun H, Wang G, Peng Y, Zeng Y, Zhu QN, Li TL, Cai JQ, Zhou HH, Zhu YS. H19 lncRNA mediates 17beta-estradiol-induced cell proliferation in MCF-7 breast cancer cells. Oncol Rep. 2015;33:3045–52.
Zhu J, Zou Z, Nie P, Kou X, Wu B, Wang S, Song Z, He J. Downregulation of microRNA-27b-3p enhances tamoxifen resistance in breast cancer by increasing NR5A2 and CREB1 expression. Cell Death Dis. 2016;7:e2454.
Martin EC, Conger AK, Yan TJ, Hoang VT, Miller DF, Buechlein A, Rusch DB, Nephew KP, Collins-Burow BM, Burow ME. MicroRNA-335-5p and -3p synergize to inhibit estrogen receptor alpha expression and promote tamoxifen resistance. FEBS Lett. 2017;591:382–92.
Muluhngwi P, Krishna A, Vittitow SL, Napier JT, Richardson KM, Ellis M, Mott JL, Klinge CM. Tamoxifen differentially regulates miR-29b-1 and miR-29a expression depending on endocrine-sensitivity in breast cancer cells. Cancer Lett. 2016;388:230–8.
Karimi G, Alamolhida N, Behravan J, Mosaffa F. MiR 221/222 as new players in tamoxifen resistance. Curr Pharm Des. 2016;
Rask L, Balslev E, Sokilde R, Hogdall E, Flyger H, Eriksen J, Litman T. Differential expression of miR-139, miR-486 and miR-21 in breast cancer patients sub-classified according to lymph node status. Cell Oncol. 2014;37:215–27.
Schulte JH, Schowe B, Mestdagh P, Kaderali L, Kalaghatgi P, Schlierf S, Vermeulen J, Brockmeyer B, Pajtler K, Thor T, et al. Accurate prediction of neuroblastoma outcome based on miRNA expression profiles. Int J Cancer. 2010;127:2374–85.
Bray I, Tivnan A, Bryan K, Foley NH, Watters KM, Tracey L, Davidoff AM, Stallings RL. MicroRNA-542-5p as a novel tumor suppressor in neuroblastoma. Cancer Lett. 2011;303:56–64.
Wang B, Li J, Sun M, Sun L, Zhang X. miRNA expression in breast cancer varies with lymph node metastasis and other clinicopathologic features. IUBMB Life. 2014;66:371–7.
Castilla MA, Moreno-Bueno G, Romero-Perez L, Van De Vijver K, Biscuola M, Lopez-Garcia MA, Prat J, Matias-Guiu X, Cano A, Oliva E, Palacios J. Micro-RNA signature of the epithelial-mesenchymal transition in endometrial carcinosarcoma. J Pathol. 2011;223:72–80.
Zhang W, Shi S, Jiang J, Li X, Lu H, Ren F. LncRNA MEG3 inhibits cell epithelial-mesenchymal transition by sponging miR-421 targeting E-cadherin in breast cancer. Biomedicine & pharmacotherapy = Biomedecine & pharmacotherapie. 2017;91:312–9.
Pan Y, Jiao G, Wang C, Yang J, Yang W. MicroRNA-421 inhibits breast cancer metastasis by targeting metastasis associated 1. Biomedicine & pharmacotherapy = Biomedecine & pharmacotherapie. 2016;83:1398–406.
Fischer KR, Durrans A, Lee S, Sheng J, Li F, Wong ST, Choi H, El Rayes T, Ryu S, Troeger J, et al. Epithelial-to-mesenchymal transition is not required for lung metastasis but contributes to chemoresistance. Nature. 2015;527:472–6.
Rane JK, Ylipaa A, Adamson R, Mann VM, Simms MS, Collins AT, Visakorpi T, Nykter M, Maitland NJ. Construction of therapeutically relevant human prostate epithelial fate map by utilising miRNA and mRNA microarray expression data. Br J Cancer. 2015;113:611–5.
Faraonio R, Salerno P, Passaro F, Sedia C, Iaccio A, Bellelli R, Nappi TC, Comegna M, Romano S, Salvatore G, et al. A set of miRNAs participates in the cellular senescence program in human diploid fibroblasts. Cell Death Differ. 2012;19:713–21.
Tobin LA, Robert C, Nagaria P, Chumsri S, Twaddell W, Ioffe OB, Greco GE, Brodie AH, Tomkinson AE, Rassool FV. Targeting abnormal DNA repair in therapy-resistant breast cancers. Molecular cancer research : MCR. 2012;10:96–107.
Berindan-Neagoe I, Monroig Pdel C, Pasculli B, Calin GA. MicroRNAome genome: a treasure for cancer diagnosis and therapy. CA Cancer J Clin. 2014;64:311–36.
Zhang Y, Li M, Wei L, Zhu L, Hu S, Wu S, Ma S, Gao Y. Differential protein expression in perfusates from metastasized rat livers. Proteome Sci. 2013;11:37.
Akekawatchai C, Roytrakul S, Kittisenachai S, Isarankura-Na-Ayudhya P, Jitrapakdee S. Protein profiles associated with Anoikis resistance of metastatic MDA-MB-231 breast cancer cells. Asian Pacific journal of cancer prevention : APJCP. 2016;17:581–90.
Tilli TM, Carels N, Tuszynski JA, Pasdar M. Validation of a network-based strategy for the optimization of combinatorial target selection in breast cancer therapy: siRNA knockdown of network targets in MDA-MB-231 cells as an in vitro model for inhibition of tumor development. Oncotarget. 2016;7:63189–203.
Graham DB, Bell MP, McCausland MM, Huntoon CJ, van Deursen J, Faubion WA, Crotty S, McKean DJ. Ly9 (CD229)-deficient mice exhibit T cell defects yet do not share several phenotypic characteristics associated with SLAM- and SAP-deficient mice. J Immunol. 2006;176:291–300.
Hagberg N, Theorell J, Schlums H, Eloranta ML, Bryceson YT, Ronnblom L. Systemic lupus erythematosus immune complexes increase the expression of SLAM family members CD319 (CRACC) and CD229 (LY-9) on plasmacytoid dendritic cells and CD319 on CD56(dim) NK cells. J Immunol. 2013;191:2989–98.
Kong LY, Xue M, Zhang QC, CF S. In vivo and in vitro effects of microRNA-27a on proliferation, migration and invasion of breast cancer cells through targeting of SFRP1 gene via Wnt/beta-catenin signaling pathway. Oncotarget. 2017;8:15507–19.
Cooper SJ, von Roemeling CA, Kang KH, Marlow LA, Grebe SK, Menefee ME, Tun HW, Colon-Otero G, Perez EA, Copland JA. Reexpression of tumor suppressor, sFRP1, leads to antitumor synergy of combined HDAC and methyltransferase inhibitors in chemoresistant cancers. Mol Cancer Ther. 2012;11:2105–15.
Atschekzei F, Hennenlotter J, Janisch S, Grosshennig A, Trankenschuh W, Waalkes S, Peters I, Dork T, Merseburger AS, Stenzl A, et al. SFRP1 CpG island methylation locus is associated with renal cell cancer susceptibility and disease recurrence. Epigenetics. 2012;7:447–57.
Acknowledgements
We would like to thank Dr. Xi Li for expert assistance.
Funding
This study was funded by National Natured Science Foundation of China (grant number No.81503563).
Availability of data and materials
All data produced by the study are disclosed in the manuscript and the additional files.
Author information
Authors and Affiliations
Contributions
Participated in research design: Q-NZ, YG. Performed data analysis: Q-NZ, YG. Wrote or contributed to the writing of the manuscript: Q-NZ, HR and YG. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Data sets of tamoxifen resistant breast cancer samples in the GEO. (XLSX 10Â kb)
Additional file 2:
Identification of DEGs in GSE66607. (XLSX 14 kb)
Additional file 3:
Identification of DEGs in GSE 67916. (XLSX 26 kb)
Additional file 4:
Rlationship between miRNA and overall survival rate for the overlap miRNAs. (XLSX 9Â kb)
Additional file 5:
Predicted targets of has-miR-542-5p. (XLSX 32 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
Zhu, QN., Renaud, H. & Guo, Y. Bioinformatics-based identification of miR-542-5p as a predictive biomarker in breast cancer therapy. Hereditas 155, 17 (2018). https://doi.org/10.1186/s41065-018-0055-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s41065-018-0055-7
Keywords
- Breast cancer
- Tamoxifen resistance
- Bioinformatics
- miRNA-542-5p