Skip to main content

Identification of biomarkers of chromophobe renal cell carcinoma by weighted gene co-expression network analysis



Chromophobe renal cell carcinoma (ChRCC) is the second common subtype of non-clear cell renal cell carcinoma (nccRCC), which accounting for 4–5% of renal cell carcinoma (RCC). However, there is no effective bio-marker to predict clinical outcomes of this malignant disease. Bioinformatic methods may provide a feasible potential to solve this problem.


In this study, differentially expressed genes (DEGs) of ChRCC samples on The Cancer Genome Atlas database were filtered out to construct co-expression modules by weighted gene co-expression network analysis and the key module were identified by calculating module-trait correlations. Functional analysis was performed on the key module and candidate hub genes were screened out by co-expression and MCODE analysis. Afterwards, real hub genes were filter out in an independent dataset GSE15641 and validated by survival analysis.


Overall 2215 DEGs were screened out to construct eight co-expression modules. Brown module was identified as the key module for the highest correlations with pathologic stage, neoplasm status and survival status. 29 candidate hub genes were identified. GO and KEGG analysis demonstrated most candidate genes were enriched in mitotic cell cycle. Three real hub genes (SKA1, ERCC6L, GTSE-1) were selected out after mapping candidate genes to GSE15641 and two of them (SKA1, ERCC6L) were significantly related to overall survivals of ChRCC patients.


In summary, our findings identified molecular markers correlated with progression and prognosis of ChRCC, which might provide new implications for improving risk evaluation, therapeutic intervention, and prognosis prediction in ChRCC patients.


Renal cell carcinoma (RCC) is a heterogenous disease, which is composed of ccRCC and nccRCC [1]. Over the past few years, targeted therapies have significantly improved overall survival (OS) and relapse free survival (RFS) of patients with ccRCC [2]. However, ascribing to relatively low incidence (25–30%) and rare clinical trails of nccRCC, the optimal targeted therapies for nccRCC patients still remain uncertain [3]. ChRCC, taking up 4–5% of RCC, is the second common subtype of nccRCC. Although the tumor grade or stage of ChRCC is relatively low, there is no significant difference between patients with localized ChRCC and ccRCC in 5-year cancer-specific survival rates (P = 0.98) [4]. Due to the poor outcomes of ChRCC, it’s urgent to identify novel molecular biomarkers to evaluate the prognosis of ChRCC patients, which might help to assess the malignancy and provide therapeutic potential for this disease.

WGCNA is a method commonly used to explore the complex relationships between genes and phenotypes. This method is able to transform gene expression data into co-expression modules and provide insights into signaling networks that may be responsible for phenotypic traits of interests [5,6,7]. WGCNA is widely used in various biological processes, such as cancer, neuroscience and genetic data analysis, which is quite helpful for the identification of potential biomarkers or therapeutic targets [8,9,10,11]. Not only can it analyses mRNA level of tumor samples, but also work on microRNA or lncRNA datasets of neoplasms to find candidate biomarkers for prognosis and treatment [12, 13].

In this study, WGCNA method was firstly used to analyze clinic traits and gene expression data of ChRCC samples provided by TCGA database to identify key genes associated with tumor prognosis and progression. Our findings may be very beneficial to assess malignant potential of ChRCC and offer therapeutic methods to this neoplasm.


Data sources and data preprocessing

