INTRODUCTION

Several observational or retrospective studies have previously been conducted to explore the possible association between lung cancer and human papillomavirus (HPV) infection. However, there may be inconsistencies in the data and conclusions due to differences in study design and HPV testing methods. There are currently no studies that provide conclusive evidence to support the involvement of HPV in the occurrence and development of lung cancer. Therefore, the relationship between HPV and lung cancer remains controversial and uncertain. This study aimed to explore whether HPV infection is causally related to lung cancer risk by systematically performing a two-way Two-Sample Mendelian Randomization (TSMR) analysis.

METHODS

In the International Lung Cancer Consortium (ILCCO) genome-wide association study dataset, we included 11,348 lung cancer (LUCA) cases, including 3275 squamous cell carcinoma (LUSC) cases, 3442 adenocarcinoma (LUAD) cases, and 15,861 cases of control. Using genetic variants associated with the HPV E7 protein as instrumental variables, we summarized statistics associated with HPV infection in the MRC IEU OpenGWAS database, which included the HPV-16 E7 protein and the HPV-18 E7 protein. Two-sample Mendelian randomization (MR) results are expressed as odds ratios (OR) and 95% confidence intervals (CI).

RESULTS

Based on a comprehensive analysis of genome-wide association study (GWAS) data from public databases, we mainly used inverse-variance weighted (IVW) to estimate causal relationships, while using MR-Egger, weighted median, simple mode, and weighted mode, and other four methods as supplements. Two-sample MR Analysis revealed no causal relationship between exposure factors (HPV-16 E7 protein and HPV-18 E7 protein) and outcome factors (lung cancer (LUCA) and its subtypes squamous cell carcinoma (LUSC) and adenocarcinoma (LUAD)) in forward MR Analysis using the IVW approach.HPV-16 E7 protein and LUCA and its subtypes LUSC and LUAD by IVW method results: [OR] = 1.002; 95% [CI]: 0.961 - 1.045; p = 0.920; [OR] = 1.023; 95% [CI]: 0.966 - 1.084; p = 0.438; [OR] = 0.994; 95% [CI]: 0.927 - 1.066; p = 0.872); HPV-18 E7 protein and LUCA and its subtypes LUSC and LUAD by IVW method results: [OR] = 0.965; 95% [CI]: 0.914 - 1.019; p = 0.197; [OR] = 0.933; 95% [CI]: 0.834 - 1.043; p = 0.222; [OR] = 1.028; 95% [CI]: 0.945 - 1.118; p = 0.524. It was observed through reverse MR that LUCA and its subtypes LUSC and LUAD were used as exposure factors, and HPV infection (HPV-16 E7 protein and HPV-18 E7 protein) was used as the outcome factors, the results of the IVW method are also invalid.LUCA and HPV-16 E7 protein and HPV-18 E7 protein by IVW method results: [OR] = 1.036; 95% [CI]: 0.761 - 1.411; p = 0.82; [OR] = 1.318; 95% [CI]: 0.949 - 1.830; p = 0.099; LUSC and HPV-16 E7 protein and HPV-18 E7 protein by IVW method results: [OR] = 1.123; 95% [CI]0.847 - 1.489; p = 0.421; [OR] = 0.931; 95% [CI]: 0.660 - 1.313; p = 0.682; LUAD and HPV-16 E7 protein and HPV-18 E7 protein by IVW method results: [OR] = 1.182; 95% [CI] 0.983 - 1.421; p = 0.075; [OR] = 1.017; 95% [CI]: 0.817 - 1.267; p = 0.877.Our results indicate that there is no causal relationship between genetically predicted HPV infection and LUCA and its subtypes LUSC and LUAD. In addition, in the reverse MR analysis, we did not observe a significant causal relationship between LUCA and its subtypes LUSC and LUAD on HPV infection.

CONCLUSIONS

Our findings do not support a genetic association between HPV infection and lung cancer.