Replication of rubella vaccine population genetic studies: validation of HLA genotype and humoral response associations

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.

Vaccinomics: Paving the Way for Personalized Immunization

Vaccines are one of the most important medical advancements in human history. They have been successfully used to control and limit the spread of many of the lethal diseases that have plagued us, such as smallpox and polio. Previous vaccine design methodologies were based on the model of "isolate-inactivateinject", which amounts to giving the same vaccine dose to everyone susceptible to infection. In recent years, the importance of how the host genetic background alters vaccine response necessitated the introduction of vaccinomics, which is aimed at studying the variability of vaccine efficacy by associating genetic variability and immune response to vaccination. Despite the rapid developments in variant screening, data obtained from association studies is often inconclusive and cannot be used to guide the new generation of vaccines. This review aims to compile the polymorphisms in HLA and immune system genes and examine the link with their immune response to vaccination. The compiled data can be used to guide the development of new strategies for vaccination for vulnerable groups. Overall, the highly polymorphic HLA locus had the highest correlation with vaccine response variability for most of the studied vaccines, and it was linked to variation in multiple stages of the immune response to the vaccines for both humoral and cellular immunity. Designing new vaccine technologies and immunization regiments to accommodate for this variability is an important step for reaching a vaccinomics-based approach to vaccination.

Genome-wide determinants of cellular immune responses to mumps vaccine

BACKGROUND

We have previously described genetic polymorphisms in candidate genes that are associated with inter-individual variations in antibody responses to mumps vaccination. To expand upon our previous work, we performed a genome-wide association study (GWAS) to discover host genetic variants associated with mumps vaccine-induced cellular immune responses.

METHODS

We performed a GWAS of mumps-specific immune response outcomes (11 secreted cytokines/chemokines) in a cohort of 1,406 subjects.

RESULTS

Among the 11 cytokine/chemokines we studied, four (IFN-γ, IL-2, IL-1β, and TNFα) demonstrated GWAS signals reaching genome-wide significance (p < 5 × 10-8). A genomic region (encoding Sialic acid-binding immunoglobulin-type lectins/SIGLEC) located on chromosome 19q13 (p < 5 × 10-8) was associated with both IL-1β and TNFα responses. The SIGLEC5/SIGLEC14 region contained 11 statistically significant single nucleotide polymorphisms (SNPs), including the intronic SIGLEC5 rs872629 (p = 1.3E-11) and rs1106476 (p = 1.32E-11) whose alternate alleles were significantly associated with decreased levels of mumps-specific IL-1β (rs872629, p = 1.77E-09; rs1106476, p = 1.78E-09) and TNFα (rs872629, p = 1.3E-11; rs1106476, p = 1.32E-11) production.

CONCLUSIONS

Our results suggest that SNPs in the SIGLEC5/SIGLEC14 genes play a role in cellular and inflammatory immune responses to mumps vaccination. These findings motivate further research into the functional roles of SIGLEC genes in the regulation of mumps vaccine-induced immunity.

Genetics and measles, mumps and rubella vaccine response in childhood and adolescence-A systematic review

Measles, mumps and rubella (MMR) are contagious infectious diseases that can be prevented by immunization. However, MMR infections can occur in previously immunized individuals. The vaccine response is, among other factors, influenced by the combined effects of many genes. This systematic review investigates the genetic influence on measles, mumps and rubella antibody responses after childhood vaccination. In accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA), systematic literature searches were conducted in the medical databases PubMed, EMBASE and PsycINFO. Search strings were adjusted for each database. Citations were included if they measured and compared the immune response with immunogenetics after vaccination with a vaccine containing one or more of the following components: measles, mumps and/or rubella, MMR. The measure of vaccine response studied was antibodies after vaccination. Forty-eight articles were included in the final analysis. The results suggest that genetic determinants, including host genes, and single nucleotide polymorphisms in immune-related genes influence the MMR antibody responses after vaccination. Specifically, replicated associations were found between HLA, CD46, RARB, IRF9, EIF2AK2, cytokine genes and MMR vaccine-induced humoral immune responses. This knowledge can be useful in understanding and predicting immune responses and may have implications for future vaccine strategies.

Genome-Wide Determinants of Cellular Immune Responses to Mumps Vaccine

Background

We have previously described genetic polymorphisms in candidate genes that are associated with inter-individual variations in antibody responses to mumps vaccination. To expand upon our previous work, we performed a genome-wide association study (GWAS) to discover host genetic variants associated with mumps vaccine-induced cellular immune responses.

Methods

We performed a GWAS of mumps-specific immune response outcomes (11 secreted cytokines/chemokines) in a cohort of 1,406 subjects.

Results

Among the 11 cytokine/chemokines we studied, four (IFN-γ, IL-2, IL-1β, and TNFα) demonstrated GWAS signals reaching genome-wide significance (p<5 x 10 -8 ). A genomic region (encoding Sialic acid-binding immunoglobulin-type lectins/SIGLEC) located on chromosome 19q13 (p<5×10 -8 ) was associated with both IL-1β and TNFα responses. The SIGLEC5/SIGLEC14 region contained 11 statistically significant single nucleotide polymorphisms (SNPs), including the intronic SIGLEC5 rs872629 (p=1.3E-11) and rs1106476 (p=1.32E-11) whose alternate alleles were significantly associated with decreased levels of mumps-specific IL-1β (rs872629, p=1.77E-09; rs1106476, p=1.78E-09) and TNFα (rs872629, p=1.3E-11; rs1106476, p=1.32E-11) production.

Conclusions

Our results suggest that SNPs in the SIGLEC5/SIGLEC14 genes play a role in cellular and inflammatory immune responses to mumps vaccination. These findings motivate further research into the functional roles of SIGLEC genes in the regulation of mumps vaccine-induced immunity.

Vaccinomics: A scoping review

BACKGROUND

This scoping review summarizes a key aspect of vaccinomics by collating known associations between heterogeneity in human genetics and vaccine immunogenicity and safety.

METHODS

We searched PubMed for articles in English using terms covering vaccines routinely recommended to the general US population, their effects, and genetics/genomics. Included studies were controlled and demonstrated statistically significant associations with vaccine immunogenicity or safety. Studies of Pandemrix®, an influenza vaccine previously used in Europe, were also included, due to its widely publicized genetically mediated association with narcolepsy.

FINDINGS

