Integrated analysis of the potential roles of miRNA‑mRNA networks in triple negative breast cancer
- Authors:
- Published online on: June 9, 2017 https://doi.org/10.3892/mmr.2017.6750
- Pages: 1139-1146
-
Copyright: © Zhu et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Among all breast cancer types, 15–20% are pathologically diagnosed as triple negative breast cancer (TNBC). A typical immunohistochemical feature of TNBC tumors is that they are negative for the estrogen, progesterone and human epidermal growth factor 2 receptors (1). TNBC is highly invasive with common symptoms including breast lumps and nipple discharge. Epidemiological studies indicate that TNBC generally develops among pre-menopausal women, and particularly among young African-American women (2). It has been reported that TNBC is frequently diagnosed at stage III, and is usually correlated with increased metastasis, disease recurrence and cancer cell proliferation and invasion (3). To date, no specific treatment guidelines have been developed for TNBC. At present, the common methods for the treatment of breast cancer include surgery, radiotherapy, chemotherapy and traditional Chinese medicine. Therefore, patients with TNBC are currently treated with methods common to all breast cancer types. However, the rapid development of metastases and high recurrence rates lead to poor patient prognoses and high mortality rates (4,5).
MicroRNAs (miRNA/miR) are endogenous non-coding small RNA sequences that are 20–24 nucleotides in length. miRNAs regulate tumor cell metabolism, proliferation, invasion and metastasis (3). Each miRNA is able to regulate multiple target genes, while different miRNAs are able to regulate the same target gene. It is thought that approximately one-third of human genes are regulated by miRNAs (6). Previous studies indicate that ~70% of miRNAs in mammalian cells are located in miRNA transcription units (7). In addition, the majority of miRNAs are located in introns, which are highly conserved in different species (8). The high conservation of miRNAs among different species suggests common and important functions. Therefore, investigating the mechanisms and functions of miRNAs may be useful for the diagnosis and treatment of human cancers. To date, a number of previous studies have attempted to elucidate the mechanisms and functions of miRNAs in the diagnosis and treatment of human cancer; however, only a limited number have focused on understanding the function of miRNAs in the pathogenesis of TNBC (9,10). In addition, the mechanisms and signaling pathways regulated by miRNAs and their associated target genes require further investigation. With the development of gene expression microarrays and bioinformatics analysis technologies, an integrated and detailed investigation of the aberrant expression of miRNA and mRNA sequences in patients with TNBC was performed in the current study. The aim was to identify differentially expressed miRNA and mRNA sequences in TNBC tissues compared with paired normal adjacent tissues using published data from the Gene Expression Omnibus (GEO). In addition, bioinformatics analysis was performed to investigate the function and signaling pathways of the identified miRNA and target mRNA sequences.
Materials and methods
Microarray data
Two expression microarrays, which included miRNAs and mRNAs, were obtained from the GEO database (www.ncbi.nlm.nih.gov/geo). The accession numbers of these arrays were GSE61723 and GSE61724, which were based on the GPL16686 [HuGene-2_0-st, Affymetrix Human Gene 2.0 ST Array, transcript (gene) version] and GPL6244 [HuGene-1_0-st, Affymetrix Human Gene 1.0 ST Array, transcript (gene) version] platforms, respectively. The GSE61723 test dataset consisted of 33 TNBC tissue samples and 17 normal adjacent tissue samples. The GSE61724 validation dataset consisted of 16 TNBC tissue samples and 4 normal adjacent tissue samples. More detailed information regarding the tumor samples have been reported in a previous study (11).
Data processing and differential expression analysis
Series matrix files were employed and analyzed in the present study. As all microarray data consisted of preprocessed normalized data, a fold-change of ≥2.0 and a p-value of ≤0.05 were used as the thresholds for screening differentially-expressed miRNAs and mRNAs, using R software (version 3.3.2; www.r-project.org) and the edgeR package. Hierarchical clustering analysis was performed using Multiple Experiment Viewer software (version, 4.90; www.mev.tm4.org), and the Pearson correlation distance metric and average linkage method.
Prediction of miRNA target genes
The target genes of differentially expressed miRNAs were predicted using TargetScan software (version 7.0; www.targetscan.org) (12). In order to reduce the probability of identifying incorrect target genes, potential target mRNAs were selected using the intersection set for the differentially expressed mRNAs identified from the GSE61723 dataset, and the predicted target genes were identified using TargetScan software. The intersection set signifies that the genes existed in multiple data sets. The miRNA-mRNA interaction regulatory network was constructed using Cytoscape software (version 3.3.0; www.cytoscape.org) (13).
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis
Enrichment analyses for GO processes and KEGG signaling pathways were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID; version 6.8) software program (www.david.ncifcrf.gov) (14,15). Biological processes, cellular components and molecular functions were the main GO terms used for functional analysis of the potential targets of miRNAs (16). GO terms and KEGG pathways were selected using a threshold of P≤0.05 and a count of ≥2.
Results
Identification of differentially expressed miRNAs and mRNAs
In order to gain an improved understanding of the regulatory mechanisms underlying miRNA-mRNA interactions in TNBC, comprehensive microarray expression profiles were employed to identify differentially expressed miRNAs and mRNAs. By comparing TNBC tissues with normal adjacent tissues, five miRNAs (upregulated miR-558 and downregulated miR-320d-1, miR-548v, miR-99a and miR-21 expression) and 97 mRNAs were identified with a threshold of P≤0.05 and a fold-change of ≥2.0 (Tables I and II). According to the GenBank database (www.ncbi.nlm.nih.gov/genbank), LOC100507381 was removed due to inconsistencies with standard genome annotation processing. In addition, the gene UBE2C was permanently removed from the NCBI database, as there is not enough evidence to prove its presence in subsequent studies. Therefore, 95 differentially expressed mRNAs were identified between the 33 TNBC tissue and 17 normal adjacent tissue samples. Unsupervised hierarchical clustering analyses were performed using a fold-change threshold of ≥2.0. The results identified five differentially expressed miRNAs and 95 differentially expressed mRNAs (Figs. 1 and 2).
Table I.miRNAs displaying significantly altered expression levels in triple negative breast cancer samples from GSE61723. |
Table II.mRNAs displaying significantly altered expression patterns in triple negative breast cancer samples from GSE61723. |
Validation of miRNA expression
To confirm the five differentially expressed miRNAs identified in TNBC samples from the GSE61723 array, microarray data from the GSE61724 array was obtained from the GEO database, which included 16 TNBC tissues and 4 normal adjacent tissue samples. Two miRNAs (miRNA-99a and miRNA-21) were selected for further validation from the five significantly dysregulated miRNAs. miRNA-99a (P=1.91×10-6) and miRNA-21 (P=1.84×10-4) were observed to be differentially expressed in TNBC samples when compared with normal adjacent samples in the GSE61723 array. Therefore, the results were consistent with the miRNA expression levels in the GSE61723 dataset, which demonstrated that the five dysregulated miRNAs were able to accurately distinguish between TNBC tissues and normal tissues. Data from GSE61723 and GSE61724 are not identical and certain miRNAs in GSE61723 were unable to be found in GSE61724; therefore, validation was only performed in two miRNAs.
Prediction of miRNA target genes and functional analysis
In order to determine the putative functions and target genes of the five dysregulated miRNAs identified in TNBC samples, the intersection set of sequences from the 95 differentially expressed mRNAs and the predicted target genes were identified using TargetScan software. A total of 49 target mRNAs were included for further analysis. The of miRNA-mRNA interaction regulatory network is shown in Fig. 3. The results indicated that miRNA-320d-1 targeted the smallest number of genes, whereas the remaining four miRNAs targeted at least 10 (Fig. 3). In addition, carboxypeptidase and the small nuclear ribonucleoprotein D1 polypeptide were identified as the hub nodes connected to three miRNAs, which suggests that they may serve a key role in cancer pathogenesis (Fig. 3). Analysis of the GO processes and KEGG pathways involving the intersectional 49 target genes was performed using the DAVID tool, with thresholds of ≥2 gene counts and P-values of ≤0.05. The results identified 31 GO terms and three signaling pathways (Table III). The 10 most significant GO terms of the predicted target genes are indicated in Fig. 4. According to the biological function analysis of GO processes, a large proportion of the predicted targets of differentially expressed miRNAs encode cellular components in extracellular regions, and are involved in cell migration, proliferation, motion, adhesion and regeneration (Table III and Fig. 4). These biological functions are considered to be crucial for human tumorigenesis and metastasis (17). In addition, KEGG pathway analysis demonstrated that cytokine-cytokine receptor interactions, chemokine signaling and hematopoietic cell lineage pathways, which are associated with cell growth, proliferation and transformation, might be regulated by the predicted target genes (Table III).
Discussion
The mortality rates of patients with TNBC are high, as the vast majority of TNBC tumors are infiltrating ductal carcinomas, which often lead to visceral metastases (18). TNBC displays characteristic biological and clinical pathologic features, including rapid progression, poor prognosis and high recurrence rates (19). Surgery and chemotherapy are the major therapeutic strategies used to treat patients with TNBC. However, few studies have investigated the pathogenesis of TNBC. Therefore, there is an urgent need to identify novel therapeutic targets and treatments for patients with TNBC.
Expression microarrays have been widely applied for the diagnosis of various cancers (20,21). An important application of microarray experiments is the analysis of differentially expressed miRNAs or mRNAs between normal and tumor samples. The differentially expressed miRNA or mRNA sequences may be used for further functional, pathway or bioinformatics analyses. Microarray analysis is generally expensive; however, it produces a large quantity of data. Therefore, integrated mining of microarray data from public databases, such as the GEO database, is popular among researchers (22–24).
In the present study, five differentially expressed miRNAs (upregulated miR-558 and downregulated miR-320d-1, miR-548v, miR-99a and miR-21 expression), and 97 mRNAs were identified in TNBC tissue samples when compared with normal adjacent tissues using the GSE61723 array. GSE61724 microarray data was subsequently used to validate the five differentially expressed miRNAs, and the results demonstrated that two of the identified miRNA sequences (miRNA-99a and miRNA-21) were differentially expressed in the GSE61724 dataset. A total of 49 predicted target mRNAs of the 5 miRNAs were selected for further analysis according to the intersection set of the 97 differentially expressed mRNAs and the predicted target genes using TargetScan software. In addition, bioinformatics analyses were employed to investigate the functions and signaling pathways of the miRNA-mRNA network, in order to further investigate the possible mechanisms underlying the development of TNBC. This may provide relevant targets for the early diagnosis and treatment of patients with TNBC.
GO analysis is divided into three sections (molecular function, biological process and cellular component), which are widely used in the organizational and functional annotation of genes (25). Previous studies have primarily focused on applying GO analysis for miRNA or mRNA microarray results (26–28). In the current study, 31 GO terms were identified using the DAVID software tool with threshold values of ≥2 counts and P≤0.05. The results revealed that the predicted target genes encode regulatory factors and proteins located in extracellular regions, and may be involved in a number of biological processes, including regulation of cell migration, cell motion, cell proliferation, and cell adhesion. These processes are closely associated with the occurrence and metastasis of human tumor cells (29).
Signaling pathways are the way of facilitating the effect of extracellular signaling molecules within cells. In the present study, three main signaling pathways were demonstrated to be regulated by 49 target genes of the five differentially expressed miRNAs. These pathways included cytokine-cytokine receptor interaction, chemokine and hematopoietic cell lineage pathways. This suggests that these pathways may serve a significant role in tumor progression, prevention and survival rates (30–32). These identified signaling pathways, may provide promising novel therapeutic targets for the treatment of patients with TNBC.
The present study has a number of limitations. Firstly, due to the lack of microarray data comparing TNBC tumor tissues with normal adjacent control tissues, the results of only two GSE series were employed in the current study. Secondly, due to the small number of differentially expressed miRNAs identified, analysis of the associated functions and signaling pathways may be incomplete. In future studies, a more comprehensive analysis will be performed using a greater number of microarray datasets.
In conclusion, five differentially expressed miRNAs and 49 potential target genes were identified in TNBC samples when compared with normal adjacent control tissues. The results of integrated GO function and KEGG signaling pathway analyses suggested that the differentially expressed miRNAs and their potential target mRNAs may influence the pathogenesis of TNBC. In addition, the results suggest that these sequences may be promising biomarkers for the early diagnosis and targeted therapy of patients with TNBC.
References
Perou CM, Sørlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al: Molecular portraits of human breast tumours. Nature. 406:747–752. 2000. View Article : Google Scholar : PubMed/NCBI | |
Pierobon M and Frankenfeld CL: Obesity as a risk factor for triple-negative breast cancers: A systematic review and meta-analysis. Breast Cancer Res Treat. 137:307–314. 2013. View Article : Google Scholar : PubMed/NCBI | |
Foulkes WD, Smith IE and Reis-Filho JS: Triple-negative breast cancer. N Engl J Med. 363:1938–1948. 2010. View Article : Google Scholar : PubMed/NCBI | |
Lehmann BD, Bauer JA, Chen X, Sanders ME, Chakravarthy AB, Shyr Y and Pietenpol JA: Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. 121:2750–2767. 2011. View Article : Google Scholar : PubMed/NCBI | |
Blows FM, Driver KE, Schmidt MK, Broeks A, van Leeuwen FE, Wesseling J, Cheang MC, Gelmon K, Nielsen TO, Blomqvist C, et al: Subtyping of breast cancer by immunohistochemistry to investigate a relationship between subtype and short and long term survival: A collaborative analysis of data for 10,159 cases from 12 studies. PLoS Med. 7:e10002792010. View Article : Google Scholar : PubMed/NCBI | |
Cui Q, Yu Z, Pan Y, Purisima EO and Wang E: MicroRNAs preferentially target the genes with high transcriptional regulation complexity. Biochem Biophys Res Commun. 352:733–738. 2007. View Article : Google Scholar : PubMed/NCBI | |
Kim VN and Nam JW: Genomics of microRNA. Trends Genet. 22:165–173. 2006. View Article : Google Scholar : PubMed/NCBI | |
Rodriguez A, Griffiths-Jones S, Ashurst JL and Bradley A: Identification of mammalian microRNA host genes and transcription units. Genome Res. 14:1902–1910. 2004. View Article : Google Scholar : PubMed/NCBI | |
Yang L, Feng Y, Qi P, Xu S and Zhou Y: Mechanism of serum miR-21 in the pathogenesis of familial and triple negative breast cancer. J Biol Regul Homeost Agents. 30:1041–1045. 2016.PubMed/NCBI | |
Zhang G, Liu Z, Han Y, Wang X and Yang Z: Overexpression of miR-509 increases apoptosis and inhibits invasion via suppression of tumor necrosis factor-α in triple-negative breast cancer Hs578T cells. Oncol Res. 24:233–238. 2016. View Article : Google Scholar : PubMed/NCBI | |
Mathe A, Wong-Brown M, Morten B, Forbes JF, Braye SG, Avery-Kiejda KA and Scott RJ: Novel genes associated with lymph node metastasis in triple negative breast cancer. Sci Rep. 5:158322015. View Article : Google Scholar : PubMed/NCBI | |
Agarwal V, Bell GW, Nam JW and Bartel DP: Predicting effective microRNA target sites in mammalian mRNAs. Elife. 4:2015. View Article : Google Scholar : | |
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B and Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498–2504. 2003. View Article : Google Scholar : PubMed/NCBI | |
da W Huang, Sherman BT and Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 4:44–57. 2009.PubMed/NCBI | |
da W Huang, Sherman BT and Lempicki RA: Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37:1–13. 2009. View Article : Google Scholar : PubMed/NCBI | |
Sam LT, Mendonca EA, Li J, Blake J, Friedman C and Lussier YA: PhenoGO: An integrated resource for the multiscale mining of clinical and biological data. BMC bioinformatics. 10 Suppl 2:S82009. View Article : Google Scholar : PubMed/NCBI | |
Labi V and Erlacher M: How cell death shapes cancer. Cell Death Dis. 6:e16752015. View Article : Google Scholar : PubMed/NCBI | |
Zhang C, Wang S, Israel HP, Yan SX, Horowitz DP, Crockford S, Gidea-Addeo D, Chao KS Clifford, Kalinsky K and Connolly EP: Higher locoregional recurrence rate for triple-negative breast cancer following neoadjuvant chemotherapy, surgery and radiotherapy. Springerplus. 4:3862015. View Article : Google Scholar : PubMed/NCBI | |
Kim S, Park HS, Kim JY, Ryu J, Park S and Kim SI: Comparisons of oncologic outcomes between triple-negative breast cancer (TNBC) and Non-TNBC among patients treated with breast-conserving therapy. Yonsei Med J. 57:1192–1198. 2016. View Article : Google Scholar : PubMed/NCBI | |
Shi KQ, Lin Z, Chen XJ, Song M, Wang YQ, Cai YJ, Yang NB, Zheng MH, Dong JZ, Zhang L and Chen YP: Hepatocellular carcinoma associated microRNA expression signature: Integrated bioinformatics analysis, experimental validation and clinical significance. Oncotarget. 6:25093–25108. 2015. View Article : Google Scholar : PubMed/NCBI | |
Dai M, Chen X, Mo S, Li J, Huang Z, Huang S, Xu J, He B, Zou Y, Chen J and Dai S: Meta-signature LncRNAs serve as novel biomarkers for colorectal cancer: Integrated bioinformatics analysis, experimental validation and diagnostic evaluation. Sci Rep. 7:465722017. View Article : Google Scholar : PubMed/NCBI | |
An N, Shi X, Zhang Y, Lv N, Feng L, Di X, Han N, Wang G, Cheng S and Zhang K: Discovery of a novel immune gene signature with profound prognostic value in colorectal cancer: A model of cooperativity disorientation created in the process from development to cancer. PloS One. 10:e01371712015. View Article : Google Scholar : PubMed/NCBI | |
Volinia S and Croce CM: Prognostic microRNA/mRNA signature from the integrated analysis of patients with invasive breast cancer. Proc Natl Acad Sci USA. 110:7413–7417. 2013. View Article : Google Scholar : PubMed/NCBI | |
Zhang X, Peng Y, Jin Z, Huang W, Cheng Y, Liu Y, Feng X, Yang M, Huang Y, Zhao Z, et al: Integrated miRNA profiling and bioinformatics analyses reveal potential causative miRNAs in gastric adenocarcinoma. Oncotarget. 6:32878–32889. 2015.PubMed/NCBI | |
Lovering RC, Camon EB, Blake JA and Diehl AD: Access to immunology through the gene ontology. Immunology. 125:154–160. 2008. View Article : Google Scholar : PubMed/NCBI | |
Xiong H, Li Q, Liu S, Wang F, Xiong Z, Chen J, Chen H, Yang Y, Tan X, Luo Q, et al: Integrated microRNA and mRNA transcriptome sequencing reveals the potential roles of miRNAs in stage I endometrioid endometrial carcinoma. PLoS One. 9:e1101632014. View Article : Google Scholar : PubMed/NCBI | |
Li C, Nong Q, Solanki MK, Liang Q, Xie J, Liu X and Li Y, Wang W, Yang L and Li Y: Differential expression profiles and pathways of genes in sugarcane leaf at elongation stage in response to drought stress. Sci Rep. 6:256982016. View Article : Google Scholar : PubMed/NCBI | |
Ma G, Luo Y, Zhu H, Luo Y, Korhonen PK, Young ND, Gasser RB and Zhou R: MicroRNAs of Toxocara canis and their predicted functional roles. Parasit Vectors. 9:2292016. View Article : Google Scholar : PubMed/NCBI | |
Ha M and Kim VN: Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 15:509–524. 2014. View Article : Google Scholar : PubMed/NCBI | |
Lapkina-Gendler L, Rotem I, Pasmanik-Chor M, Gurwitz D, Sarfstein R, Laron Z and Werner H: Identification of signaling pathways associated with cancer protection in Laron syndrome. Endocr Relat Cancer. 23:399–410. 2016. View Article : Google Scholar : PubMed/NCBI | |
Shi KQ, Lin Z, Chen XJ, Song M, Wang YQ, Cai YJ, Yang NB, Zheng MH, Dong JZ, Zhang L and Chen YP: Hepatocellular carcinoma associated microRNA expression signature: Integrated bioinformatics analysis, experimental validation and clinical significance. Oncotarget. 6:25093–25108. 2015. View Article : Google Scholar : PubMed/NCBI | |
Geraldo MV and Kimura ET: Integrated analysis of thyroid cancer public datasets reveals role of post-transcriptional regulation on tumor progression by targeting of immune system mediators. PLoS One. 10:e01417262015. View Article : Google Scholar : PubMed/NCBI |