Abstract
Hypoxia is a universal feature in the tumor microenvironment (TME). Nonetheless, the heterogeneous hypoxia patterns of TME have still not been elucidated in hepatocellular carcinoma (HCC). Using consensus clustering algorithm and public datasets, we identified heterogeneous hypoxia subtypes. We also revealed the specific biological and clinical characteristics via bioinformatic methods. The principal component analysis algorithm was employed to develop a hypoxia-associated risk score (HARS). We identified the two hypoxia subtypes: low hypoxia pattern (C1) and high hypoxia pattern (C2). C1 was less sensitive to immunotherapy compared to C2, consistent with the lack of immune cells and immune checkpoints (ICPs) in C1, whereas C2 was the opposite. C2 displayed worse prognosis and higher sensitivity to obatoclax relative to C1, while C1 was more sensitive to sorafenib. The two subtypes also demonstrated subtype-specific genomic variations including mutation, copy number alteration, and methylation. Moreover, we developed and validated a risk signature: HARS, which had excellent performance for predicting prognosis and immunotherapy. We revealed two hypoxia subtypes with distinct biological and clinical characteristics in HCC, which enhanced the understanding of hypoxia pattern. The risk signature was a promising biomarker for predicting prognosis and immunotherapy.