Of the 2,300 articles manually screened, 214 were included for data extraction. Six included articles examined genetic influences on vaccine safety; the rest examined vaccine immunogenicity. Hepatitis B vaccine immunogenicity was reported in 92 articles and associated with 277 genetic determinants across 117 genes. Thirty-three articles identified 291 genetic determinants across 118 genes associated with measles vaccine immunogenicity, 22 articles identified 311 genetic determinants across 110 genes associated with rubella vaccine immunogenicity, and 25 articles identified 48 genetic determinants across 34 genes associated with influenza vaccine immunogenicity. Other vaccines had fewer than 10 studies each identifying genetic determinants of their immunogenicity. Genetic associations were reported with 4 adverse events following influenza vaccination (narcolepsy, GBS, GCA/PMR, high temperature) and 2 adverse events following measles vaccination (fever, febrile seizure).

CONCLUSION

This scoping review identified numerous genetic associations with vaccine immunogenicity and several genetic associations with vaccine safety. Most associations were only reported in one study. This illustrates both the potential of and need for investment in vaccinomics. Current research in this field is focused on systems and genetic-based studies designed to identify risk signatures for serious vaccine reactions or diminished vaccine immunogenicity. Such research could bolster our ability to develop safer and more effective vaccines.

Genetic predisposition to adverse events in Chinese children aged 3-24 months after diphtheria, tetanus, acellular pertussis and haemophilus influenzae type b combined vaccination

BACKGROUND

Post-vaccination safety is a major public health concern. The genetic predisposition on immune response has not been clearly identified. Clarifying whether individual genetic predisposition plays a role on adverse events (AEs) is critical for the prevention of AEs.

METHODS

From July 2019 to June 2020, we performed a case-control study among children aged 3-24 months in seven Chinese provinces. Each child received a combination vaccination against diphtheria, tetanus, acellular pertussis, and Haemophilus influenzae type b (DTaP-Hib). Through daily telephone follow-up, we collected AEs within seven days. Oral swab samples were collected to investigate the effects of single nucleotide polymorphisms (SNPs) on the risk of AEs.

RESULTS

304 participants were included in the study. In univariate analysis, we discovered three protective SNPs (rs452204, OR = 0.67, P = 0.0352; rs9282763 and rs839, OR = 0.64, P = 0.0256) and one risk SNP (rs9610, OR = 2.20, P = 0.0397). In multivariate analysis, the effects of rs452204 and rs839 were found to be stable. The interaction between rs452204 and rs9610 was observed (OR = 7.25, 95% CI: 1.44-36.58, P = 0.0165).

CONCLUSION

Genetic predisposition was associated with the risk of AEs after DTaP-Hib vaccination, emphasizing the potential application in the prevention of AEs.

Current Challenges in Vaccinology

The development of vaccines, which prime the immune system to respond to future infections, has led to global declines in morbidity and mortality from dreadful infectious communicable diseases. However, many pathogens of public health importance are highly complex and/or rapidly evolving, posing unique challenges to vaccine development. Several of these challenges include an incomplete understanding of how immunity develops, host and pathogen genetic variability, and an increased societal skepticism regarding vaccine safety. In particular, new high-dimensional omics technologies, aided by bioinformatics, are driving new vaccine development (vaccinomics). Informed by recent insights into pathogen biology, host genetic diversity, and immunology, the increasing use of genomic approaches is leading to new models and understanding of host immune system responses that may provide solutions in the rapid development of novel vaccine candidates.

Polymorphisms in the Wilms Tumor Gene Are Associated With Interindividual Variations in Rubella Virus-Specific Cellular Immunity After Measles-Mumps-Rubella II Vaccination

Rubella vaccination induces widely variable immune responses in vaccine recipients. While rubella vaccination is effective at inducing immunity to rubella infection in most subjects, up to 5% of individuals do not achieve or maintain long-term protective immunity. To expand upon our previous work identifying genetic polymorphisms that are associated with these interindividual differences in humoral immunity to rubella virus, we performed a genome-wide association study in a large cohort of 1843 subjects to discover single-nucleotide polymorphisms (SNPs) associated with rubella virus-specific cellular immune responses. We identified SNPs in the Wilms tumor protein gene (WT1) that were significantly associated (P < 5 × 10-8) with interindividual variations in rubella-specific interleukin 6 secretion from subjects' peripheral blood mononuclear cells postvaccination. No SNPs were found to be significantly associated with variations in rubella-specific interferon-γ secretion. Our findings demonstrate that genetic polymorphisms in the WT1 gene in subjects of European ancestry are associated with interindividual differences in rubella virus-specific cellular immunity after measles-mumps-rubella II vaccination.

Genome-wide compound heterozygosity analysis highlighted 4 novel susceptibility loci for congenital heart disease in Chinese population

Genome-wide association studies (GWASs) have achieved great success in deciphering the genetic cause of congenital heart disease (CHD). However, the heritability of CHD remains to be clarified, and numerous genetic factors responsible for occurrence of CHD are yet unclear. In this study, we performed a genome-wide search for relaxed forms of compound heterozygosity (CH) in association with CHD using our existing GWAS data including 2265 individuals (957 CHD cases and 1308 controls). CollapsABEL was used to iteratively test the association between the CH genotype and the CHD phenotype in a sliding window manner. We highlighted 17 genetic loci showing suggestive CH-like associations with CHD (P < 5 × 10^-8 ), among which 4 genetic loci had expression quantitative trait loci (eQTL) effects in blood (P_eQTL  < 0.01). After conditional association analysis, each loci had only 1 independently effective signal reaching the significance threshold (rs2071477/rs3129299 at 6p21.32, P = 2.47 × 10^-10 ; rs10773097/rs2880921 at 12q24.31, P = 3.30 × 10^-8 ; rs73032040/rs7259476 at 19q13.11, P = 1.14 × 10^-8 ; rs10416386/rs4239517 at 19q13.31, P = 1.15 × 10^-9 ), together explained 7.83% of the CHD variance. Among these 4 associated loci, outstanding candidates for CHD-associated genes included UBC, CFM2, ZNF302, LYPD3 and CADM4. Although replication studies with larger sample size are warranted, the first CH GWAS of CHD may extend our current knowledge of the genetic contributions to CHD in the Han Chinese population.

A study of vaccine-induced immune pressure on breakthrough infections in the Phambili phase 2b HIV-1 vaccine efficacy trial

INTRODUCTION

The Merck Adenovirus-5 Gag/Pol/Nef HIV-1 subtype-B vaccine evaluated in predominately subtype B epidemic regions (Step Study), while not preventing infection, exerted vaccine-induced immune pressure on HIV-1 breakthrough infections. Here we investigated if the same vaccine exerted immune pressure when tested in the Phambili Phase 2b study in a subtype C epidemic.

