Expression Profiles of HOXC6 Predict the Survival of Glioblastoma Patients and Correlate with Cell Cycle

The goal of this study was to investigate the homeobox (HOX) gene expression status and its prognostic value in glioblastoma multiforme (GBM) and to uncover the biological processes related to its expression. The prognostic value of HOX genes in GBM was systematically investigated by a genome-wide analysis of HOX gene expression profiles in GBM patient samples in The Cancer Genome Atlas (TCGA) project (microarray dataset) and validation datasets. Using the differentially expressed gene (DEG) analysis and a Cox regression model, we discovered that the HOXC6 could stratify patients into significantly different survival (p = 0.0012, log-rank test) groups in the training cohort. TCGA RNA-seq and GSE16011 datasets were used for validation. Multivariate Cox and stratification analysis indicated that HOXC6 was an independent prognostic factor after adjusting for other clinical covariates. Bioinformatic analysis suggested that the HOXC6 might be involved in the cell cycle-related biological processes and pathways that are well established in the context of glioblastoma tumorigenesis. We further explored the bioinformatic implications by gene set enrichment analysis (GSEA). Tumor cell biology experiments verified the role of HOXC6 in proliferation and cell cycle progression. In conclusion, HOXC6 might be a candidate biomarker gene for individual treatment optimization of glioblastoma. HOXC6 expression has a significant prognostic value and is related to the cell cycle process in glioblastoma.