The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10^-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRS_ZEBRA and GRS_MCV were associated in opposite directions with risk of IDH wild-type gliomas (OR_ZEBRA = 0.91, p = 0.0099/OR_MCV = 1.11, p = 0.0054). GRS_EBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1^∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10^-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.