Identification and validation of a gap junction protein related signature for predicting the prognosis of renal clear cell carcinoma.

BACKGROUND

Gap junction proteins (GJPs) are a class of channel proteins that are closely related to cell communication and tumor development. The objective of this study was to screen out GJPs related prognostic signatures (GRPS) associated with clear cell renal cell carcinoma (ccRCC).

MATERIALS AND METHODS

GJPs microarray data for ccRCC patients were obtained from The Gene Expression Omnibus (GEO) database, along with RNA sequencing data for tumor and paired normal tissues from The Cancer Genome Atlas (TCGA) database. In the TCGA database, least absolute shrinkage and selection Operator (LASSO) and Cox regression models were used to identify GJPs with independent prognostic effects as GRPS in ccRCC patients. According to the GRPS expression and regression coefficient from the multivariate Cox regression model, the risk score (RS) of each ccRCC patient was calculated, to construct the RS prognostic model to predict survival. Overall survival (OS) and progression-free survival (PFS) analyses; gene pan-cancer analysis; single gene survival analysis; gene joint effect analysis; functional enrichment analysis; tumor microenvironment (TME) analysis; tumor mutational burden (TMB) analysis; and drug sensitivity analysis were used to explore the biological function, mechanism of action and clinical significance of GRPS in ccRCC. Further verification of the genetic signature was performed with data from the GEO database. Finally, the cytofunctional experiments were used to verify the biological significance of GRPS associated GJPs in ccRCC cell lines.

RESULTS

GJA5 and GJB1, which are GRPS markers of ccRCC patients, were identified through LASSO and Cox regression models. Low expression of GJA5 and GJB1 is associated with poor patient prognosis. Patients with high-RS had significantly shorter OS and PFS than patients with low-RS (< 0.001). The risk of death for individuals with high-RS was 1.695 times greater than that for those with low-RS (HR = 1.695, 95%CI= 1.439-1.996, < 0.001). Receiver Operating Characteristic (ROC) curve showed the great predictive power of the RS prognostic model for the survival rate of patients. The area under curve (AUC) values for predicting 1-year, 3-year and 5-year survival rates were 0.740, 0.781 and 0.771, respectively. The clinical column chart was also reliable for predicting the survival rate of patients, with AUC values of 0.859, 0.846 and 0.796 for predicting 1-year, 3-year and 5-year survival, respectively. The GRPS was associated with immune cell infiltration, the TME, the TMB, and sensitivity to chemotherapy drugs. Further experiments showed that knockdown of GJA5 or GJB1 could promote the proliferation, migration and epithelial-mesenchymal transition (EMT) and inhibit apoptosis of ccRCC cells.

CONCLUSION

GJA5 and GJB1 could be potential biological markers for predicting survival in patients with ccRCC.