Abstract
The tumor immune microenvironment of hepatocellular carcinoma (HCC) is complex, yet the causal relationship between immune cell subpopulations and HCC risk remains incompletely elucidated. This study aims to systematically evaluate the causal association between immune cell subpopulations and HCC using Mendelian randomization (MR) analysis, and to validate the biological mechanisms underlying these associations through multi-omics data. Bidirectional two-sample MR analysis was performed to examine causal relationships between 731 immune cell subpopulations and HCC. Inverse-variance weighting (IVW) served as the primary analysis method, with robustness validation through Bayesian weighted MR (BWMR) and machine learning algorithms. Therefore, for significantly associated immune subpopulations, independent analyses of gene expression, prognosis, and tumor immune microenvironment were conducted using HCC data from the Cancer Genome Atlas (TCGA) LIHC cohort. MR analysis and validation identified 21 immune cell subpopulations with significant causal associations to HCC risk. Among these, 12 were identified as risk factors, and 9 as protective factors. Validation in the TCGA cohort revealed that risk-associated immune subpopulations were predominantly enriched for markers of T cell exhaustion and immunosuppressive microenvironments, whereas protective subpopulations likely represented a distinct regulatory B cell subset whose function was associated with the anti-inflammatory factor interleukin-10. This study genetically confirms that specific immune cell functional subpopulations constitute causal risk factors for HCC. These subpopulations exert their effects by shaping distinct tumor immune microenvironments. These findings provide novel mechanisms for understanding the immunopathogenesis of HCC and identify potential targets for developing novel immune intervention strategies.