Comprehensive analysis of the lncRNA-miRNA-mRNA regulatory network for bladder cancer.

Background: Long non-coding RNAs (lncRNAs) are essential regulators for various human cancers. However, these lncRNAs need to be further classified for cancer. In the present study, we identified novel competing endogenous RNA (ceRNA) network for bladder cancer (BC) and explored the gene functions of the ceRNA regulatory network.

Methods: Differential gene expression analysis were performed on The Cancer Genome Atlas Urothelial Bladder Carcinoma (TCGA-BLCA) datasets to identify differentially expressed messenger RNAs (mRNAs), lncRNAs, and microRNAs (miRNAs). Based on the competing endogenous RNA (ceRNA) hypothesis, a lncRNA-miRNA-mRNA network was constructed using the StarBase database and visualization by Cytoscape software. Functional enrichment analyses of Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were performed via R package ClusterProfiler. The protein-protein interaction network was constructed by STRING database and visualization by Cytoscape. Finally, we used CIBERSORT and the TIMER database to analyze the immune infiltrations for BC.

Results: The regulatory network was constructed via TCGA BLCA cohort. The differential expressions of lncRNA, miRNA, and mRNA were 186, 200, and 2,661, respectively. There were 106 lncRNA, miRNA, and mRNA included in the ceRNA network. In this network, Calcium Voltage-gated Channel Auxiliary Subunit Alpha2delta1 (CACNA2D1, P<0.001), domain containing engulfment adaptor1 (GULP1, P=0.001), latent transforming growth factor beta binding protein 1 (LTBP1, P=0.006), myosin light chain kinase (MYLK, P=0.001), serpin family E member 2 (SERPINE2, P=0.002), spectrin beta non-erythrocytic 2 (SPTBN2, P=0.047), and hsa-miR-590-3p (P<0.001) significantly affected the prognosis of BC patients. Functional enrichment analyses showed that the biological functions included negative regulation of protein phosphorylation, cell morphogenesis, and sensory organ morphogenesis. Important cancer pathways of KEGG included parathyroid hormone synthesis secretion action, the notch signaling pathway, MAPK signaling pathway, the Rap1 signaling pathway, signaling pathways regulating the pluripotency of stem cells, and the transforming growth factor-β signaling pathway. Our findings demonstrated that the ceRNA network has important biological functions and a significant influence on the prognosis of BC.

Conclusions: The lncRNA-miRNA-mRNA network constructed in the present study could provide useful insight into the underlying tumorigenesis of BC, and can determine new molecular biomarkers for the diagnosis and therapeutical treatment of BC.