MATERIALS AND METHODS

A sieve analysis, which compares breakthrough viruses from placebo and vaccine arms, was performed on 277 near full-length genomes generated from 23 vaccine and 20 placebo recipients. Vaccine coverage was estimated by computing the percentage of 9-mers that were exact matches to the vaccine insert.

RESULTS

There was significantly greater protein distances from the vaccine immunogen sequence in Gag (p=0.045) and Nef (p=0.021) in viruses infecting vaccine recipients compared to placebo recipients. Twenty-seven putative sites of vaccine-induced pressure were identified (p<0.05) in Gag (n=10), Pol (n=7) and Nef (n=10), although they did not remain significant after adjustment for multiple comparisons. We found the epitope sieve effect in Step was driven by HLA A∗02:01; an allele which was found in low frequency in Phambili participants compared to Step participants. Furthermore, the coverage of the vaccine against subtype C Phambili viruses was 31%, 46% and 14% for Gag, Pol and Nef, respectively, compared to subtype B Step virus coverage of 56%, 61% and 26%, respectively.

DISCUSSION

This study presents evidence of sieve effects in Gag and Nef; however could not confirm effects on specific amino acid sites. We propose that this weaker signal of vaccine immune pressure detected in the Phambili study compared to the Step study may have been influenced by differences in host genetics (HLA allele frequency) and reduced impact of vaccine-induced immune responses due to mismatch between the viral subtype in the vaccine and infecting subtypes.

Searching for the human genetic factors standing in the way of universally effective vaccines

Vaccines have revolutionized modern public health. The effectiveness of some vaccines is limited by the variation in response observed between individuals and across populations. There is compelling evidence that a significant proportion of this variability can be attributed to human genetic variation, especially for those vaccines administered in early life. Identifying and understanding the determinants of this variation could have a far-reaching influence upon future methods of vaccine design and deployment. In this review, we summarize the genetic studies that have been undertaken attempting to identify the genetic determinants of response heterogeneity for the vaccines against hepatitis B, measles and rubella. We offer a critical appraisal of these studies and make a series of suggestions about how modern genetic techniques, including genome-wide association studies, could be used to characterize the genetic architecture of vaccine response heterogeneity. We conclude by suggesting how the findings from such studies could be translated to improve vaccine effectiveness and target vaccination in a more cost-effective manner.

Variability in Humoral Immunity to Measles Vaccine: New Developments

Despite the existence of an effective measles vaccine, resurgence in measles cases in the USA and across Europe has occurred, including in individuals vaccinated with two doses of the vaccine. Host genetic factors result in inter-individual variation in measles vaccine-induced antibodies, and play a role in vaccine failure. Studies have identified HLA (human leukocyte antigen) and non-HLA genetic influences that individually or jointly contribute to the observed variability in the humoral response to vaccination among healthy individuals. In this exciting era, new high-dimensional approaches and techniques including vaccinomics, systems biology, GWAS, epitope prediction and sophisticated bioinformatics/statistical algorithms provide powerful tools to investigate immune response mechanisms to the measles vaccine. These might predict, on an individual basis, outcomes of acquired immunity post measles vaccination.

Overall conceptual framework for studying the genetics of autoimmune diseases following vaccination: a regulatory perspective

The US Vaccine Adverse Event Reporting System contains case reports of autoimmune diseases (ADs) occurring following vaccinations. ADs are rare and occur in unvaccinated people, making the potential association between vaccines and ADs challenging to evaluate. Developing mechanistic pathways that link genes, immune mediators, vaccine components and ADs would be helpful for hypothesis generation, enhancing theories of biologic plausibility and grouping rare autoimmune adverse events to increase the ability to detect and evaluate safety signals. Here, we propose a conceptual framework for investigating the genetics of ADs as safety signals following vaccination, potentially contributing to the identification of relevant biomarkers. We also discuss a study design that incorporates genetic information into postmarket clinical evaluation of autoimmune adverse events following vaccination.

Single nucleotide polymorphisms/haplotypes associated with multiple rubella-specific immune response outcomes post-MMR immunization in healthy children

The observed heterogeneity in rubella-specific immune response phenotypes post-MMR vaccination is thought to be explained, in part, by inter-individual genetic variation. In this study, single nucleotide polymorphisms (SNPs) and multiple haplotypes in several candidate genes were analyzed for associations with more than one rubella-specific immune response outcome, including secreted IFN-γ, secreted IL-6, and neutralizing antibody titers. Overall, we identified 23 SNPs in 10 different genes that were significantly associated with at least two rubella-specific immune responses. Of these SNPs, we detected eight in the PVRL3 gene, five in the PVRL1 gene, one in the TRIM22 gene, two in the IL10RB gene, two in the TLR4 gene, and five in other genes (PVR, ADAR, ZFP57, MX1, and BTN2A1/BTN3A3). The PVRL3 gene haplotype GACGGGGGCAGCAAAAAGAAGAGGAAAGAACAA was significantly associated with both higher IFN-γ secretion (t-statistic 4.43, p < 0.0001) and higher neutralizing antibody titers (t-statistic 3.14, p = 0.002). Our results suggest that there is evidence of multigenic associations among identified gene SNPs and that polymorphisms in these candidate genes contribute to the overall observed differences between individuals in response to live rubella virus vaccine. These results will aid our understanding of mechanisms behind rubella-specific immune response to MMR vaccine and influence the development of vaccines in the future.

Rubella

Rubella remains an important pathogen worldwide, with roughly 100,000 cases of congenital rubella syndrome estimated to occur every year. Rubella-containing vaccine is highly effective and safe and, as a result, endemic rubella transmission has been interrupted in the Americas since 2009. Incomplete rubella vaccination programmes result in continued disease transmission, as evidenced by recent large outbreaks in Japan and elsewhere. In this Seminar, we provide present results regarding rubella control, elimination, and eradication policies, and a brief review of new laboratory diagnostics. Additionally, we provide novel information about rubella-containing vaccine immunogenetics and review the emerging evidence of interindividual variability in humoral and cell-mediated innate and adaptive immune responses to rubella-containing vaccine and their association with haplotypes and single-nucleotide polymorphisms across the human genome.

Genetics and vaccines in the era of personalized medicine

