Genome-wide Polygenic Risk Scores Predict Risk of Glioma and Molecular Subtypes.

BACKGROUND: Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to efficiently capture genetic risk using available data. METHODS: We applied a method based on continuous shrinkage priors (PRS-CS) to model the joint effects of over 1 million common variants on disease risk and compared this to an approach (PRS-CT) that only selects a limited set of independent variants that reach genome-wide significance (P<5×10-8). PRS models were trained using GWAS stratified by histological (10,346 cases, 14,687 controls) and molecular subtype (2,632 cases, 2,445 controls), and validated in two independent cohorts. RESULTS: PRS-CS was generally more predictive than PRS-CT with a median increase in explained variance (R2) of 24% (interquartile range=11-30%) across glioma subtypes. Improvements were pronounced for glioblastoma (GBM), with PRS-CS yielding larger odds ratios (OR) per standard deviation (OR=1.93, P=2.0×10-54 vs. OR=1.83, P=9.4×10-50) and higher explained variance (R2=2.82% vs. R2=2.56%). Individuals in the 80th percentile of the PRS-CS distribution had significantly higher risk of GBM (0.107%) at age 60 compared to those with average PRS (0.046%, P=2.4×10-12). Lifetime absolute risk reached 1.18% for glioma and 0.76% for IDH wildtype tumors for individuals in the 95th PRS percentile. PRS-CS augmented the classification of IDH mutation status in cases when added to demographic factors (AUC=0.839 vs. AUC=0.895, PΔAUC=6.8×10-9). CONCLUSIONS: Genome-wide PRS has potential to enhance the detection of high-risk individuals and help distinguish between prognostic glioma subtypes.

Association of immunoglobulin E levels with glioma risk and survival.

Background: Previous epidemiological studies have reported an association of serum immunoglobulin E (IgE) levels with reduced glioma risk, but the association between IgE and glioma prognosis is not well characterized. This study aimed to examine how sex, tumor subtype, and IgE class modulate the association of serum IgE levels with glioma risk and survival. Methods: We conducted a case-control study using participants from the University of California, San Francisco Adult Glioma Study (1997-2010). Serum IgE levels for total, respiratory and food allergy were measured in adults diagnosed with glioma (n=1,696) and cancer-free controls (n=1,135) matched based on age, sex, and race/ethnicity. Logistic regression was adjusted for patient demographics to assess the association between IgE levels and glioma risk. Multivariable Cox regression adjusted for patient-specific and tumor-specific factors compared survival between the elevated and normal IgE groups. Results: Elevated total IgE was associated with reduced risk of IDH wildtype (OR=0.65, 95% CI: 0.54-0.78) and IDH mutant glioma (OR=0.65, 95% CI: 0.50-0.85). In multivariable Cox regression, elevated respiratory IgE was associated with improved survival for IDH wildtype glioma (HR=0.78, 95% CI: 0.67-0.91). The reduction in mortality risk was more pronounced in females (HR=0.71, 95% CI: 0.53-0.96) than in males (HR=0.80, 95% CI: 0.66-0.97), with improvements in median survival of 6.2 months (P<.001) and 1.6 months (P=0.003), respectively. Conclusion: Elevated serum IgE was associated with improved prognosis for IDH wildtype glioma, with a more pronounced protective effect in females. These results suggest a possible sexual dimorphism and antitumor activity of IgE-mediated immune responses.

Genetic determinants of IgG antibody response to COVID-19 vaccination.

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRß1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

Genome-wide Polygenic Risk Scores Predict Risk of Glioma and Molecular Subtypes.

