Abstract
INTRODUCTION: More and more studies have focused on the associations between human papillomavirus (HPV) infection and pan-cancers. However, current evidence is largely based on retrospective studies, which are susceptible to confounding factors and do not enable the establishment of causal relationships. METHODS: A bidirectional two-sample Mendelian randomization (MR) design was employed to thoroughly evaluate the causal relationships between HPV and 12 site-specific cancers except cervical cancer. Single nucleoside polymers (SNPs) with strong evidence from genome-wide association studies (GWAS) were selected from HPV exposure datasets and used as instrumental variables (IVs) in this study. For the MR analysis results, MR-Egger's intercept P test, MR-PRESSO global test, Cochran's Q test and a leave-one-out test were applied for sensitivity analysis. Using HPVTIMER, we also performed immune infiltration analyses in head and neck squamous cell carcinoma (HNSCC), oropharyngeal squamous cell carcinoma (OPSCC) and vulval squamous cell carcinoma (VSCC) to evaluate the tumor-immune microenvironment. RESULTS: Based on the evidence of MR analysis, our study conclusively identified HPV16 as a risk factor implicated in the development of bladder cancer, colorectal cancer, and breast cancer, while HPV18 was identified as a risk factor for prostate cancer, ovarian cancer, lung cancer and breast cancer. The MR results also showed that HPV16 may be a protective factor for prostate cancer, anal cancer, lung cancer and oropharyngeal cancer, while HPV18 may be a protective factor for vaginal cancer. CONCLUSION: An HPV infection may modulate the immune microenvironment and therefore has a potential inhibitory effect on the development of certain cancers. These conclusions provided new insights into the potential mechanisms of carcinogenesis and needed further research for validation.