Vaccines represent the most successful and sustainable tactic to prevent and counteract infection. A vaccine generally improves immunity to a particular disease upon administration by inducing specific protective and efficient immune responses in all of the receiving population. The main known factors influencing the observed heterogeneity for immune re-sponses induced by vaccines are gender, age, co-morbidity, immune system, and genetic background. This review is mainly focused on the genetic status effect to vaccine immune responses and how this could contribute to the development of novel vaccine candidates that could be better directed and predicted relative to the genetic history of an individual and/or population. The text offers a brief history of vaccinology as a field, a description of the genetic status of the most relevant and studied genes and their functionality and correlation with exposure to specific vaccines; followed by an inside look into autoimmunity as a concern when designing vaccines as well as perspectives and conclusions looking towards an era of personalized and predictive vaccinology instead of a one size fits all approach.

Adversomics: a new paradigm for vaccine safety and design

Despite the enormous population benefits of routine vaccination, vaccine adverse events (AEs) and reactions, whether real or perceived, have posed one of the greatest barriers to vaccine acceptance--and thus to infectious disease prevention--worldwide. A truly integrated clinical, translational, and basic science approach is required to understand the mechanisms behind vaccine AEs, predict them, and then apply this knowledge to new vaccine design approaches that decrease, or avoid, these events. The term 'adversomics' was first introduced in 2009 and refers to the study of vaccine adverse reactions using immunogenomics and systems biology approaches. In this review, we present the current state of adversomics research, review known associations and mechanisms of vaccine AEs/reactions, and outline a plan for the further development of this emerging research field.

Polymorphisms in HLA-DPB1 are associated with differences in rubella virus-specific humoral immunity after vaccination

Vaccination with live attenuated rubella virus induces a strong immune response in most individuals. However, small numbers of subjects never reach or maintain protective antibody levels, and there is a high degree of variability in immune response. We have previously described genetic polymorphisms in HLA and other candidate genes that are associated with interindividual differences in humoral immunity to rubella virus. To expand our previous work, we performed a genome-wide association study (GWAS) to discover single-nucleotide polymorphisms (SNPs) associated with rubella virus-specific neutralizing antibodies. We identified rs2064479 in the HLA-DPB1 genetic region as being significantly associated with humoral immune response variations after rubella vaccination (P = 8.62 × 10(-8)). All other significant SNPs in this GWAS were located near the HLA-DPB1 gene (P ≤ 1 × 10(-7)). These findings demonstrate that polymorphisms in HLA-DPB1 are strongly associated with interindividual differences in neutralizing antibody levels to rubella vaccination and represent a validation of our previous HLA work.

HLA alleles associated with the adaptive immune response to smallpox vaccine: a replication study

We previously reported HLA allelic associations with vaccinia virus (VACV)-induced adaptive immune responses in a cohort of healthy individuals (n = 1,071 subjects) after a single dose of the licensed smallpox (Dryvax) vaccine. This study demonstrated that specific HLA alleles were significantly associated with VACV-induced neutralizing antibody (NA) titers (HLA-B*13:02, *38:02, *44:03, *48:01, and HLA-DQB1*03:02, *06:04) and cytokine (HLA-DRB1*01:03, *03:01, *10:01, *13:01, *15:01) immune responses. We undertook an independent study of 1,053 healthy individuals and examined associations between HLA alleles and measures of adaptive immunity after a single dose of Dryvax-derived ACAM2000 vaccine to evaluate previously discovered HLA allelic associations from the Dryvax study and determine if these associations are replicated with ACAM2000. Females had significantly higher NA titers than male subjects in both study cohorts [median ID50 discovery cohort 159 (93, 256) vs. 125 (75, 186), p < 0.001; replication cohort 144 (82, 204) vs. 110 (61, 189), p = 0.024]. The association between the DQB1*03:02 allele (median ID50 discovery cohort 152, p = 0.015; replication cohort 134, p = 0.010) and higher NA titers was replicated. Two HLA associations of comparable magnitudes were consistently found between DRB1*04:03 and DRB1*08:01 alleles and IFN-γ ELISPOT responses. The association between the DRB1*15:01 allele with IFN-γ secretion was also replicated (median pg/mL discovery cohort 182, p = 0.052; replication cohort 203, p = 0.014). Our results suggest that smallpox vaccine-induced adaptive immune responses are significantly influenced by HLA gene polymorphisms. These data provide information for functional studies and design of novel candidate smallpox vaccines.

Genetic polymorphisms associated with rubella virus-specific cellular immunity following MMR vaccination

Rubella virus causes a relatively benign disease in most cases, although infection during pregnancy can result in serious birth defects. An effective vaccine has been available since the early 1970s and outbreaks typically do not occur among highly vaccinated (≥2 doses) populations. Nevertheless, considerable inter-individual variation in immune response to rubella immunization does exist, with single-dose seroconversion rates ~95 %. Understanding the mechanisms behind this variability may provide important insights into rubella immunity. In the current study, we examined associations between single nucleotide polymorphisms (SNPs) in selected cytokine, cytokine receptor, and innate/antiviral genes and immune responses following rubella vaccination in order to understand genetic influences on vaccine response. Our approach consisted of a discovery cohort of 887 subjects aged 11-22 at the time of enrollment and a replication cohort of 542 older adolescents and young adults (age 18-40). Our data indicate that SNPs near the butyrophilin genes (BTN3A3/BTN2A1) and cytokine receptors (IL10RB/IFNAR1) are associated with variations in IFNγ secretion and that multiple SNPs in the PVR gene, as well as SNPs located in the ADAR gene, exhibit significant associations with rubella virus-specific IL-6 secretion. This information may be useful, not only in furthering our understanding immune responses to rubella vaccine, but also in identifying key pathways for targeted adjuvant use to boost immunity in those with weak or absent immunity following vaccination.

Analysis of HLA A*02 association with vaccine efficacy in the RV144 HIV-1 vaccine trial