Background: Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to find efficient ways of capturing genetic risk factors using available germline data. Methods: We developed a novel PRS (PRS-CS) that uses continuous shrinkage priors to model the joint effects of over 1 million polymorphisms on disease risk and compared it to an approach (PRS-CT) that selects a limited set of independent variants that reach genome-wide significance (P<5×10-8). PRS models were trained using GWAS results stratified by histological (10,346 cases, 14,687 controls) and molecular subtype (2,632 cases, 2,445 controls), and validated in two independent cohorts. Results: PRS-CS was consistently more predictive than PRS-CT across glioma subtypes with an average increase in explained variance (R2) of 21%. Improvements were particularly pronounced for glioblastoma tumors, with PRS-CS yielding larger effect sizes (odds ratio (OR)=1.93, P=2.0×10-54 vs. OR=1.83, P=9.4×10-50) and higher explained variance (R2=2.82% vs. R2=2.56%). Individuals in the 95th percentile of the PRS-CS distribution had a 3-fold higher lifetime absolute risk of IDH mutant (0.63%) and IDH wildtype (0.76%) glioma relative to individuals with average PRS. PRS-CS also showed high classification accuracy for IDH mutation status among cases (AUC=0.895). Conclusions: Our novel genome-wide PRS may improve the identification of high-risk individuals and help distinguish between prognostic glioma subtypes, increasing the potential clinical utility of germline genetics in glioma patient management.

Genetic predisposition to altered blood cell homeostasis is associated with glioma risk and survival.

Glioma is a highly fatal brain tumor comprised of molecular subtypes with distinct clinical trajectories. Observational studies have suggested that variability in immune response may play a role in glioma etiology. However, their findings have been inconsistent and susceptible to reverse causation due to treatment effects and the immunosuppressive nature of glioma. We applied genetic variants associated (p<5×10-8) with blood cell traits to a meta-analysis of 3418 glioma cases and 8156 controls. Genetically predicted increase in the platelet to lymphocyte ratio (PLR) was associated with an increased risk of glioma (odds ratio (OR)=1.25, p=0.005), especially in IDH-mutant (IDHmut OR=1.38, p=0.007) and IDHmut 1p/19q non-codeleted (IDHmut-noncodel OR=1.53, p=0.004) tumors. However, reduced glioma risk was observed for higher counts of lymphocytes (IDHmut-noncodel OR=0.70, p=0.004) and neutrophils (IDHmut OR=0.69, p=0.019; IDHmut-noncodel OR=0.60, p=0.009), which may reflect genetic predisposition to enhanced immune-surveillance. In contrast to susceptibility, there was no association with survival in IDHmut-noncodel; however, in IDHmut 1p/19q co-deleted tumors, we observed higher mortality with increasing genetically predicted counts of lymphocytes (hazard ratio (HR)=1.65, 95% CI: 1.24-2.20), neutrophils (HR=1.49, 1.13-1.97), and eosinophils (HR=1.59, 1.18-2.14). Polygenic scores for blood cell traits were also associated with tumor immune microenvironment features, with heterogeneity by IDH status observed for 17 signatures related to interferon signaling, PD-1 expression, and T-cell/Cytotoxic responses. In summary, we identified novel, immune-mediated susceptibility mechanisms for glioma with potential disease management implications.

Pediatric brain tumors: Origins, epidemiology, and classification - The 2022 Brain Tumor Epidemiology Consortium meeting report.

The Brain Tumor Epidemiology Consortium (BTEC) is an international organization that fosters collaboration among scientists focused on understanding the epidemiology of brain tumors with interests ranging from the etiology of brain tumor development and outcomes to the control of morbidity and mortality. The 2022 annual BTEC meeting with the theme "Pediatric Brain Tumors: Origins, Epidemiology, and Classification" was held in Lyon, France on June 20 - 22, 2022. Scientists from North America and Europe presented recent research and progress in the field. The meeting content is summarized in this report.

Antibodies to varicella-zoster virus and three other herpesviruses and survival in adults with glioma.

