Abstract
Infection by high-risk human papillomavirus is known to exacerbate cervical cancer development. The host immune response is crucial in disease regression. Large-scale genetic association studies for cervical cancer have identified few susceptibility variants, mainly at the human leukocyte antigen locus on chromosome 6. We hypothesized that the host immune response modifies cervical cancer risk and performed three genome-wide association analyses for HPV16, HPV18 and HPV16/18 seropositivity in 7814, 7924, and 7924 samples from the UK Biobank, followed by validation genotyping in the German Cervigen case-control series of cervical cancer and dysplasia. In GWAS analyses, we identified two loci associated with HPV16 seropositivity (6p21.32 and 15q26.2), two loci associated with HPV18 seropositivity (5q31.2 and 14q24.3), and one locus for HPV16 and/or HPV18 seropositivity (at 6p21.32). MAGMA gene-based analysis identified HLA-DQA1 and HLA-DQB1 as genome-wide significant (GWS) genes. In validation genotyping, the genome-wide significant lead variant at 6p21.32, rs9272293 associated with overall cervical disease (OR = 0.86, p = 0.004, 95% CI = 0.78-0.95, n = 3710) and HPV16 positive invasive cancer (OR = 0.73, p = 0.005, 95% CI = 0.59-0.91, n = 1431). This variant was found to be a robust eQTL for HLA-DRB1, HLA-DQB1-AS1, C4B, HLA-DRB5, HLA-DRB6, HLA-DQB1, and HLA-DPB1 in a series of cervical epithelial tissue samples. We additionally genotyped twenty-four HPV seropositivity variants below the GWS threshold out of which eleven variants were found to be associated with cervical disease in our cohort, suggesting that further seropositivity variants may determine cervical disease outcome. Our study identifies novel genomic risk loci that associate with HPV type-specific cervical cancer and dysplasia risk and provides evidence for candidate genes at one of the risk loci.