The RV144 HIV-1 vaccine trial demonstrated partial efficacy of 31% against HIV-1 infection. Studies into possible correlates of protection found that antibodies specific to the V1 and V2 (V1/V2) region of envelope correlated inversely with infection risk and that viruses isolated from trial participants contained genetic signatures of vaccine-induced pressure in the V1/V2 region. We explored the hypothesis that the genetic signatures in V1 and V2 could be partly attributed to selection by vaccine-primed T cells. We performed a T-cell-based sieve analysis of breakthrough viruses in the RV144 trial and found evidence of predicted HLA binding escape that was greater in vaccine versus placebo recipients. The predicted escape depended on class I HLA A*02- and A*11-restricted epitopes in the MN strain rgp120 vaccine immunogen. Though we hypothesized that this was indicative of postacquisition selection pressure, we also found that vaccine efficacy (VE) was greater in A*02-positive (A*02(+)) participants than in A*02(-) participants (VE = 54% versus 3%, P = 0.05). Vaccine efficacy against viruses with a lysine residue at site 169, important to antibody binding and implicated in vaccine-induced immune pressure, was also greater in A*02(+) participants (VE = 74% versus 15%, P = 0.02). Additionally, a reanalysis of vaccine-induced immune responses that focused on those that were shown to correlate with infection risk suggested that the humoral responses may have differed in A*02(+) participants. These exploratory and hypothesis-generating analyses indicate there may be an association between a class I HLA allele and vaccine efficacy, highlighting the importance of considering HLA alleles and host immune genetics in HIV vaccine trials.

IMPORTANCE

The RV144 trial was the first to show efficacy against HIV-1 infection. Subsequently, much effort has been directed toward understanding the mechanisms of protection. Here, we conducted a T-cell-based sieve analysis, which compared the genetic sequences of viruses isolated from infected vaccine and placebo recipients. Though we hypothesized that the observed sieve effect indicated postacquisition T-cell selection, we also found that vaccine efficacy was greater for participants who expressed HLA A*02, an allele implicated in the sieve analysis. Though HLA alleles have been associated with disease progression and viral load in HIV-1 infection, these data are the first to suggest the association of a class I HLA allele and vaccine efficacy. While these statistical analyses do not provide mechanistic evidence of protection in RV144, they generate testable hypotheses for the HIV vaccine community and they highlight the importance of assessing the impact of host immune genetics in vaccine-induced immunity and protection. (This study has been registered at ClinicalTrials.gov under registration no. NCT00223080.).

Associations between single nucleotide polymorphisms in cellular viral receptors and attachment factor-related genes and humoral immunity to rubella vaccination

Background

Viral attachment and cell entry host factors are important for viral replication, pathogenesis, and the generation and sustenance of immune responses after infection and/or vaccination, and are plausible genetic regulators of vaccine-induced immunity.

Methods

Using a tag-SNP approach in candidate gene study, we assessed the role of selected cell surface receptor genes, attachment factor-related genes, along with other immune genes in the genetic control of immune response variations after live rubella vaccination in two independent study cohorts.

Results

Our analysis revealed evidence for multiple associations between genetic variants in the PVR, PVRL2, CD209/DC-SIGN, RARB, MOG, IL6 and other immune function-related genes and rubella-specific neutralizing antibodies after vaccination (meta p-value <0.05).

Conclusion

Our results indicate that multiple SNPs from genes involved in cell adhesion, viral attachment, and viral entry, as well as others in genes involved in signaling and/or immune response regulation, play a role in modulating humoral immune responses following live rubella vaccination.

HLA genotypes and rubella vaccine immune response: additional evidence

Recent population-based studies have demonstrated the genetic heritability of rubella vaccine response and assessed that the HLA system may explain about 20% of the inter-individual variance in humoral immune response to this vaccine. Our earlier studies compared HLA allelic associations with rubella vaccine-specific antibodies between two smaller cohorts of healthy Rochester, MN, children (346 and 396 subjects) after two doses of rubella-containing vaccine. This study found that specific HLA alleles were consistently associated with rubella-specific antibody titers (B*27:05, DPA1*02:01, and DPB1*04:01 alleles). The current study examined HLA associations in an independent larger cohort of 1012 healthy San Diego, CA, subjects (age 19-40 years) after rubella vaccine in order to replicate our previous findings in the Rochester subjects. Two HLA associations of comparable magnitudes were consistently observed between B*27:05 (median NT50 Rochester cohort 48.9, p=0.067; San Diego cohort 54.8, p=0.047) and DPB1*04:01 (median NT50 Rochester cohort 61.6, p<0.001; San Diego cohort 70.8, p=0.084) alleles and rubella virus-neutralizing antibody titers. Additional HLA alleles resulted in consistent effects on IL-6 production in both cohorts, but did not meet criteria for statistical significance. Our data suggest these HLA alleles play a role in rubella vaccine-induced immunity and provide the basis for future studies that may explain the mechanism(s) by which these HLA polymorphisms affect immune responses to rubella vaccine.

Genome-wide SNP associations with rubella-specific cytokine responses in measles-mumps-rubella vaccine recipients

Genetic polymorphisms are known to affect responses to both viral infection and vaccination. Our previous work has described genetic polymorphisms significantly associated with variations in immune response to rubella vaccine from multiple gene families with known immune function, including HLA, cytokine and cytokine receptor genes, and in genes controlling innate and adaptive immunity. In this study, we assessed cellular immune responses (IFNγ and IL-6) in a cohort of healthy younger individuals and performed genome-wide SNP analysis on these same individuals. Here, we report the first genome-wide association study focused on immune responses following rubella vaccination. Our results indicate that rs16928280 in protein tyrosine phosphatase delta (PTPRD) and a collection of SNPs in ACO1 (encoding an iron regulatory protein) are associated with interindividual variations in IFNγ response to rubella virus stimulation. In contrast, we did not identify any significant genetic associations with rubella-specific IL-6 response. These genetic regions may influence rubella vaccine-induced IFNγ responses and warrant further studies in additional cohorts in order to confirm these findings.

The personal touch: strategies toward personalized vaccines and predicting immune responses to them

The impact of vaccines on public health and wellbeing has been profound. Smallpox has been eradicated, polio is nearing eradication, and multiple diseases have been eliminated from certain areas of the world. Unfortunately, we now face diseases such as hepatitis C, malaria or tuberculosis, as well as new and re-emerging pathogens for which we lack effective vaccines. Empirical approaches to vaccine development have been successful in the past, but may not be up to the current infectious disease challenges facing us. New, directed approaches to vaccine design, development, and testing need to be developed. Ideally these approaches will capitalize on cutting-edge technologies, advanced analytical and modeling strategies, and up-to-date knowledge of both pathogen and host. These approaches will pay particular attention to the causes of inter-individual variation in vaccine response in order to develop new vaccines tailored to the unique needs of individuals and communities within the population.

High-throughput assay optimization and statistical interpolation of rubella-specific neutralizing antibody titers

