Identification of Hub Genes Associated with Progression and Prognosis of Bladder Cancer by Integrated Bioinformatics Analysis.

BACKGROUND

Bladder cancer (BLCA) is an extremely common carcinoma of the urinary system that has a high incidence of relapse. Although intensive studies have investigated its pathology in the past decades, there are significant knowledge gaps regarding the characterization of the molecular processes underlying the progression of disease and consequently its prognosis. The purpose of current research was to identify significant genes that could serve as prognostic and progression biomarkers.

METHODS

Gene expression profiles were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Differential gene expression analysis (DGEA) and weighted gene co-expression network analysis (WGCNA) were conducted to recognize differential co-expression genes. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to explore gene function. Moreover, protein-protein interactions (PPI) network, overall survival (OS) and disease-free survival (DFS) were used to identify survival-related hub genes. Additionally, associations between these gene's expression and clinical parameters were determined. Finally, the Human Protein Atlas (HPA) database and qRT-PCR were used to validate gene's expression.

RESULTS

About 124 differential co-expression genes were identified. These genes were mainly enriched in muscle system process and muscle contraction (biological process, BP), contractile fiber, myofibril, sarcomere, focal adhesion and cell-substrate junction (cellular component, CC) and actin binding (molecular function, MF) in GO enrichment analysis, while enriched in vascular smooth muscle contraction, focal adhesion, cardiac muscle contraction, hypertrophic cardiomyopathy, dilated cardiomyopathy and regulation of actin cytoskeleton in KEGG analysis. Furthermore, five survival-related hub genes (MYH11, ACTA2, CALD1, TPM1, MYLK) were identified via OS and DFS. In addition, these survival-related gene's expression was correlated with grade, stage and TNM stage. Finally, all survival-related hub genes were determined to be down-regulated in BLCA tissues by qRT-PCR and HPA databases.

CONCLUSIONS

Our current study verified five new key genes in BLCA, which may participate in the prognosis and progression and serve as novel biomarkers of BLCA.