BACKGROUND: Lifetime exposure to the varicella-zoster virus (VZV) has been consistently inversely associated with glioma risk, however, the relationship of VZV with survival in adults with glioma has not been investigated. In this study, we analyzed the survival of adults with glioma in relation to their antibody measurements to 4 common herpes viral infections, including VZV, measured post-diagnosis. METHODS: We analyzed IgG antibody measurements to VZV, cytomegalovirus (CMV), herpes simplex virus 1/2 (HSV), and Epstein-Barr virus (EBV) collected from 1378 adults with glioma diagnosed between 1991 and 2010. Blood was obtained a median of 3 months after surgery. Associations of patient IgG levels with overall survival were estimated using Cox models adjusted for age, sex, self-reported race, surgery type, dexamethasone usage at blood draw, and tumor grade. Models were stratified by recruitment series and meta-analyzed to account for time-dependent treatment effects. RESULTS: VZV antibody seropositivity was associated with improved survival outcomes in adults with glioma (Hazard ratio, HR = 0.70, 95% Confidence Interval 0.54-0.90, P = .006). Amongst cases who were seropositive for VZV antibodies, survival was significantly improved for those above the 25th percentile of continuous reactivity measurements versus those below (HR = 0.76, 0.66-0.88, P = .0003). Antibody seropositivity to EBV was separately associated with improved survival (HR = 0.71, 0.53-0.96, P = .028). Antibody positivity to 2 other common viruses (CMV, HSV) was not associated with altered survival. CONCLUSIONS: Low levels of VZV or EBV antibodies are associated with poorer survival outcomes for adults with glioma. Differential immune response rather than viral exposure may explain these findings.

Accurate Identification of Mucinous Pancreatic Cystic Lesions Using Small-Volume Analytes.

INTRODUCTION: The accurate identification of mucinous pancreatic cystic lesions (PCLs) is paramount for cancer risk stratification. Cyst fluid carcinoembryonic antigen (CEA), the only routinely used test, requires high volumes and has low sensitivity. We aimed to compare the performance of two investigational small-volume biomarkers, glucose and the protease gastricsin, to CEA for PCL classification. METHODS: We obtained cyst fluid samples from 81 patients with pathologically confirmed PCLs from four institutions between 2003 and 2016. Gastricsin activity was measured using an internally quenched fluorescent substrate. Glucose levels were measured with a standard glucometer. CEA levels were obtained from the medical record. Models using Classification and Regression Trees were created to predict mucinous status. Model performance was evaluated using nested cross-validation. RESULTS: Gastricsin activity, CEA, and glucose levels from patients with mucinous (n = 50) and nonmucinous (n = 31) PCLs were analyzed. Area under the curve (AUC) was similar for individual classifiers (gastricsin volume normalized [GVN] 0.88; gastricsin protein concentration normalized [GPN] 0.95; glucose 0.83; CEA 0.84). The combination of two classifiers did not significantly improve AUC, with CEA + GVN (0.88) performing similarly to CEA + GPN (0.95), GVN + glucose (0.87), GPN + glucose (0.95), and CEA + glucose (0.84). The three-analyte combination performed similarly to single and dual classifiers (GPN + glucose + CEA AUC 0.95; GVN + glucose + CEA AUC 0.87). After multiple comparison corrections, there were no significant differences between the individual, dual, and triple classifiers. CONCLUSIONS: Gastricsin and glucose performed similarly to CEA and required <5% of the volume required for CEA; these classifiers may be useful in patients with limited cyst fluid. Future multicenter prospective studies are needed to validate and compare these novel small-volume biomarkers.

Pre-surgery immune profiles of adult glioma patients.

INTRODUCTION: Although immunosuppression is a known characteristic of glioma, no previous large studies have reported peripheral blood immune cell profiles prior to patient surgery and chemoradiation. This report describes blood immune cell characteristics and associated variables prior to surgery among typical glioma patients seen at a large University practice. METHODS: We analyzed pre-surgery blood samples from 139 glioma patients diagnosed with a new or recurrent grade II/III glioma (LrGG, n = 64) or new glioblastoma (GBM, n = 75) and 454 control participants without glioma. Relative cell fractions of CD4, CD8, B-cells, Natural Killer cells, monocytes, and neutrophils, were estimated via a validated deconvolution algorithm from blood DNA methylation measures from Illumina EPIC arrays. RESULTS: Dexamethasone use at time of blood draw varied by glioma type being highest among patients with IDH wild-type (wt) GBM (75%) and lowest for those with oligodendroglioma (14%). Compared to controls, glioma patients showed statistically significant lower cell fractions for all immune cell subsets except for neutrophils which were higher (all p-values < 0.001), in part because of the higher prevalence of dexamethasone use at time of blood draw for IDHwt GBM. Patients who were taking dexamethasone were more likely to have a low CD4 count (< 200, < 500), increased neutrophils, low absolute lymphocyte counts, higher total cell count and higher NLR. CONCLUSION: We show that pre-surgery blood immune profiles vary by glioma subtype, age, and more critically, by use of dexamethasone. Our results highlight the importance of considering dexamethasone exposures in all studies of immune profiles and of obtaining immune measures prior to use of dexamethasone, if possible.

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). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA 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.

