Uncovering the potential functions of lymph node metastasis-associated aberrant methylation differentially expressed genes and their association with the immune infiltration and prognosis in bladder urothelial carcinoma

Abstract

Background Bladder urothelial carcinoma (BLCA) is a malignant tumor of the urinary system. This study aimed to explore the potential role of lymph node metastasis-associated aberrant methylation differentially expressed genes (DEGs) in BLCA. Methods CHAMP and limma packages were used to identify lymph node metastasis-associated aberrant methylation DEGs. Univariate Cox analysis and Lasso analysis were performed to identify the signature genes, and multivariate Cox analysis was used to construct the risk score. Subsequently, the molecular characteristics of the signature genes and the relationship between risk score and prognosis, clinical characteristics and immune cell infiltration were analyzed. The signature gene AKAP7 was selected for functional verification. Results A novel risk score model was constructed based on 12 signature genes. The risk score had a good ability to predict overall survival (OS). The nomogram constructed based on age, N stage and risk score had a higher value in predicting the prognosis of patients. It was also found that stromal activation in TIME may inhibit the antitumor effects of immune cells. Functional enrichment analysis revealed that ECM receptor interaction and focal adhesion were two important pathways involved in the regulation of BLCA. Immunohistochemistry showed that AKAP7 may be associated with the occurrence, clinical stages and grades, and lymph node metastasis of BLCA. In vitro cell experiments showed that the migration and invasion ability of EJ cells was significantly inhibited after AKAP7 overexpression, while the migration and invasion ability of T24 cells was significantly promoted after AKAP7 knockdown. Conclusion The risk score model based on lymph node metastasis-associated aberrant methylation DEGs has a good ability to predict OS and is an independent prognostic factor for BLCA. It was also found that stromal activation in TIME may inhibit the antitumor effects of immune cells. This implicates aberrant methylation modifications as an important factor contributing to the heterogeneity and complexity of individual tumor microenvironments. Functional enrichment analysis revealed that ECM receptor interaction and focal adhesion were two important pathways involved in the regulation of BLCA, which contributed to the exploration of the pathological mechanism of BLCA. In addition, immunohistochemistry showed that AKAP7 may be associated with the occurrence, progression and lymph node metastasis of BLCA. In vitro cell experiments showed that AKAP7 could also inhibit the migration and invasion of cancer cells.