Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 11 September 2018
Sec. Cancer Genetics
This article is part of the Research Topic Accomplishments, Collaborative Projects and Future Initiatives in Breast Cancer Genetic Predisposition View all 13 articles

Prognostic Genes of Breast Cancer Identified by Gene Co-expression Network Analysis

\r\nJianing TangJianing Tang1Deguang KongDeguang Kong2Qiuxia CuiQiuxia Cui1Kun WangKun Wang3Dan ZhangDan Zhang3Yan Gong*Yan Gong4*Gaosong Wu*Gaosong Wu1*
  • 1Department of Thyroid and Breast Surgery, Zhongnan Hospital of Wuhan University, Wuhan, China
  • 2Department of General Surgery, Zhongnan Hospital of Wuhan University, Wuhan, China
  • 3Department of Thyroid and Breast Surgery, Tongji Hospital, Huazhong University of Science and Technology, Wuhan, China
  • 4Department of Biological Repositories, Zhongnan Hospital of Wuhan University, Wuhan, China

Breast cancer is one of the most common malignancies. The molecular mechanisms of its pathogenesis are still to be investigated. The aim of this study was to identify the potential genes associated with the progression of breast cancer. Weighted gene co-expression network analysis (WGCNA) was used to construct free-scale gene co-expression networks to explore the associations between gene sets and clinical features, and to identify candidate biomarkers. The gene expression profiles of GSE1561 were selected from the Gene Expression Omnibus (GEO) database. RNA-seq data and clinical information of breast cancer from TCGA were used for validation. A total of 18 modules were identified via the average linkage hierarchical clustering. In the significant module (R2 = 0.48), 42 network hub genes were identified. Based on the Cancer Genome Atlas (TCGA) data, 5 hub genes (CCNB2, FBXO5, KIF4A, MCM10, and TPX2) were correlated with poor prognosis. Receiver operating characteristic (ROC) curve validated that the mRNA levels of these 5 genes exhibited excellent diagnostic efficiency for normal and tumor tissues. In addition, the protein levels of these 5 genes were also significantly higher in tumor tissues compared with normal tissues. Among them, CCNB2, KIF4A, and TPX2 were further upregulated in advanced tumor stage. In conclusion, 5 candidate biomarkers were identified for further basic and clinical research on breast cancer with co-expression network analysis.

Introduction

Breast cancer is the most frequently diagnosed malignancy and the second leading cause of cancer death in females worldwide, accounting for 30% of cancer diagnoses and 14% of cancer death. In 2017, it was estimated that nearly 252,710 new cases were diagnosed in the United States, with ~40,610 deaths (1). Therapeutic strategies of breast cancer have been markedly improved. A number of treatments such as surgery, chemotherapy, radiotherapy, hormone therapy, and targeted therapy are available for breast cancer (2). However, the patients with distant metastases were usually diagnosed with a late stage and nearly incurable (3). Moreover, 30% patients diagnosed with early stage were easy to recur in distant organs even after surgery of removing the primary tumor (4). The classification of breast cancer affects treatment decision and prognosis: hormone-based therapy for ER+ patients; targeted therapy for HER2+ patients; and poorly differentiated cancer often has the worse prognosis (57).

Inheritance plays an important role in the development of breast cancer. BRCA1 and BRCA2 are 2 biomarkers which are currently used clinically to assess the familial breast cancer risk. BRCA-associated breast cancer has relatively distinct pathologic characteristics. Up to 20% women with triple-negative breast cancer present BRCA mutations, while BRCA mutations occur less common in general population (8, 9). HER2 expression was found to be upregulated in over 30% patients with breast cancer (10). Previous data suggested that high HER2 levels not only indicated prognostic value, but also affected treatment decisions. Lapatinib and trastuzumab presented dramatically therapeutic effects in patients with HER2-positive breast cancer (11, 12). Expression levels of hormone receptors (ER/PR) predicted the efficacy of endocrine therapies, and their upregulation was often associated with a favorable prognosis (13). Ki-67 was reported to be associated with disease-free survival (14). High CXCR4 levels were associated with lymph node metastasis and distant metastasis (15). Despite the substantial improvements in the treatment of breast cancer, to date, the ability to treat the advanced ones is still limited due to the lack of precise molecular targets for breast cancer (16). Therefore, it is important to explore the molecule mechanisms involved in the occurrence and development of breast cancer. More novel candidate genes are needed to improve the early diagnosis and treatment decisions.