The Epidemiology of Central Nervous System Tumors.

This article reviews the current epidemiology of central nervous system tumors. Population-level basic epidemiology, nationally and internationally, and current understanding of germline genetic risk are discussed, with a focus on known and well-studied risk factors related to the etiology of central nervous system tumors.

Epidemiology of Brain and Other CNS Tumors.

PURPOSE OF REVIEW: Brain and other central nervous system (CNS) tumors, while rare, cause significant morbidity and mortality across all ages. This article summarizes the current state of the knowledge on the epidemiology of brain and other CNS tumors. RECENT FINDINGS: For childhood and adolescent brain and other CNS tumors, high birth weight, non-chromosomal structural birth defects and higher socioeconomic position were shown to be risk factors. For adults, increased leukocyte telomere length, proportion of European ancestry, higher socioeconomic position, and HLA haplotypes increase risk of malignant brain tumors, while immune factors decrease risk. Although no risk factor accounting for a large proportion of brain and other CNS tumors has been discovered, the use of high throughput "omics" approaches and improved detection/measurement of environmental exposures will help us refine our current understanding of these factors and discover novel risk factors for this disease.

Genomic epidemiology of SARS-CoV-2 in the UAE reveals novel virus mutation, patterns of co-infection and tissue specific host immune response.

To unravel the source of SARS-CoV-2 introduction and the pattern of its spreading and evolution in the United Arab Emirates, we conducted meta-transcriptome sequencing of 1067 nasopharyngeal swab samples collected between May 9th and Jun 29th, 2020 during the first peak of the local COVID-19 epidemic. We identified global clade distribution and eleven novel genetic variants that were almost absent in the rest of the world and that defined five subclades specific to the UAE viral population. Cross-settlement human-to-human transmission was related to the local business activity. Perhaps surprisingly, at least 5% of the population were co-infected by SARS-CoV-2 of multiple clades within the same host. We also discovered an enrichment of cytosine-to-uracil mutation among the viral population collected from the nasopharynx, that is different from the adenosine-to-inosine change previously reported in the bronchoalveolar lavage fluid samples and a previously unidentified upregulation of APOBEC4 expression in nasopharynx among infected patients, indicating the innate immune host response mediated by ADAR and APOBEC gene families could be tissue-specific. The genomic epidemiological and molecular biological knowledge reported here provides new insights for the SARS-CoV-2 evolution and transmission and points out future direction on host-pathogen interaction investigation.

Socioeconomic status and childhood central nervous system tumors in California.

PURPOSE: Childhood central nervous system (CNS) tumors are the leading cause of cancer mortality in children. Previous studies have suggested that some childhood cancers, including primary CNS tumors, may be associated with higher socioeconomic status (SES). METHODS: We linked data from the California Cancer Registry to California birth records for children (age 0-19 years) diagnosed with primary CNS tumors during 1988-2011 and analyzed multiple measures of parental SES around the birth of their children and subsequent risk for childhood CNS tumors. Our SES measures included birth record-derived parental education and insurance utilization. For a subset of subjects born between 1997 and 2007, we geocoded addresses and examined census-derived median household income and educational level. RESULTS: We analyzed data for 3,022 children with primary CNS tumors and 10,791 matched controls. We found consistent evidence across multiple measures that lower estimates of SES are associated with a reduced risk of CNS tumors. In tumor subgroup analyses, this relationship was most consistent in astrocytomas and ependymomas, with varying findings for embryonal tumors. CONCLUSION: Higher parental SES appears to be a risk factor for childhood CNS tumors in California. Further research is needed to determine specific exposures that may explain this increased risk.

The landscape of host genetic factors involved in immune response to common viral infections.

