Skip to main content
Download PDF
  • Primary research
  • Open access
  • Published:

Identification of biological targets of therapeutic intervention for clear cell renal cell carcinoma based on bioinformatics approach

Cancer Cell International volume 16, Article number: 16 (2016) Cite this article

Abstract

Background

We aimed to discover the potential microRNA (miRNA) targets and to explore the underlying molecular mechanisms of clear cell renal cell carcinoma (ccRCC).

Methods

Microarray data of GSE16441 was downloaded from Gene Expression Omnibus database. Differentially expressed genes (DEGs) and differentially expressed miRNAs between ccRCC tumors and matched non-tumor samples were analyzed. Target genes of differentially expressed miRNAs were screened. Besides, functional enrichment analysis of DEGs was performed, followed by protein–protein interaction (PPI) network construction and sub-module analysis. Finally, the integrated miRNA-DEGs network was constructed.

Results

A total of 1758 up- and 2465 down-regulated DEGs were identified. Moreover, 15 up- and 12 down-regulated differentially expressed miRNAs were screened. The up-regulated DEGs were significantly enriched in pathways such as cell adhesion molecules and focal adhesion. Besides, the down-regulated DEGs were enriched in oxidative phosphorylation, and citrate cycle (TCA cycle). Moreover, eight sub-modules of PPI network were obtained. Totally, eight down-regulated miRNAs were identified to significantly regulate the DEGs and miRNA-200c that could regulate collagen, type V, alpha 2 (COL5A2) as well as COL5A3 was found to be the most significant. Additionally, 10 up-regulated miRNAs were identified to be significantly associated with the DEGs. Thereinto, miRNA-15a that could regulate ATPase, H+ transporting, lysosomal 21 kDa, V0 subunit b (ATP6V0B) and miRNA-155 were found to be the most significant.

Conclusions

miRNA-200c that could regulate COL5A2 and COL5A3, miRNA-15a that could regulate ATP6V0B and miRNA-155 may play key roles in ccRCC progression. These miRNAs may be potential targets for ccRCC treatment.

Background

Renal cell carcinoma (RCC) represents the most lethal genitourinary malignancy and about 90 % of all kidney tumors are RCC [1]. Clear cell RCC (ccRCC) is the most common subtype of RCC, accounting for about 75 % of cases [2]. Because of being generally asymptomatic, one-third of ccRCC are already metastatic at initial diagnosis, leading to a high mortality rate up to 95 % [3]. However, the benefit of surgical, chemotherapeutic and radiological approaches for ccRCC is limited [3, 4]. Thus, it is of great importance to gain a better understanding of the pathogenesis of ccRCC, which may facilitate the development of effective biomarkers and novel, targeted therapeutic strategies.

Major efforts have been carried out to explore the underlying etiology and molecular mechanisms of ccRCC. Cigarette smoking, obesity, hypertension, and acquired cystic kidney disease were shown to be associated with renal cancer [5]. There is ample evidence indicating that the most frequent event in ccRCC is inactivation or mutations in the von-Hippel Lindau (VHL) tumor suppressor gene, which induces the expression of the hypoxia inducible factors (HIF)1α and HIF2α [6, 7]. Two of the other most commonly mutated genes in ccRCC are PBRM1 (about 40–50 %) and BAP1 (10–15 %) [8]. In addition, extracellular signal-regulated kinase 5 (ERK5) is demonstrated to mediate angiogenesis, proliferation, and tumor aggressiveness in ccRCC [9]. By contrast, microRNAs (miRNAs), a group of small non-coding RNAs (21–25 nucleotides) that regulate expression of target genes at the post-transcriptional level, are crucial for gene-regulatory networks which play an important role in different biological processes, including metabolic disorder, tissue injury, and oxidative stress [10, 11]. Aberrant alterations in miRNA expression have been shown to be related to human malignancies [12, 13]. Numerous miRNAs have also been identified in ccRCC in which some were demonstrated to function as oncogenes, while others were identified to play tumor suppressor roles [14, 15]. For instance, Mathew et al. have demonstrated that the suppression of miR-30c-2-3p as well as miR-30a-3p increases HIF2α levels in ccRCC, promoting angiogenesis, cellular proliferation, and tumor growth [16]. In addition, Liu et al. have suggested miR-23b* was up-regulated in renal cancer as an important regulator of proline oxidase [17]. Additionally, a number of array platform-based studies recently indicate that a considerable number of miRNAs are dysregulated in ccRCC [18]. Although the results of previous studies may not be consistent, all the data indicate that dysregulated miRNAs may play pivotal roles in the pathogenesis of ccRCC. However, the present knowledge is still poorly understood.