Rubella remains a social and economic burden due to the high incidence of congenital rubella syndrome (CRS) in some countries. For this reason, an accurate and efficient high-throughput measure of antibody response to vaccination is an important tool. In order to measure rubella-specific neutralizing antibodies in a large cohort of vaccinated individuals, a high-throughput immunocolorimetric system was developed. Statistical interpolation models were applied to the resulting titers to refine quantitative estimates of neutralizing antibody titers relative to the assayed neutralizing antibody dilutions. This assay, including the statistical methods developed, can be used to assess the neutralizing humoral immune response to rubella virus and may be adaptable for assessing the response to other viral vaccines and infectious agents.

Characterization of humoral and cellular immunity to rubella vaccine in four distinct cohorts

Although vaccination campaigns have significantly reduced the global burden of rubella disease, there are still regional outbreaks and cases of congenital rubella syndrome. Rubella vaccination elicits a strong humoral as well as cellular response. The relationship between these two measures in response to rubella vaccine is poorly understood. We have previously reported no correlation between rubella-virus-specific cytokine secretion and IgG antibody levels after rubella vaccination. In the current study, we extend our previous work to report correlations between secreted cytokines and functional neutralizing antibodies after rubella vaccination in four distinct cohorts. There was evidence of significant differences (p < 0.05) in rubella-virus-specific humoral and cellular responses between cohorts. When investigating relationships between rubella-vaccine-specific humoral and cellular immunity, we observed a significant correlation between neutralizing antibodies and IFN-γ (r(s) = 0.21, p = 0.0004). We also observed correlations in subjects with extreme humoral immune phenotypes and IFN-γ levels in two of the four cohorts (r(s) = 0.32, p = 0.01; r(s) = 0.36, p = 0.01, respectively). These findings indicate that there is a high level of heterogeneity in rubella-specific immune responses between study populations. We believe that the novel correlation discovered between IFN-γ and neutralizing antibody titers will give future insight into the functional mechanisms of immunity induced by rubella virus and other live viral vaccines.

Genetic variants within the MHC region are associated with immune responsiveness to childhood vaccinations

The influence of genetic variability within the major histocompatibility complex (MHC) region on variations in immune responses to childhood vaccination was investigated. The study group consisted of 135 healthy infants who had been immunized with hepatitis B (HBV), 7-valent pneumococcal conjugate (PCV7), and diphtheria, tetanus, acellular pertussis (DTaP) vaccines according to standard childhood immunization schedules. Genotype analysis was performed on genomic DNA using Illumina Goldengate MHC panels (Mapping and Exon Centric). At the 1 year post vaccination check-up total, isotypic, and antigen-specific serum antibody levels were measured using multiplex immunoassays. A number of single nucleotide polymorphisms (SNPs) within MHC Class I and II genes were found to be associated with variations in the vaccine specific antibody responses and serum levels of immunoglobulins (IgG, IgM) and IgG isotypes (IgG1, IgG4) (all at p<0.001). Linkage disequilibrium patterns and functional annotations showed that significant SNPs were strongly correlated with other functional regulatory SNPs. These SNPs were found to regulate the expression of a group of genes involved in antigen processing and presentation including HLA-A, HLA-C, HLA-G, HLA-H, HLA-DRA, HLA-DRB1, HLA-DRB5, HLA-DQA1, HLA-DQB1, HLA-DOB, and TAP-2. The results suggest that genetic variations within particular MHC genes can influence immune response to common childhood vaccinations, which in turn may influence vaccine efficacy.

Human leukocyte antigens and cellular immune responses to anthrax vaccine adsorbed

Interindividual variations in vaccine-induced immune responses are in part due to host genetic polymorphisms in the human leukocyte antigen (HLA) and other gene families. This study examined associations between HLA genotypes, haplotypes, and homozygosity and protective antigen (PA)-specific cellular immune responses in healthy subjects following immunization with Anthrax Vaccine Adsorbed (AVA). While limited associations were observed between individual HLA alleles or haplotypes and variable lymphocyte proliferative (LP) responses to AVA, analyses of homozygosity supported the hypothesis of a "heterozygote advantage." Individuals who were homozygous for any HLA locus demonstrated significantly lower PA-specific LP than subjects who were heterozygous at all eight loci (median stimulation indices [SI], 1.84 versus 2.95, P = 0.009). Similarly, we found that class I (HLA-A) and class II (HLA-DQA1 and HLA-DQB1) homozygosity was significantly associated with an overall decrease in LP compared with heterozygosity at those three loci. Specifically, individuals who were homozygous at these loci had significantly lower PA-specific LP than subjects heterozygous for HLA-A (median SI, 1.48 versus 2.13, P = 0.005), HLA-DQA1 (median SI, 1.75 versus 2.11, P = 0.007), and HLA-DQB1 (median SI, 1.48 versus 2.13, P = 0.002) loci, respectively. Finally, homozygosity at an increasing number (≥ 4) of HLA loci was significantly correlated with a reduction in LP response (P < 0.001) in a dose-dependent manner. Additional studies are needed to reproduce these findings and determine whether HLA-heterozygous individuals generate stronger cellular immune response to other virulence factors (Bacillus anthracis LF and EF) than HLA-homozygous subjects.

Vaccinomics, adversomics, and the immune response network theory: individualized vaccinology in the 21st century

Vaccines, like drugs and medical procedures, are increasingly amenable to individualization or personalization, often based on novel data resulting from high throughput "omics" technologies. As a result of these technologies, 21st century vaccinology will increasingly see the abandonment of a "one size fits all" approach to vaccine dosing and delivery, as well as the abandonment of the empiric "isolate-inactivate-inject" paradigm for vaccine development. In this review, we discuss the immune response network theory and its application to the new field of vaccinomics and adversomics, and illustrate how vaccinomics can lead to new vaccine candidates, new understandings of how vaccines stimulate immune responses, new biomarkers for vaccine response, and facilitate the understanding of what genetic and other factors might be responsible for rare side effects due to vaccines. Perhaps most exciting will be the ability, at a systems biology level, to integrate increasingly complex high throughput data into descriptive and predictive equations for immune responses to vaccines. Herein, we discuss the above with a view toward the future of vaccinology.

Characterizing vaccine responses using host genomic and transcriptomic analysis

Vaccines have had a profound influence on human health with no other health intervention rivaling their impact on the morbidity and mortality associated with infectious disease. However, the magnitude and persistence of vaccine immunity varies considerably between individuals, a phenomenon that is not well understood. Recent studies have used contemporary technologies to correlate variation in the genome and transcriptome to immunological measures of vaccine responsiveness. These approaches have provided fresh insight into the intrinsic factors determining the potency and duration of vaccine-induced immunity. The fundamental challenge will be to translate these findings into innovative and pragmatic strategies to develop new and more effective vaccines.

