- Research article
- Open Access
Exploring microRNA target genes and identifying hub genes in bladder cancer based on bioinformatic analysis
BMC Urology volume 21, Article number: 90 (2021)
Bladder cancer (BC) is the second most frequent malignancy of the urinary system. The aim of this study was to identify key microRNAs (miRNAs) and hub genes associated with BC as well as analyse their targeted relationships.
According to the microRNA dataset GSE112264 and gene microarray dataset GSE52519, differentially expressed microRNAs (DEMs) and differentially expressed genes (DEGs) were obtained using the R limma software package. The FunRich software database was used to predict the miRNA-targeted genes. The overlapping common genes (OCGs) between miRNA-targeted genes and DEGs were screened to construct the PPI network. Then, gene ontology (GO) analysis was performed through the “cluster Profiler” and “org.Hs.eg.db” R packages. The differential expression analysis and hierarchical clustering of these hub genes were analysed through the GEPIA and UCSC Cancer Genomics Browser databases, respectively. KEGG pathway enrichment analyses of hub genes were performed through gene set enrichment analysis (GSEA).
A total of 12 DEMs and 10 hub genes were identified. Differential expression analysis of the hub genes using the GEPIA database was consistent with the results for the UCSC Cancer Genomics Browser database. The results indicated that these hub genes were oncogenes, but VCL, TPM2, and TPM1 were tumour suppressor genes. The GSEA also showed that hub genes were most enriched in those pathways that were closely associated with tumour proliferation and apoptosis.
In this study, we built a miRNA-mRNA regulatory targeted network, which explores an understanding of the pathogenesis of cancer development and provides key evidence for novel targeted treatments for BC.
Bladder cancer (BC) is a common urogenital malignancy and is associated with a high recurrence rate and high mortality . Despite great progress in the treatment of BC, its recurrence rate is still high, and its prognosis is poor. Drug resistance and tumour metastasis often contribute to the poor prognosis of bladder cancer. Numerous studies have shown that miRNAs play an important role in these processes. It has been reported that miR-22-3p could enhance chemoresistance by targeting NET1. miRNA-124-3p suppressed cell migration and invasion by targeting ITGA3 signalling in bladder cancer , and the miR-626/EYA4 axis was associated with cancer proliferation and metastasis . Accumulating evidence suggests that various miRNAs contribute to cancer initiation, development, and survival prognosis and have potential applications in cancer therapy in numerous malignancies, including BC [4,5,6].
MicroRNAs (miRNAs) are a series of small endogenous single-stranded noncoding RNAs that are usually approximately 21–25 nucleotides in length. They regulate gene expression post-transcriptionally by binding to the 3’-untranslated region of target mRNAs . Many studies have shown that miRNAs regulate tumour proliferation and apoptosis mainly by regulating mRNA expression [8,9,10]. miR-133 was reported to play an essential role in muscle development by regulating many genes and has been identified as an important factor in cancer development . miRNA-217 could accelerate the proliferation and migration of BC by inhibiting KMT2D . Thus, there is a regulatory network between miRNAs and mRNAs that can regulate the occurrence, differentiation and development of tumours.
In this study, we screened differentially expressed microRNAs (DEMs) and differentially expressed genes (DEGs) between cancer tissues and normal tissues via bioinformatic analyses. The overlapping common genes (OCGs) between DEGs and DEM-targeted genes were extracted, classified, and extensively analysed. We constructed a miRNA-targeted gene network and investigated a protein–protein interaction (PPI) network along with its hub genes and significant module. Finally, we determined which hub genes and miRNAs might promote the progression of BC, which might help researchers examine molecular mechanisms involved in disease diagnosis, pathogenesis, and prognosis, thus providing information on precise gene therapy for BC research.
Data acquisition and processing
We downloaded the gene expression profiles GSE112264 (miRNA) and GSE52519 (mRNA) from the National Center for Biotechnology Information (NCBI) GEO database (https://www.ncbi.nlm.nih.gov/geo). A flowchart was created to illustrate our study programme (Fig. 1). The dataset GSE112264 series contained 50 BC samples and 41 normal controls. The GSE52519 series included 9 cancer tissues and 3 normal controls.
Identification of DEMs and DEGs
The DEMs and DEGs were obtained from the GEO database using the R language 3.6.1 version limma package. |log fold change (FC)|> 1 and p-value < 0.05 were considered statistically significant for the exploration of DEGs and DEMs.
PPI network construction and analysis
We predicted the miRNA-targeted genes from the FunRich software database [13, 14] and identified OCGs between miRNA-targeted genes and DEGs. The PPI network was built according to these OCGs using the STRING database (http://string-db.org). Finally, we used Cytoscape 3.7.2  and its plug-in cytoHubba to visualize the PPI network.
GO enrichment analyses of the OCGs
We predicted the potential biological functions of those network genes by Gene Ontology (GO) analysis through the “cluster Profiler” and “org.Hs.eg.db” R packages. A p-value < 0.05 was considered statistically significant.
Hub gene screening and analysis
We used Cytoscape 3.7.2 and its plug-in cytoHubba to find hub genes. Hierarchical clustering of hub genes was built using the UCSC Cancer Genomics Browser (http://genome-cancer.ucsc.edu). The differential expression analysis of hub genes was performed via the GEPIA database (http://gepia.cancer-pku.cn/). To further investigate the potential functions of hub genes, we performed gene set enrichment analysis (GSEA) and multiple GSEA.
All statistical analyses were performed using R version 3.6.1, and a p-value < 0.05 was accepted for statistical significance.
Identification of DEMs and DEGs
Heat map analyses were performed to identify DEGs and DEMs. There were 406 differentially expressed miRNAs, including 127 upregulated and 279 downregulated miRNAs, from the microRNA dataset GSE112264 (Fig. 2). In addition, 370 upregulated genes and 253 downregulated genes were screened from the gene microarray dataset GSE52519.
Overlapping common gene analysis between DEM-targeted genes and DEGs
We screened 5218 DEM-targeted genes according to the FunRich database. After determining the overlap between DEM-targeted genes and DEGs through a Venn diagram, we found 166 OCGs (Fig. 3a) and used the STRING database to visualize the protein relationship. Then, we built a protein–protein interaction (PPI) network (Fig. 3b) and calculated the protein node degrees through cytoHubba (Fig. 3c).
GO enrichment analyses of OCGs
We predicted the potential biological functions of these genes through GO enrichment analysis. GO analysis results were listed in Table 1. The top ten terms of GO annotation in biological process (BP), cellular component (CC) and molecular function (MF) were displayed. The BP-enriched functions were mainly associated with muscle organ development and muscle cell proliferation (Fig. 4a, b). The CC biological functions were mainly enriched in focal adhesion and cell substrate adherent junctions (Fig. 4c, d). The MF biological functions were mainly enriched in actin binding and actin filament binding (Fig. 4e, f).
The network between DEMs and hub genes
Ten hub genes were identified according to the node degrees (Table 2) through Cytoscape and its plug-in cytoHubba (Fig. 5a). Then, the network between DEMs and hub genes was visualized through Cytoscape software (Fig. 5b). According to the miRNA-mRNA targeted network, the network included hsa-miR-103a-3p, hsa-miR-107, hsa-miR-130a-3p, hsa-miR-133b, hsa-miR-142-3p, hsa-miR-193b-3p, hsa-miR-22-3p, hsa-miR-29b-3p, hsa-miR-367-3p, hsa-miR-4295, hsa-miR-4500 and hsa-miR-503-5p, which had targeted relationships with the hub genes. The top 10 hub genes included vascular endothelial growth factor A (VEGFA), ribonucleotide reductase regulatory subunit M2 (RRM2), H2A histone family member X (H2AFX), vinculin (VCL), cell division cycle associated 8 (CDCA8), checkpoint kinase 1 (CHEK1), aurora kinase A (AURKA), centrosomal protein 55 (CEP55), tropomyosin 2 (TPM2), and tropomyosin 1 (TPM1), and their target miRNAs are shown in Table 3.
Hub gene expression and functional analysis
Hierarchical clustering analysis indicated that the hub genes were significantly differentially expressed between the normal and tumour tissues (Fig. 6). GEPIA database analysis indicated that those hub genes were significantly upregulated in the tumour tissues; however, VCL, TPM2, and TPM1 were significantly downregulated (Fig. 7). This result was consistent with hierarchical clustering analysis. Therefore, the results indicated that these hub genes were oncogenes, but VCL, TPM2, and TPM1 were tumour suppressor genes. To further investigate the potential functions of the hub genes, we performed GSEA. The high enrichment plot of those hub genes was most enriched in the “cell cycle” and “nucleotide excision repair” pathways. VEGFA was most enriched in bladder cancer, and VCL was highly related to vascular smooth muscle contraction (Fig. 8). Multiple GSEA also showed that these hub genes with lower expression levels were most enriched in the “vascular smooth muscle contraction”, “snare interactions in vesicular transport” and “abc transporters” pathways. The hub genes with higher expression levels were most enriched in the “cell cycle”, “nucleotide excision repair” and “cell apoptosis” pathways (Fig. 9). These pathways were closely associated with tumour proliferation and apoptosis.
miRNAs regulate the progression of tumours by regulating their target genes, and some miRNAs have been identified as being involved in several types of cancer [16, 17]. Therefore, it is of great significance to study relationships between miRNAs and targeted genes in BC. In our study, we screened the DEGs and DEMs through GSE52519 and GSE112264 expression profiles, respectively. There were 127 upregulated and 279 downregulated miRNAs that were differentially expressed in BC patients. We found that hsa-miR-103a-3p, hsa-miR-107, hsa-miR-130a-3p, hsa-miR-133b, hsa-miR-142-3p, hsa-miR-193b-3p, hsa-miR-22-3p, hsa-miR-29b-3p, hsa-miR-367-3p, hsa-miR-4295, hsa-miR-4500 and hsa-miR-503-5p were potential core miRNAs that play an important role in the development and prognosis of BC.
hsa-miR-29b-3p was the most significantly expressed in the tumour tissues and regulated two hub genes, VEGFA and TPM1. hsa-miR-29b-3p plays an important role in many pathophysiological processes, such as cell proliferation, apoptosis and metastasis [18,19,20], and these processes have been reported in a series of tumour types [21,22,23]. hsa-miR-29b-3p could control tumour proliferation, invasion and angiogenesis by regulating VEGFA . Chou et al.  reported that the hsa-miR‐29b-VEGFA axis was closely related to metastasis and deterioration in breast cancer. Szczyrba et al.  showed that induced VEGF protein expression regulated by hsa-miR‐29b was involved in cancer development. It was shown that hsa-miR-103a-3p played an important role in bladder functions and was significantly correlated with patient survival [26, 27]. Moreover, some studies showed that miR-103a-3p and hsa-miR-107 were significantly upregulated in tumour tissues of BC, and they were positively associated with tumour stages . hsa-miR-133b was also associated with the prognosis of patients , and hsa-miR-142-3p was upregulated in bladder urothelial carcinoma patients . miR-193b-3p could downregulate the cell cycle to inhibit cell proliferation . Promoting miR-4295 transcription could reduce p63α translation and enhance urothelial transformation for tumorigenesis . As a result, these miRNAs have been reported to be closely correlated with the occurrence and development of BC.
We selected 10 hub genes according to node degrees, such as VEGFA, VCL, TPM2, TPM1, RRM2, H2AFX, CHEK1, CEP55, CDCA8 and AURKA. In the present study, VEGFA, RRM2, H2AFX, CHEK1, CEP55, CDCA8 and AURKA were significantly upregulated; however, VCL, TPM2 and TPM1 were significantly downregulated in cancer samples of BC. Furthermore, GEPIA database analysis showed that the result was consistent with the hierarchical clustering analysis and indicated that VCL, TPM2 and TPM1 were tumour suppressor genes and the others were oncogenic genes. Our study showed that VEGFA was the core of these hub genes, which further validated its important role in tumour deterioration and development. VEGFA is a member of the VEGF growth factor family, which can promote the proliferation and migration of vascular endothelial cells. VEGFA is significantly upregulated in many cancers and is essential for tumour proliferation, invasion and metastasis [34, 35]. VCL is a cytoskeletal protein associated with cell–cell and cell–matrix junctions, where it is thought to function as one of several interacting proteins involved in anchoring F-actin to the membrane. Recent research has shown that morphological and mechanical stability are directly related to the actin filament organization used by cancer cells to adapt to altered laminin-rich microenvironments . VCL might affect cancer development through the focal adhesion pathway . TPM2 encodes beta-tropomyosin, a member of the actin filament binding protein family that is poorly expressed in high-grade, relapsed, and metastatic prostate tumours and may be a potential prognostic biomarker in prostate cancer . TPM1 is also a member of the tropomyosin family of highly conserved, widely distributed actin-binding proteins involved in the contractile system of striated and smooth muscles and the cytoskeleton of non-muscle cells. The upregulation of TPM1 inhibited cell proliferation, induced apoptosis, and inhibited cancer growth in BC cells . Moreover, the migratory and invasion ability of cancer cells was enhanced by the destruction of stress fibres and by TPM1-mediated relevant adhesive structures [40, 41]. Additionally, Thorsen et al.  also revealed that the mRNA expression level and protein expression level of TPM1 were downregulated in BC. These results supported our idea of TPM1 as a tumour suppressor. To further investigate the potential functions of these hub genes, we performed GSEA and multiple GSEA. The results showed that these hub genes were most enriched in those pathways that were closely associated with tumour proliferation and apoptosis.
In conclusion, we constructed a miRNA-mRNA regulatory network and showed the connections between hub genes and miRNAs in the regulatory network. This study indicated the potential mechanisms of the miRNA-mRNA regulatory network in the physiological and pathological processes of BC, and it might provide insight into therapy and improve prognosis in the future. Since the number of gene chips included in this study is not large, there may be some shortcomings and contradictions, such as data bias and differences in detection methods. Thus, we are planning to further perform a series of experiments to verify the results of this study and explore the molecular mechanism of the regulatory network.
Availability of data and materials
Publicly available datasets were analysed in this study. This data can be found here: Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo).
Differentially expressed microRNAs
Differentially expressed genes
Overlapping common genes
Kyoto Encyclopedia of Genes and Genomes
Gene Set Enrichment Analysis
Vascular endothelial growth factor A
Ribonucleotide reductase regulatory subunit M2
H2A histone family member X
Cell division cycle associated 8
Checkpoint kinase 1
Aurora kinase A
Centrosomal protein 55
Gene expression omnibus
Martinez Rodriguez RH, Buisan Rueda O, Ibarz L. Bladder cancer: present and future. Med Clin (Barc). 2017;149(10):449–55.
Wang JR, Liu B, Zhou L, Huang YX. MicroRNA-124-3p suppresses cell migration and invasion by targeting ITGA3 signaling in bladder cancer. Cancer Biomark. 2019;24(2):159–72.
Dong W, Bi JM, Liu HW, Yan D, He QQ, Zhou QH, Wang Q, Xie RH, Su YJ, Yang MH, et al. Circular RNA ACVR2A suppresses bladder cancer cells proliferation and metastasis through miR-626/EYA4 axis. Mol Cancer. 2019;18:1–16.
Babaei K, Shams S, Keymoradzadeh A, Vahidi S, Hamami P, Khaksar R, Norollahi SE, Samadani AA. An insight of microRNAs performance in carcinogenesis and tumorigenesis; an overview of cancer therapy. Life Sci. 2020;240:117077.
Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6(11):857–66.
Enokida H, Yoshino H, Matsushita R, Nakagawa M. The role of microRNAs in bladder cancer. Investig Clin Urol. 2016;57(Suppl 1):S60-76.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Amuran GG, Eyuboglu IP, Tinay I, Akkiprik M. New Insights in Bladder Cancer Diagnosis: Urinary miRNAs and Proteins. Med Sci (Basel). 2018;6(4):113.
Dong F, Xu T, Shen Y, Zhong S, Chen S, Ding Q, Shen Z. Dysregulation of miRNAs in bladder cancer: altered expression with aberrant biogenesis procedure. Oncotarget. 2017;8(16):27547–68.
Kutwin P, Konecki T, Borkowska EM, Traczyk-Borszynska M, Jablonowski Z. Urine miRNA as a potential biomarker for bladder cancer detection: a meta-analysis. Cent European J Urol. 2018;71(2):177–85.
Yu H, Lu Y, Li Z, Wang Q. microRNA-133: expression, function and therapeutic potential in muscle diseases and cancer. Curr Drug Targets. 2014;15(9):817–28.
Hou G, Xu W, Jin Y, Wu J, Pan Y, Zhou F. MiRNA-217 accelerates the proliferation and migration of bladder cancer via inhibiting KMT2D. Biochem Biophys Res Commun. 2019;519(4):747–53.
Pathan M, Keerthikumar S, Chisanga D, Alessandro R, Ang CS, Askenase P, Batagov AO, Benito-Martin A, Camussi G, Clayton A, et al. A novel community driven software for functional enrichment analysis of extracellular vesicles data. J Extracell Vesicles. 2017;6(1):1321455.
Benito-Martin A, Peinado H. FunRich proteomics software analysis, let the fun begin! Proteomics. 2015;15(15):2555–6.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Ying Y, Li J, Xie H, Yan H, Jin K, He L, Ma X, Wu J, Xu X, Fang J, et al. CCND1, NOP14 and DNMT3B are involved in miR-502-5p-mediated inhibition of cell migration and proliferation in bladder cancer. Cell Prolif. 2020;53(2):e12751.
Wang F, Wu H, Fan M, Yu R, Zhang Y, Liu J, Zhou X, Cai Y, Huang S, Hu Z, et al. Sodium butyrate inhibits migration and induces AMPK-mTOR pathway-dependent autophagy and ROS-mediated apoptosis via the miR-139-5p/Bmi-1 axis in human bladder cancer cells. FASEB J. 2020;34(3):4266–82.
Li Y, Wang F, Xu J, Ye F, Shen Y, Zhou J, Lu W, Wan X, Ma D, Xie X. Progressive miRNA expression profiles in cervical carcinogenesis and identification of HPV-related target genes for miR-29. J Pathol. 2011;224(4):484–95.
Martinez I, Cazalla D, Almstead LL, Steitz JA, DiMaio D. miR-29 and miR-30 regulate B-Myb expression during cellular senescence. Proc Natl Acad Sci U S A. 2011;108(2):522–7.
Ugalde AP, Ramsay AJ, de la Rosa J, Varela I, Marino G, Cadinanos J, Lu J, Freije JM, Lopez-Otin C. Aging and chronic DNA damage response activate a regulatory pathway involving miR-29 and p53. EMBO J. 2011;30(11):2219–32.
Espinosa-Parrilla Y, Munoz X, Bonet C, Garcia N, Vencesla A, Yiannakouris N, Naccarati A, Sieri S, Panico S, Huerta JM, et al. Genetic association of gastric cancer with miRNA clusters including the cancer-related genes MIR29, MIR25, MIR93 and MIR106: results from the EPIC-EURGAST study. Int J Cancer. 2014;135(9):2065–76.
Sandhu R, Rivenbark AG, Mackler RM, Livasy CA, Coleman WB. Dysregulation of microRNA expression drives aberrant DNA hypermethylation in basal-like breast cancer. Int J Oncol. 2014;44(2):563–72.
Szczyrba J, Nolte E, Hart M, Doll C, Wach S, Taubert H, Keck B, Kremmer E, Stohr R, Hartmann A, et al. Identification of ZNF217, hnRNP-K, VEGF-A and IPO7 as targets for microRNAs that are downregulated in prostate carcinoma. Int J Cancer. 2013;132(4):775–84.
Zhao X, Liu Y, Li Z, Zheng S, Wang Z, Li W, Bi Z, Li L, Jiang Y, Luo Y, et al. Linc00511 acts as a competing endogenous RNA to regulate VEGFA expression through sponging hsa-miR-29b-3p in pancreatic ductal adenocarcinoma. J Cell Mol Med. 2018;22(1):655–67.
Chou J, Lin JH, Brenot A, Kim JW, Provot S, Werb Z. GATA3 suppresses metastasis and modulates the tumour microenvironment by regulating microRNA-29b expression. Nat Cell Biol. 2013;15(2):201–13.
Jiang X, Du L, Duan W, Wang R, Yan K, Wang L, Li J, Zheng G, Zhang X, Yang Y, et al. Serum microRNA expression signatures as novel noninvasive biomarkers for prediction and prognosis of muscle-invasive bladder cancer. Oncotarget. 2016;7(24):36733–42.
Gheinani AH, Kiss B, Moltzahn F, Keller I, Bruggmann R, Rehrauer H, Fournier CA, Burkhard FC, Monastyrskaya K. Characterization of miRNA-regulated networks, hubs of signaling, and biomarkers in obstruction-induced bladder dysfunction. JCI Insight. 2017;2(2):e89560.
Zhong Z, Lv M, Chen J. Screening differential circular RNA expression profiles reveals the regulatory role of circTCF25-miR-103a-3p/miR-107-CDK6 pathway in bladder carcinoma. Sci Rep. 2016;6:30919.
Yu QF, Liu P, Li ZY, Zhang CF, Chen SQ, Li ZH, Zhang GY, Li JC. MiR-103/107 induces tumorigenicity in bladder cancer cell by suppressing PTEN. Eur Rev Med Pharmacol Sci. 2018;22(24):8616–23.
Shi Z, Kadeer A, Wang M, Wen B, Li M, Huang J, Gao Y, Liu E, Liu D, Jia D, et al. The deregulation of miR-133b is associated with poor prognosis in bladder cancer. Pathol Res Pract. 2019;215(2):354–7.
Zhang ZC, Huang Y, Wang XJ, Wang M, Ma LL. Expression of circulating microRNAs in patients with bladder urothelial carcinoma. Beijing Da Xue Xue Bao Yi Xue Ban. 2013;45(4):532–6.
Lin SR, Yeh HC, Wang WJ, Ke HL, Lin HH, Hsu WC, Chao SY, Hour TC, Wu WJ, Pu YS, et al. MiR-193b mediates CEBPD-induced cisplatin sensitization through targeting ETS1 and cyclin D1 in human urothelial carcinoma cells. J Cell Biochem. 2017;118(6):1563–73.
Jin H, Xu J, Guo X, Huang H, Li J, Peng M, Zhu J, Tian Z, Wu XR, Tang MS, et al. XIAP RING domain mediates miR-4295 expression and subsequently inhibiting p63alpha protein translation and promoting transformation of bladder epithelial cells. Oncotarget. 2016;7(35):56540–57.
Paez-Ribes M, Allen E, Hudock J, Takeda T, Okuyama H, Vinals F, Inoue M, Bergers G, Hanahan D, Casanovas O. Antiangiogenic therapy elicits malignant progression of tumors to increased local invasion and distant metastasis. Cancer Cell. 2009;15(3):220–31.
Yu X, Ye F. Role of Angiopoietins in Development of Cancer and Neoplasia Associated with Viral Infection. Cells. 2020;9(2):457.
Lekka M, Pabijan J, Orzechowska B. Morphological and mechanical stability of bladder cancer cells in response to substrate rigidity. Biochim Biophys Acta Gen Subj. 2019;1863(6):1006–14.
Wu Z, Wang S, Jiang F, Li Q, Wang C, Wang H, Zhang W, Xue P, Wang SL. Mass spectrometric detection combined with bioinformatic analysis identified possible protein markers and key pathways associated with bladder cancer. Gene. 2017;626:407–13.
Varisli L. Identification of new genes downregulated in prostate cancer and investigation of their effects on prognosis. Genet Test Mol Biomark. 2013;17(7):562–6.
Liu G, Zhao X, Zhou J, Cheng X, Ye Z, Ji Z. Long non-coding RNA MEG3 suppresses the development of bladder urothelial carcinoma by regulating miR-96 and TPM1. Cancer Biol Ther. 2018;19(11):1039–56.
Redwood C, Robinson P. Alpha-tropomyosin mutations in inherited cardiomyopathies. J Muscle Res Cell Motil. 2013;34(3–4):285–94.
Lin J, Shen J, Yue H, Cao Z. miRNA1835p1 promotes the migration and invasion of gastric cancer AGS cells by targeting TPM1. Oncol Rep. 2019;42(6):2371–81.
Thorsen K, Sorensen KD, Brems-Eskildsen AS, Modin C, Gaustadnes M, Hein AM, Kruhoffer M, Laurberg S, Borre M, Wang K, et al. Alternative splicing in colon, bladder, and prostate cancer identified by exon array analysis. Mol Cell Proteomics. 2008;7(7):1214–24.
The authors wish to thank the National Center for Biotechnology Information (NCBI) GEO database (https://www.ncbi.nlm.nih.gov/geo).
Ethics approval and consent to participate
Consent to publish
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wu, H., Jiang, W., Ji, G. et al. Exploring microRNA target genes and identifying hub genes in bladder cancer based on bioinformatic analysis. BMC Urol 21, 90 (2021). https://doi.org/10.1186/s12894-021-00857-w
- Bladder cancer
- MiRNA-mRNA network
- Hub genes