Gene expression data and patient clinic traits of ChRCC were downloaded from The Cancer Genome Atlas (TCGA) database ( Annotation information of microarray was used to match probes with corresponding genes. The average value was calculated out for those genes corresponding to more than one probes, while probes matched with more than one gene were eliminated.

Screening for differentially expressed genes

The “DEseq2” R package was used to screen DEGs between ChRCC samples and paired normal tissues in the expression data. The DEGs threshold was set at a |log2FoldChange| > 0.585 and adj.P.value < 0.05. After DEGs were screen out, flashClust tool package in R language was used to perform cluster analysis of ChRCC samples.

Co-expression network construction and module analysis

Firstly, expression data profile of DEGs was tested by to check if they were the good samples and good genes. Afterwards, power value was screened out by WGCNA [14] algorithm which is implemented in the R software package ( Scale independence and average connectivity degree of modules with different power value was tested by gradient method (the power value ranging from 1 to 20). The appropriate power value was determined when the degree of independence was above 0.85 and average connectivity degree is relatively higher [8]. Once the power value was determined, the scale-free gene co-expression networks were constructed by WGCNA algorithm. Besides, the corresponding genes information of each module was extracted out. The minimum number of genes in each module was set as 50 for the high reliability of the results.

Interaction analysis of co-expression modules of ChRCC

After co-expression modules were identified by WGCNA algorithm. Heatmap was painted to describe the strength of the interactions between different modules by heatmap function embedded in R language.

Construct module-trait relationships of ChRCC and key module identification

The correlation between module eigengenes (MEs) and phenotype (clinic traits) was used to evaluated module-trait associations. MEs were considered as the major component in the principal component analysis for each gene module. We calculated the correlation between MEs and clinical trait to identify the relevant module. Gene significance (GS), which was defined as the log10 transformation of the P value (GS = lgP) in the linear regression between gene expression and clinical information, was calculated to evaluate correlation strength. Modules with the highest correlation coefficients among all modules were usually considered as the key module and selected for further analysis.

Functional enrichment analysis of the key module

The information of genes in key modules was upload to Enrichr online database to perform GO and KEGG pathway analysis [15, 16]. Analysis results were extracted out under the condition of P < 0.05 after correction [17]. The top 10 GO terms were visualized if there were more than ten terms, so as to KEGG pathways.

Selection of candidate hub genes

Hub genes, with highly intramodular connectivity, have been shown to be functionally significant. In this study, candidate hub genes were defined by module connectivity, measured by module membership (MM) > 0.8 and clinical trait relationship, measured by significance of the Pearson’s correlation (GS. Pathologic stage > 0.2) [18]. Furthermore, we uploaded all genes in the hub module to the STRING database, choosing confidence > 0.4 to construct protein–protein interaction (PPI). After setting degree cut-off = 5, node score cut-off = 0.2, k-core = 2, and max. depth = 100, the most significant sub-module was selected by using plug-in MCODE [9]. Genes both in co-expression network and MCODE sub-module were regarded as candidate hub genes for further analysis.

Hub genes identification and validation

Another independent dataset GSE15641 containing six ChRCC samples and 23 normal tissues was extracted and performed on GEO2R analysis ( |log2FoldChange| > 1 and P.Value < 0.05 as criterion, DEGs were identified. Candidate hub genes were mapped to DEGs in GSE15641 to find the real hub genes. For further validation, real hub genes were upload to online TCGA database to perform survival analysis ( ChRCC patients were divided into two groups according to median expression of each hub gene (high vs. low) and Kaplan–Meier survival analysis was conducted.

Statistical analyses

Two-tailed Student’s t-test was used for significance of differences between groups. Statistical analyses were performed with GraphPad Prism 6.02. Statistical significance was set at probability values of P < 0.05.


DEGs screening

The gene expression profile of 66 ChRCC samples were download from TCGA database. A total of 2215 DEGs (1741 up-regulated and 384 down-regulated) between ChRCC samples and normal tissues were screened out under the threshold of |log2FoldChange| > 0.585 and adj.P.Value < 0.05. These 2215 DEGs were then selected for subsequent analysis.

Construction of co-expression modules of ChRCC

All DEGs were included for constructing co-expression modules by WGCNA algorithm. The “flashClust tool” package was used to perform the cluster and trait analysis and results were shown in (Fig. 1). First of all, the appropriate power value was screened out (Additional file 1: Fig S1). When the power value was equal to 5, scale free networks were constructed with independence degree up to 0.85 and the relatively higher average connectivity (mean connectivity = 6.04). As a result, eight distinct gene co-expression modules were identified by the appropriate power value (5) in ChRCC after excluding the grey module. Constructed modules painted with different colors and cluster trees of DEGs were shown in (Fig. 2). Gene numbers and module names were shown in Table 1. Interactions between eight co-expression modules were subsequently analyzed. The heatmap demonstrates the Topological Overlap Matrix (TOM) among all genes in the analysis. Light color represents low overlap and darker red color represents higher overlap. As a result, each module showed independent validation to each other (Fig. 3).

Fig. 1
figure 1

Clustering dendrogram of 66 tumor samples and heatmaps of clinical traits. The clustering was based on the expression data of DEGs between ChRCC samples and non-tumor samples. The color intensity in heatmaps was proportional to older age, higher pathological stage. In neoplasm and survival status, white means patient neoplasm free or live, while red means patient with tumor or dead

Fig. 2
figure 2

Clustering dendrogram of DEGs, with dissimilarity based on topological overlap, together with assigned module colors

Table 1 The number of genes in eight constructed modules
Fig. 3
figure 3

Visualizing the gene network using a heatmap plot. The heatmap depicts the Topological Overlap Matrix (TOM) among all genes in the analysis. Light color represents low overlap and progressively darker red color represents higher overlap. The gene dendrogram and module assignment are also shown along the left side and the top

Identification of key modules corresponding to clinic traits

The correlations between module eigengene and clinic traits were shown in (Fig. 4). The brown module was selected as key module for taking up top three highest correlations with progression and survival of ChRCC (R2 = 0.48 and P = 5e−05 with pathologic stage, R2 = 0.41 and P = 6e−04 with neoplasm status, R2 = 0.52 and P = 8e−06 with living status). Afterwards, we plotted scatter plots of GS vs. MM (module membership) in the brown modules with clinic traits respectively (Additional file 2: Fig S2).

Fig. 4
figure 4

Module-trait heatmap. Each row corresponds to a module eigengene, column to a trait. Each cell contains the corresponding correlation and P-value

Functional enrichment analysis of genes in the key module

To obtain further insight into the function of genes in the hub module, all genes in the brown module were uploaded to Enrichr online database to conduct GO and KEGG pathway analysis. According to P-value of each term, top 10 biological process and KEGG pathways were extracted out (Additional file 3: Table S1 and Additional file 4: Table S2) and visualized. GO analysis results showed genes in brown module significantly enriched in mitotic cell cycle transition, mitotic spindle assembly, mitotic spindle organization, regulation of cell cycle process, etc. KEGG pathway analysis revealed cell cycle, oocyte meiosis, progesterone-mediated oocyte maturation pathways were significantly enriched (Fig. 5). All these results implied that dysfunctional mitotic cell cycle may contribute to tumorigeneses of ChRCC.

Fig. 5
figure 5

Functional analysis. a Gene ontology analysis of genes in brown module. b KEGG pathway enrichment analysis of genes in brown module. The x-axis shows the − log10 (P-value) of each term and the y-axis shows the GO and KEGG pathway terms

Selection of candidate hub genes

Under the condition of MM > 0.8 and GS. Pathologic stage > 0.2, 39 genes in brown module were taken out (Fig. 6a). PPI network were constructed under the cutoff of confidence > 0.4, 32 genes were filtered out after MCODE process (Fig. 6b). 29 candidate genes identified both in co-expression network and MCODE sub-module were show in Table 2.

Fig. 6
figure 6

Candidate hub gene identification. a Scatter plot of module eigengenes in brown module. The vertical line represents cutoff of module membership = 0.8, the horizontal line represents cutoff of gene significances in pathologic stage = 0.2. Genes on upper right were selected out. b PPI network of genes in MCODE sub-module

Table 2 Candidate hub genes identified by co-expression and MCODE methods

Identification and validation of hub genes

1794 DEGs (884 up-regulated, 910 down-regulate) were identified in GSE15641. All candidate genes were mapped to DEGs in GSE15641, four genes (SKA1, ERCC6L, POLQ, GTSE1) were found to be differentially expressed both in TCGA and GSE15641 datasets. While expressions of POLQ in two datasets were opposite (up-regulated in TCGA, down-regulated in GSE15641). Thus, SKA1, ERCC6L and GTSE1 were regarded as the real hub genes in ChRCC by WGCNA analysis. Expressions of real hub genes between ChRCC samples and normal tissues were shown in (Fig. 7). Survival analysis of these genes showed over-expressions of SKA1 and ERCC6L are significantly related to shorter overall survival time (Fig. 8).

Fig. 7
figure 7

Expression levels of real hub genes. a Expression levels of real hub genes between tumor samples and normal tissues in TCGA database. b Expression levels of real hub genes between tumor samples and normal tissues in GSE15641 dataset

Fig. 8
figure 8

Validation of real hub genes by survival analysis between high expression group (red line) and low expression group (blue line). a SKA1, b ERCC6L, c GTSE1


ChRCC is an uncommon, but not rare, malignant disease representing ~ 5% of histologic spectrum of cancer arising from kidney [19]. However, patients with ChRCC have a relatively low risk of tumor progression and cancer-specific death [20], the clinical outcomes of target therapies between ChRCC and ccRCC were not significantly different in metastatic disease [21]. To the best of our knowledge, there are no effective treatments for patients with metastatic ChRCC or any preventions for tumor recurrence. Thus, it’s important to get a better understanding of molecular mechanism of ChRCC and identify potential biomarkers to evaluate the biological behavior of this malignancy.

WGCNA has many prominent advantages over other methods since the analysis explores associations between co-expression modules and clinic traits and the results had much higher reliability and biological significance [22]. In this study, a total of eight co-expression modules were constructed by the 2215 DEGs from the 66 human ChRCC samples by WGCNA method, which was used to detect the relationship between ChRCC transcriptome and clinic traits. We calculated the correlations between co-expression modules and pathologic stage, neoplasm status and living status. Brown module with the highest correlation with clinic traits were regarded as the key module to explore primary cause for disease progression.

For further analysis, genes in brown module with high module membership (> 0.8) and significance of correlations with pathologic stage (> 0.2) were filtered out. Meanwhile, MCODE analysis was performed on brown module to find a panel of genes with high connective degrees [9]. 29 candidate hub genes were identified by overlapping results of co-expression analysis and MCODE analysis. GO and KEGG analysis showed that dysregulation of cell cycle may be the underlying mechanism of tumorigeneses of ChRCC.

29 candidate hub genes were subsequently validated in GSE15641 dataset, four candidate genes (SKA1, ERCC6L, POLQ, GTSE1) differentially expressed both in TCGA and GSE15641 database were considered as real hub genes. After excluding POLQ for conflicting expression level in two databases, two of three hub genes (SKA1, ERCC6L) were found to be significantly correlated with shorter overall survival time under their over expressions.

SKA1 was a microtubule-binding subcomplex of the outer kinetochore which is essential for proper chromosome segregation. It played an important role on tumorigenesis in multiple malignancies [23, 24]. SKA1 over-expression led to cancerization in human prostate epithelial cells via the induction of centriole over-duplication [25]. Depletion of SKA1 inhibited cell proliferation in gastric cancer by blocking cell cycle in S phase [26]. What’s more, SKA1 was proved to be an oncogene in ccRCC which could be down-regulated by antitumor miR-10a-5p transfection [27]. ERCC6L gene, also named PLK1-interacting checkpoint helicase (PICH), was a member of the SNF2 protein family (SWI/SNF catalytic subunit SNF2). It was a mitotic target and substrate of polo-like kinases (PLKs) which regulates multiple processes in mammalian cell mitosis [28]. Downregulation of ERCC6L decreased cell viability in RCC cell lines by blocking mitogen-activated protein kinase (MAPK) signaling pathway and interactions with protein PLK1 [29]. Protein encoded by GTSE-1 modulated cell migrations in an EB1-dependent manner. Up-regulation of GTSE1 expression could be associated with increased invasive potential in breast cancer [30]. Depletion of GTSE-1 enhanced mitotic centromere-associated kinesin (MCAK) activity in mitotic cells, leading to chromosomal instability (CIN) which was presented in most solid tumors [31]. Silencing GTSE-1 expression inhibits proliferation and invasion of hepatocellular carcinoma cells [32]. Based on many studies about three hub genes, we could find their important roles in tumorigenesis and metastasis. Targeting these biomarkers and other candidate hub genes participating in cell cycles could provide therapeutic potentials for ChRCC.

Some limitations exiting in our studies should be mentioned. The most vital genes out of 29 candidate hub genes couldn’t been filtered out under restrictions of the bioinformatics methods. Multi-central and large sample studies on ChRCC expression profiles were needed for validation of our findings. However, this couldn’t be fulfilled due to limited online databases for ChRCC.


Our study used weighted gene co-expression analysis to construct a gene co-expression network, identify and validate the key module and hub genes associated with the progression and prognosis of ChRCC. Functional analysis of hub genes reveal mitotic cell cycle dysregulation could be the underlying mechanism of ChRCC. Three real hub genes including SKA1, ERCC6L and GTSE-1 were identified. SKA1 and ERCC6L were validated in association with poor prognosis of ChRCC. However, further molecular biological experiments are needed to confirm the function of these biomarkers in ChRCC.



chromophobe renal cell carcinoma


non-clear cell renal cell carcinoma


renal cell carcinoma


The Cancer Genome Atlas


weighted gene co-expression network analysis

GO analysis:

gene ontology analysis


clear cell renal cell carcinoma


overall survival


relapse free survival


gene significance


module membership


  1. Barata PC, Rini BI. Treatment of renal cell carcinoma: current status and future directions. CA Cancer J Clin. 2017;67(6):507–24.

    Article  Google Scholar 

  2. Jonasch E, Gao JJ, Rathmell WK. Renal cell carcinoma. BMJ. 2014;349:g4797.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bellmunt J, Dutcher J. Targeted therapies and the treatment of non-clear cell renal cell carcinoma. Ann Oncol. 2013;24(7):1730–40.

    CAS  Article  PubMed  Google Scholar 

  4. Lee WK, Byun SS, Kim HH, Rha KH, Hwang TK, Sung GT, Lee W, Lim JS, Jeong YB, Kwon TG. Characteristics and prognosis of chromophobe non-metastatic renal cell carcinoma: a multicenter study. Int J Urol. 2010;17(11):898–904.

    Article  PubMed  Google Scholar 

  5. Horvath S, Zhang B, Carlson M, Lu KV, Zhu S, Felciano RM, Laurance MF, Zhao W, Qi S, Chen Z, Lee Y, Scheck AC, Liau LM, Wu H, Geschwind DH, Febbo PG, Kornblum HI, Cloughesy TF, Nelson SF, Mischel PS. Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target. Proc Natl Acad Sci USA. 2006;103(46):17402–7.

    CAS  Article  PubMed  Google Scholar 

  6. Mao Q, Zhang L, Zhang Y, Dong G, Yang Y, Xia W, Chen B, Ma W, Hu J, Jiang F, Xu L. A network-based signature to predict the survival of non-smoking lung adenocarcinoma. Cancer Manag Res. 2018;10:2683–93.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Yuan L, Qian G, Chen L, Wu CL, Dan HC, Xiao Y, Wang X. Co-expression network analysis of biomarkers for adrenocortical carcinoma. Front Genet. 2018;9:328.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Chen L, Yuan L, Wang Y, Wang G, Zhu Y, Cao R, Qian G, Xie C, Liu X, Xiao Y, Wang X. Co-expression network analysis identified FCER1G in association with progression and prognosis in human clear cell renal cell carcinoma. Int J Biol Sci. 2017;13(11):1361–72.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Yuan LS, Zeng G, Chen L, Wang G, Wang XL, Cao XY, Lu MX, Liu XF, Qian GF, Xiao Y, Wang XH. Identification of key genes and pathways in human clear cell renal cell carcinoma (ccRCC) by co-expression analysis. Int J Biol Sci. 2018;14(3):266–79.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. Spiers H, Hannon E, Schalkwyk LC, Smith R, Wong CCY, O’Donovan MC, Bray NJ, Mill J. Methylomic trajectories across human fetal brain development. Genome Res. 2015;25(3):338–52.

    CAS  Article  Google Scholar 

  11. Liu Q, Jiang C, Xu J, Zhao MT, Van Bortle K, Cheng X, Wang GW, Chang HY, Wu JC, Snyder MP. Genome-wide temporal profiling of transcriptome and open chromatin of early cardiomyocyte differentiation derived from hiPSCs and hESCs. Circ Res. 2017;121(4):376–91.

    CAS  Article  Google Scholar 

  12. Giulietti M, Occhipinti G, Principato G, Piva F. Identification of candidate miRNA biomarkers for pancreatic ductal adenocarcinoma by weighted gene co-expression network analysis. Cell Oncol. 2017;40(2):181–92.

    CAS  Article  Google Scholar 

  13. Giulietti M, Righetti A, Principato G, Piva F. LncRNA co-expression network analysis reveals novel biomarkers for pancreatic cancer. Carcinogenesis. 2018.

    Article  PubMed  Google Scholar 

  14. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559.

    CAS  Article  Google Scholar 

  15. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, Gundersen GW, Ma’ayan A. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR, Ma’ayan A. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013;14:128.

    Article  Google Scholar 

  17. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. Chen L, Yuan L, Qian K, Qian G, Zhu Y, Wu CL, Dan HC, Xiao Y, Wang X. Identification of biomarkers associated with pathological stage and prognosis of clear cell renal cell carcinoma by co-expression network analysis. Front Physiol. 2018;9:399.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Reuter VE. The pathology of renal epithelial neoplasms. Semin Oncol. 2006;33(5):534–43.

    Article  Google Scholar 

  20. Volpe A, Novara G, Antonelli A, Bertini R, Billia M, Carmignani G, Cunico SC, Longo N, Martignoni G, Minervini A, Mirone V, Simonato A, Terrone C, Zattoni F, Ficarra V, LUN-PF A. Chromophobe renal cell carcinoma (RCC): oncological outcomes and prognostic factors in a large multicentre series. BJU Int. 2012;110(1):76–83.

    Article  Google Scholar 

  21. Keizman D, Sarid D, Lee JL, Sella A, Gottfried M, Hammers H, Eisenberger MA, Carducci MA, Sinibaldi V, Neiman V, Rosenbaum E, Peer A, Neumann A, Mermershtain W, Rouvinov K, Berger R, Yildiz I, Grp TO. Outcome of patients with metastatic chromophobe renal cell carcinoma treated with sunitinib. Oncologist. 2016;21(10):1212–7.

    Article  Google Scholar 

  22. Chou WC, Cheng AL, Brotto M, Chuang CY. Visual gene-network analysis reveals the cancer gene co-expression in human endometrial cancer. BMC Genom. 2014;15:300.

    CAS  Article  Google Scholar 

  23. Tian F, Xing X, Xu F, Cheng W, Zhang Z, Gao J, Ge J, Xie H. Downregulation of SKA1 Gene expression inhibits cell growth in human bladder cancer. Cancer Biother Radiopharm. 2015;30(7):271–7.

    CAS  Article  PubMed  Google Scholar 

  24. Qin X, Yuan B, Xu X, Huang H, Liu Y. Effects of short interfering RNA-mediated gene silencing of SKA1 on proliferation of hepatocellular carcinoma cells. Scand J Gastroenterol. 2013;48(11):1324–32.

    CAS  Article  PubMed  Google Scholar 

  25. Li J, Xuan JW, Khatamianfar V, Valiyeva F, Moussa M, Sadek A, Yang BB, Dong BJ, Huang YR, Gao WQ. SKA1 over-expression promotes centriole over-duplication, centrosome amplification and prostate tumourigenesis. J Pathol. 2014;234(2):178–89.

    CAS  Article  PubMed  Google Scholar 

  26. Sun W, Yao L, Jiang B, Guo L, Wang Q. Spindle and kinetochore-associated protein 1 is overexpressed in gastric cancer and modulates cell growth. Mol Cell Biochem. 2014;391(1–2):167–74.

    CAS  Article  PubMed  Google Scholar 

  27. Arai T, Okato A, Kojima S, Idichi T, Koshizuka K, Kurozumi A, Kato M, Yamazaki K, Ishida Y, Naya Y, Ichikawa T, Seki N. Regulation of spindle and kinetochore-associated protein 1 by antitumor miR-10a-5p in renal cell carcinoma. Cancer Sci. 2017;108(10):2088–101.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Baumann C, Korner R, Hofmann K, Nigg EA. PICH, a centromere-associated SNF2 family ATPase, is regulated by Plk1 and required for the spindle checkpoint. Cell. 2007;128(1):101–14.

    CAS  Article  PubMed  Google Scholar 

  29. Zhang G, Yu Z, Fu S, Lv C, Dong Q, Fu C, Kong C, Zeng Y. ERCC6L that is up-regulated in high grade of renal cell carcinoma enhances cell viability in vitro and promotes tumor growth in vivo potentially through modulating MAPK signalling pathway. Cancer Gene Ther. 2018.

    Article  PubMed  Google Scholar 

  30. Scolz M, Widlund PO, Piazza S, Bublik DR, Reber S, Peche LY, Ciani Y, Hubner N, Isokane M, Monte M, Ellenberg J, Hyman AA, Schneider C, Bird AW. GTSE1 is a microtubule plus-end tracking protein that regulates EB1-dependent cell migration. PLoS ONE. 2012;7(12):e51259.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Bendre S, Rondelet A, Hall C, Schmidt N, Lin YC, Brouhard GJ, Bird AW. GTSE1 tunes microtubule stability for chromosome alignment and segregation by inhibiting the microtubule depolymerase MCAK. J Cell Biol. 2016;215(5):631–47.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. Guo L, Zhang S, Zhang B, Chen W, Li X, Zhang W, Zhou C, Zhang J, Ren N, Ye Q. Silencing GTSE-1 expression inhibits proliferation and invasion of hepatocellular carcinoma cells. Cell Biol Toxicol. 2016;32(4):263–74.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

JZ designed the project for this study. XY analyzed datasets and wrote the initial edition of this manuscript. JW revised the manuscript and assisted with formatting and language editing. All authors read and approved the final manuscript.


Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used during the current study are available from TCGA database ( and GEO database (

Consent for publication

Not applicable.

Ethics approval and consent to participate

For this type of study, ethics approval is not applicable.


No funding was received.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jin Zhang.

Additional files

Additional file 1: Fig S1.

Analysis of network topology for various soft thresholding powers. The left panel shows the scale-free fit index (y-axis) as a function of the soft-thresholding power (x-axis). The right panel displays the mean connectivity (degree, y-axis) as a function of the soft-thresholding power (x-axis).

Additional file 2: Fig S2.

(a) Scatterplot of Gene Significance (GS) for pathologic stage vs. Module Membership (MM) in the brown module. (b) scatterplot of Gene Significance (GS) for neoplasm status vs. Module Membership (MM) in the brown module. (c) scatterplot of Gene Significance (GS) for survival status vs. Module Membership (MM) in the brown module.

Additional file 3: Table S1.

GO enrichment analysis in brown module.

Additional file 4: Table S2.

KEGG pathways enrichment in brown module.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yin, X., Wang, J. & Zhang, J. Identification of biomarkers of chromophobe renal cell carcinoma by weighted gene co-expression network analysis. Cancer Cell Int 18, 206 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Chromophobe renal cell carcinoma
  • Weighted gene co-expression network analysis (WGCNA)
  • Biomarker