Genome-wide characterization of transcriptional patterns in high and low antibody responders to rubella vaccination

Immune responses to current rubella vaccines demonstrate significant inter-individual variability. We performed mRNA-Seq profiling on PBMCs from high and low antibody responders to rubella vaccination to delineate transcriptional differences upon viral stimulation. Generalized linear models were used to assess the per gene fold change (FC) for stimulated versus unstimulated samples or the interaction between outcome and stimulation. Model results were evaluated by both FC and p-value. Pathway analysis and self-contained gene set tests were performed for assessment of gene group effects.

Of 17,566 detected genes, we identified 1,080 highly significant differentially expressed genes upon viral stimulation (p<1.00E⁻¹⁵, FDR<1.00E⁻¹⁴), including various immune function and inflammation-related genes, genes involved in cell signaling, cell regulation and transcription, and genes with unknown function. Analysis by immune outcome and stimulation status identified 27 genes (p≤0.0006 and FDR≤0.30) that responded differently to viral stimulation in high vs. low antibody responders, including major histocompatibility complex (MHC) class I genes (HLA-A, HLA-B and B2M with p = 0.0001, p = 0.0005 and p = 0.0002, respectively), and two genes related to innate immunity and inflammation (EMR3 and MEFV with p = 1.46E⁻⁰⁸ and p = 0.0004, respectively). Pathway and gene set analysis also revealed transcriptional differences in antigen presentation and innate/inflammatory gene sets and pathways between high and low responders. Using mRNA-Seq genome-wide transcriptional profiling, we identified antigen presentation and innate/inflammatory genes that may assist in explaining rubella vaccine-induced immune response variations. Such information may provide new scientific insights into vaccine-induced immunity useful in rational vaccine development and immune response monitoring.

Consistency of HLA associations between two independent measles vaccine cohorts: a replication study

Associations between HLA genotypes and measles vaccine humoral and cellular immune responses were examined to better understand immunogenetic drivers of vaccine response. Two independent study cohorts of healthy schoolchildren were examined: cohort one, 346 children between 12 and 18 years of age; and cohort two, 388 children between 11 and 19 years of age. All received two age-appropriate doses of measles-containing vaccine. The purpose of this study was to identify and replicate associations between HLA genes and immune responses following measles vaccination found in our first cohort. Associations of comparable magnitudes and with similar p-values were observed between B*3503 (1st cohort p=0.01; 2nd cohort p=0.07), DQA1*0201 (1st cohort p=0.03; 2nd cohort p=0.03), DQB1*0303 (1st cohort p=0.10; 2 cohort p=0.02), DQB1*0602 (1st cohort p=0.07; 2nd cohort p=0.10), and DRB1*0701 (1st cohort p=0.03; 2nd cohort p=0.07) alleles and measles-specific antibody levels. Suggestive, yet consistent, associations were observed between the B7 (1st cohort p=0.01; 2nd cohort p=0.08) supertype and higher measles antibody levels in both cohorts. Also, in both cohorts, the B*0801 and DRB1*0301 alleles, C*0802 and DPA1*0202 alleles, and DRB1*1303 alleles displayed consistent associations with variations in IFN-γ, IL-2 and IL-10 secretion, respectively. This study emphasizes the importance of replicating HLA associations with measles vaccine-induced humoral and cellular immune responses and increases confidence in the results. These data will inform strategies for functional studies and novel vaccine development, including epitope-based measles vaccines. This is the first HLA association replication study with measles vaccine-specific immune responses to date.

Vaccinomics: current findings, challenges and novel approaches for vaccine development

Recent years have witnessed a growing interest in a field of vaccinology that we have named vaccinomics. The overall idea behind vaccinomics is to identify genetic and other mechanisms and pathways that determine immune responses, and thereby provide new candidate vaccine approaches. Considerable data show that host genetic polymorphisms act as important determinants of innate and adaptive immunity to vaccines. This review highlights examples of the role of immunogenetics and immunogenomics in understanding immune responses to vaccination, which are highly variable across the population. The influence of HLA genes, non-HLA, and innate genes in inter-individual variations in immune responses to viral vaccines are examined using population-based gene/SNP association studies. The ability to understand relationships between immune response gene variants and vaccine-specific immunity may assist in designing new vaccines. At the same time, application of state-of-the-art next-generation sequencing technology (and bioinformatics) is desired to provide new genetic information and its relationship to the immune response.

Vaccinomics and a new paradigm for the development of preventive vaccines against viral infections

In this article we define vaccinomics as the integration of immunogenetics and immunogenomics with systems biology and immune profiling. Vaccinomics is based on the use of cutting edge, high-dimensional (so called "omics") assays and novel bioinformatics approaches to the development of next-generation vaccines and the expansion of our capabilities in individualized medicine. Vaccinomics will allow us to move beyond the empiric "isolate, inactivate, and inject" approach characterizing past vaccine development efforts, and toward a more detailed molecular and systemic understanding of the carefully choreographed series of biological processes involved in developing viral vaccine-induced "immunity." This enhanced understanding will then be applied to overcome the obstacles to the creation of effective vaccines to protect against pathogens, particularly hypervariable viruses, with the greatest current impact on public health. Here we provide an overview of how vaccinomics will inform vaccine science, the development of new vaccines and/or clinically relevant biomarkers or surrogates of protection, vaccine response heterogeneity, and our understanding of immunosenescence.

Genome-wide association study identifies HLA-DP as a susceptibility gene for pediatric asthma in Asian populations

Asthma is a complex phenotype influenced by genetic and environmental factors. We conducted a genome-wide association study (GWAS) with 938 Japanese pediatric asthma patients and 2,376 controls. Single-nucleotide polymorphisms (SNPs) showing strong associations (P<1×10⁻⁸) in GWAS were further genotyped in an independent Japanese samples (818 cases and 1,032 controls) and in Korean samples (835 cases and 421 controls). SNP rs987870, located between HLA-DPA1 and HLA-DPB1, was consistently associated with pediatric asthma in 3 independent populations (P_combined = 2.3×10⁻¹⁰, odds ratio [OR] = 1.40). HLA-DP allele analysis showed that DPA1*0201 and DPB1*0901, which were in strong linkage disequilibrium, were strongly associated with pediatric asthma (DPA1*0201: P = 5.5×10⁻¹⁰, OR = 1.52, and DPB1*0901: P = 2.0×10⁻⁷, OR = 1.49). Our findings show that genetic variants in the HLA-DP locus are associated with the risk of pediatric asthma in Asian populations.

