Evaluation of the Role of hsa-mir-124 in Predicting Clinical Outcome in Breast Invasive Carcinoma Based on Bioinformatics Analysis

Feng, Tongbao; Zhang, Ping; Sun, Yingxin; Han, Xiu; Tong, Jichun; Hua, Zichun

doi:https://doi.org/10.1155/2020/1839205

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Abbreviations Data Availability Ethical Approval Consent Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1839205 | https://doi.org/10.1155/2020/1839205

Evaluation of the Role of hsa-mir-124 in Predicting Clinical Outcome in Breast Invasive Carcinoma Based on Bioinformatics Analysis

Tongbao Feng,^1,2Ping Zhang,²Yingxin Sun,²Xiu Han,²Jichun Tong,³and Zichun Hua^1,4

Academic Editor: Hakan Bermek

Received14 Nov 2019

Revised10 Jan 2020

Accepted28 Jan 2020

Published04 Mar 2020

Abstract

Purpose. Breast invasive carcinoma (BRCA) is the most common malignant tumor. MiR-124 plays a tumor-suppressive role in human cancer. However, the clinical significance of miR-124 in BRCA remains unclear. The aim of this study was to evaluate the association of hsa-mir-124 expression and the clinicopathological characteristics in BRCA using database analysis. Methods. The clinical data and expression profiles of hsa-mir-124 were obtained from the cancer genome atlas for BRCA (TCGA_BRCA). Then, the prognostic value of hsa-mir-124 in BRCA was investigated using the Cox Regression test, and the association of hsa-mir-124 and pathology TNM stages and pathologic stages were measured by the Kruskal–Wallis test and Wilcox. test. In addition, the association of hsa-mir-124 and tumor molecular phenotypes was performed using the Chi-Square test. Results. We found that the overall survival of patients with high expression of hsa-mir-124-1 and hsa-mir-124-2 was better than that of patients with low expression of hsa-mir-124-1 and hsa-mir-124-2. And the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 was mainly enriched in T1/T2 stages, NO/N1 stages, and M0 stages. Then, the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 was negatively associated with tumor lymph node metastasis. Moreover, the expression of hsa-mir-124 was associated with tumor molecular phenotype in breast invasive carcinoma. Conclusion. Our findings indicated that hsa-mir-124 expressions were associated with overall survival, TNM stages, pathologic characteristics, and tumor molecular phenotype in BRCA via TCGA_BRCA database, providing a new biomarker and a potential therapeutic target for BRCA patients.

1. Introduction

Breast tumor is one of the most common malignant tumors and the leading cause of cancer-related deaths among women worldwide [1, 2]. Using the classification scheme based on the World Health Organization (WHO) criteria, various histologic types of breast invasive carcinoma (BRCA) can be diagnosed by applying cytomorphologic criteria. Invasive ductal carcinoma (IDC) and invasive lobular carcinoma (ILC) account for about 90% of BRCA, of which IDC is the most common type, making up about 80% [3–5]. With IDC, cancer cells start in a milk duct, break through the walls, and invade breast tissue [4]. Previous studies reported that patients with BRCA had the worst prognosis of the breast cancer [6–8]. Therefore, it is urgently needed to identify new molecular markers to guide the diagnosis, treatment, and prognosis of BRCA.

Emerging studies have implicated that miR-124 plays the tumor-suppressive role in the development of the human cancer and acts as a diagnostic prognostic biomarker for cancer [9–13]. For instance, downregulation of miR-124 was an independent prognostic factor in patients with colorectal cancer [11, 14]. Moreover, miR-124-3p acted as a potential marker and suppresses tumor growth in gastric cancer [9]. In humans, three subtypes of hsa-mir-124 have been identified, namely, hsa-mir-124-1 (8p23.1), hsa-mir-124-2 (8q12.3), and hsa-mir-124-3 (20q13.33), with different chromosomal locations. However, the alterations in three subtypes of hsa-mir-124 profiles in different clinicopathological characteristics of BRCA are not yet fully understood.

In the present study, we investigated the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in breast invasive carcinoma by using bioinformatics analysis of the clinical parameters and survival data in TCGA_BRCA database. We found that the subtypes of hsa-mir-124 expressions were associated with overall survival and TNM stages in BRCA. In addition, we observed that the subtypes of hsa-mir-124 expressions were also associated with pathologic characteristics and tumor molecular phenotype in BRCA. These findings suggested that the alterations in hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 profiles could provide a new biomarker and a potential therapeutic target for BRCA patients.

2. Materials and Methods

2.1. Data Acquisition