BACKGROUND: Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities. METHODS: We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort. RESULTS: Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRß1 at positions 11, 13, 71, and 74 for Epstein-Barr virus (EBV), Varicella zoster virus (VZV), human herpesvirus 7, (HHV7), and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P < 5.0 × 10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, and CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P = 5.0 × 10-15 (MCV), NTN5: P = 1.1 × 10-9 (BKV), and P2RY13: P = 1.1 × 10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases, from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions. CONCLUSIONS: Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.

Pediatric glioma and medulloblastoma risk and population demographics: a Poisson regression analysis.

BACKGROUND: The incidence of pediatric brain tumors varies by race and ethnicity, but these relationships may be confounded by socioeconomic status (SES). In this study, the Surveillance, Epidemiology, and End Results Program (SEER) database was evaluated for associations between race/ethnicity and pediatric glioma and medulloblastoma risk with adjustment for SES. METHODS: Pediatric glioma and medulloblastoma cases from the SEER database (years: 2000-2016) were included. Differences in incidence rates by ethnicity, sex, age, and SES-related factors were evaluated by calculation of age-adjusted incidence rates (AAIRs) and annual percent change (APC). SES-related factors (percentage without less than high school graduation, median household income, and percentage foreign-born) were derived from the census at the county-level (year: 2000). Multivariable Poisson regression models with adjustment for selected covariates were constructed to evaluate risk factors. RESULTS: The highest AAIRs of pediatric glioma were observed among non-Hispanic Whites (AAIR: 2.91 per 100 000, 95%-CI: 2.84-2.99). An increasing incidence of pediatric glioma by calendar time was observed among non-Hispanic Whites and non-Hispanic Blacks (APC: 0.97%, 95%-CI: 0.28-1.68 and APC: 1.59%, 95%-CI: 0.03-3.18, respectively). Hispanic and non-Hispanic Black race/ethnicity was associated with lower risk when compared with non-Hispanic White (incidence rate ratios [IRRs]: 0.66, 95%-CI: 0.63-0.70; and 0.69, 95%-CI: 0.65-0.74, respectively). For medulloblastoma, the highest AAIR was observed for non-Hispanic Whites with a positive APC (1.52%, 95%-CI: 0.15-2.91). Hispanics and non-Hispanic Blacks had statistically significant lower IRRs compared with non-Hispanic Whites (IRRs: 0.83, 95%-CI: 0.73-0.94; and 0.72, 95%-CI: 0.59-0.87, respectively). CONCLUSION: Non-Hispanic White race/ethnicity was associated with higher pediatric glioma and medulloblastoma IRRs in models with adjustments for SES.

The landscape of host genetic factors involved in immune response to common viral infections.

INTRODUCTION: Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and therapeutic opportunities. METHODS: We conducted a comprehensive study including genome-wide and transcriptome-wide association analyses to identify genetic loci associated with immunoglobulin G antibody response to 28 antigens for 16 viruses using serological data from 7924 European ancestry participants in the UK Biobank cohort. RESULTS: Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. We identified specific amino acid (AA) residues that are associated with seroreactivity, the strongest associations presented in a range of AA positions within DRß1 at positions 11, 13, 71, and 74 for Epstein-Barr Virus (EBV), Varicella Zoster Virus (VZV), Human Herpes virus 7, (HHV7) and Merkel cell polyomavirus (MCV). Genome-wide association analyses discovered 7 novel genetic loci outside the HLA associated with viral antibody response (P<5.×10-8), including FUT2 (19q13.33) for human polyomavirus BK (BKV), STING1 (5q31.2) for MCV, as well as CXCR5 (11q23.3) and TBKBP1 (17q21.32) for HHV7. Transcriptome-wide association analyses identified 114 genes associated with response to viral infection, 12 outside of the HLA region, including ECSCR: P=5.0×10-15 (MCV), NTN5: P=1.1×10-9 (BKV), and P2RY13: P=1.1×10-8 EBV nuclear antigen. We also demonstrated pleiotropy between viral response genes and complex diseases; from autoimmune disorders to cancer to neurodegenerative and psychiatric conditions. CONCLUSIONS: Our study confirms the importance of the HLA region in host response to viral infection and elucidates novel genetic determinants beyond the HLA that contribute to host-virus interaction.