- 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)
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.
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).
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).
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.
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.
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.
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
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.
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.
Escudier B, Kataja V. Renal cell carcinoma: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2010;21:v137–9.
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.
Cohen HT, Mcgovern FJ. Renal-cell carcinoma. N Engl J Med. 2005;353:2477–90.
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.
Kaelin WG. Treatment of kidney cancer. Cancer. 2009;115:2262–72.
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.
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.
Rottiers V, Näär AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13:239–50.
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.
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.
Allison DB, Cui X, Page GP, Sabripour M. Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006;7:55–65.
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.
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.
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.
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.
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.
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.
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.
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.
Smyth GK, Speed T. Normalization of cDNA microarray data. Methods. 2003;31:265–73.
Smyth GK. Limma: linear models for microarray data, in bioinformatics and computational biology solutions using R and Bioconductor. Springer. 2005. p. 397–420.
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.
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.
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.
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.
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.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–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.
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.
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.
Grignon DJ, Che M. Clear cell renal cell carcinoma. Clin Lab Med. 2005;25:305–16.
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.
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.
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.
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.
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.
Marian B, Danner M. Skin tumor promotion is associated with increased type V collagen content in the dermis. Carcinogenesis. 1987;8:151–4.
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.
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.
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.
Zheng J. Energy metabolism of cancer: glycolysis versus oxidative phosphorylation (Review). Oncol Lett. 2012;4:1151–7.
Semenza GL. Oxygen-dependent regulation of mitochondrial respiration by hypoxia-inducible factor 1. Biochem J. 2007;405:1–9.
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.
Li P, Jiao J, Gao G, Prabhakar BS. Control of mitochondrial activity by miRNAs. J Cell Biochem. 2012;113:1104–10.
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.
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.
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.
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.
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.
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.
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.
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.
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
Corresponding author
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12935-016-0291-8