Linkedomics (http://www.linkedomics.org/) [15] was used to analyze the expression of hsa-mir-124 via a nonparametric test (attribute-dependent) in TCGA_BRCA. The expression of hsa-mir-124 was estimated by Human TCGA_BRCA miRNASeq HS miR 01/28/2016 Gene Firehose RPKM log2 and shown as (RPM, log 2 (Val + 1)) on the website. The clinical information (January 2016) regarding TCGA_BRCA was obtained from Human TCGA_BRCA Clinical 01/28/2016 Clinical Firehose on public databases. The clinicopathological data of 294 BRCA patients were used to analyze the association between the expression of hsa-mir-124-1 and hsa-mir-124-2 and the overall survival of BRCA patients; 742 BRCA patients were tested to analyze the association between the expression of hsa-mir-124-3 and the overall survival of BRCA patients. Then, the clinicopathological data of 754 patients, 744 patients, and 603 patients were enrolled to analyze the association of hsa-mir-124 expression and TNM stages, respectively, in BRCA patients. In addition, the clinicopathological data of 747 patients and 678 patients were enrolled to analyze the association of the expression of hsa-mir-124 and pathologic characteristics in BRCA. Moreover, the clinicopathological data of 43 patients and 44 patients were enrolled to analyze the association of the expression of hsa-mir-124 and tumor molecular phenotype in BRCA. And the demographic information and exclusion criteria for BRCA patients were according to the TCGA database.

2.2. Statistics Analysis

All the patients were divided into high and low expression groups according to the median expression of hsa-mir-124. The Cox regression test was used to investigate the prognostic value of hsa-mir-124 in BRCA. The Kruskal–Wallis test and Wilcox test were used to analyse the association of hsa-mir-124 and pathology TNM stages and pathologic stages in TCGA_BRCA. The chi-square test was performed to reveal the association of hsa-mir-124 and tumor molecular phenotypes in TCGA_BRCA. Differences with were considered statistically significant.

3. Results

3.1. The Expression of hsa-mir-124 Was Associated with Overall Survival in Breast Invasive Carcinoma

To investigate the prognostic value of hsa-mir-124 in BRCA, we firstly performed the Cox regression test by deep sequencing data in TCGA_BRCA. The results showed that the overall survival of patients with high expression of hsa-mir-124-1 and hsa-mir-124-2 was better than that of patients with low expression of hsa-mir-124-1 and hsa-mir-124-2 (Figures 1(a) and 1(b), N = 294, ). In addition, there was no statistically significant difference between the overall survival of patients with high expression of hsa-mir-124-3 and that of patients with low expression of hsa-mir-124-3 (Figure 1(c), N = 742, ). These results suggested that the expression of hsa-mir-124 was associated with overall survival in breast invasive carcinoma patients.

(a)

(b)

(c)

3.2. The Expression of hsa-mir-124 Was Associated with TNM Stages in Breast Invasive Carcinoma

To assess the pathological association between hsa-mir-124 and BRCA patients, we used the Kruskal–Wallis test and Wilcox test to analyse the association of hsa-mir-124 and pathology TNM stages in TCGA_BRCA. We found that the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in BRCA was mainly enriched in T1/T2 stages (Figures 2(a)–2(c), N = 754), N0/N1 stages (Figures 2(d)–2(f), N = 744), and M0 stages (Figures 2(g)–2(i), N = 603). These results indicated that hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 primarily expressed in the early stages of BRCA development.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.3. The Expression of hsa-mir-124 Was Associated with Pathologic Characteristics in Breast Invasive Carcinoma

To confirm the pathological association between hsa-mir-124 and BRCA patients, we performed the Kruskal–Wallis test to detect the association of hsa-mir-124 and pathologic stages in TCGA_BRCA. The data showed that the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in patients with BRCA was mostly concentrated in stage I/II (Figures 3(a)–3(c), N = 747). Moreover, we also performed spearman correlation to analysis the association of hsa-mir-124 and tumor lymph node metastasis in TCGA_BRCA. We found that the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 was negatively associated with tumor lymph node metastasis in BRCA (Figures 3(d)–3(f), N = 678). These findings suggested that hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 inhibited tumor growth in the BRCA development.

(a)

(b)

(c)

(d)

(e)

(f)

3.4. The Expression of hsa-mir-124 Was Associated with Tumor Molecular Phenotype in Breast Invasive Carcinoma

To further evaluate the role of hsa-mir-124 in the progression of BRCA, we performed the chi-square test to reveal the association of hsa-mir-124 and tumor molecular phenotype in TCGA-BRCA. The results showed that the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in patients with estrogen receptor (ER) positive was higher than that in patients with ER negative (Figures 4(a)–4(c), N = 43). In addition, the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in patients with progesterone receptor (PR) negative was higher than that in patients with PR positive (Figures 4(d)–4(f), N = 44). Moreover, the expression of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 in patients with human epidermal growth factor receptor (HER2) negative was higher than that in patients with HER2 positive (Figures 4(g)–4(i), N = 43). These results indicated that hsa-mir-124 was associated with the molecular phenotype of BRCA tumors, which could be used to predict the molecular phenotype of tumors by detecting the expression of hsa-mir-124 in BRCA patients.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

4. Discussion

As described in our research, we showed that noncoding RNAs were involved in breast cancer metastasis and miR-124 inhibited cell proliferation in breast cancer [16–18]. Consistent with our findings, miR-124 was decreased in colorectal tissue samples and was correlated with adverse clinical indicators and poor patient survival time [19–21]. Downregulation of plasma miR-124 expression was a predictive biomarker for prognosis of glioma [10, 22]. A large number of previous studies have suggested that miR-124 could be used as a molecular marker for tumor progression in cervical cancer [23, 24], breast cancer [12, 25, 26], lung cancer [27, 28], and osteosarcoma [29]. Numerous clinical studies have shown that BRCA is the most common type of breast cancer and has poor prognosis survival, and the pathogenesis of breast cancer involves a variety of factors and genes [4, 30]. Therefore, exploring the molecular diagnostic markers for BRCA could reduce cancer-related mortality and morbidity. In our present study, we investigated the significance of hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 expression in breast invasive carcinoma by using the TCGA_BRCA database. We found that the subtypes of hsa-mir-124 expressions were associated with overall survival and TNM stages in BRCA. Moreover, we observed that the subtypes of hsa-mir-124 expressions were also associated with pathologic characteristics and tumor molecular phenotype in BRCA.

Based on the abovementioned findings, we concluded that the abnormal expression of hsa-mir-124 in BRCA was a favorable prognostic marker associated with less aggressive tumor characteristics. However, there were still some concerns that need to be expressed. Firstly, our conclusions were based on analysis and summary of the TCGA database, and we did not conduct further experimental analysis and verification. Secondly, our findings showed that hsa-mir-124 expressions were associated with overall survival, TNM stages, pathologic characteristics, and tumor molecular phenotype in BRCA, which were basically congruent with those of previous studies. However, this study was unable to determine the practical significance of hsa-mir-124 abnormal expression as a favorable prognosis marker and the relationship between existing known prognostic markers in BRCA. Thirdly, there was currently no effective targeted therapy for hsa-mir-124, which was questionable to evaluate the benefit of the aberrant expression of hsa-mir-124 for prognostic significance in BRCA. Therefore, in future studies, we need to further research the sensitivity and specificity of hsa-mir-124 in diagnosis of BRCA and explore the feasibility of targeted hsa-mir-124 in BRCA clinical treatment.

In summary, our present study indicated that hsa-mir-124 expressions were associated with overall survival, TNM stages, pathologic characteristics, and tumor molecular phenotype in BRCA via TCGA-BRCA database. These findings suggested that the alterations in hsa-mir-124-1, hsa-mir-124-2, and hsa-mir-124-3 profiles could predict clinical outcome in breast invasive carcinoma, which provided a molecular foundation to further research the development and progression of BRCA, and ultimately provided vital insights that may lead to personalized treatment of BRCA patients.

Abbreviations

BRCA:	Breast invasive carcinoma.
IDC:	Invasive ductal carcinoma.
ILC:	Invasive lobular carcinoma.
TCGA:	The Cancer Genome Atlas.
TNM:	Stage, the primary tumor site and the regional lymph node involvement and the presence or otherwise of distant metastatic spread.
ER:	Estrogen receptor.
PR:	Progesterone receptor.
HER2:	Human epidermal growth factor receptor.

Data Availability

The data used to support the findings of this study areincluded within the article.

Ethical Approval

The patients’ information involved in our research was obtained from The Cancer Genome Atlas (TCGA). All the patients and treatments complied with the principles laid down in the Declaration of Helsinki in 1964 and its later amendments or comparable ethical standards.

Informed consent was confirmed by all the patients who participated in the TCGA_Breast invasive carcinoma project.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Z.H. was responsible for original idea and study design; T.F. was involved in analysis and interpretation of data, statistical analyses, and writing of the manuscript; P.Z., X.H., Y.S, and J.T. were involved in analysis and interpretation of data and statistical analysis; Y.S. and J.T. were responsible for technical and material support and statistical analyses; and X.H and P.Z. were responsible for data collection and technical support. All authors have read and approved the final version of this manuscript.

Acknowledgments

This work was supported by grants from the Jiangsu Natural Science Foundation (20191156 to T.F), National Natural Science Foundation of China (31601156 to T.F), Changzhou Science and Technology Program (CJ20180062 to T.F), Jiangsu Province “333 Talent Project” (BRA2018170 to P.Z), and Changzhou Municipal Health Commission Youth Talents Science and Technology Project (QN201926 to X.H). The authors would like to thank LinkedOmics for free database obtained.

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
R.-M. Feng, Y.-N. Zong, S.-M. Cao, and R.-H. Xu, “Current cancer situation in China: good or bad news from the 2018 Global Cancer Statistics?” Cancer Communications, vol. 39, no. 1, pp. 1–12, 2019.
View at: Publisher Site | Google Scholar
H. Bozkurt, I. B. Karakaya, E. Aktas, and O. Irkorucu, “Coexistence of phylloides tumour and invasive ductal cancer in the breast,” Nigerian Journal of Clinical Practice, vol. 22, no. 8, pp. 1169–1171, 2019.
View at: Google Scholar
M. Alsaleem, M. S. Toss, C. Joseph et al., “The molecular mechanisms underlying reduced E-cadherin expression in invasive ductal carcinoma of the breast: high throughput analysis of large cohorts,” Modern Pathology, vol. 32, no. 7, pp. 967–976, 2019.
View at: Publisher Site | Google Scholar
D. Flores-Díaz, C. Arce, L. Flores-Luna et al., “Impact of invasive lobular carcinoma on long-term outcomes in Mexican breast cancer patients,” Breast Cancer Research and Treatment, vol. 176, no. 1, pp. 243–249, 2019.
View at: Publisher Site | Google Scholar
C. D. Khillare, S. G. Sinai Khandeparkar, A. R. Joshi, M. M. Kulkarni, B. P. Gogate, and S. Battin, “Immunohistochemical expression of vimentin in invasive breast carcinoma and its correlation with clinicopathological parameters,” Nigerian Medical Journal, vol. 60, no. 1, pp. 17–21, 2019.
View at: Publisher Site | Google Scholar
Y. He, X. Li, Y. Meng et al., “A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma,” Journal of Cellular Biochemistry, vol. 120, no. 10, pp. 16692–16702, 2019.
View at: Publisher Site | Google Scholar
T. Du, S. Islam, C. Apostolatos, A. Apostolatos, and M. Acevedo-Duncan, “Analysis of PKC-ζ protein levels in normal and malignant breast tissue subtypes,” Oncology Letters, vol. 17, no. 2, pp. 1537–1546, 2019.
View at: Publisher Site | Google Scholar
F. Liu, H. Hu, J. Zhao et al., “miR-124-3p acts as a potential marker and suppresses tumor growth in gastric cancer,” Biomedical Reports, vol. 9, no. 2, pp. 147–155, 2018.
View at: Publisher Site | Google Scholar
X. Jia, X. Wang, X. Guo et al., “MicroRNA-124: an emerging therapeutic target in cancer,” Cancer Medicine, vol. 8, no. 12, pp. 5638–5650, 2019.
View at: Publisher Site | Google Scholar
Z. Qin and X. Liu, “miR-124, a potential therapeutic target in colorectal cancer,” OncoTargets and Therapy, vol. 12, pp. 749–751, 2019.
View at: Publisher Site | Google Scholar
G. Yan, Y. Li, L. Zhan et al., “Decreased miR-124-3p promoted breast cancer proliferation and metastasis by targeting MGAT5,” American journal of cancer research., vol. 9, no. 3, pp. 585–596, 2019.
View at: Google Scholar
Z. Cao, L. Xu, S. Zhao, and X. Zhu, “The functions of microRNA-124 on bladder cancer,” OncoTargets and Therapy, vol. 12, pp. 3429–3439, 2019.
View at: Publisher Site | Google Scholar
C. Wei, L. Lei, H. Hui, and Z. Tao, “MicroRNA-124 regulates TRAF6 expression and functions as an independent prognostic factor in colorectal cancer,” Oncology Letters, vol. 18, no. 1, pp. 856–863, 2019.
View at: Publisher Site | Google Scholar
S. V. Vasaikar, P. Straub, J. Wang, and B. Zhang, “LinkedOmics: analyzing multi-omics data within and across 32 cancer types,” Nucleic Acids Research, vol. 46, pp. D956–D963, 2018.
View at: Publisher Site | Google Scholar
T. Feng, P. Zhang, Y. Sun et al., “High throughput sequencing identifies breast cancer-secreted exosomal LncRNAs initiating pulmonary pre-metastatic niche formation,” Gene, vol. 710, pp. 258–264, 2019.
View at: Publisher Site | Google Scholar
T. Feng, F. Shao, Q. Wu et al., “miR-124 downregulation leads to breast cancer progression via LncRNA-MALAT1 regulation and CDK₄/E2F1 signal activation,” Oncotarget, vol. 7, no. 13, pp. 16205–16216, 2016.
View at: Publisher Site | Google Scholar
T. Feng, D. Xu, C. Tu et al., “miR-124 inhibits cell proliferation in breast cancer through downregulation of CDK₄,” Tumor Biology, vol. 36, no. 8, pp. 5987–5997, 2015.
View at: Publisher Site | Google Scholar
H.-Q. Ying, H.-X. Peng, B.-S. He et al., “miR-608, pre-miR-124-1 and pre-miR26a-1 polymorphisms modify susceptibility and recurrence-free survival in surgically resected CRC individuals,” Oncotarget, vol. 7, no. 46, pp. 75865–75873, 2016.
View at: Publisher Site | Google Scholar
K. Liu, H. Yao, S. Lei et al., “The miR-124-p63 feedback loop modulates colorectal cancer growth,” Oncotarget, vol. 8, no. 17, pp. 29101–29115, 2017.
View at: Publisher Site | Google Scholar
M. L. Slattery, A. J. Pellatt, F. Y. Lee et al., “Infrequently expressed miRNAs influence survival after diagnosis with colorectal cancer,” Oncotarget, vol. 8, no. 48, pp. 83845–83859, 2017.
View at: Publisher Site | Google Scholar
D. Hu and Y. Zhang, “Circular RNA HIPK3 promotes glioma progression by binding to miR-124-3p,” Gene, vol. 690, pp. 81–89, 2019.
View at: Publisher Site | Google Scholar
S. M. Wilting, R. A. van Boerdonk, F. E. Henken et al., “Methylation-mediated silencing and tumour suppressive function of hsa-miR-124 in cervical cancer,” Molecular Cancer, vol. 9, no. 1, p. 167, 2010.
View at: Publisher Site | Google Scholar
C. D. Rogeri, H. C. S. Silveira, R. L. Causin et al., “Methylation of the hsa-miR-124, SOX₁, TERT, and LMX₁A genes as biomarkers for precursor lesions in cervical cancer,” Gynecologic Oncology, vol. 150, no. 3, pp. 545–551, 2018.
View at: Publisher Site | Google Scholar
A. Sueta, Y. Yamamoto, M. Tomiguchi, T. Takeshita, M. Yamamoto-Ibusuki, and H. Iwase, “Differential expression of exosomal miRNAs between breast cancer patients with and without recurrence,” Oncotarget, vol. 8, no. 41, pp. 69934–69944, 2017.
View at: Publisher Site | Google Scholar
H. Danesh, M. Hashemi, F. Bizhani, S. M. Hashemi, and G. Bahari, “Association study of miR-100, miR-124-1, miR-218-2, miR-301b, miR-605, and miR-4293 polymorphisms and the risk of breast cancer in a sample of Iranian population,” Gene, vol. 647, pp. 73–78, 2018.
View at: Publisher Site | Google Scholar
M. M. Qi, F. Ge, X. J. Chen, C. Tang, and J. Ma, “MiR-124 changes the sensitivity of lung cancer cells to cisplatin through targeting STAT3,” European Review for Medical and Pharmacological Sciences., vol. 23, no. 12, pp. 5242–5250, 2019.
View at: Google Scholar
H. Li, X. Guo, Q. Li et al., “Long non-coding RNA 1308 promotes cell invasion by regulating the miR-124/ADAM 15 axis in non-small-cell lung cancer cells,” Cancer Management and Research, vol. 10, pp. 6599–6609, 2018.
View at: Publisher Site | Google Scholar
D. Cheng, X. Qiu, M. Zhuang, C. Zhu, H. Zou, and Z. Liu, “MicroRNAs with prognostic significance in osteosarcoma: a systemic review and meta-analysis,” Oncotarget, vol. 8, no. 46, pp. 81062–81074, 2017.
View at: Publisher Site | Google Scholar
Y. Li and Y. Wang, “Bioinformatics analysis of gene expression data for the identification of critical genes in breast invasive carcinoma,” Molecular Medicine Reports, vol. 16, no. 6, pp. 8657–8664, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Tongbao Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

566

Downloads

1225

Citations