ABSTRACT: Background and Objectives: Bladder cancer (BC) is a complex tumor associated with high recurrence and mortality. To discover key molecular changes in BC, we analyzed next-generation sequencing data of BC and surrounding tissue samples from clinical specimens. Methods. Gene expression profiling datasets of bladder cancer were analyzed online. The Database for Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/) was used to perform Gene Ontology (GO) functional and KEGG pathway enrichment analyses. Molecular Complex Detection (MCODE) in Cytoscape software (Cytoscape_v3.6.1) was applied to identify hub genes. Protein expression and survival data were downloaded from OncoLnc (http://www.oncolnc.org/). Gene expression data were obtained from the ONCOMINE website (https://www.oncomine.org/). Results. We identified 4211 differentially expressed genes (DEGs) by analysis of surrounding tissue vs. cancer tissue (SC analysis) and 410 DEGs by analysis of cancer tissue vs. recurrent tissue cluster (CR analysis). GO function analysis revealed enrichment of DEGs in genes related to the cytoplasm and nucleoplasm for both clusters, and KEGG pathway analysis showed enrichment of DEGs in the PI3K-Akt signaling pathway. We defined the 20 genes with the highest degree of connectivity as the hub genes. Cox regression revealed CCNB1, ESPL1, CENPM, BLM, and ASPM were related to overall survival. The expression levels of CCNB1, ESPL1, CENPM, BLM, and ASPM were 4.795-, 5.028-, 8.691-, 2.083-, and 3.725-fold higher in BC than the levels in normal tissues, respectively. Conclusions. The results suggested that the functions of CCNB1, ESPL1, CENPM, BLM, and ASPM may contribute to BC development and the functions of CCNB1, ESPL1, CENPM, and BLM may also contribute to BC recurrence.