Identification of a hypoxia-related gene prognostic signature in colorectal cancer based on bulk and single-cell RNA-seq

Abstract Colorectal cancer (CRC) is the most common and fatal tumor in the gastrointestinal system. Its incidence and mortality rate have increased in recent years. Hypoxia, a persistent physiological tumor feature, plays a vital role in CRC tumorigenesis, metastasis, and tumor microenvironment (TME). Therefore, we constructed a hypoxia-related gene (HRG) nomogram to predict overall survival (OS) and explored the role of HRGs in the CRC TME. The Cancer Genome Atlas (TCGA) dataset was used as the training set, and two Gene Expression Omnibus datasets (GSE39582 and GSE103479) were used as the testing sets. HRGs were identified using the Gene Set Enrichment Analysis (GSEA) database. An HRG prognostic model was constructed in the training set using the least absolute shrinkage and selection operator regression algorithm and validated in the testing sets. Then, we analyzed tumor-infiltrating cells (TICs) using the cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm. Furthermore, single-cell next-generation RNA sequencing (RNA-seq) was used to investigate HRG expression in different TICs in the GSE139555 dataset. Finally, reverse transcription polymerase chain reactions (RT-PCR) were used to validate HRG mRNA expression in ten pairs of CRC normal and cancer tissue samples. A six HRG prognostic signature was constructed, with a superior OS prediction ability in CRC patients (area under the receiver operating characteristic curve (AUC) at one year: 0.693, AUC at three years: 0.712, and AUC at five years: 0.780). GSEA enrichment analysis identified six pathways enriched in the high-risk group. The TIC analysis indicated that the high-risk group had lower T-cell expression and higher neutrophil expression than the low-risk group. Furthermore, immune-related genes had an inseparable relationship with the HRG prognostic signature. Based on single-cell RNA-seq data, we found elevated hexokinase 1 (HK1) and glucose-6-phosphate isomerase (GPI) gene expression in natural killer (NK) and CD8+ T cells. RT-PCR in ten CRC normal-tumor tissue pairs showed that expression of the signature’s six HRGs varied differently in cancerous and paracancerous tissues. The constructed HRG signature successfully predicted the OS of whole-stage CRC patients. In addition, we showed that the signature’s six HRGs were closely associated with the TME in CRC, where hypoxia inhibits the antitumor function of T cells.