Dynamic expression changes between non-muscle-invasive bladder cancer and muscle-invasive bladder cancer

Despite elaborate characterization of the risk factors, bladder cancer is still a major epidemiological problem whose incidence continues to rise each year. We aim to investigate the dynamic expression changes between non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC). The gene expression profile GSE13507 was obtained from the Gene Expression Omnibus, and the R package was used to identify gene expression signatures (GESs) between NMIBC and MIBC. Gene ontology enrichment analysis was performed for GES function analysis. We used miRTarBase and TargetScan to identify the differentially regulated microRNAs, and TfactS to identify transcription factors between NMIBC and MIBC. Bionet was used to identify the differentially expressed subnetwork. A total of 802 upregulated NMIBC GESs and 668 downregulated MIBC GESs were identified. Functional enrichment analysis revealed that the MIBC GESs were majorly involved in cell cycle and inflammatory response. miR-29c and miR-9 were regarded as key microRNAs in MIBC. SMAD3 in MIBC and SMAD5 and SMAD7 in NMIBC were potential activated transcription factors. In addition, a subnetwork that was considered to capture the differences between MIBC and NMIBC was identified, of which GRB2 and UBC were the hub nodes. Some key microRNAs, activated transcription factors and hub nodes have been identified in this study, which may be used as potential biomarkers or targets for the diagnosis, treatment and detection of bladder cancer at different stages.