Co-expression analysis is a powerful technique to construct free-scale gene co-expression networks. The weighted gene co-expression network analysis (WGCNA) was widely used to analyze large-scale data sets and to find modules of highly correlated genes. WGCNA was successfully used to explore the associations between gene sets and clinical features, and to identify candidate biomarkers (17). Thus, we described the correlation patterns among genes through a systematic biology method based on WGCNA and identified novel biomarkers associated with breast cancer prognosis.

Materials and Methods

Data Procession

A workflow of this study was indicated in Figure 1. The gene expression profiles of GSE1561 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE1561) submitted by Richard Iggo et al. was downloaded from the Gene Expression Omnibus (GEO) database. The GSE1561 was an expression profiling based on GPL96 platform (Affymetrix Human Genome U133A Array) and contained 49 samples. Most patients had 2 trucut biopsies taken, and both biopsies were analyzed from 2 tumors to test the reproducibility of the technique. Repeat amplifications and duplicate biopsies clustered together suggested that biological variation was greater than technical variation in this data set. The results of immunohistochemistry (IHC) also suggested the high quality of this data set (18). Robust Multi-array Average (RMA) algorithm in affy package within Bioconductor (http://www.bioconductor.org) in R was used to preprocess the gene expression profile data. After background correction, quantile normalization and probe summarization, the data set with 12,413 genes was further processed, and the top 50% most variant genes by analysis of variance (6,206 genes) were selected for WGCNA analysis.

FIGURE 1
www.frontiersin.org

Figure 1. Flow chart of data preparation, processing, analysis, and validation.

Co-expression Network Construction

After validation, the expression data profile of these 6,206 genes were constructed to a gene co-expression network using WGCNA package in R (Supplementary Data Sheet 1) (17). The analysis was performed as described previously (17).

The adjacency matrix aij which calculated the connection strength between each pair of nodes was calculated as follows:

sij=|cor(xi,xj)|aij=Sijβ

Where Xi and Xj were vectors of expression value for gene i and j, sij represented the Pearson's correlation coefficient of gene i and gene j, aij encoded the network connection strength between gene i and gene j. In the presented study, the power of β = 9 (scale free R2 = 0.95) was selected as the soft-thresholding parameter to ensure a scale-free network. In the co-expression network, genes with high absolute correlations were clustered into the same module. WGCNA method not only considers the association between the 2 connected genes, but also takes associated genes into account. Modules were also identified via hierarchical clustering of the weighting coefficient matrix. To further identify functional modules in the co-expression network with these 6,206 genes, the topological overlap measure (TOM) representing the overlap in shared neighbors, was calculated using the adjacency matrix.

TOMi,j=K=1NAi,k · Ak,j+Ai,jmin (Ki,Kj)+1Ai,j

Where A is the weighted adjacency matrix given by Aij=|cor(xi,xj)|β and β= 9 is the soft thresholding power. According to the TOM-based dissimilarity measure with a minimum size (gene group) of 30 for the gene dendrogram, average linkage hierarchical clustering was conducted, and genes with similar expression profiles were classified into the same gene modules using the DynamicTreeCut algorithm.

Identification of Clinical Significant Modules

Two approaches were used to identify modules associated with clinical information of breast cancer. First, module eigengenes (MEs) were defined as the first principal component of each gene module and the expression of MEs was considered as a representative of all genes in a given module. The correlation between MEs and clinical trait was calculated to identify the clinical significant module. In addition, the gene significance (GS) was defined as mediated p-value of each gene (GS = lgP) in the linear regression between gene expression and the clinical traits. Then, the module significance (MS) were defined as the average GS of all the genes involved in the module.MS was measured to incorporate clinical information into the co-expression network. Module significance (MS) was defined as the average absolute gene significance measured for all genes in a given module.

Gene Ontology and Pathway Enrichment Analysis

DAVID (http://david.abcc.ncifcrf.gov/) is a database for annotation, visualization and integrated discovery. Gene Ontology (GO) and KEGG pathway analysis of differentially expressed mRNAs were carried out using DAVID (version 6.8) online tools: functional annotation. The ontology contains three categories: biological process (BP), molecular function (MF), and cellular component (CC). Enriched GO terms and KEGG pathways were identified according to the cut-off criterion of adjusted P < 0.001.

Hub Gene Identification and Validation

The connectivity of genes was measured by absolute value of the Pearson's correlation. Genes with high within-module connectivity were considered as hub genes of the modules (cor.geneModuleMembership > 0.8). Hub genes inside a given module tended to have a strong correlation with certain clinical trait, which was measured by absolute value of the Pearson's correlation (cor.geneTraitSignificance > 0.2). To validate the hub genes, the clinical information and RNA sequencing data of breast cancer were obtained from the Cancer Genome Atlas Project database (TCGA, https://cancergenome.nih.gov/). The mRNA sequencing data was normalized using edgeR package in R language. The Human Protein Atlas (http://www.proteinatlas.org) was also used to validate the immunohistochemistry of candidate hub genes. The direct link to these images in the human protein atlas are as follows: http://www.proteinatlas.org/ENSG00000112029-FBXO5/tissue/breast#img (FBXO5 in normal tissue); http://www.proteinatlas.org/ENSG00000112029-FBXO5/pathology/tissue/breast$+$cancer#img (FBXO5 in tumor tissue); http://www.proteinatlas.org/ENSG00000157456-CCNB2/tissue/breast#img (CCNB2 in normal tissue); http://www.proteinatlas.org/ENSG00000157456-CCNB2/pathology/tissue/breast$+$cancer#img (CCNB2 in tumor tissue); http://www.proteinatlas.org/ENSG00000090889-KIF4A/tissue/breast#img (CCNB2 in normal tissue); http://www.proteinatlas.org/ENSG00000090889-KIF4A/pathology/tissue/breast$+$cancer#img (CCNB2 in tumor tissue); http://www.proteinatlas.org/ENSG00000065328-MCM10/tissue/breast#img (MCM10 in normal tissue); http://www.proteinatlas.org/ENSG00000065328-MCM10/pathology/tissue/breast$+$cancer#img (MCM10 in tumor tissue); http://www.proteinatlas.org/ENSG00000088325-TPX2/tissue/breast#img (TPX2 in normal tissue); http://www.proteinatlas.org/ENSG00000088325-TPX2/pathology/tissue/breast$+$cancer#img (TPX2 in tumor tissue). Survival analysis of hub genes were performed using Kaplan Meier-plotter (www.kmplot.com) (19).

Results

Weighted Co-expression Network Construction and Key Modules Identification

The samples of GSE1561 were clustered using average linkage method and Pearson's correlation method (Figure 2). The co-expression analysis was carried out to construct the co-expression network. In this study, the power of β = 9 (scale free R2 = 0.95) was selected as the soft-thresholding parameter to ensure a scale-free network (Figure 3). A total of 18 modules were identified via the average linkage hierarchical clustering. Blue module was found to have the highest association with tumor grade (Figure 4), and this module was selected as the clinical significant module for further analysis.

FIGURE 2
www.frontiersin.org

Figure 2. Clustering dendrogram of 49 samples.

FIGURE 3
www.frontiersin.org

Figure 3. Determination of soft-thresholding power in the WGCNA. (A) Analysis of the scale-free fit index for various soft-thresholding powers (β). (B) Analysis of the mean connectivity for various soft-thresholding powers. (C) Checking the scale free topology when β = 9.

FIGURE 4
www.frontiersin.org

Figure 4. Identification of modules associated with the clinical traits of breast cancer. (A) Dendrogram of all differentially expressed genes clustered based on a dissimilarity measure (1-TOM). (B) Heatmap of the correlation between module eigengenes and clinical traits of breast cancer. (C) Distribution of average gene significance and errors in the modules associated with tumor grades of breast cancer.

Gene Ontology and Pathway Enrichment Analysis

The genes in the clinical significant module were categorized into 3 functional groups (BP, CC, and MF). Clinical significant module genes in the BP group were mainly enriched in cell division, DNA replication, sister chromatid cohesion, mitotic nuclear division, and DNA replication initiation; The genes in the MF group were mainly enriched in protein binding, poly(A) RNA binding, RNA binding, and ATP binding; the genes in the CC group were significantly enriched in nucleoplasm, nucleus, nucleolus, cytosol, and cytoplasm (Figure 5). According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, our results demonstrated that these genes were mainly involved in cell cycle, DNA replication, spliceosome, ribosome biogenesis in eukaryotes and RNA transport. These results indicated that the clinical significant module genes were mainly involved in mitotic cell cycle process.

FIGURE 5
www.frontiersin.org

Figure 5. Gene ontology and pathway enrichment analysis of blue module genes. (A) Biological process analysis. (B) Cellular component analysis. (C) Molecular function analysis. (D) KEGG pathway analysis.

Identification and Validation of Hub Genes

Based the cut-off criteria (|MM| > 0.8 and |GS| > 0.2), 42 genes with high connectivity in the clinical significant module were identified as hub genes. Among them, CCNB2, FBXO5, KIF4A, MCM10, and TPX2 were negatively associated with the overall survival and relapse free survival (Figures 6, 7). Moreover, based on the TCGA data, the expression levels of these 5 genes were significantly higher in tumor tissues, especially in the triple negative breast cancers. The expression of CCNB2, KIF4A, and TPX2 were upregulated in the advanced tumor stages. ROC curve indicated that CCNB2, FBXO5, KIF4A, MCM10, and TPX2 exhibited excellent diagnostic efficiency for normal and tumor tissues (Figures 8, 9). In addition, the protein levels of these 5 genes were significantly higher in tumor tissues compared with normal tissues based on the Human Protein Atlas database (Figure 10). Since these 5 genes were all hub genes in the clinical significant module, they might have a tendency to co-express. Our results of correlation analysis demonstrated a strong correlation of mRNA expression levels between KIF4A and TPX2 (Supplementary Data Sheet 2).

FIGURE 6
www.frontiersin.org

Figure 6. Overall survival of the five hub genes in breast cancer based on Kaplan Meier-plotter. The patients were stratified into high-level group and low-level group according to median expression. (A) CCNB2. (B) FBXO5. (C) KIF4A. (D) MCM10. (E) TPX2.

FIGURE 7
www.frontiersin.org

Figure 7. Relapse free survival analysis of the five hub genes in breast cancer based on Kaplan Meier-plotter. The patients were stratified into high-level group and low-level group according to median expression (A) CCNB2. (B) FBXO5. (C) KIF4A. (D) MCM10. (E) TPX2.

FIGURE 8
www.frontiersin.org

Figure 8. Validation of CCNB2, FBXO5, KIF4A, MCM10, and TPX2. (A) The correlation of CCNB2 (A), FBXO5 (B), KIF4A (C), MCM10 (D), and TPX2 (E) expression with breast cancer molecular subtypes. (F) The correlation of CCNB2 expression with pathological stage. (G) The correlation of KIF4A expression with pathological stage. (H) The correlation of TPX2 expression with pathological stage.*P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. One-way analysis of variance (ANOVA) was used to evaluate the statistical significance of differences.

FIGURE 9
www.frontiersin.org

Figure 9. Gene expression levels of CCNB2, FBXO5, KIF4A, MCM10, and TPX2 between normal breast and tumor samples. The mRNA levels of CCNB2 (A), CCNB2 (B), FBXO5 (C), KIF4A (D), and TPX2 (E). ROC curve of CCNB2 (F), FBXO5 (G), KIF4A (H), MCM10 (I), and TPX2 (J). (A–E) *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Two-tailed Student's t-tests was used to evaluate the statistical significance of differences.

FIGURE 10
www.frontiersin.org

Figure 10. Immunohistochemistry of the five hub genes based on the Human Protein Atlas. (A) Protein levels of FBXO5 in normal tissue (staining: medium; intensity: moderate; quantity: >75%). (B) Protein levels of FBXO5 in tumor tissue (staining: high; intensity: strong; quantity: >75%). (C) Protein levels of CCNB2 in normal tissue (staining: low; intensity: moderate; quantity: <25%). (D) Protein levels of CCNB2 in tumor tissue (staining: medium; intensity: strong; quantity: <25%). (E) Protein levels of KIF4A in normal tissue (staining: low; intensity: weak; quantity: 25–75%). (F) Protein levels of KIF4A in tumor tissue (staining: high; intensity: strong; quantity: >75%). (G) Proteins level of MCM10 in normal tissue (staining: not detected; intensity: weak; quantity: <25%). (H) Protein levels of MCM10 in tumor tissue (staining: low; intensity: moderate; quantity: <25%). (I) Protein levels of TPX2 in normal tissue (staining: medium; intensity: strong; quantity: <25%). (J) Protein levels of TPX2 in tumor tissue (staining: medium; intensity: strong; quantity: <25%).

Discussion

Breast cancer seriously endangers female health, and it is easy to recur even after combined therapy. Although the treatment of breast cancer was improved during the last decades, the ability to treat the advanced ones is still limited due to the lack of precise molecular targets for breast cancer. Therefore, it is important to explore the molecule mechanisms involved in the occurrence and development of breast cancer. Better biomarkers for cancer specific prognosis and progression are highly demanded. In the presented study, we used gene expression datasets from GEO database to screen potential biomarkers related to the progression and prognosis of breast cancer. We also obtained the clinical information and RNA sequencing data of breast cancer from TCGA database for validation.

WGCNA was performed to explore gene co-expression modules associated with progression of breast cancer. A total of 6,206 most variant genes were used to construct co-expression network and 18 modules were identified. Blue module was found to have the highest association with tumor grades and 42 genes with high connectivity were screened out from the module. Among them, CCNB2, FBXO5, KIF4A, MCM10, and TPX2 were negatively associated with the overall survival (Figure 6).

CCNB2, also known as cyclin B2, is a member of cyclin family. CCNB2 was reported to regulate cell cycle by activating CDC2 kinase in eukaryotes, and inhibition of CCNB2 induced cell cycle arrest. CCNB2 was overexpressed in multiple tumors, including bladder cancer, uterine corpus endometrial carcinoma, prostate cancer, and gastric cancer (2023). In addition, compared with normal controls, the levels of serum circulating CCNB2 are higher in digestive tract cancer and lung cancer patients, and they are found to be significantly associated with tumor stage and metastasis status (24). In invasive breast carcinoma, cytoplasmic CCNB2 protein levels were significantly correlated with a poor disease specific survival. CCNB2 expression level was reported to be an independent prognostic factor for the disease specific survival of breast cancer (25). Our results indicated that CCNB2 was upregulated in breast cancer tissues compared to normal tissues, and that its expression was significantly associated with molecular subtypes of breast cancer and tumor stages (Figure 8). The underlying mechanisms of CCNB2 on tumor progression need to be further clarified.

F-Box Protein 5 (FBXO5) is a key cell cycle regulatory gene which regulates the progression to S phase and mitosis by inhibiting the anaphase promoting complex (APC). FBXO5 is overexpressed in various solid tumors. In the G0 and early G1 phases, the expression of FBXO5 is low, while in the S phase it is upregulated. In ovarian clear cell carcinoma, FBXO5 accumulation was related to mitotic errors with centrosome overduplication and abnormal spindle formation. These findings demonstrated that it might be involved in human cell cycle disorders and genomic stability to promote tumor growth (2628). In breast carcinoma tissues, FBXO5 induced proliferation through the PI3K/Akt pathway. Overexpression of FBXO5 was reported to correlate with poor prognosis. In addition, PI3K inhibitor reduced FBXO5 expression (29).

The protein encoded by Kinesin family member 4A (KIF4A) was reported to be involved in the intracellular transport of membranous organelles and chromosome integrity during mitosis. In patients with colorectal cancer, KIF4A was upregulated, and downregulation of KIF4A reduced cell proliferation in colorectal cancer cells (30). In hepatocellular carcinoma (HCC) patients, KIF4A overexpression was associated with poorer overall and disease-free survival. In HCC cells, higher levels of KIF4A dramatically increased cellular clonogenic abilities and proliferation, while KIF4A depletion caused a significant augmentation of apoptosis (31). In breast cancer, high KIF4A levels were associated with poor relapse-free survival of ER-positive patients. In tamoxifen-resistant and sensitive breast cancer cells, KIF4A knockdown significantly impeded cellular proliferation and induced apoptosis (32).

Mini-chromosome maintenance complex component 10 (MCM10) is one of the highly conserved mini-chromosome maintenance proteins. MCM10 is bound to chromatin through the interaction with MCM2-7, and plays crucial roles both in initiation and elongation during eukaryotic genome replication (33). For urothelial carcinoma, high MCM10 levels were significantly correlated with advanced tumors stages, vascular invasion, and nodal status. MCM10 overexpression also predicted poor disease-specific survival and inferior metastasis-free survival (34). In our analysis of GSE1561, MCM10 was one of the hub genes in the blue module which was significantly associated with tumor grade (Figure 3). In the validation dataset of TCGA, our results indicated that MCM10 was significantly upregulated in breast tumor tissues, and even higher in the triple negative breast cancer (Figures 8, 9).

Targeting protein for Xenopus kinesin-like protein 2 (TPX2) plays a critical role in chromosome segregation machinery during mitosis (35). It was reported to be overexpressed in multiple tumors: lung cancer, kidney renal clear cell carcinoma, hepatocellular Carcinoma, prostate cancer, and breast cancer (36). TPX2 activates PI3K/Akt pathway and upregulates matrix metalloproteinases (MMP) family members in colon cancer. Previous studies showed that TPX2 expression promoted proliferation, migration, and invasion of liver cancer and breast cancer cells via upregulating expressions of MMP2 and MMP9 (37, 38). In patients with HCC, overexpression TPX2 was correlated with worse prognosis. In addition, knockdown TPX2 in HCC cells strongly reduced cellular proliferation, induced apoptosis and inhibited EMT (39).

Co-expression analysis is a powerful technique for multigene analysis of large-scale data sets. In cancer research, co-expression analyses revealed the mRNA and microRNA expression network in multiple cancers. In the present study, we used WGCNA to construct a gene co-expression network, to measure the relationships between genes and modules, and to explore the relationships between modules and clinical traits. We also screened out a clinical significant module which was associated with the progression of breast cancer. KEGG pathway analysis demonstrated that this module was mostly involved in cell cycle. In addition, 5 hub genes, CCNB2, FBXO5, KIF4A, MCM10, and TPX2 were identified and validated to be associated with the progression and worse prognosis of breast cancer. Our results provided valuable indication for basic and clinical research on breast cancer. The underlying concept of gene co-expression analysis is guilt-by-association. The groups of genes known as co-expression modules were found to maintain a consistent expression relationship independent of phenotype, and might share a common biological role. Similar to the limitations of most other data mining methods, our results of WGCNA can be biased or invalid when dealing with technical artifacts or tissue contaminations (6). To increase the credibility of WGCNA results, TCGA RNA-seq data and IHC data from the Human Protein Atlas database were used for validation. While due to the limitation of the database, the related IHC of each sample can't be found, tumor and normal samples were from different patients.

Author Contributions

JT, YG, and GW reviewed relevant literature and drafted the manuscript. DK, KW, DZ, and QC conducted all statistical analyses. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2018.00374/full#supplementary-material

References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. (2011) 61:69–90. doi: 10.3322/caac.20107

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Shi H, Zhang L, Qu Y, Hou L, Wang L, Zheng M. Prognostic genes of breast cancer revealed by gene co-expression network analysis. Oncol Lett. (2017) 14:4535–42. doi: 10.3892/ol.2017.6779

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Redig AJ, McAllister SS. Breast cancer as a systemic disease: a view of metastasis. J Intern Med. (2013) 274:113–26. doi: 10.1111/joim.12084

PubMed Abstract | CrossRef Full Text | Google Scholar

4. McAllister SS, Gifford AM, Greiner AL, Kelleher SP, Saelzler MP, Ince TA, et al. Systemic endocrine instigation of indolent tumor growth requires osteopontin. Cell (2008) 133:994–1005. doi: 10.1016/j.cell.2008.04.045

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zhang Y, Lv F, Yang Y, Qian X, Lang R, Fan Y, et al. Clinicopathological features and prognosis of metaplastic breast carcinoma: experience of a major Chinese cancer center. PLoS ONE (2015) 10:e0131409. doi: 10.1371/journal.pone.0131409

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Clarke C, Madden SF, Doolan P, Aherne ST, Joyce H, O'Driscoll L, et al. Correlating transcriptional networks to breast cancer survival: a large-scale coexpression analysis. Carcinogenesis (2013) 34:2300–8. doi: 10.1093/carcin/bgt208

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Krishnamurti U, Silverman JF. HER2 in breast cancer: a review and update. Adv Anat Pathol. (2014) 21:100–7. doi: 10.1097/PAP.0000000000000015

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Trainer AH, Thompson E, James PA. BRCA and beyond: a genome-first approach to familial breast cancer risk assessment. Discov Med. (2011) 12:433–43.

PubMed Abstract | Google Scholar

9. Weitzel JN. The genetics of breast cancer: what the surgical oncologist needs to know. Surgical Oncol Clin N Am. (2015) 24:705–32. doi: 10.1016/j.soc.2015.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL. Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science (1987) 235:177–82.

PubMed Abstract | Google Scholar

11. Romond EH, Perez EA, Bryant J, Suman VJ, Geyer CE Jr, Davidson NE, et al. Trastuzumab plus adjuvant chemotherapy for operable HER2-positive breast cancer. N Engl J Med. (2005) 353:1673–84. doi: 10.1056/NEJMoa052122

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Geyer CE, Forster J, Lindquist D, Chan S, Romieu CG, Pienkowski T, et al. Lapatinib plus capecitabine for HER2-positive advanced breast cancer. N Engl J Med. (2006) 355:2733–43. doi: 10.1056/NEJMoa064320

PubMed Abstract | CrossRef Full Text | Google Scholar

13. M Braden A, V Stankowski R, M Engel J, A Onitilo A. Breast cancer biomarkers: risk assessment, diagnosis, prognosis, prediction of treatment efficacy and toxicity, and recurrence. Curr Pharm Design (2014) 20:4879–98. doi: 10.2174/1381612819666131125145517

CrossRef Full Text | Google Scholar

14. Kontzoglou K, Palla V, Karaolanis G, Karaiskos I, Alexiou I, Pateras I, et al. Correlation between Ki67 and breast cancer prognosis. Oncology (2013) 84:219–25. doi: 10.1159/000346475

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhao Q, Song W, He DY, Li Y. Identification of key gene modules and pathways of human breast cancer by co-expression analysis. Breast Cancer (2017) 25:213–23. doi: 10.1007/s12282-017-0817-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Meng L, Xu Y, Xu C, Zhang W. Biomarker discovery to improve prediction of breast cancer survival: using gene expression profiling, meta-analysis, and tissue validation. Oncotargets Ther. (2016) 9:6177–85. doi: 10.2147/OTT.S113855

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics (2008) 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Farmer P, Bonnefoi H, Becette V, Tubiana-Hulin M, Fumoleau P, Larsimont D, et al. Identification of molecular apocrine breast tumours by microarray analysis. Oncogene (2005) 24:4660–71. doi: 10.1038/sj.onc.1208561

PubMed Abstract | CrossRef Full Text

19. Lanczky A, Nagy A, Bottai G, Munkácsy G, Szabó A, Santarpia L, et al. 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. doi: 10.1007/s10549-016-4013-7

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Shen L, Liu M, Liu W, Cui J, Li C. Bioinformatics analysis of RNA sequencing data reveals multiple key genes in uterine corpus endometrial carcinoma. Oncol Lett. (2018) 15:205–12. doi: 10.3892/ol.2017.7346

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Han Y, Jin X, Zhou H, Liu B. Identification of key genes associated with bladder cancer using gene expression profiles. Oncol Lett. (2018) 15:297–303. doi: 10.3892/ol.2017.7310

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Huang CG, Li FX, Pan S, Xu CB, Dai JQ, Zhao XH. Identification of genes associated with castrationresistant prostate cancer by gene expression profile analysis. Mol Med Rep. (2017) 16:6803–13. doi: 10.3892/mmr.2017.7488

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Shi Q, Wang W, Jia Z, Chen P, Ma K, Zhou C. ISL1, a novel regulator of CCNB1, CCNB2 and c-MYC genes, promotes gastric cancer cell proliferation and tumor growth. Oncotarget (2016) 7:36489–500. doi: 10.18632/oncotarget.9269

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Mo ML, Chen Z, Li J, Li HL, Sheng Q, Ma HY, et al. Use of serum circulating CCNB2 in cancer surveillance. Int J Biol Markers (2010) 25:236–42. doi: 10.5301/JBM.2010.6088

CrossRef Full Text | Google Scholar

25. Shubbar E, Kovacs A, Hajizadeh S, Parris E, Nemes S, Gunnarsdóttir K, et al. Elevated cyclin B2 expression in invasive breast carcinoma is associated with unfavorable clinical outcome. BMC Cancer (2013) 13:1. doi: 10.1186/1471-2407-13-1

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Guan C, Zhang J, Zhang J, Shi H, Ni R. Enhanced expression of early mitotic inhibitor-1 predicts a poor prognosis in esophageal squamous cell carcinoma patients. Oncol Lett. (2016) 12:114–20. doi: 10.3892/ol.2016.4611

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Vaidyanathan S, Cato K, Tang L, Pavey S, Haass NK, Gabrielli BG, et al. In vivo overexpression of Emi1 promotes chromosome instability and tumorigenesis. Oncogene (2016) 35:5446–55. doi: 10.1038/onc.2016.94

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gutgemann I, Lehman NL, Jackson PK, Longacre TA. Emi1 protein accumulation implicates misregulation of the anaphase promoting complex/cyclosome pathway in ovarian clear cell carcinoma. Mod Pathol. (2008) 21:445–54. doi: 10.1038/modpathol.3801022

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Liu X, Wang H, Ma J, Xu J, Sheng C, Yang S, et al. The expression and prognosis of Emi1 and Skp2 in breast carcinoma: associated with PI3K/Akt pathway and cell proliferation. Med Oncol. (2013) 30:735. doi: 10.1007/s12032-013-0735-0

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Matsumoto Y, Saito M, Saito K, Saito K, Kanke Y, Watanabe Y, et al. Enhanced expression of KIF4A in colorectal cancer is associated with lymph node metastasis. Oncol Lett. (2018) 15:2188–94. doi: 10.3892/ol.2017.7555

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Huang Y, Wang H, Lian Y, Wu X, Zhou L, Wang J, et al. Upregulation of kinesin family member 4A enhanced cell proliferation via activation of Akt signaling and predicted a poor prognosis in hepatocellular carcinoma. Cell Death Dis. (2018) 9:141. doi: 10.1038/s41419-017-0114-4

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zou JX, Duan Z, Wang J, Sokolov A, Xu J, Chen CZ, et al. Kinesin family deregulation coordinated by bromodomain protein ANCCA and histone methyltransferase MLL for breast cancer cell growth, survival, and tamoxifen resistance. Mol Cancer Res. (2014) 12:539–49. doi: 10.1158/1541-7786.MCR-13-0459

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Izumi M, Mizuno T, Yanagi KI, Sugimura K, Okumura K, Imamoto N, et al. The Mcm2-7-interacting domain of human mini-chromosome maintenance 10 (Mcm10) protein is important for stable chromatin association and origin firing. J Biol Chem. (2017) 292:13008–21. doi: 10.1074/jbc.M117.779371

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li WM, Huang CN, Ke HL, Li CC, Wei YC, Yeh HC, et al. MCM10 overexpression implicates adverse prognosis in urothelial carcinoma. Oncotarget (2016) 7:77777–92. doi: 10.18632/oncotarget.12795

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Pan HW, Su HH, Hsu CW, Huang GJ, Wu TT. Targeted TPX2 increases chromosome missegregation and suppresses tumor cell growth in human prostate cancer. Oncotargets Ther. (2017) 10:3531–43. doi: 10.2147/OTT.S136491

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Geiger TR, Ha NH, Faraji F, Michael HT, Rodriguez L, Walker RC, et al. Functional analysis of prognostic gene expression network genes in metastatic breast cancer models. PLoS ONE (2014) 9:e111813. doi: 10.1371/journal.pone.0111813

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Liu Q, Tu K, Zhang H, Zheng X, Yao Y, Liu Q. TPX2 as a novel prognostic biomarker for hepatocellular carcinoma. Hepatol Res. (2015) 45:906–18. doi: 10.1111/hepr.12428

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Yang Y, Li DP, Shen N, Yu XC, Li JB, Song Q, et al. TPX2 promotes migration and invasion of human breast cancer cells. Asian Pac J Trop Med. (2015) 8:1064–70. doi: 10.1016/j.apjtm.2015.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liang B, Jia C, Huang Y, He H, Li J, Liao HE, et al. TPX2 level correlates with hepatocellular carcinoma cell proliferation, apoptosis, and EMT. Digest Dis Sci. (2015) 60:2360–72. doi: 10.1007/s10620-015-3730-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: breast cancer, weighted gene co-expression network analysis (WGCNA), prognosis, GEO, TCGA

Citation: Tang J, Kong D, Cui Q, Wang K, Zhang D, Gong Y and Wu G (2018) Prognostic Genes of Breast Cancer Identified by Gene Co-expression Network Analysis. Front. Oncol. 8:374. doi: 10.3389/fonc.2018.00374

Received: 21 April 2018; Accepted: 21 August 2018;
Published: 11 September 2018.

Edited by:

Luis G. Carvajal-Carmona, University of California, Davis, United States

Reviewed by:

Parvin Mehdipour, Tehran University of Medical Sciences, Iran
Tracy A. O'Mara, QIMR Berghofer Medical Research Institute, Australia

Copyright © 2018 Tang, Kong, Cui, Wang, Zhang, Gong and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yan Gong, yan.gong@whu.edu.cn
Gaosong Wu, wugaosongtj@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.