Identification of HOXB9 to predict prognosis of endometrial cancer based on comprehensive bioinformatics analysis

Abstract

Abstract Background The HOXB9 gene, which plays a key role in embryonic development, is also involved in the regulation of various human cancers. However, the potential relationship between HOXB9 and endometrial cancer (EC) has not yet been comprehensively analyzed and fully understood. Methods We used multiple bioinformatics tools to explore the role of HOXB9 in EC. Results The expression of HOXB9 was significantly upregulated in pan-cancer, including EC (P < 0.05). Quantitative real time polymerase chain reaction (qRT-PCR) experiment confirmed the high expression of HOXB9 in EC from clinical samples (P < 0.001). Double validated by Enrichr and Metascape, HOXB9 showed a strong correlation with HOX family, suggesting that HOX family may also involve in the development of EC (P < 0.05). Enrichment analysis revealed HOXB9 is mainly associated with cellular process, developmental process, P53 signaling pathway, etc. At the single-cell level, the clusters of cells ranked were glandular and luminal cells c-24, glandular and luminal cells c-9, endothelial cells c-15, compared with the other cells. At the genetic level, promoter methylation levels of HOXB9 were significantly higher in tumors than in normal tissues. Furthermore, variations of HOXB9 were closely associated with overall survival (OS) and recurrence free survival (RFS) in EC patients (P < 0.05). The agreement between univariate and multivariate Cox regression indicated that the results were more reliable. Stages III and IV, G2 and G3, tumor invasion ≥ 50%, mixed or serous histological type, age > 60 years, and high expression of HOXB9 were risk factors strongly associated with OS in EC patients (P < 0.05). Therefore, six factors were incorporated to construct a nomogram for survival prediction. Finally, we used the Kaplan-Meier (KM) curve, receiver operating characteristic (ROC) curve, and time-dependent ROC to assess predictive power of HOXB9. KM curve showed EC patients overexpressing HOXB9 had a worse OS. AUC of diagnostic ROC was 0.880. AUCs of time-dependent ROC were 0.602, 0.591, and 0.706 for 1-year, 5-year, and 10-year survival probabilities (P < 0.001). Conclusions Our study provids new insights into the diagnosis and prognosis of HOXB9 in EC and constructs a model that can accurately predict the prognosis of EC.