In the current study, we employed the bioinformatics methods to identify the differentially expressed genes (DEGs) and differentially expressed miRNAs between ccRCC tumors and corresponding non-tumor samples. Target genes of differentially expressed miRNAs were screened. Then, functional enrichment analysis of DEGs was performed, followed by protein–protein interaction (PPI) network construction and sub-module analysis. Finally, an integrated miRNA-DEG network was constructed. We intended to give a systematic perspective on understanding the molecular mechanism and to investigate more potential therapeutic targets for ccRCC.

Methods

Data collection

The expression profile of GSE16441, deposited by Liu et al. [19] was downloaded from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database [20] (http://www.ncbi.nlm.nih.gov/geo/). The GSE16441 dataset is composed of gene expression microarray data and miRNA microarray data. There were 34 samples (GSM413237-GSM413270) including 17 ccRCC tumors and 17 adjacent normal kidney tissue samples in the gene expression microarray data. The platform is Agilent-014850 Whole Human Genome Microarray 4x44 K G4112F (Probe Name version). Moreover, 34 samples from 17 ccRCC tumors and 17 corresponding non-tumor samples were obtained using a platform of Agilent Human miRNA Microarray Rel12.0 for miRNA expression analysis. In this study, all the samples were used for the follow-up analysis.

Data preprocessing and identification of DEGs and differentially expressed miRNAs

Data in CEL format were downloaded and converted into expression values, followed by background correction using Normexp algorithm [21]. Then, data were normalized by the Loess method [22]. Subsequently, the linear models for microarray data (Limma) package [23] from Bioconductor was applied to screen DEGs between RCC tumors and 17 corresponding non-tumor samples. The associated P values were adjusted to false discovery rates (FDR) based on the Benjamini and Hochberg method [24]. FDR <0.05 and |log FC| >1 were selected as the cutoff criterion for DEGs screening.

Similarly, with respect to the screening of differentially expressed miRNAs, Normexp method was used to perform the background correction. Then, data were normalized by the Quantile method [25]. Subsequently, Limma package from Bioconductor was applied to identify differentially expressed miRNAs between ccRCC tumors and 17 corresponding non-tumor samples. FDR <0.05 and |log FC|> 1 were defined as the cutoff criterion for the screening of differentially expressed miRNAs.

The analysis of target genes of differentially expressed miRNAs

The genes targeted by the miRNAs were screened based on seven databases (miRanda, MirTarget2, PicTar, PITA, TargetScan, miRecords as well as MirWalk). The combination of the seven databases was the criteria for target gene databases. Then, these target genes were further processed to investigate the more potential miRNA-DEG pairs by using the Lasso regression method, which was proposed to identify miRNA–mRNA targeting relationships with considerable advantages in both sensitivity and specificity [26]. The criteria was set as select.score >90 and FDR <0.05.

Gene ontology (GO) and pathway enrichment analysis of DEGs

GO analysis has been used as functional annotation of large-scale genes frequently [27]. Database for Annotation, Visualization and Integrated Discovery (DAVID) [28] is a tool providing a comprehensive set of functional annotation. In this study, DAVID was utilized for GO enrichment analysis to identify the significantly enriched biological process (BP) terms of DEGs, with P value <0.05.

Kyoto Encyclopedia of Genes and Genomes (KEGG) database includes all kinds of biochemistry pathways [29]. In this study, DAVID was also applied for KEGG pathway enrichment analysis of DEGs. P < 0.05 was chosen as the cutoff value.

PPI network construction and significant sub-module extraction

The Search Tool for the Retrieval of Interacting Genes (STRING) database offers both experimental and predicted interaction information [30]. In this study, STRING database was applied to identify functional interactions between DEGs and a PPI network was then constructed, which was visualized with Cytoscape software [31]. The combined score >0.9 was selected as the threshold.

In order to further analyze the more significant interactions from the PPI network, functional sub-modules were extracted using ClusterOne [32].

Construction of integrated miRNA-DEG network

Based on the results of the identified potential target genes of differentially expressed miRNAs from databases and results of the significant sub-modules of PPI network, the integrated miRNA-DEG network was constructed.

Results

Identification of DEGs and differentially expressed miRNAs

Totally, 1758 up-regulated DEGs (corresponding to 1804 up-regulated transcripts), and 2465 down-regulated DEGs (corresponding to 2566 down-regulated transcripts) were screened out in ccRCC samples compared with matched non-tumor samples. Moreover, a total of 15 up-regulated and 12 down-regulated differentially expressed miRNAs were identified in ccRCC samples compared with matched control samples.

The analysis of target genes of miRNAs

According to the criteria of select.score >90 and FDR <0.05, target genes of the up- and down-regulated miRNAs were obtained. Number of target genes for the up- and down-regulated miRNAs is listed in Tables 1, 2 respectively.

Table 1 The number of target genes of up-regulated miRNAs
Full size table
Table 2 The number of target genes of down-regulated miRNAs
Full size table

GO and pathway enrichment analysis of DEGs

The top 5 GO BP terms of up-regulated DEGs and down-regulated DEGs were shown in Table 3. From the results, we found that the up-regulated DEGs were significantly enriched in immune response (P = 9.10E−18), cell activation (P = 3.82E−13), positive regulation of immune system process (P = 2.70E−12), regulation of cell activation (P = 5.70E−12) and leukocyte activation (P = 1.86E−11). While the down-regulated DEGs mainly participated in oxidation reduction (P = 2.56E−16), generation of precursor metabolites and energy (P = 2.93E−15), ion transport (P = 4.19E−13), hydrogen transport (P = 1.38E−12) and proton transport (P = 3.91E−12).

Table 3 The top five enriched gene ontology (GO) terms functional enrichment of the differentially expressed genes
Full size table

The pathways that were significantly enriched by the up- and down-regulated genes were shown in Tables 4, 5, respectively. The results showed that the up-regulated DEGs were significantly associated with several pathways such as primary immunodeficiency (P = 1.63E−06), cell adhesion molecules (CAMs) (P = 6.66E−06), natural killer cell mediated cytotoxicity (P = 2.13E−04), cytokine–cytokine receptor interaction (P = 3.52E−04). By contrast, the down-regulated DEGs were significantly enriched in oxidative phosphorylation (P = 1.83E−13), valine, leucine and isoleucine degradation (P = 1.31E−07) and citrate cycle (TCA cycle) (P = 4.63E−07).

Table 4 Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of up-regulated differentially expressed genes
Full size table
Table 5 Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of down-regulated differentially expressed genes
Full size table

PPI network construction and sub-modules mining

Totally, four most significant sub-modules of up-regulated DEGs (sub-module a, b, c, and d) were extracted from PPI network as shown in Fig. 1. Sub-module a including 28 nodes and 162 edges was enriched in cell cycle. Sub-module b containing 9 nodes and 36 edges was related with ECM-receptor interaction. Sub-module c with 16 nodes and 54 edges was enriched in chemokine signaling pathway. Besides, 19 nodes combined with 89 edges made up sub-module d which was enriched in T cell receptor signaling pathway.

Fig. 1
figure 1

The four significant sub-modules of up-regulated genes (sub-module a, b, c, d) in PPI network. The red oval nodes denote the up-regulated genes

Full size image

Additionally, a total of four significant sub-modules of down-regulated genes (sub-module a, b, c, d) were achieved as shown in Fig. 2. Sub-module a including 14 nodes and 80 edges was enriched in oxidative phosphorylation. Sub-module b with 13 nodes and 78 edges was involved in oxidative phosphorylation. Sub-module c with 15 nodes and 56 edges was enriched in neuroactive ligand-receptor interaction. Moreover, 15 nodes combined with 56 edges made up sub-module d which was related with neuroactive ligand-receptor interaction.

Fig. 2
figure 2

The four significant sub-modules of down-regulated genes (sub-module a, b, c, d) in PPI network. The green oval nodes denote the down-regulated genes

Full size image

Construction of integrated miRNA-DEG network

The miRNA-DEG network was constructed by integrating significant DEGs in sub-modules and potential miRNAs-DEG pairs. Up-regulated DEGs in the four identified sub-modules and their related miRNAs (down-regulated) were shown in Fig. 3. Moreover, down-regulated DEGs in the four identified sub-modules and their related miRNAs (up-regulated) were shown in Fig. 4.

Fig. 3
figure 3

The network of up-regulated DEGs in four sub-modules (a, b, c, d) and their related miRNAs. The red oval nodes are the up-regulated DEGs, and green triangle nodes are the down-regulated miRNAs

Full size image
Fig. 4
figure 4

The network of down-regulated DEGs in four sub-modules (a, b, c, d) and their related miRNAs. The green oval nodes are the down-regulated DEGs, and red triangle nodes are the up-regulated miRNAs

Full size image

Discussion

ccRCC is the most common histological subtype of RCC that occurs in adults and associated with worse prognosis [33]. In this study, we applied bioinformatics method to predict the potential miRNA targets for the treatment of ccRCC progression. Our results suggested that 1758 up- and 2465 down-regulated DEGs were screened out in ccRCC samples. Moreover, a total of 15 up- and 12 down-regulated differentially expressed miRNAs were identified. The up-regulated DEGs were enriched in significant pathways such as CAMs and focal adhesion. Besides, the down-regulated DEGs were significantly associated with oxidative phosphorylation, and TCA cycle. Several significant differentially expressed miRNAs were identified and miRNA-200 family was found to be the most significant.

miRNA-200 family includes miRNA-200a, miRNA-200b, miRNA-200c, miRNA-429, and miRNA-141 [34]. In the present work, miR-200a, miR-200b, miR-200c and miR-429 were identified from the network of DEGs and their related miRNAs. Additionally, these miRNAs were down-regulated in ccRCC. As previously reported, the members of the miRNA-200 family (especially miR-200c and miR-141) play an outstanding role as metastasis suppressor genes via inhibiting the expression of zinc finger E-box binding homeobox 1 (ZEB1) [35, 36]. Under-expression of miRNA-200 family members is correlated with renal cancer [37]. Moreover, the elevation of collagens and fibronectin in obstructed kidneys can be repressed by the injection of miR-200b [38]. Previous reports exhibited that the increased level of type V collagen has been detected in human breast cancer and in mouse skin tumors [39, 40]. In the current study, we found that under-expression of miR-200c targeted and up-regulated the level of collagen, type V, alpha 2 (COL5A2) and COL5A3. Besides, we found that the up-regulated DEGs were significantly associated with CAMs and focal adhesion (Table 4). Amounting evidence has demonstrated that collagens and fibronectin contributes to cell adhesion dynamics and cell migration which are significantly concerned with tumor metastasis [41, 42]. Taken together, we infer that miRNA-200 family may regulate several genes such as COL5A2 and play a critical role in the progression in ccRCC through participating in cell adhesion and migration. The members of the miRNA-200 family, especially miRNA-200c, may constitute novel therapeutic targets in ccRCC and further experimental verifications are needed.

Out of the highly up-regulated miRNAs, miRNA-15a was a member of the miRNA-15 precursor family including miRNA-15a, miRNA-15b, miRNA-16-1, miRNA-16-2, miRNA-195 and miRNA-497. Metabolic activities for energy in normal cells depend on mitochondrial oxidative phosphorylation (OXPHOS) primarily, but OXPHOS capacity is decreased in various cancer cells [43]. The process of OXPHOS needs oxygen. Oxygen deficiency leads to the inhibition of OXPHOS, mainly mediated by the HIF-1 [44]. HIF is a master regulator of RCC metabolism [45]. Moreover, the miRNA-15a can regulate oxygen consumption and adenosine triphosphate (ATP) production via targeting uncoupling-protein 2 [46]. The low ATP synthase often observed in ccRCC [47]. In the present study, the up-regulated miRNA-15a was identified to regulate the ATPase, H+ transporting, lysosomal 21 kDa, V0 subunit b (ATP6V0B) gene which was significantly enriched in the pathway OXPHOS. Moreover, pathway enrichment analysis of down-regulated genes showed these genes were significantly linked with oxidative phosphorylation (Table 5). Hence, we infer that miRNA-15a may play a crucial role in the pathogenesis of ccRCC via being involved in oxidative phosphorylation. miRNA-15a may provide a novel therapeutic target in ccRCC and further experiments are needed to verify this finding.

Another up-regulated miRNA in human cancers is miRNA-155, which is described as an oncogene [48]. Previous studies demonstrated that miRNA-155 was up-regulated in various kinds of human malignancy, including breast cancer, non-small cell lung cancer and lymphomas [49–51]. Recent studies have also exhibited the up-regulation of miRNA-155 levels in renal cancers [52, 53]. In accordance with the previous studies, our results found that miRNA-155 was up-regulated in RCC samples. Moreover, the inhibition of miRNA-155 expression induced apoptosis, suppressed proliferation and migration in renal cancer cells [54]. Collectively, our results further confirmed that therapy targeting miRNA-155 inhibition may be an effective approach for ccRCC treatment.

However, our study has several limitations. For instance, analysis of the other data in similar topic with larger samples may be beneficial to cross-check of the results in our study. Besides, there was a lack of experimental verifications. We intend to perform further experimental verifications in our future studies using different methods as the other study performed [14, 19], such as immunohistochemistry, and quantitative RT-PCR.

Conclusions

In summary, the identified DEGs and their related miRNAs, especially miRNA-200c and its target genes COL5A2 as well as COL5A3, miRNA-15a and its target genes ATP6V0B and miRNA-155 may play key roles in the progression of ccRCC, and these may be useful biomarkers for predicting tumor metastasis and therapeutic targets for the treatment of ccRCC. However, further experiments are needed to validate the effects and mechanisms of miRNA-200c, miRNA-15a and miRNA-155 in ccRCC.

Abbreviations

ATP:

adenosine triphosphate

CAMs:

cell adhesion molecules

ccRCC:

clear cell renal cell carcinoma

DEGs:

differentially expressed genes

ERK5:

extracellular signal-regulated kinase 5

GEO:

gene expression omnibus

NCBI:

National Center of Biotechnology Information

OXPHOS:

oxidative phosphorylation

PBRM1:

polybromo 1

PPI:

protein–protein interaction

VHL:

von-Hippel Lindau

ZEB1:

zinc finger E-box binding homeobox 1

References

  1. Figlin RA. Personalized immunotherapy (AGS-003) when combined with sunitinib for the treatment of metastatic renal cell carcinoma. Expert Opin Biol Ther. 2015;15:1241–8.

    Article  CAS  PubMed  Google Scholar 

  2. Srinivasan R, Ricketts CJ, Sourbier C, Linehan WM. New strategies in renal cell carcinoma: targeting the genetic and metabolic basis of disease. Clin Cancer Res. 2015;21:10–7.

    Article  CAS  PubMed  Google Scholar 

  3. Escudier B, Kataja V. Renal cell carcinoma: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2010;21:v137–9.

    Article  PubMed  Google Scholar 

  4. Coppin C, Kollmannsberger C, Le L, Porzsolt F, Wilt TJ. Targeted therapy for advanced renal cell cancer (RCC): a Cochrane systematic review of published randomised trials. BJU Int. 2011;108:1556–63.

    Article  CAS  PubMed  Google Scholar 

  5. Cohen HT, Mcgovern FJ. Renal-cell carcinoma. N Engl J Med. 2005;353:2477–90.

    Article  CAS  PubMed  Google Scholar 

  6. Beroukhim R, Brunet J-P, Di Napoli A, Mertz KD, Seeley A, Pires MM, Linhart D, Worrell RA, Moch H, Rubin MA. Patterns of gene expression and copy-number alterations in von-hippel lindau disease-associated and sporadic clear cell carcinoma of the kidney. Cancer Res. 2009;69:4674–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Kaelin WG. Treatment of kidney cancer. Cancer. 2009;115:2262–72.

    Article  CAS  PubMed  Google Scholar 

  8. Joseph RW, Kapur P, Serie DJ, Parasramka M, Ho TH, Cheville JC, Frenkel E, Parker AS, Brugarolas J. Clear cell renal cell carcinoma subtypes identified by BAP1 and PBRM1 expression. J Urol. 2016;195:180–7.

    Article  CAS  PubMed  Google Scholar 

  9. Zubac DP, Bostad L, Kihl B, Seidal T, Wentzel-Larsen T, Haukaas SA. The expression of thrombospondin-1 and p53 in clear cell renal cell carcinoma: its relationship to angiogenesis, cell proliferation and cancer specific survival. J Urol. 2009;182:2144–9.

    Article  PubMed  Google Scholar 

  10. Rottiers V, Näär AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13:239–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Karthikeyan Chandrasekaran DSK. Sepramaniam S, Armugam A, Wintour EM, Bertram JF, Jeyaseelan K: role of microRNAs in kidney homeostasis and disease. Kidney Int. 2012;81:617–27.

    Article  Google Scholar 

  12. Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, Levy R. Prediction of survival in diffuse large-B-cell lymphoma based on the expression of six genes. N Engl J Med. 2004;350:1828–37.

    Article  CAS  PubMed  Google Scholar 

  13. Allison DB, Cui X, Page GP, Sabripour M. Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006;7:55–65.

    Article  CAS  PubMed  Google Scholar 

  14. Xiang W, He J, Huang C, Chen L, Tao D, Wu X, Wang M, Luo G, Xiao X, Zeng F. miR-106b-5p targets tumor suppressor gene SETD2 to inactive its function in clear cell renal cell carcinoma. Oncotarget. 2015;6:4066.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Yi Z, Fu Y, Zhao S, Zhang X, Ma C. Differential expression of miRNA patterns in renal cell carcinoma and nontumorous tissues. J Cancer Res Clin Oncol. 2010;136:855–62.

    Article  CAS  PubMed  Google Scholar 

  16. Mathew LK, Lee SS, Skuli N, Rao S, Keith B, Nathanson KL, Lal P, Simon MC. Restricted expression of miR-30c-2-3p and miR-30a-3p in clear cell renal cell carcinomas enhances HIF2α activity. Cancer Discov. 2014;4:53–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Liu W, Zabirnyk O, Wang H, Shiao Y, Nickerson M, Khalil S, Anderson L, Perantoni A, Phang J. miR-23b* targets proline oxidase, a novel tumor suppressor protein in renal cancer. Oncogene. 2010;29:4914–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Osanto S, Qin Y, Buermans HP, Berkers J, Lerut E, Goeman JJ, Van Poppel H. Genome-wide microRNA expression analysis of clear cell renal cell carcinoma by next generation deep sequencing. PLoS One. 2012;7:e38298.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Liu H, Brannon AR, Reddy AR, Alexe G, Seiler MW, Arreola A, Oza JH, Yao M, Juan D, Liou LS. Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma. BMC Syst Biol. 2010;4:51.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  PubMed  Google Scholar 

  22. Smyth GK, Speed T. Normalization of cDNA microarray data. Methods. 2003;31:265–73.

    Article  CAS  PubMed  Google Scholar 

  23. Smyth GK. Limma: linear models for microarray data, in bioinformatics and computational biology solutions using R and Bioconductor. Springer. 2005. p. 397–420.

  24. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125:279–84.

    Article  CAS  PubMed  Google Scholar 

  25. Bolstad BM, Irizarry RA, Strand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–93.

    Article  CAS  PubMed  Google Scholar 

  26. Lu Y, Zhou Y, Qu W, Deng M, Zhang C. A Lasso regression model for the construction of microRNA-target regulatory networks. Bioinformatics. 2011;27:2406–13.

    Article  CAS  PubMed  Google Scholar 

  27. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4:P3.

    Article  PubMed  Google Scholar 

  29. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, Von Mering C, Jensen LJ. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41:D808–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kohl M, Wiese S, Warscheid B, Cytoscape: software for visualization and analysis of biological networks, in Data Mining in Proteomics. Springer. 2011. p. 291–303.

  32. Maraziotis IA, Dimitrakopoulou K, Bezerianos A. An in silico method for detecting overlapping functional modules from composite biological networks. BMC Syst Biol. 2008;2:93.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Grignon DJ, Che M. Clear cell renal cell carcinoma. Clin Lab Med. 2005;25:305–16.

    Article  PubMed  Google Scholar 

  34. Howe EN, Cochrane DR, Richer JK. The miR-200 and miR-221/222 microRNA families: opposing effects on epithelial identity. J Mammary Gland Biol Neoplasia. 2012;17:65–77.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Hurteau GJ, Carlson JA, Spivack SD, Brock GJ. Overexpression of the microRNA hsa-miR-200c leads to reduced expression of transcription factor 8 and increased expression of E-cadherin. Cancer Res. 2007;67:7972–6.

    Article  CAS  PubMed  Google Scholar 

  36. Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, Vadas MA, Khew-Goodall Y, Goodall GJ. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol. 2008;10:593–601.

    Article  CAS  PubMed  Google Scholar 

  37. Nakada C, Matsuura K, Tsukamoto Y, Tanigawa M, Yoshimoto T, Narimatsu T, Nguyen L, Hijiya N, Uchida T, Sato F. Genome-wide microRNA expression profiling in renal cell carcinoma: significant down-regulation of miR-141 and miR-200c. J Pathol. 2008;216:418–27.

    Article  CAS  PubMed  Google Scholar 

  38. Oba S, Kumano S, Suzuki E, Nishimatsu H, Takahashi M, Takamori H, Kasuya M, Ogawa Y, Sato K, Kimura K. miR-200b precursor can ameliorate renal tubulointerstitial fibrosis. PLoS One. 2010;5:e13614.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Marian B, Danner M. Skin tumor promotion is associated with increased type V collagen content in the dermis. Carcinogenesis. 1987;8:151–4.

    Article  CAS  PubMed  Google Scholar 

  40. Barsky SH, Rao C, Grotendorst GR, Liotta LA. Increased content of type V collagen in desmoplasia of human breast carcinoma. Am J Pathol. 1982;108:276.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Grimaldi C, Pisanti S, Laezza C, Malfitano AM, Santoro A, Vitale M, Caruso MG, Notarnicola M, Iacuzzo I, Portella G. Anandamide inhibits adhesion and migration of breast cancer cells. Exp Cell Res. 2006;312:363–73.

    Article  CAS  PubMed  Google Scholar 

  42. Von Sengbusch A, Gassmann P, Fisch KM, Enns A, Nicolson GL, Haier J. Focal adhesion kinase regulates metastatic adhesion of carcinoma cells within liver sinusoids. Am J Pathol. 2005;166:585–96.

    Article  Google Scholar 

  43. Zheng J. Energy metabolism of cancer: glycolysis versus oxidative phosphorylation (Review). Oncol Lett. 2012;4:1151–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Semenza GL. Oxygen-dependent regulation of mitochondrial respiration by hypoxia-inducible factor 1. Biochem J. 2007;405:1–9.

    Article  CAS  PubMed  Google Scholar 

  45. Sudarshan S, Karam JA, Brugarolas J, Thompson RH, Uzzo R, Rini B, Margulis V, Patard JJ, Escudier B, Linehan WM. Metabolism of kidney cancer: from the lab to clinical practice. Eur Urol. 2013;63:244–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Li P, Jiao J, Gao G, Prabhakar BS. Control of mitochondrial activity by miRNAs. J Cell Biochem. 2012;113:1104–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Simonnet H, Alazard N, Pfeiffer K, Gallou C, Beroud C, Demont J, Bouvier R, Schagger H, Godinot C. Low mitochondrial respiratory chain content correlates with tumor aggressiveness in renal cell carcinoma. Carcinogenesis. 2002;23:759–68.

    Article  CAS  PubMed  Google Scholar 

  48. Tam W, Hughes SH, Hayward WS, Besmer P. Avian bic, a gene isolated from a common retroviral site in avian leukosis virus-induced lymphomas that encodes a noncoding RNA, cooperates with c-myc in lymphomagenesis and erythroleukemogenesis. J Virol. 2002;76:4275–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Eis PS, Tam W, Sun L, Chadburn A, Li Z, Gomez MF, Lund E, Dahlberg JE. Accumulation of miR-155 and BIC RNA in human B cell lymphomas. Proc Natl Acad Sci. 2005;102:3627–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Iorio MV, Ferracin M, Liu C-G, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M. MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005;65:7065–70.

    Article  CAS  PubMed  Google Scholar 

  51. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006;9:189–98.

    Article  CAS  PubMed  Google Scholar 

  52. Neal CS, Michael MZ, Rawlings LH, Van Der Hoek MB, Gleadle JM. The VHL-dependent regulation of microRNAs in renal cancer. BMC Med. 2010;8:64.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Juan D, Alexe G, Antes T, Liu H, Madabhushi A, Delisi C, Ganesan S, Bhanot G, Liou LS. Identification of a microRNA panel for clear-cell kidney cancer. Urology. 2010;75:835–41.

    Article  PubMed  Google Scholar 

  54. Li S, Chen T, Zhong Z, Wang Y, Li Y, Zhao X. microRNA-155 silencing inhibits proliferation and migration and induces apoptosis by upregulating BACH1 in renal cancer cells. Mol Med Rep. 2012;5:949–54.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

YC and LT participated in the design of this study. WL and YC performed the statistical analysis. DD carried out the study, together with WW, and collected important background information. RA drafted the manuscript. HC and CL conceived of this study, and participated in the design and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This project was supported by research grants from Heilongjiang Provincial Health and Family Planning Commission issued in 2014 (2014-344).

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

  1. Department of Urology, Harbin Medical University Cancer Hospital, Harbin, 150040, Heilongjiang Province, China

    Yongsheng Chen, Lichen Teng, Yan Cao, Dexin Ding, Wentao Wang, Hui Chen & Changfu Li

  2. Intensive Care Unit (ICU) Department, The Second Affiliated Hospital of Harbin Medical University, Harbin, 150086, Province Heilongjiang, China

    Wenhua Liu

  3. Department of Urology, The First Affiliated Hospital of Harbin Medical University, No.31 Youzheng Street, Harbin, 150001, Province Heilongjiang, China

    Ruihua An

Authors
  1. Yongsheng Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lichen Teng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Wenhua Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yan Cao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Dexin Ding
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Wentao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Hui Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Changfu Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Ruihua An
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Ruihua An.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, Y., Teng, L., Liu, W. et al. Identification of biological targets of therapeutic intervention for clear cell renal cell carcinoma based on bioinformatics approach. Cancer Cell Int 16, 16 (2016). https://doi.org/10.1186/s12935-016-0291-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-016-0291-8

Keywords

Cancer Cell International

ISSN: 1475-2867

Contact us