- Primary research
- Open Access
Identification and analysis of long non-coding RNA related miRNA sponge regulatory network in bladder urothelial carcinoma
Cancer Cell International volume 19, Article number: 327 (2019)
The aim of this study was to investigate the regulatory network of lncRNAs as competing endogenous RNAs (ceRNA) in bladder urothelial carcinoma (BUC) based on gene expression data derived from The Cancer Genome Atlas (TCGA).
Materials and methods
RNA sequence profiles and clinical information from 414 BUC tissues and 19 non-tumor adjacent tissues were downloaded from TCGA. Differentially expressed RNAs derived from BUC and non-tumor adjacent samples were identified using the R package “edgeR”. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed using the “clusterProfiler” package. Gene ontology and protein–protein interaction (PPI) networks were analyzed for the differentially expressed mRNAs using the “STRING” database. The network for the dysregulated lncRNA associated ceRNAs was then constructed for BUC using miRcode, miRTarBase, miRDB, and TargetScan. Cox regression analysis was performed to identify independent prognostic RNAs associated with BUC overall survival (OS). Survival analysis for the independent prognostic RNAs within the ceRNA network was calculated using Kaplan–Meier curves.
Based on our analysis, a total of 666, 1819 and 157 differentially expressed lncRNAs, mRNAs and miRNAs were identified respectively. The ceRNA network was then constructed and contained 59 lncRNAs, 23 DEmiRNAs, and 52 DEmRNAs. In total, 5 lncRNAs (HCG22, ADAMTS9-AS1, ADAMTS9-AS2, AC078778.1, and AC112721.1), 2 miRNAs (hsa-mir-145 and hsa-mir-141) and 6 mRNAs (ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3) were found to be related to OS. Two lncRNAs (ADAMTS9-AS1 and ADAMTS9-AS2) and 4 mRNA (DUSP2, JUN, MAP1B, and TMEM100) were validated using GEPIA. Thirty key hub genes were identified using the ranking method of degree. KEGG analysis demonstrated that the majority of the DEmRNAs were involved in pathways associated with cancer.
Our findings provide an understanding of the important role of lncRNA–related ceRNAs in BUC. Additional experimental and clinical validations are required to support our findings.
Bladder cancer (BC) is the most common urologic cancer. Approximately 429,800 newly diagnosed cases and 165,100 deaths are recorded worldwide every year . The main pathological bladder cancer types are bladder urothelial carcinoma (BUC), bladder adenocarcinoma and bladder squamous cell carcinoma. The most common type is BUC and accounts for more than 90% of all BCs . Unfortunately, the high recurrence rate is a characteristic of BUC [3, 4]. Currently, the principal treatment strategy for BUC consists of surgery and adjuvant combination chemotherapy. However, chemotherapy resistance reduces the sensitivity of BUC to chemotherapeutic drugs and frequently results in treatment failure resulting in BUC clinical management being a major challenge . In addition, only a limited number of biomarkers are available for diagnosing BUC compared to other cancers. Hence, identifying sensitive and specific BUC biomarkers, as well as therapeutic targets for BUC are critically needed.
Long non-coding RNAs (lncRNAs) are a subtype of ncRNAs with transcript lengths over 200 nucleotides and have recently attracted increased attention . lncRNAs were initially regarded as transcriptional noise without the capacity to encode proteins . However, growing evidence has demonstrated that lncRNAs may play crucial biological roles in a variety of biological processes that are associated with carcinogenesis and cancer metastasis . With regards to bladder cancer, several studies have suggested that lncRNAs may function as oncogenes or tumor suppressors and may affect overall patient survival and mortality [9, 10]. To date, only a few lncRNAs have been verified experimentally, but their roles in regulating gene expression remains to be deciphered.
Considerable efforts have been made to demonstrate how lncRNAs exert their diverse biological functions in human malignant tumors. Rapid progression has been made to elucidate the role of lncRNAs in miRNA function. miRNAs are endogenous single-stranded RNA with lengths between 20 and 25 nucleotides that do not encode proteins. They repress gene expression by complement binding to their target mRNA sequences (i.e. microRNA response element, MRE) . In 2011, Salmena et al.  proposed the competing endogenous RNA (ceRNA) hypothesis, which states that mRNAs, transcribed pseudogenes, and lncRNAs could act as natural miRNA “sponges” and inhibit miRNA function by competing with the binding of one or more MREs in complex and comprehensive regulatory networks, leading to pathogenic conditions. The ceRNA regulation theory has been proven to be involved in bladder cancer initiation and progression in several studies [13, 14]. Similarly, lncRNAs acting as ceRNAs has also been reported in other cancers [15,16,17]. Recently, Kouhsar et al.  constructed a ceRNA network related to the staging of Non-Muscle Invasive Bladder Cancer (Ta and T1) derived from public data sources. They identified several biomarkers associated with tumor stage. We hypothesized that lncRNAs may function as ceRNAs during BUC initiation and progression. Understanding how lncRNAs function as ceRNAs will be important in deciphering BUC carcinogenesis.
In the present study, we aimed to decipher the regulatory ceRNA network of lncRNAs–miRNAs–mRNAs in BUC by analyzing gene expression data. This was performed using bioinformatics prediction and correlation analyses. In addition, using clinical trials and survival analyses, we identified potential prognostic genes.
Materials and methods
RNA sequence data from 406 BUC patients were retrieved from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/) in 2018. LncRNA, miRNA, and mRNAseq data were obtained using the Data Transfer Tool (provided by GDC Apps). Patient clinical information was also downloaded using the Data Transfer Tool. Sequencing data derived from the Illumina HiSeq RNAseq and Illumina HiSeq miRNAseq platforms were publicly available. This study met the publication guidelines stated by TCGA (https://cancergenome.nih.Gov/publications/publicationguidelines). All data used in the study were obtained from TCGA, and hence ethics approval and informed consent were not required.
Differential expression analysis
The mRNAseq and lncRNAseq data derived from 414 BUC tissue samples and 19 non-tumor adjacent tissue samples were downloaded from TCGA. The BUC miRNAseq data were derived from 418 BUC tissue samples and 19 non-tumor adjacent tissue samples. For tumor and non-tumor group comparison, differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs) and lncRNAs (DElncRNAs) were identified using the “edgeR” package (http://bioconductor.org/packages/release/bioc/html/edgeR.html) with a cut-off criteria of |log2 (fold change [FC])| > 2.0 and adjusted P value < 0.01 . Differentially expressed lncRNAs (DElncRNAs) were defined and annotated using the Encyclopedia of DNA Elements (ENCODE), which included 15,877 human lncRNAs. All P-values used the False discovery rate (FDR) to correct for statistical significance of multiple testing (Benjamini–Hochberg method) . FDR significance level was set at 0.05.
KEGG enrichment analysis of DEmRNAs
KEGG enrichment analysis was performed using the “clusterProfiler” package in R software based on the retrieved DEmRNAs and visualized using the Cytoscape v 3.5.1 software.
Protein–protein interaction (PPI) network
To understand the interactions of the DEmRNAs, we constructed a PPI network using the Search Tool for the Retrieval of Interacting Genes (STRING, http://string.embl.de/). Combined scores greater than 0.4 were considered statistically significant. The PPI network was visualized using the Cytoscape v 3.5.1 software. Subsequently, the top 30 mRNAs were identified using the ranking method of degree. In addition, gene ontology (GO) enrichment analysis was performed using STRING to functionally annotate the DEmRNAs in BUC.
In order to investigate the role of the differentially expressed RNAs in the ceRNA network, a dysregulated lncRNA–miRNA–mRNA ceRNA network was constructed and visualized using the Cytoscape v 3.5.1 software. LncRNA–miRNA interactions were predicted using miRcode (http://www.mircode.org/). miRNA–targeted mRNAs were predicted using the miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/), miRDB (http://www.mirdb.org/), and TargetScan (http://www.targetscan.org/vert_71/). Each regulatory pair of miRNAs and mRNAs were verified using quantitative real-time PCR, western blot, reporter assays, microarrays, and next-generation sequencing data derived from miRTarBase.
Independent prognostic factors for overall survival (OS)
Independent prognostic differentially expressed RNAs for OS were identified by univariate and multivariate cox regression analyses. The risk score (RS) was estimated using the following formula:
where Exp(i) denotes the expression value for RNA; n represents the number of RNA modules, and R(i) denotes the estimated regression coefficient of RNA. If the RS for a given sample was greater than the mean RS of all samples, the latter was regarded as a high-risk sample, otherwise, it was regarded as a low-risk sample. Kaplan–Meier method was used to evaluate the survival curves in the high- and low-risk groups. Additionally, the sensitivity and specificity were accessed using the receiver operating characteristic (ROC) curves and area under the ROC curves (AUC values).
Survival curves were plotted using the “survival” package in R for the independent prognostic RNAs that were identified. Long-rank test was used to evaluate statistical significance and P < 0.05 was considered statistically significant.
Validating the prognostic value of lncRNAs and mRNAs in the ceRNA network using GEPIA
Independent prognostic lncRNAs and mRNAs were validated using Gene Expression Profiling Interactive Analysis (GEPIA,http://gepia.cancerpku.cn/index.htm), which was based on RNA sequencing data from 9736 tumors and 8587 normal samples in the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) dataset project . In addition, the correlation between lncRNAs and mRNAs was confirmed using the Pearson correlation statistic. The correlation of two RNAs was considered significant when R was greater than 0.4 and the P-value was less than 0.05.
Identifying DElncRNAs, DEmRNAs, and DEmiRNAs
A total of 666 DElncRNAs, 1819 DEmRNAs and 157 DEmiRNAs were identified with |log2FC| > 2.0 and adjusted P-value < 0.01 using the “edgeR” package. Heat maps with complete linkage clustering of differentially expressed RNAs were performed using the “gplots” package (Additional file 1: Figure S1A–C). The results identified; 420 (63.1%) up-regulated and 246 (36.9%) down-regulated DElncRNAs, 1030 (56.6%) up-regulated and 789 (43.4%) down-regulated DEmRNAs, and 131 (83.4%) up-regulated and 26 (16.6%) down-regulated DEmiRNAs.
Pathway enrichment analysis of DEmRNAs
To investigate the mechanisms associated with BUC tumorigenesis, 1819 DEmRNAs were used for KEGG enrichment analysis. The threshold was set at P < 0.01. The top 15 significantly enriched pathways are presented in Table 1. The results showed that the majority of DEmRNAs were enriched for “neuroactive ligand–receptor interaction”, “viral carcinogenesis”, “protein digestion and absorption”, “ECM–receptor interaction”, and “cAMP signaling pathway”. In addition, a network of pathways and DEmRNAs was constructed based on the above KEGG analysis and visualized using Cytoscape v 3.5.1 (Additional file 2: Figure S2). Interestingly, several DEmRNAs were related to more than one pathway, such as SLC8A1, CAMK2B, and CACNA2D1. CAMK2B was enriched for Circadian entrainment, calcium signaling pathway, adrenergic signaling in cardiomyocytes, cAMP signaling pathway, oxytocin signaling pathway, and insulin secretion. These DEmRNAs may be associated with BUC carcinogenesis and progression.
Functional enrichment analysis of DEmRNAs
Functional enrichment of the top 600 DEmRNAs were analyzed using STRING (300 significantly up-regulated and 300 significantly down-regulated DEmRNAs based on log fold change (log FC)). The GO function for the various genes were divided into biological processes (BP), cellular function (CF) and molecular component (MC). The GO results are presented in Fig. 1. The results from the GO functional enrichment analysis suggested that DEmRNAs were significantly enriched in multicellular organisms (ontology: BP), binding (ontology: CF) and extracellular region (MC). These results demonstrated that the enriched DEmRNAs were associated with BUC proliferation and migration.
Protein–protein interaction (PPI) network
The interrelationship between the top 600 significant DEmRNAs (300 significantly up-regulated and 300 significantly down-regulated DEmRNAs based on log fold change (log FC)) were retrieved from the STRING database to construct the PPI network (Fig. 2). The PPI network consisted of 418 nodes and 1937 edges. The nodes denoted DEmRNAs, while the edges denoted interactions among the DEmRNAs. Additionally, the top 30 mRNAs in the PPI network were analyzed using the ranking method of degree (Additional file 3: Figure S3). KEGG pathway enrichment analysis for the 30 mRNAs identified with a high degree was performed using the “clusterProfiler” package in the R software with a P < 0.05 as the cut-off criteria. The results demonstrated that these mRNAs were enriched for pathways related to “Systemic lupus erythematosus”, “Alcoholism”, “viral carcinogenesis” and “transcriptional dysregulation in cancer” (Additional file 4: Figure S4A). Furthermore, we identified 21 (HIST2H2BF, HIST1H2BO, HIST1H2BH, HIST1H2BE, HIST1H2BF, HIST1H4E, HIST1H2BM, HIST1H4D, HIST1H2BI, HIST1H4C, HIST1H4B, HIST1H4A, HIST1H2BB, HIST1H4F, HIST1H4L, IL6, HIST1H3G, HIST1H3B, HIST1H3C, HIST1H3F and HIST1H3I) of these 30 DEmRNAs that were enriched for cancer-related pathways: “viral carcinogenesis” or “transcriptional dysregulation in cancer”. These 21 DEmRNAs were classified as hub genes. In addition, a network linking the pathways and mRNAs were constructed and visualized using Cytoscape v 3.5.1 (Additional file 4: Figure S4B).
Construction of the ceRNA network for BUC
To better understand the role of the identified differentially expressed RNAs in BUC, a dysregulated ceRNA network based on DElncRNA–DEmiRNA–DEmRNA interactions was constructed using Cytoscape v 3.5.1 (Fig. 3). The regulatory relationship between DEmiRNAs and DElncRNAs pairs were retrieved from miRcode (http://www.mircode.org/). All the DEmiRNAs and DElncRNAs involved in the ceRNA network meet the cut-off criteria (|log2FC| > 2.0 and adjusted P-value < 0.01) mentioned above. We found 259 lncRNA–miRNA interaction pairs containing 59 lncRNAs and 23 DEmiRNAs from the miRcode database. Subsequently, we searched for mRNAs that were targeted by the 23 DEmiRNAs using miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/), miRDB (http://www.mirdb.org/) and TargetScan (http://www.targetscan.org/vert_71/). All the DEmRNAs meet the cut-off criteria mentioned above (|log2FC| > 2.0 and adjusted P-value < 0.01). Finally, 52 DEmRNAs including 70 miRNA–mRNA interaction pairs were included in our ceRNA network.
Independent prognostic factors for overall survival
59 DElncRNAs, 23 DEmiRNAs and 52 DEmRNAs in the ceRNA network were included in univariate cox regression analysis. 17 DElncRNAs, 6 DEmiRNAs and 19 DEmRNAs (P < 0.05 in univariate cox regression analysis) were identified for multivariable cox regression analysis. Multivariable analysis suggested that 7 lncRNAs (HCG22, ADAMTS9-AS1, ADAMTS9-AS2, AC078778.1, AC112721.1, LINC00525 and NAV2-AS2), 3 DEmiRNAs (hsa-mir-145, hsa-mir-141 and hsa-mir-373), and 7 DEmRNAs (ZEB1, TMEM100, MAP1B, DUSP2, JUN, AIFM3 and MEST) were closely related to OS in BUC patients (P < 0.05) (Table 2 and heat maps in Additional file 5: Figure S5A–C). Based on risk scores (RS) of the independent prognostic RNAs identified above, patients included in the present study were assigned into high- and low-risk groups. The low-risk group had a significantly better prognosis compared to the high-risk group (Fig. 4a for lncRNA, B for miRNA and C for mRNA). The results showed that the 5-year OS for the low-risk group was 56.4%, 55.0%, and 55.8% respectively, while it was 26.5%, 29.1% and 27.8% for the high-risk group, respectively. To determine the prognostic power of the 7 lncRNAs, 3 DEmiRNAs, and 7 DEmRNAs identified above, time-dependent receiver operating characteristic (ROC) curve analysis was performed and the area under the curve (AUC values) was estimated. AUC values for the independent prognostic DElncRNAs, DEmiRNAs, and DEmRNAs were 0.707, 0.624 and 0.681, respectively (Fig. 5, a for lncRNA, b for miRNA, C for mRNA), indicating good specificity and sensitivity.
Kaplan–Meier curve analysis was performed to determine the OS for the independent prognostic RNAs. One patient was lost during follow-up and was excluded from the survival analysis. Five DElncRNAs were significantly related to OS, of which, four DElncRNAs (HCG22, ADAMTS9-AS1, ADAMTS9-AS2, and AC112721.1) were negatively related to OS (Fig. 6b–e), while AC078778.1 was positively related to OS (log-rank P < 0.05) (Fig. 6a). In addition, Kaplan–Meier curve analysis for the three DEmiRNAs and seven DEmRNAs showed that two DEmiRNAs (hsa-mir-141 and hsa-mir-145) and 6 DEmRNAs (ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3) were significantly related to OS (log-rank P < 0.05) (Figs. 7a, b, 8a–f).
Prognostic value of lncRNAs and mRNAs in the ceRNA network analyzed using GEPIA
GEPIA was used to validate the expression levels and prognostic value of the 5 independent lncRNAs. The expression levels of HCG22, ADAMTS9-AS1 and ADAMTS9-AS2 were negative or nearly negative, while AC112721.1 expression levels were positive in BUC tissues. These findings were concordant with previous results using TCGA analysis. However, only ADAMTS9-AS1 and ADAMTS9-AS2 expression levels were down-regulated with statistical significance in BUC. The expression levels of these 4 lncRNAs are shown in Fig. 9. However, AC078778.1 was not found in GEPIA, while the overall survival of the 2 lncRNAs, ADAMTS9-AS1 and ADAMTS9-AS2, were estimated using GEPIA. The results are shown in Fig. 10 and were similar to our previous results. Lower expression of ADAMTS9-AS1 and ADAMTS9-AS2 were associated with a good prognosis in BUC patients.
Similar to the lncRNAs, 6 independent prognostic mRNAs (ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3) in the ceRNA network were also analyzed using GEPIA. The six mRNA expression levels are shown in Fig. 11a–f. ZEB1, TMEM100, MAP1B, DUSP2 and JUN were down-regulated with statistical significance, while AIFM3 was up-regulated without statistical significance in BUC. The correlation with overall survival of these 5 mRNAs, ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3 was also estimated using GEPIA. The results are shown in Fig. 12a–e. Four mRNAs (TMEM100, MAP1B, DUSP2, and JUN) except ZEB1 (P = 0.073) had statistical significance with OS in BUC patients. Lower expression of JUN, MAP1B, and TMEM100 was associated with a good prognosis, while lower expression of DUSP2 was associated with a poor prognosis.
Correlation of the independent prognostic factors; lncRNAs and mRNAs
With regards to the lncRNA and mRNA independent prognostic factors involved in ceRNA network, we observed that ADAMTS9-AS1 interacted with ZEB1 through hsa-mir-96, and ADAMTS9-AS2 interacted with three DEmRNAs (TMEM100, DUSP2, and ZEB1) through seven different DEmiRNAs (hsa-mir-96, hsa-mir-372, hsa-mir-183, hsa-mir-200a, hsa-mir-141, hsa-mir-373, and hsa-mir-205). Pearson correlation analysis was performed to verify the correlation of the independent prognostic lncRNA and mRNA factors. The results are shown in Fig. 13 (a for ADAMTS9-AS2 and TMEM100, B for ADAMTS9-AS2 and ZEB1, C for ADAMTS9-AS1 and ZEB1). We found that ADAMTS9-AS1 had a strong positive correlation with ZEB1, and ADAMTS9-AS2 had a strong positive correlation with ZEB1 and TMEM100. ADAMTS9-AS1 and ADAMTS9-AS2 interacted with ZEB1 and TMEM100 during BUC development. At present, no studies have demonstrated the relationship between ADAMTS9-AS1 and ZEB1, or ADAMTS9-AS2 and ZEB1 or TMEM100 in cancer.
In this study, five DElncRNA factors (HCG22, ADAMTS9-AS1, ADAMTS9-AS2, AC078778.1, and AC112721.1) in the ceRNAs network were identified as independent prognostic factors for OS in BUC patients. HCG22 expression levels have been reported to be down-regulated in oral cancer and its low expression was associated with poor survival in a recent study based on TCGA data analysis . Lu et al.  investigated HCG22 expression levels in 20 oral cavity and oropharyngeal squamous cell carcinoma (OSCC) samples and 10 control samples by qRT-PCR. They demonstrated that HCG22 was downregulated in OSCC tissues compared to controls, while no association was observed between HCG22 expression levels and overall survival. ADAMTS9-AS1 (ADAMTS9 antisense transcript), is a novel lncRNA without any functional annotation but could interact with two RNA-binding proteins, DCGR8 and FUS . lncRNAs play an important role in cancer mainly via their associations with RNA-binding g proteins, which include HOTTIP, MaLAT1, H19, and HOTAIR. They participate in several biological pathways involved in cell differentiation and proliferation, apoptosis and tumorigenesis by interacting with RNA-binding proteins in hepatocellular carcinoma . We hypothesized that ADAMTS9-AS1 may play a role in the development of cancer. Wang et al.  found that ADAMTS9-1 and ADAMTS9-2 expression levels were decreased in malignant epithelial ovarian cancer tissues compared to normal ovary tissues and benign ovarian cysts using lncRNA and mRNA microarray analysis. These results were confirmed using 8 normal ovarian, 17 benign ovarian cysts and 15 malignant epithelial ovarian cancer samples by qPCR assays. Low ADAMTS9-AS2 levels were found to be a significant independent predictor of poor survival in glioma patients . Liu et al.  suggested that lncRNA ADAMTS9-AS2 could suppress cancer progression by inhibiting miR-223-3p and activating TGFBR3. Additionally, increased ADAMTS9-AS2 levels could reduce lung cancer tumor size and improve OS. However, to date, no studies have been performed to determine the role of AC078778.1 and AC112721.1 in cancer.
miRNAs are involved in multiple roles during carcinogenesis. In this study, we found two independent prognostic DEmiRNA factors (hsa-mir-141 and hsa-mir-145) that were involved in the ceRNA network. miR-141 was found to be up-regulated in malignant bladder tissue samples compared to healthy tissues and was a favorable prognostic biomarker . microRNA-141 (hsa-mir-141) has been shown to exert a regulatory role during epithelial to mesenchymal transition process, and its expression levels have been associated with tumorigenicity and invasiveness in several human cancers. hsa-mir-141 has been shown to be associated with the development of certain epithelial cancer cell types, including prostate , colorectal  and breast cancer . Huang et al.  demonstrated that miR-141 could inhibit gastric cancer cell proliferation and tumor growth, while low miR-141 levels were associated with poor prognosis. Another study on gastric cancer suggested that miR-141 could play an important anti-tumor role by interacting with MEG3 and targeting E2F3 during gastric cancer pathogenesis and may be a therapeutic target. miR-145 has been frequently observed to be down-regulated in cancers and restoration of miR-145 levels suppressed cancer cell invasion by reversing the EMT phenotype . Tan et al.  demonstrated that TUG1 promoted bladder cancer cell metastasis and radio-resistance by negatively regulating miR-145 expression.
In the present study, 6 prognostic DEmRNA factors (ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3) were involved in the ceRNAs network and functioned as independent prognostic factors for OS in BUC patients. Several studies have demonstrated that ZEB1 was significantly overexpressed in bladder cancer tissues compared to normal healthy adjacent tissues . Li et al.  reported that ZEB1 was significantly overexpressed in bladder cancers compared to normal tissues, and played a crucial role during VM formation, and was closely associated with invasion, metastasis and poor prognosis of malignant tumors [38, 39]. Transmembrane protein 100 (TMEM100), located at 17q32, was first identified as a transcript in the mouse genome. Han et al.  found that TMEM100 could function as a tumor suppressor by inhibiting the growth and metastasis of non-small-cell lung cancer via the inhibition of the TNF pathway. Low TMEM100 expression levels were associated with poor prognosis. Similar results have been reported for hepatocellular carcinoma . MAP1B, which encodes for the microtubule-associated protein 1B (MAP1B), is one of the main cytoskeletal proteins. Several studies have demonstrated that MAP1B plays an important role in a number of cellular processes, including synaptic transmission, autophagy, and cancer [41,42,43]. DUSP2 is a member of the class 1 DUSP family of proteins and is localized in the nucleus. DUSP2 levels are significantly decreased in bladder cancer and low expression of DUSP2 is correlated with poor prognosis . JUN (C-Jun, AP-1 transcription factor subunit) is specifically phosphorylated by JNK and plays a central role in the AP-1 complex. It is involved in cellular DNA damage response by regulating the expression of several genes [45, 46]. c-Jun is a proto-oncogene and is involved in transformation and tumor development [47, 48]. AIFM3 (apoptosis-inducing factor mitochondria associated 3) is a gene with homology to apoptosis-inducing factor (AIF). AIF induces apoptosis in a caspase-dependent manner. AIFM3 was observed to be highly expressed in breast cancer tissues and associated with shorter overall survival and disease-free survival .
BUC specific DEmRNA pathways were assessed using KEGG pathway analysis. Our results demonstrated that the majority of DEmRNAs were frequently enriched for cancer-related pathways. Of these, the neuroactive ligand-receptor interaction signaling pathway has been reported to be associated with the progression of renal cell carcinoma in bioinformatics studies . In addition, previous studies have suggested that cAMP-related signaling could control apoptosis induction and cell growth [51, 52], while another study demonstrated that cAMP was an inhibitor of cell cycle progression and apoptosis in gastric cancer cells . Zhang et al.  showed that the ECM–receptor interaction pathway played a significant role in tumor progression and metastasis. We found several overlapping DEmRNAs that were involved in multiple pathways, such as SLC8A1, CAMK2B, and CACNA2D1, and have been demonstrated to play important roles in cancer pathogenesis. Muñoz reported that lower levels of SLC8A1, which was at least partly mediated by miR-223, was associated with reduced calcium and apoptosis levels in penile carcinoma . Based on an integrative meta-analysis, CAMK2B was found to be associated with the development of cancer cachexia . Recently, Feng et al.  assayed 15 colorectal cancer tissues and 10 paracancerous tissues using microarrays and found that CAMK2B was involved in the progression of Fusobacterium nucleatum-induced colorectal cancer. Another study also reported that CAMK2B played an important role in glioblastoma multiforme using bioinformatics analysis on publicly available datasets . High expression levels of CACNA2D1 in epithelial ovarian cancers were significantly correlated with histological subtypes, advanced FIGO stages and tumor differentiation . In addition, the majority of PLB4 hotspot mutations are gain-of function mutations that have been demonstrated to be involved in uveal melanoma tumorigenesis by activating the same signaling pathway .
In terms of GO functional enrichment analysis, we found that DEmRNAs were significantly enriched for the regulation of multiple processes, such as binding, protein binding, receptor binding, ion channel activity, and endopeptidase activity. Interestingly, several studies have reported that these functions overlapped in different cancers. Several studies have demonstrated that knocking down ezrin and P65 expression induces tumor metastasis in different cancers [61, 62]. Furthermore, Tang et al.  demonstrated that ezrin and P65 were physically associated with one another. We hypothesize that the interaction between ezrin and P65 is associated with the activation of the NF-κB pathway leading to breast cancer metastasis. A previous study suggested that Estrogen receptor β could increase the levels of miR-92a by binding to the estrogen-response-element (ERE) leading to a decrease in DAB2IP tumor suppressor expression to ultimately promote bladder cancer growth and invasion . Numerous studies have demonstrated that ion channels play an important role in tumorigenesis and progression, such as inducing neo-angiogenesis , apoptosis resistance , proliferative potential  as well as cell migration and invasiveness [68, 69]. In terms of endopeptidase activity, Zhu et al. . demonstrated that asparaginyl endopeptidase (AEP) was highly expressed in tissues and ascites of patients with epithelial ovarian cancer and promoted tumor growth and progression both in vivo and in vitro.
We then constructed a PPI network to identify hub DEmRNAs. Proteins that corresponded to genes were used to build the PPI network, and the top 30 DEmRNAs with a high degree were selected. KEGG pathway analysis of these 30 DEmRNAs was performed using the “clusterProfiler” package. The results showed that 21 of these 30 DEmRNAs were enriched for “viral carcinogenesis” and “transcriptional dysregulation in cancer”. The 21 identified DEmRNAs may play an important role in cancer. For example, IL6, HIST1H3C, and HIST1H3G, which were classified with high degrees are present in pathways associated with transcriptional dysregulation in cancer. In addition, these three genes have been previously reported to be closely associated with tumorigenesis and development [71,72,73].
In terms of correlation analysis, a positive correlation between ZEB1 and ADAMTS9-AS1-AS2 and TMEM100 was observed using the Pearson correlation statistic based on the expression levels of these genes. The correlation between these four genes may play an important role in the initiation and progression of BUC. Only a few studies on their interactions have been published in public databases, such as Pubmed and Embase.
Several limitations of the present study should be stated. First, the number of normal bladder tissues (19 samples) was limited and may have compromised the reliability of our results. Second, BUC patient information from TCGA was not validated using experimental procedures. Third, we only investigated the ceRNAs network associated with lncRNAs, miRNAs, and mRNAs, and did not include other regulatory models. Finally, several novel lncRNAs with significant clinical value needs to be investigated to determine their functional role and underlying mechanism during BUC carcinogenesis. We verified the independent prognostic lncRNAs and mRNAs using GEPIA and this supported our findings. Our results provide a better understanding of lncRNA-related ceRNAs and its important role in BUC. Additional experimental and clinical studies are needed to validate our findings.
Five independent prognostic DElncRNAs (HCG22, ADAMTS9-AS1, ADAMTS9-AS2, AC078778.1, and AC112721.1) in the ceRNA network, two DEmiRNAs (hsa-mir-141 and hsa-mir-145) and six DEmRNAs (ZEB1, TMEM100, MAP1B, DUSP2, JUN, and AIFM3) were identified to be closely associated with BUC pathogenesis. Two independent prognostic lncRNAs (ADAMTS9-AS1 and ADAMTS9-AS2) and four independent prognostic mRNAs (DUSP2, JUN, MAP1B, and TMEM100) were validated using GEPIA. ADAMTS9-AS1 and ADAMTS9-AS2 interacting with ZEB1 and TMEM100 may play significant roles in BUC development. Key hub DEmRNAs and their relevant pathways may play central or significant roles in BUC tumorigenesis and progression.
Availability of data and materials
The authors declare that the data supporting the findings of this study are available within the article. The R script, which was used to generate figures and reproduce key findings in this study, was stored as Additional file 6.
bladder urothelial carcinoma
long non-coding RNA
The Cancer Genome Atlas
microRNA response element
competing endogenous RNA
differentially expressed mRNAs
differentially expressed lncRNAs
false discovery rate
Gene Expression Profiling Interactive Analysis
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Martinez Rodriguez RH, Buisan Rueda O, Ibarz L. Bladder cancer: present and future. Med Clin (Barc). 2017;149(10):449–55.
Zhao F, Lin T, He W, Han J, Zhu D, Hu K, Li W, Zheng Z, Huang J, Xie W. Knockdown of a novel lincRNA AATBC suppresses proliferation and induces apoptosis in bladder cancer. Oncotarget. 2015;6(2):1064–78.
Huang M, Zhong Z, Lv M, Shu J, Tian Q, Chen J. Comprehensive analysis of differentially expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in bladder carcinoma. Oncotarget. 2016;7(30):47186–200.
Anghel RM, Gales LN, Trifanescu OG. Outcome of urinary bladder cancer after combined therapies. J Med Life. 2016;9(1):95–100.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62.
Yang X, Meng T. Long noncoding RNA in preeclampsia: transcriptional noise or innovative indicators? Biomed Res Int. 2019;2019:5437621.
Bhan A, Soleimani M, Mandal SS. Long noncoding RNA and cancer: a new paradigm. Cancer Res. 2017;77(15):3965–81.
Qin Z, Wang Y, Tang J, Zhang L, Li R, Xue J, Han P, Wang W, Qin C, Xing Q, et al. High LINC01605 expression predicts poor prognosis and promotes tumor progression via upregulation of MMP9 in bladder cancer. Biosci Rep. 2018. https://doi.org/10.1042/BSR20180562.
Li Z, Hong S, Liu Z. LncRNA LINC00641 predicts prognosis and inhibits bladder cancer progression through miR-197-3p/KLF10/PTEN/PI3K/AKT cascade. Biochem Biophys Res Commun. 2018;503:1825–9.
Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17(5):272–83.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Yu C, Longfei L, Long W, Feng Z, Chen J, Chao L, Peihua L, Xiongbing Z, Hequn C. LncRNA PVT1 regulates VEGFC through inhibiting miR-128 in bladder cancer cells. J Cell Physiol. 2018;234:1346–53.
Jiao D, Li Z, Zhu M, Wang Y, Wu G, Han X. LncRNA MALAT1 promotes tumor growth and metastasis by targeting miR-124/foxq1 in bladder transitional cell carcinoma (BTCC). Am J Cancer Res. 2018;8(4):748–60.
Wang Y, Zhang R, Cheng G, Xu R, Han X. Long non-coding RNA HOXA-AS2 promotes migration and invasion by acting as a ceRNA of miR-520c-3p in osteosarcoma cells. Cell Cycle. 2018;37:6316–26.
Huang X, Xie X, Liu P, Yang L, Chen B, Song C, Tang H, Xie X. Adam12 and lnc015192 act as ceRNAs in breast cancer by regulating miR-34a. Oncogene. 2018;17:1637–48.
Xiong DD, Li ZY, Liang L, He RQ, Ma FC, Luo DZ, Hu XH, Chen G. The LncRNA NEAT1 accelerates lung adenocarcinoma deterioration and binds to Mir-193a-3p as a competitive endogenous RNA. Cell Physiol Biochem. 2018;48(3):905–18.
Kouhsar M, Azimzadeh Jamalkandi S, Moeini A, Masoudi-Nejad A. Detection of novel biomarkers for early detection of Non-Muscle-Invasive Bladder Cancer using Competing Endogenous RNA network analysis. Sci Rep. 2019;9(1):8434.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc: Ser B (Methodol). 1995;57(1):289–300.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Nohata N, Abba MC, Gutkind JS. Unraveling the oral cancer lncRNAome: identification of novel lncRNAs associated with malignant progression and HPV infection. Oral Oncol. 2016;59:58–66.
Feng L, Houck JR, Lohavanichbutr P, Chen C. Transcriptome analysis reveals differentially expressed lncRNAs between oral squamous cell carcinoma and healthy oral mucosa. Oncotarget. 2017;8(19):31521–31.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.
Mohamadkhani A. Long noncoding RNAs in interaction With RNA binding proteins in hepatocellular carcinoma. Hepat Mon. 2014;14(5):e18794.
Wang H, Fu Z, Dai C, Cao J, Liu X, Xu J, Lv M, Gu Y, Zhang J, Hua X, et al. LncRNAs expression profiling in normal ovary, benign ovarian cyst and malignant epithelial ovarian cancer. Sci Rep. 2016;6:38983.
Yao J, Zhou B, Zhang J, Geng P, Liu K, Zhu Y, Zhu W. A new tumor suppressor LncRNA ADAMTS9-AS2 is regulated by DNMT1 and inhibits migration of glioma cells. Tumour Biol. 2014;35(8):7935–44.
Basu A, Domb AJ. Recent advances in polyanhydride based biomaterials. Adv Mater. 2018;30:e1706815.
Wang XL, Xie HY, Zhu CD, Zhu XF, Cao GX, Chen XH, Xu HF. Increased miR-141 expression is associated with diagnosis and favorable prognosis of patients with bladder cancer. Tumour Biol. 2015;36(2):877–83.
Lynch SM, O’Neill KM, McKenna MM, Walsh CP, McKenna DJ. Regulation of miR-200c and miR-141 by methylation in prostate cancer. Prostate. 2016;76(13):1146–59.
Ding L, Yu LL, Han N, Zhang BT. miR-141 promotes colon cancer cell proliferation by inhibiting MAP2K4. Oncol Lett. 2017;13(3):1665–71.
Zhang G, Zhang W, Li B, Stringer-Reasor E, Chu C, Sun L, Bae S, Chen D, Wei S, Jiao K, et al. MicroRNA-200c and microRNA-141 are regulated by a FOXP3-KAT2B axis and associated with tumor metastasis in breast cancer. Breast Cancer Res. 2017;19(1):73.
Huang M, Wu L, Qin Y, Li Z, Luo S, Qin H, Yang Y, Chen J. Anti-proliferative role and prognostic implication of miR-141 in gastric cancer. Am J Transl Res. 2016;8(8):3549–57.
Ren D, Wang M, Guo W, Huang S, Wang Z, Zhao X, Du H, Song L, Peng X. Double-negative feedback loop between ZEB2 and miR-145 regulates epithelial–mesenchymal transition and stem cell properties in prostate cancer cells. Cell Tissue Res. 2014;358(3):763–78.
Tan J, Qiu K, Li M, Liang Y. Double-negative feedback loop between long non-coding RNA TUG1 and miR-145 promotes epithelial to mesenchymal transition and radioresistance in human bladder cancer cells. FEBS Lett. 2015;589(20 Pt B):3175–81.
Mahdavinezhad A, Yadegarazari R, Mousavi-Bahar SH, Poorolajal J, Jafari M, Amirzargar MA, Effatpanah H, Saidijam M. Evaluation of zinc finger E-box binding homeobox 1 and transforming growth factor-beta2 expression in bladder cancer tissue in comparison with healthy adjacent tissue. Investig Clin Urol. 2017;58(2):140–5.
Li B, Mao X, Wang H, Su G, Mo C, Cao K, Qiu S. Vasculogenic mimicry in bladder cancer and its association with the aberrant expression of ZEB1. Oncol Lett. 2018;15(4):5193–200.
Liu Q, Qiao L, Liang N, Xie J, Zhang J, Deng G, Luo H, Zhang J. The relationship between vasculogenic mimicry and epithelial–mesenchymal transitions. J Cell Mol Med. 2016;20(9):1761–9.
Wang H, Lin H, Pan J, Mo C, Zhang F, Huang B, Wang Z, Chen X, Zhuang J, Wang D, et al. Vasculogenic mimicry in prostate cancer: the roles of EphA2 and PI3K. J Cancer. 2016;7(9):1114–24.
Ou D, Yang H, Hua D, Xiao S, Yang L. Novel roles of TMEM100: inhibition metastasis and proliferation of hepatocellular carcinoma. Oncotarget. 2015;6(19):17379–90.
Palenzuela R, Gutierrez Y, Draffin JE, Lario A, Benoist M, Esteban JA. MAP1B light chain modulates synaptic transmission via AMPA receptor intracellular trapping. J Neurosci. 2017;37(41):9945–63.
Auburger G, Gispert S, Torres-Odio S, Jendrach M, Brehm N, Canet-Pons J, Key J, Sen NE. SerThr-PhosphoProteome of brain from aged PINK1-KO + A53T-SNCA mice reveals pT1928-MAP1B and pS3781-ANK2 deficits, as hub between autophagy and synapse changes. Int J Mol Sci. 2019;20:13.
Tessema M, Yingling CM, Picchi MA, Wu G, Liu Y, Weissfeld JL, Siegfried JM, Tesfaigzi Y, Belinsky SA. Epigenetic repression of CCDC37 and MAP1B links chronic obstructive pulmonary disease to lung cancer. J Thor Oncol. 2015;10(8):1181–8.
Yin H, He W, Li Y, Xu N, Zhu X, Lin Y, Gou X. Loss of DUSP2 predicts a poor prognosis in patients with bladder cancer. Hum Pathol. 2019;85:152–61.
Chien MH, Lee TH, Lee WJ, Yeh YH, Li TK, Wang PC, Chen JJ, Chow JM, Lin YW, Hsiao M, et al. Trichodermin induces c-Jun N-terminal kinase-dependent apoptosis caused by mitotic arrest and DNA damage in human p53-mutated pancreatic cancer cells and xenografts. Cancer Lett. 2017;388:249–61.
Sharmila R, Sindhu G. Evaluate the antigenotoxicity and anticancer role of beta-sitosterol by determining oxidative DNA damage and the expression of phosphorylated mitogen-activated protein kinases’, C-fos, C-jun, and endothelial growth factor receptor. Pharmacognosy Mag. 2017;13(49):95–101.
Gao GY, Ma J, Lu P, Jiang X, Chang C. Ophiopogonin B induces the autophagy and apoptosis of colon cancer cells by activating JNK/c-Jun signaling pathway. Biomed Pharmacother. 2018;108:1208–15.
He C, Zhang H, Wang B, He J, Ge G. SDF-1/CXCR4 axis promotes the growth and sphere formation of hypoxic breast cancer SP cells by c-Jun/ABCG2 pathway. Biochem Biophys Res Commun. 2018;505(2):593–9.
Zheng A, Zhang L, Song X, Wang Y, Wei M, Jin F. Clinical implications of a novel prognostic factor AIFM3 in breast cancer patients. BMC Cancer. 2019;19(1):451.
Liu X, Wang J, Sun G. Identification of key genes and pathways in renal cell carcinoma through expression profiling data. Kidney Blood Press Res. 2015;40(3):288–97.
Wiley SZ, Sriram K, Liang W, Chang SE, French R, McCann T, Sicklick J, Nishihara H, Lowy AM, Insel PA. GPR68, a proton-sensing GPCR, mediates interaction of cancer-associated fibroblasts and cancer cells. FASEB J. 2018;32(3):1170–83.
Choudhary P, Gutteridge A, Impey E, Storer RI, Owen RM, Whiting PJ, Bictash M, Benn CL. Targeting the cAMP and transforming growth factor-beta pathway increases proliferation to promote re-epithelialization of human stem cell-derived retinal pigment epithelium. Stem Cells Transl Med. 2016;5(7):925–37.
Abudoureyimu A, Muhemaitibake A. Arsenic trioxide regulates gastric cancer cell apoptosis by mediating cAMP. Eur Rev Med Pharmacol Sci. 2017;21(3):612–7.
Zhang HJ, Tao J, Sheng L, Hu X, Rong RM, Xu M, Zhu TY. Twist2 promotes kidney cancer cell proliferation and invasion by regulating ITGA6 and CD44 expression in the ECM-receptor interaction pathway. Onco Targets Ther. 2016;9:1801–12.
Munoz JJ, Drigo SA, Barros-Filho MC, Marchi FA, Scapulatempo-Neto C, Pessoa GS, Guimaraes GC, Trindade Filho JC, Lopes A, Arruda MA, et al. Down-Regulation of SLC8A1 as a putative apoptosis evasion mechanism by modulation of calcium levels in penile carcinoma. J Urol. 2015;194(1):245–51.
Freire PP, Fernandez GJ, Cury SS, de Moraes D, Oliveira JS, de Oliveira G, Dal-Pai-Silva M, Dos Reis PP, Carvalho RF. The pathway to cancer cachexia: MicroRNA-regulated networks in muscle wasting based on integrative meta-analysis. Int J Mol Sci. 2019;20(8):e1962.
Feng YY, Zeng DZ, Tong YN, Lu XX, Dun GD, Tang B, Zhang ZJ, Ye XL, Li Q, Xie JP, et al. Alteration of microRNA-4474/4717 expression and CREB-binding protein in human colorectal cancer tissues infected with Fusobacterium nucleatum. PLoS ONE. 2019;14(4):e0215088.
Xiong DD, Xu WQ, He RQ, Dang YW, Chen G, Luo DZ. In silico analysis identified miRNAbased therapeutic agents against glioblastoma multiforme. Oncol Rep. 2019;41(4):2194–208.
Yu D, Holm R, Goscinski MA, Trope CG, Nesland JM, Suo Z. Prognostic and clinicopathological significance of Cacna2d1 expression in epithelial ovarian cancers: a retrospective study. Am J Cancer Res. 2016;6(9):2088–97.
Johansson P, Aoude LG, Wadt K, Glasson WJ, Warrier SK, Hewitt AW, Kiilgaard JF, Heegaard S, Isaacs T, Franchina M, et al. Deep sequencing of uveal melanoma identifies a recurrent mutation in PLCB4. Oncotarget. 2016;7(4):4624–31.
Wu B, Xie J, Du Z, Wu J, Zhang P, Xu L, Li E. PPI network analysis of mRNA expression profile of ezrin knockdown in esophageal squamous cell carcinoma. Biomed Res Int. 2014;2014:651954.
Qu Y, Zhang X, Wu R. Knockdown of NF-kappaB p65 subunit expression suppresses growth of nude mouse lung tumor cell xenografts by activation of Bax apoptotic pathway. Neoplasma. 2015;62(1):34–40.
Tang R, Li FX, Shao WF, Wen QS, Yu XR, Xiong JB. Protein-protein interaction between ezrin and p65 in human breast cancer cells. Genetics Mol Res. 2016. https://doi.org/10.4238/gmr.15028334.
Ou Z, Wang Y, Chen J, Tao L, Zuo L, Sahasrabudhe D, Joseph J, Wang L, Yeh S. Estrogen receptor beta promotes bladder cancer growth and invasion via alteration of miR-92a/DAB2IP signals. Exp Mol Med. 2018;50(11):152.
Munaron L. Systems biology of ion channels and transporters in tumor angiogenesis: an omics view. Biochim Biophys Acta. 2015;1848(10 Pt B):2647–56.
Hoffmann EK, Lambert IH. Ion channels and transporters in the development of drug resistance in cancer cells. Philos Trans R Soc Lond B Biol Sci. 2014;369(1638):20130109.
Pardo LA, Stuhmer W. The roles of K(+) channels in cancer. Nat Rev Cancer. 2014;14(1):39–48.
Litan A, Langhans SA. Cancer as a channelopathy: ion channels and pumps in tumor development and progression. Front Cell Neurosci. 2015;9:86.
Djamgoz MB, Onkal R. Persistent current blockers of voltage-gated sodium channels: a clinical opportunity for controlling metastatic disease. Recent Pat Anti-Cancer Drug Discovery. 2013;8(1):66–84.
Zhu Q, Tang M, Wang X. The expression of asparaginyl endopeptidase promotes growth potential in epithelial ovarian cancer. Cancer Biol Ther. 2017;18(4):222–8.
Patel SA, Gooderham NJ. IL6 mediates immune and colorectal cancer cell cross-talk via miR-21 and miR-29b. Mol Cancer Res. 2015;13(11):1502–8.
Benezeder T, Tiran V, Treitler AAN, Suppan C, Rossmann C, Stoeger H, Cote RJ, Datar RH, Balic M, Dandachi N. Multigene methylation analysis of enriched circulating tumor cells associates with poor progression-free survival in metastatic breast cancer patients. Oncotarget. 2017;8(54):92483–96.
Zhang M, Luo SC. Gene expression profiling of epithelial ovarian cancer reveals key genes and pathways associated with chemotherapy resistance. Genetics Mol Res. 2016. https://doi.org/10.4238/gmr.15017496.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 2: Figure S2. The network of significant-top 15 KEGG pathways enriched in the DEmRNAs. Red nodes represent increased expression levels, while blue nodes represent decreased expression levels. Triangle nodes represent DEmRNAs; Green ellipse nodes represent enrichment pathways. Gray edges indicate mRNAs involved in the pathway.
Additional file 4: Figure S4. (A) KEEG pathways enriched using the hub DEmRNAs; (B) Network between hub DEmRNAs and KEGG pathways enriched by hub DEmRNAs. The red nodes represent increased expression levels, while the blue nodes represent decreased expression levels. Triangle nodes represent DEmRNAs; Green ellipse nodes represent enrichment pathways. Gray edges indicate mRNAs involved in the pathway.
About this article
Cite this article
Wang, J., Zhang, C., Wu, Y. et al. Identification and analysis of long non-coding RNA related miRNA sponge regulatory network in bladder urothelial carcinoma. Cancer Cell Int 19, 327 (2019) doi:10.1186/s12935-019-1052-2