Asthma is the most common chronic disorder in children, and asthma exacerbation is an important cause of childhood morbidity and hospitalization. Here, taking advantage of recent technological advances in human genetics, we performed a genome-wide association study and follow-up validation studies to identify genetic variants for asthma. By examining 6,428 Asians, we found rs987870 and HLA-DPA1*0201/DPB1*0901 were associated with pediatric asthma. The association signal was stretched in the region of HLA-DPB2, collagen, type XI, alpha 2 (COL11A2), and Retinoid X receptor beta (RXRB), but strong linkage disequilibrium in this region made it difficult to specifically identify causative variants. Interestingly, the SNP (or the HLA-DP allele) associated with pediatric asthma (Th-2 type immune diseases) in the present study confers protection against Th-1 type immune diseases, such as type 1 diabetes and rheumatoid arthritis. Therefore, the association results obtained in the present study could partially explain the inverse relationship between asthma and Th-1 type immune diseases and may lead to better understanding of Th-1/Th-2 immune diseases.

Atopy history and the genomics of wheezing after influenza vaccination in children 6-59 months of age

BACKGROUND

A multinational clinical trial compared the safety and efficacy of intranasal trivalent live attenuated influenza vaccine (LAIV) with intramuscular trivalent inactivated vaccine (TIV) in very young children prior to the 2004-5 influenza season [1]. Wheezing was noted more often in recipients of LAIV and laboratory-confirmed influenza infection was noted more often in recipients of TIV. We sought to determine whether epidemiologic or genetic factors were associated with these outcomes.

METHODS

Atopy surveys and DNA collections were performed in trial participants at two United States sites, Nashville, TN and Boston, MA. DNA samples were genotyped on Illumina Infinium 610 or 660-Quad. Standard allelic tests of association were performed.

RESULTS

At the Nashville and Boston sites, a total of 99 children completed the trial, 6 (1 TIV, 5 LAIV) developed medically attended wheezing within 42 days following vaccination, and 8 (5 TIV, 3 LAIV) developed laboratory-confirmed influenza during the season. Eighty-one surveys and 70 DNA samples were collected. Family history of asthma (p=0.001) was associated with wheezing after vaccination. Of 468,458 single nucleotide polymorphisms tested in the genome-wide association study (GWAS), none achieved genome-wide significance for either wheezing after vaccination or laboratory-confirmed influenza infection.

CONCLUSIONS

Family history of asthma appears to be a risk factor for wheezing after influenza vaccination. Given the limitations of the sample size, our pilot study demonstrated the feasibility of performing a GWAS but was not able to determine genetic polymorphisms associated with wheezing after influenza immunization.

Associations between SNPs in candidate immune-relevant genes and rubella antibody levels: a multigenic assessment

BACKGROUND

The mechanisms of immune response are structured within a highly complex regulatory system. Genetic associations with variation in the immune response to rubella vaccine have typically been assessed one locus at a time. We simultaneously assessed the associations between 726 SNPs tagging 84 candidate immune response genes and rubella-specific antibody levels. Blood samples were obtained from 714 school-aged children who had received two doses of MMR vaccine. Associations between rubella-specific antibody levels and 726 candidate tagSNPs were assessed both one SNP at a time and in a variety of multigenic analyses.

RESULTS

Single-SNP assessments identified 4 SNPs that appeared to be univariately associated with rubella antibody levels: rs2844482 (p = 0.0002) and rs2857708 (p = 0.001) in the 5'UTR of the LTA gene, rs7801617 in the 5'UTR of the IL6 gene (p = 0.0005), and rs4787947 in the 5'UTR of the IL4R gene (p = 0.002). While there was not significant evidence in favor of epistatic genetic associations among the candidate SNPs, multigenic analyses identified 29 SNPs significantly associated with rubella antibody levels when selected as a group (p = 0.017). This collection of SNPs included not only those that were significant univariately, but others that would not have been identified if only considered in isolation from the other SNPs.

CONCLUSIONS

For the first time, multigenic assessment of associations between candidate SNPs and rubella antibody levels identified a broad number of genetic associations that would not have been deemed important univariately. It is important to consider approaches like those applied here in order to better understand the full genetic complexity of response to vaccination.

Extended LTA, TNF, LST1 and HLA gene haplotypes and their association with rubella vaccine-induced immunity

Background

Recent studies have suggested the importance of HLA genes in determining immune responses following rubella vaccine. The telomeric class III region of the HLA complex harbors several genes, including lymphotoxin alpha (LTA), tumor necrosis factor (TNF) and leukocyte specific transcript -1 (LST1) genes, located between the class I B and class II DRB1 loci. Apart from HLA, little is known about the effect of this extended genetic region on HLA haplotypic backgrounds as applied to immune responses.

Methodology/Principal Findings

We examined the association between immune responses and extended class I-class II-class III haplotypes among 714 healthy children after two doses of rubella vaccination. These extended haplotypes were then compared to the HLA-only haplotypes. The most significant association was observed between haplotypes extending across the HLA class I region, ten-SNP haplotypes, and the HLA class II region (i.e. A-C-B-LTA-TNF-LST1-DRB1-DQA1-DQB1-DPA1-DPB1) and rubella-specific antibodies (global p-value of 0.03). Associations were found between both extended A*02-C*03-B*15-AAAACGGGGC-DRB1*04-DQA1*03-DQB1*03-DPA1*01-DPB1*04 (p = 0.002) and HLA-only A*02-C*03-B*15-DRB1*04-DQA1*03-DQB1*03-DPA1*01-DPB1*04 haplotypes (p = 0.009) and higher levels of rubella antibodies. The class II HLA-only haplotype DRB1*13-DQA1*01-DQB1*06-DPA1*01-DPB1*04 (p = 0.04) lacking LTA-TNF-LST1 SNPs was associated with lower rubella antibody responses. Similarly, the class I-class II HLA-only A*01-C*07-B*08-DRB1*03-DQA1*05-DQB1*02-DPA1*01-DPB1*04 haplotype was associated with increased TNF-α secretion levels (p = 0.009). In contrast, the extended AAAACGGGGC-DRB1*01-DQA1*01-DQB1*05-DPA1*01-DPB1*04 (p = 0.01) haplotype was found to trend with decreased rubella-specific IL-6 secretion levels.

Conclusions/Significance

These data suggest the importance of examining both HLA genes and genes in the class III region as part of the extended haplotypes useful in understanding genomic drivers regulating immune responses to rubella vaccine.