Genome-wide genetic associations with IFNγ response to smallpox vaccine

Identification and functionalization of thyrotropin receptor antibodies with different antigenic epitopes

One of the sensitive markers for autoimmune thyroid disease (AITD) clinical identification is thyroid-stimulating hormone receptor antibodies (TRAbs). To quickly distinguish TRAb with distinct antigenic epitopes, a straightforward and uncomplicated technique has not yet been created. The objective of this study is to search for molecular diagnostic targets for different types of AITD {Graves' disease (GD), Graves' orbitopathy (GO), GD with third-degree goiter [GD(3)], hypothyroidism combined with positive TRAb [HT(TRAb+)]} as molecular diagnostic targets. Following action on thyroid cells, differential genes (DEGs) generated by TRAb with distinct antigenic epitopes were detected and identified by RNA sequencing (RNA-Seq), bioinformatics analysis, and quantitative reverse transcription-polymerase chain reaction (RT-qPCR) in the serum of patients with AITD. Using the 5-ethynyl-2'-deoxyuridine (EdU) assay, the effect of coculturing thyroid cells with different antigenic TRAb epitopes on the cells' capacity to proliferate was investigated. Bioinformatics analysis and RT-qPCR validation identified one GD key gene alpha 2-HS glycoprotein (AHSG), two GO key genes [adrenoceptor alpha 1D (ADRA1D) and H2B clustered histone 18 (H2BC18)], two GD(3) key genes [suppressor of cytokine signaling 1 (SOCS1) and cytochrome b-245 beta (CYBB)], and one HT(TRAb+) key gene (MASP2). Correlation analysis and ROC curves showed that the abovementioned genes could be used as molecular diagnostic targets for different types of AITD. Finally, EdU results showed that TRAb inhibited thyroid cell proliferation in the HT(TRAb+) group compared with the normal control group, whereas the remaining three groups promoted thyroid cell proliferation, with a statistically significant difference (P < 0.05). We identified six key genes for different types of AITD, which have diagnostic value for different types of AITD. Meanwhile, we found that TRAbs with different antigenic epitopes in AITD have different biological functions.NEW & NOTEWORTHY We identified six molecular targets of different types of AITD [GD, GO, GD(3), and HT(TRAb+)], which have diagnostic value for different types of AITD. Meanwhile, we found that TRAb with different antigenic epitopes extracted from the sera of patients with AITD had different biological functions, which also provided a new idea for further research on the mechanism of action of TRAb with different antigenic epitopes in AITD.

Polygenic Prediction of Cellular Immune Responses to Mumps Vaccine

In this report, we provide a follow-up analysis of a previously published genome-wide association study of host genetic variants associated with inter-individual variations in cellular immune responses to mumps vaccine. Here we report the results of a polygenic score (PGS) analysis showing how common variants can predict mumps vaccine response. We found higher PGS for IFNγ, IL-2, and TNFα were predictive of higher post-vaccine IFNγ (p-value = 2e-6), IL-2 (p = 2e-7), and TNFα (p = 0.004) levels, respectively. Control of immune responses after vaccination is complex and polygenic in nature. Our results suggest that the PGS-based approach enables better capture of the combined genetic effects that contribute to mumps vaccine-induced immunity, potentially offering a more comprehensive understanding than traditional single-variant GWAS. This approach will likely have broad utility in studying genetic control of immune responses to other vaccines and to infectious diseases.

Candidate Genes for IgA Nephropathy in Pediatric Patients: Exome-Wide Association Study

IgA nephropathy (IgAN) is an autoimmune disorder which is believed to be non-monogenic. We performed an exome-wide association study of 70 children with IgAN and 637 healthy donors. The HLA allele frequencies were compared between the patients and healthy donors from the bone marrow registry of the Pirogov University. We tested 78,020 gene markers for association and performed functional enrichment analysis and transcription factor binding preference detection. We identified 333 genetic variants, employing three inheritance models. The most significant association with the disorder was observed for rs143409664 (PRAG1) in the case of the additive and dominant models (PBONF = 1.808 × 10-15 and PBONF = 1.654 × 10-15, respectively), and for rs13028230 (UBR3) in the case of the recessive model (PBONF = 1.545 × 10-9). Enrichment analysis indicated the strongly overrepresented "immune system" and "kidney development" terms. The HLA-DQA1*01:01:01G allele (p = 0.0076; OR, 2.021 [95% CI, 1.322-3.048]) was significantly the most frequent among IgAN patients. Here, we characterized, for the first time, the genetic background of Russian IgAN patients, identifying the risk alleles typical of the population. The most important signals were detected in previously undescribed loci.

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.

Virus-specific and shared gene expression signatures in immune cells after vaccination in response to influenza and vaccinia stimulation

Background

In the vaccine era, individuals receive multiple vaccines in their lifetime. Host gene expression in response to antigenic stimulation is usually virus-specific; however, identifying shared pathways of host response across a wide spectrum of vaccine pathogens can shed light on the molecular mechanisms/components which can be targeted for the development of broad/universal therapeutics and vaccines.

Method

We isolated PBMCs, monocytes, B cells, and CD8+ T cells from the peripheral blood of healthy donors, who received both seasonal influenza vaccine (within <1 year) and smallpox vaccine (within 1 - 4 years). Each of the purified cell populations was stimulated with either influenza virus or vaccinia virus. Differentially expressed genes (DEGs) relative to unstimulated controls were identified for each in vitro viral infection, as well as for both viral infections (shared DEGs). Pathway enrichment analysis was performed to associate identified DEGs with KEGG/biological pathways.

Results

We identified 2,906, 3,888, 681, and 446 DEGs in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, in response to influenza stimulation. Meanwhile, 97, 120, 20, and 10 DEGs were identified as gene signatures in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, upon vaccinia stimulation. The majority of DEGs identified in PBMCs were also found in monocytes after either viral stimulation. Of the virus-specific DEGs, 55, 63, and 9 DEGs occurred in common in PBMCs, monocytes, and B cells, respectively, while no DEGs were shared in infected CD8+ T cells after influenza and vaccinia. Gene set enrichment analysis demonstrated that these shared DEGs were over-represented in innate signaling pathways, including cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, Toll-like receptor signaling, RIG-I-like receptor signaling pathways, cytosolic DNA-sensing pathways, and natural killer cell mediated cytotoxicity.

Conclusion

Our results provide insights into virus-host interactions in different immune cells, as well as host defense mechanisms against viral stimulation. Our data also highlights the role of monocytes as a major cell population driving gene expression in ex vivo PBMCs in response to viral stimulation. The immune response signaling pathways identified in this study may provide specific targets for the development of novel virus-specific therapeutics and improved vaccines for vaccinia and influenza. Although influenza and vaccinia viruses have been selected in this study as pathogen models, this approach could be applicable to other pathogens.

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.

Machine Learning for the Identification of a Common Signature for Anti-SSA/Ro 60 Antibody Expression Across Autoimmune Diseases

OBJECTIVE

Anti-Ro autoantibodies are among the most frequently detected extractable nuclear antigen autoantibodies, mainly associated with primary Sjögren's syndrome (SS), systemic lupus erythematosus (SLE), and undifferentiated connective tissue disease (UCTD). This study was undertaken to determine if there is a common signature for all patients expressing anti-Ro 60 autoantibodies regardless of their disease phenotype.

METHODS

Using high-throughput multiomics data collected from the cross-sectional cohort in the PRECISE Systemic Autoimmune Diseases (PRECISESADS) study Innovative Medicines Initiative (IMI) project (genetic, epigenomic, and transcriptomic data, combined with flow cytometry data, multiplexed cytokines, classic serology, and clinical data), we used machine learning to assess the integrated molecular profiling of 520 anti-Ro 60+ patients compared to 511 anti-Ro 60- patients with primary SS, patients with SLE, and patients with UCTD, and 279 healthy controls.

RESULTS

The selected clinical features for RNA-Seq, DNA methylation, and genome-wide association study data allowed for a clear distinction between anti-Ro 60+ and anti-Ro 60- patients. The different features selected using machine learning from the anti-Ro 60+ patients constituted specific signatures when compared to anti-Ro 60- patients and healthy controls. Remarkably, the transcript Z score of 3 genes (ATP10A, MX1, and PARP14), presenting with overexpression associated with hypomethylation and genetic variation and independently identified using the Boruta algorithm, was clearly higher in anti-Ro 60+ patients compared to anti-Ro 60- patients regardless of disease type. Our findings demonstrated that these signatures, enriched in interferon-stimulated genes, were also found in anti-Ro 60+ patients with rheumatoid arthritis and those with systemic sclerosis and remained stable over time and were not affected by treatment.

CONCLUSION

Anti-Ro 60+ patients present with a specific inflammatory signature regardless of their disease type, suggesting that a dual therapeutic approach targeting both Ro-associated RNAs and anti-Ro 60 autoantibodies should be considered.

Proteomics study of Mycoplasma pneumoniae pneumonia reveals the Fc fragment of the IgG-binding protein as a serum biomarker and implicates potential therapeutic targets

Macrolide and corticosteroid resistance has been reported in patients with Mycoplasma pneumoniae (MP) pneumonia (MPP). MP clearance is difficult to achieve through antibiotic treatment in sensitive patients with severe MPP (SMPP). SMPP in children might progress to airway remodeling and even bronchiolitis/bronchitis obliterans. Therefore, identifying serum biomarkers that indicate MPP progression and exploring new targeted drugs for SMPP treatment require urgency. In this study, serum samples were collected from patients with general MPP (GMPP) and SMPP to conduct proteomics profiling. The Fc fragment of the IgG-binding protein (FCGBP) was identified as the most promising indicator of SMPP. Biological enrichment analysis indicated uncontrolled inflammation in SMPP. ELISA results proved that the FCGBP level in patients with SMPP was substantially higher than that in patients with GMPP. Furthermore, the FCGBP levels showed a decreasing trend in patients with GMPP but the opposite trend in patients with SMPP during disease progression. Connectivity map analyses identified 25 possible targeted drugs for SMPP treatment. Among them, a mechanistic target of rapamycin kinase (mTOR) inhibitor, which is a macrolide compound and a cell proliferation inhibitor, was the most promising candidate for targeting SMPP. To our knowledge, this study was the first proteomics-based characterization of patients with SMPP and GMPP.

Polymorphisms in STING Affect Human Innate Immune Responses to Poxviruses

We conducted a large genome-wide association study (GWAS) of the immune responses to primary smallpox vaccination in a combined cohort of 1,653 subjects. We did not observe any polymorphisms associated with standard vaccine response outcomes (e.g., neutralizing antibody, T cell ELISPOT response, or T cell cytokine production); however, we did identify a cluster of SNPs on chromosome 5 (5q31.2) that were significantly associated (p-value: 1.3 x 10-12 - 1.5x10-36) with IFNα response to in vitro poxvirus stimulation. Examination of these SNPs led to the functional testing of rs1131769, a non-synonymous SNP in TMEM173 causing an Arg-to-His change at position 232 in the STING protein-a major regulator of innate immune responses to viral infections. Our findings demonstrate differences in the ability of the two STING variants to phosphorylate the downstream intermediates TBK1 and IRF3 in response to multiple STING ligands. Further downstream in the STING pathway, we observed significantly reduced expression of type I IFNs (including IFNα) and IFN-response genes in cells carrying the H232 variant. Subsequent molecular modeling of both alleles predicted altered ligand binding characteristics between the two variants, providing a potential mechanism underlying differences in inter-individual responses to poxvirus infection. Our data indicate that possession of the H232 variant may impair STING-mediated innate immunity to poxviruses. These results clarify prior studies evaluating functional effects of genetic variants in TMEM173 and provide novel data regarding genetic control of poxvirus immunity.

TMEM203 is a binding partner and regulator of STING-mediated inflammatory signaling in macrophages

Activators of interferons have received a great deal of interest in recent decades, both due to the central role they play in host defense against a range of pathogens, as well as the now well-recognized importance of dysregulated interferon activation/signaling in the pathogenesis of a number of highly prevalent and hard-to-treat diseases, such as systemic lupus erythematosus (SLE). Therefore, novel regulators of interferon activation are being sought as they may provide better targets to treat these diseases. We report the discovery of TMEM203 as an SLE-associated gene and a regulator of ligand-dependent activation of interferon production via STING. Thus, our work could form the basis of a novel therapeutic strategy for the treatment of interferonopathies, including SLE.

Regulation of IFN signaling is critical in host recognition and response to pathogens while its dysregulation underlies the pathogenesis of several chronic diseases. STimulator of IFN Genes (STING) has been identified as a critical mediator of IFN inducing innate immune pathways, but little is known about direct coregulators of this protein. We report here that TMEM203, a conserved putative transmembrane protein, is an intracellular regulator of STING-mediated signaling. We show that TMEM203 interacts, functionally cooperates, and comigrates with STING following cell stimulation, which in turn leads to the activation of the kinase TBK1, and the IRF3 transcription factor. This induces target genes in macrophages, including IFN-β. Using Tmem203 knockout bone marrow-derived macrophages and transient knockdown of TMEM203 in human monocyte-derived macrophages, we show that TMEM203 protein is required for cGAMP-induced STING activation. Unlike STING, TMEM203 mRNA levels are elevated in T cells from patients with systemic lupus erythematosus, a disease characterized by the overexpression of type I interferons. Moreover, TMEM203 mRNA levels are associated with disease activity, as assessed by serum levels of the complement protein C3. Identification of TMEM203 sheds light into the control of STING-mediated innate immune responses, providing a potential novel mechanism for therapeutic interventions in STING-associated inflammatory diseases.

OMIC Technologies and Vaccine Development: From the Identification of Vulnerable Individuals to the Formulation of Invulnerable Vaccines

Routine vaccination is among the most effective clinical interventions to prevent diseases as it is estimated to save over 3 million lives every year. However, the full potential of global immunization programs is not realised because population coverage is still suboptimal. This is also due to the inadequate immune response and paucity of informative correlates of protection upon immunization of vulnerable individuals such as newborns, preterm infants, pregnant women, and elderly individuals as well as those patients affected by chronic and immune compromising medical conditions. In addition, these groups are undervaccinated for a number of reasons, including lack of awareness of vaccine-preventable diseases and uncertainty or misconceptions about the safety and efficacy of vaccination by parents and healthcare providers. The presence of these nonresponders/undervaccinated individuals represents a major health and economic burden to society, which will become particularly difficult to address in settings with limited public resources. This review describes innovative and experimental approaches that can help identify specific genomic profiles defining nonresponder individuals for whom specific interventions might be needed. We will provide examples that show how such information can be useful to identify novel biomarkers of safety and immunogenicity for future vaccine trials. Finally, we will discuss how system biology “OMICs” data can be used to design bioinformatic tools to predict the vaccination outcome providing genetic and molecular “signatures” of protective immune response. This strategy may soon enable identification of signatures highly predictive of vaccine safety, immunogenicity, and efficacy/protection thereby informing personalized vaccine interventions in vulnerable populations.

Personalized vaccinology: A review

At the current time, the field of vaccinology remains empirical in many respects. Vaccine development, vaccine immunogenicity, and vaccine efficacy have, for the most part, historically been driven by an empiric “isolate-inactivate-inject” paradigm. In turn, a population-level public health paradigm of “the same dose for everyone for every disease” model has been the normative thinking in regard to prevention of vaccine-preventable infectious diseases. In addition, up until recently, no vaccines had been designed specifically to overcome the immunosenescence of aging, consistent with a post-WWII mentality of developing vaccines and vaccine programs for children. It is now recognized that the current lack of knowledge concerning how immune responses to vaccines are generated is a critical barrier to understanding poor vaccine responses in the elderly and in immunoimmaturity, discovery of new correlates of vaccine immunogenicity (vaccine response biomarkers), and a directed approach to new vaccine development.

The new fields of vaccinomics and adversomics provide models that permit global profiling of the innate, humoral, and cellular immune responses integrated at a systems biology level. This has advanced the science beyond that of reductionist scientific approaches by revealing novel interactions between and within the immune system and other biological systems (beyond transcriptional level), which are critical to developing “downstream” adaptive humoral and cellular responses to infectious pathogens and vaccines. Others have applied systems level approaches to the study of antibody responses (a.k.a. “systems serology”), [1] high-dimensional cell subset immunophenotyping through CyTOF, [2,3] and vaccine induced metabolic changes [4]. In turn, this knowledge is being utilized to better understand the following: identifying who is at risk for which infections; the level of risk that exists regarding poor immunogenicity and/or serious adverse events; and the type or dose of vaccine needed to fully protect an individual. In toto, such approaches allow for a personalized approach to the practice of vaccinology, analogous to the substantial inroads that individualized medicine is playing in other fields of human health and medicine. Herein we briefly review the field of vaccinomics, adversomics, and personalized vaccinology.

Integration of Immune Cell Populations, mRNA-Seq, and CpG Methylation to Better Predict Humoral Immunity to Influenza Vaccination: Dependence of mRNA-Seq/CpG Methylation on Immune Cell Populations

The development of a humoral immune response to influenza vaccines occurs on a multisystems level. Due to the orchestration required for robust immune responses when multiple genes and their regulatory components across multiple cell types are involved, we examined an influenza vaccination cohort using multiple high-throughput technologies. In this study, we sought a more thorough understanding of how immune cell composition and gene expression relate to each other and contribute to interindividual variation in response to influenza vaccination. We first hypothesized that many of the differentially expressed (DE) genes observed after influenza vaccination result from changes in the composition of participants' peripheral blood mononuclear cells (PBMCs), which were assessed using flow cytometry. We demonstrated that DE genes in our study are correlated with changes in PBMC composition. We gathered DE genes from 128 other publically available PBMC-based vaccine studies and identified that an average of 57% correlated with specific cell subset levels in our study (permutation used to control false discovery), suggesting that the associations we have identified are likely general features of PBMC-based transcriptomics. Second, we hypothesized that more robust models of vaccine response could be generated by accounting for the interplay between PBMC composition, gene expression, and gene regulation. We employed machine learning to generate predictive models of B-cell ELISPOT response outcomes and hemagglutination inhibition (HAI) antibody titers. The top HAI and B-cell ELISPOT model achieved an area under the receiver operating curve (AUC) of 0.64 and 0.79, respectively, with linear model coefficients of determination of 0.08 and 0.28. For the B-cell ELISPOT outcomes, CpG methylation had the greatest predictive ability, highlighting potentially novel regulatory features important for immune response. B-cell ELISOT models using only PBMC composition had lower performance (AUC = 0.67), but highlighted well-known mechanisms. Our analysis demonstrated that each of the three data sets (cell composition, mRNA-Seq, and DNA methylation) may provide distinct information for the prediction of humoral immune response outcomes. We believe that these findings are important for the interpretation of current omics-based studies and set the stage for a more thorough understanding of interindividual immune responses to influenza vaccination.

Genome-wide associations of CD46 and IFI44L genetic variants with neutralizing antibody response to measles vaccine

Population-based studies have revealed 2-10% measles vaccine failure rate even after two vaccine doses. While the mechanisms behind this remain unknown, we hypothesized that host genetic factors are likely to be involved. We performed a genome-wide association study of measles specific neutralizing antibody and IFNγ ELISPOT response in a combined sample of 2872 subjects. We identified two distinct chromosome 1 regions (previously associated with MMR-related febrile seizures), associated with vaccine-induced measles neutralizing antibody titers. The 1q32 region contained 20 significant SNPs in/around the measles virus receptor-encoding CD46 gene, including the intronic rs2724384 (p value = 2.64 × 10^-09) and rs2724374 (p value = 3.16 × 10^-09) SNPs. The 1q31.1 region contained nine significant SNPs in/around IFI44L, including the intronic rs1333973 (p value = 1.41 × 10^-10) and the missense rs273259 (His73Arg, p value = 2.87 × 10^-10) SNPs. Analysis of differential exon usage with mRNA-Seq data and RT-PCR suggests the involvement of rs2724374 minor G allele in the CD46 STP region exon B skipping, resulting in shorter CD46 isoforms. Our study reveals common CD46 and IFI44L SNPs associated with measles-specific humoral immunity, and highlights the importance of alternative splicing/virus cellular receptor isoform usage as a mechanism explaining inter-individual variation in immune response after live measles vaccine.

Statistical Methods for Testing Genetic Pleiotropy

Genetic pleiotropy is when a single gene influences more than one trait. Detecting pleiotropy and understanding its causes can improve the biological understanding of a gene in multiple ways, yet current multivariate methods to evaluate pleiotropy test the null hypothesis that none of the traits are associated with a variant; departures from the null could be driven by just one associated trait. A formal test of pleiotropy should assume a null hypothesis that one or no traits are associated with a genetic variant. For the special case of two traits, one can construct this null hypothesis based on the intersection-union (IU) test, which rejects the null hypothesis only if the null hypotheses of no association for both traits are rejected. To allow for more than two traits, we developed a new likelihood-ratio test for pleiotropy. We then extended the testing framework to a sequential approach to test the null hypothesis that [Formula: see text] traits are associated, given that the null of k traits are associated was rejected. This provides a formal testing framework to determine the number of traits associated with a genetic variant, while accounting for correlations among the traits. By simulations, we illustrate the type I error rate and power of our new methods; describe how they are influenced by sample size, the number of traits, and the trait correlations; and apply the new methods to multivariate immune phenotypes in response to smallpox vaccination. Our new approach provides a quantitative assessment of pleiotropy, enhancing current analytic practice.

The Integration of Epistasis Network and Functional Interactions in a GWAS Implicates RXR Pathway Genes in the Immune Response to Smallpox Vaccine

Although many diseases and traits show large heritability, few genetic variants have been found to strongly separate phenotype groups by genotype. Complex regulatory networks of variants and expression of multiple genes lead to small individual-variant effects and difficulty replicating the effect of any single variant in an affected pathway. Interaction network modeling of GWAS identifies effects ignored by univariate models, but population differences may still cause specific genes to not replicate. Integrative network models may help detect indirect effects of variants in the underlying biological pathway. In this study, we used gene-level functional interaction information from the Integrative Multi-species Prediction (IMP) tool to reveal important genes associated with a complex phenotype through evidence from epistasis networks and pathway enrichment. We test this method for augmenting variant-based network analyses with functional interactions by applying it to a smallpox vaccine immune response GWAS. The integrative analysis spotlights the role of genes related to retinoid X receptor alpha (RXRA), which has been implicated in a previous epistasis network analysis of smallpox vaccine.

An interaction quantitative trait loci tool implicates epistatic functional variants in an apoptosis pathway in smallpox vaccine eQTL data

Expression quantitative trait loci (eQTL) studies have functionalized nucleic acid variants through the regulation of gene expression. Although most eQTL studies only examine the effects of single variants on transcription, a more complex process of variant-variant interaction (epistasis) may regulate transcription. Herein, we describe a tool called interaction QTL (iQTL) designed to efficiently detect epistatic interactions that regulate gene expression. To maximize biological relevance and minimize the computational and hypothesis testing burden, iQTL restricts interactions such that one variant is within a user-defined proximity of the transcript (cis-regulatory). We apply iQTL to a data set of 183 smallpox vaccine study participants with genome-wide association study and gene expression data from unstimulated samples and samples stimulated by inactivated vaccinia virus. While computing only 0.15% of possible interactions, we identify 11 probe sets whose expression is regulated through a variant-variant interaction. We highlight the functional epistatic interactions among apoptosis-related genes, DIABLO, TRAPPC4 and FADD, in the context of smallpox vaccination. We also use an integrative network approach to characterize these iQTL interactions in a posterior network of known prior functional interactions. iQTL is an efficient, open-source tool to analyze variant interactions in eQTL studies, providing better understanding of the function of epistasis in immune response and other complex phenotypes.

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.

Cytokine production associated with smallpox vaccine responses

Smallpox was eradicated 34 years ago due to the success of the smallpox vaccine; yet, the vaccine continues to be studied because of its importance in responding to potential biological warfare and the adverse events associated with current smallpox vaccines. Interindividual variations in vaccine response are observed and are, in part, due to genetic variation. In some cases, these varying responses lead to adverse events, which occur at a relatively high rate for the smallpox vaccine compared with other vaccines. Here, we aim to summarize the cytokine responses associated with smallpox vaccine response to date. Along with a description of each of these cytokines, we describe the genetic and adverse event data associated with cytokine responses to smallpox vaccination.

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.

The genetic regulation of infant immune responses to vaccination

A number of factors are recognized to influence immune responses to vaccinations including age, gender, the dose, and quality of the antigen used, the number of doses given, the route of administration, and the nutritional status of the recipient. Additionally, several immunogenetic studies have identified associations between polymorphisms in genes encoding immune response proteins, both innate and adaptive, and variation in responses to vaccines. Variants in the genes encoding Toll-like receptors, HLA molecules, cytokines, and cytokine receptors have associated with heterogeneity of responses to a wide range of vaccines including measles, hepatitis B, influenza A, BCG, Haemophilus influenzae type b, and certain Neisseria meningitidis serotypes, amongst others. However, the vast majority of these studies have been conducted in older children and adults and there are very few data available from studies conducted in infants. This paper reviews the evidence to date that host genes influencing vaccines responses in these older population and identifies a large gap in our understanding of the genetic regulation of responses in early life. Given the high mortality from infection in early life and the challenges of developing vaccines that generate effective immune responses in the context of the developing immune system further research on infant populations is required.

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.

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.).

Quantitative genetics in the study of virus-induced disease

While the role of viral variants has long been known to play a key role in causing variation in disease severity, it is also clear that host genetic variation plays a critical role in determining virus-induced disease responses. However, a variety of factors, including confounding environmental variables, rare genetic variants requiring extremely large cohorts, the temporal dynamics of infections, and ethical limitation on human studies, have made the identification and dissection of variant host genes and pathways difficult within human populations. This difficulty has led to the development of a variety of experimental approaches used to identify host genetic contributions to disease responses. In this chapter, we describe the history of genetic associations within the human population, the development of experimentally tractable systems, and the insights these specific approaches provide. We conclude with a discussion of recent advances that allow for the investigation of the role of complex genetic networks that underlie host responses to infection, with the goal of drawing connections to human infections. In particular, we highlight the need for robust animal models with which to directly control and assess the role of host genetics on viral infection outcomes.

Changing face of vaccination in immunocompromised hosts

Infection prevention is a key component of care and an important determinant of clinical outcomes in a diverse population of immunocompromised hosts. Vaccination remains a fundamental preventative strategy, and clear guidelines exist for the vaccination of immunocompromised individuals and close contacts. Unfortunately, adherence to such guidelines is frequently suboptimal, with consequent missed opportunities to prevent infection. Additionally, vaccination of immunocompromised individuals is known to produce responses inferior to those observed in immunocompetent hosts. Multiple factors contribute to this finding, and developing improved vaccination strategies for those at high risk of infectious complications remains a priority of care providers. Herein, we review potential factors contributing to vaccine outcomes, focusing on host immune responses, and propose a means for applying modern, innovative systems biology technology to model critical determinants of vaccination success. With influenza vaccine in solid organ transplants used as a case in point, novel means for stratifying individuals using a host "immunophenotype" are explored, and strategies for individualizing vaccine approaches tailored to safely optimize vaccine responses in those most at risk are discussed.

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.

Associations between race, sex and immune response variations to rubella vaccination in two independent cohorts

INTRODUCTION

Immune response variations after vaccination are influenced by host genetic factors and demographic variables, such as race, ethnicity and sex. The latter have not been systematically studied in regard to live rubella vaccine, but are of interest for developing next generation vaccines for diverse populations, for predicting immune responses after vaccination, and for better understanding the variables that impact immune response.

METHODS

We assessed associations between demographic variables, including race, ethnicity and sex, and rubella-specific neutralizing antibody levels and secreted cytokines (IFNγ, IL-6) in two independent cohorts (1994 subjects), using linear and linear mixed models approaches, and genetically defined racial and ethnic categorizations.

RESULTS

Our replicated findings in two independent, large, racially diverse cohorts indicate that individuals of African descent have significantly higher rubella-specific neutralizing antibody levels compared to individuals of European descent and/or Hispanic ethnicity (p<0.001).

CONCLUSION

Our study provides consistent evidence for racial/ethnic differences in humoral immune response following rubella vaccination.

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.

Race and sex-based differences in cytokine immune responses to smallpox vaccine in healthy individuals

We assessed the effects of sex, race and ethnicity on smallpox vaccine-induced immune responses in 1071 armed forces members after primary Dryvax(®) smallpox vaccination, including 790 males and 281 females; 580 Caucasians, 217 African-Americans, and 217 Hispanics. Analysis of vaccinia-specific cytokine responses revealed that Caucasians had higher total IFNγ ELISPOT responses (median 57 spot-forming units/SFUs per 200,000 cells, p=0.01) and CD8(+)IFNγ ELISPOT responses (12 SFUs, p<0.001) than African-Americans (51 and 4 SFUs, respectively) and Hispanics (47 and 8 SFUs, respectively). Similarly, Caucasians secreted higher levels of vaccinia-specific IL-2 (p=0.003) and IFNα (p<0.001) compared to other racial/ethnic groups. Males had higher total IFNγ ELISPOT responses (median 55 SFUs) compared to females (41 SFUs, p<0.001). We observed statistically significant sex-related differences in the secretion of IL-2 (p<0.001), IL-1β (p<0.001) and IL-10 (p=0.017). These data suggest that vaccinia-specific cytokine responses following primary smallpox vaccination are significantly influenced by race and sex of vaccinees.

Genetic variation in IL18R1 and IL18 genes and Inteferon γ ELISPOT response to smallpox vaccination: an unexpected relationship

BACKGROUND

Genetic association studies demonstrated a role for cytokine proteins and cytokine or cytokine receptor gene polymorphisms in smallpox vaccine-induced adaptive immunity.

METHODS

We examined the association of genetic polymorphisms with cellular (interferon [IFN] γ enzyme-linked immunospot assay [ELISPOT]) immune response to smallpox vaccine in 1076 immunized individuals.

RESULTS

The majority of significant associations were discovered between single-nucleotide polymorphisms/haplotypes in IL18R1 and IL18 genes, in which we previously reported an association with vaccinia virus-induced neutralizing antibody titers in this study cohort. A functional coding IL18R1 polymorphism (rs1035130/Phe251Phe; P = .01) was significantly associated with an allele dose-related increase in IFN-γ production and was also associated with vaccinia-specific neutralizing antibody titers. Significant associations were also found between IL18R1 haplotypes and variations in IFN-γ ELISPOT responses (global P < .0001).

CONCLUSIONS

Our data suggest the importance of variants in the IL18R1 and IL18 genetic loci for broad-based smallpox vaccine-induced adaptive immunity.

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.

NK cells and poxvirus infection

In recent years, our understanding of the role of natural killer (NK) cells in the response to viral infection has grown rapidly. Not only do we realize viruses have many immune-evasion strategies to escape NK cell responses, but that stimulation of NK cell subsets during an antiviral response occurs through receptors seemingly geared directly at viral products and that NK cells can provide a memory response to viral pathogens. Tremendous knowledge has been gained in this area through the study of herpes viruses, but appreciation for the significance of NK cells in the response to other types of viral infections is growing. The function of NK cells in defense against poxviruses has emerged over several decades beginning with the early seminal studies showing the role of NK cells and the NK gene complex in susceptibility of mouse strains to ectromelia, a poxvirus pathogen of mice. More recently, greater understanding has emerged of the molecular details of the response. Given that human diseases caused by poxviruses can be as lethal as smallpox or as benign as Molluscum contagiosum, and that vaccinia virus, the prototypic member of the pox family, persists as a mainstay of vaccine design and has potential as an oncolytic virus for tumor therapy, further research in this area remains important. This review focuses on recent advances in understanding the role of NK cells in the immune response to poxviruses, the receptors involved in activation of NK cells during poxvirus infection, and the viral evasion strategies poxviruses employ to avoid the NK response.

A genome-wide perspective of human diversity and its implications in infectious disease

Progress in genomic technologies, such as DNA arrays and next-generation sequencing, is allowing systematic characterization of the degree of human genetic variation at the scale of individual genomes. Public efforts, such as the International HapMap Project and the 1000 Genomes Project, have provided a realistic picture of the levels of genetic diversity in individuals and populations. These genomic techniques are also making it possible to evaluate the contribution of host genetic diversity to differences in susceptibility to both rare and common infectious diseases. Recent studies have revealed the power of whole-exome sequencing for dissecting the immunological mechanisms underlying the pathogenesis of severe, rare infectious diseases. Likewise, genome-wide association studies on common viral, bacterial, and parasitic infections have shed light on the host genetic basis of susceptibility to infectious diseases and, in some cases, of disease progression and drug responses.

Human genetics and respiratory syncytial virus disease: current findings and future approaches

Infection with respiratory syncytial virus (RSV) can result in a wide spectrum of pulmonary manifestations, from mild upper respiratory symptoms to severe bronchiolitis and pneumonia. Although there are several known risk factors for severe RSV disease, namely, premature birth, chronic lung disease, congenital heart disease, and T cell immunodeficiency, the majority of young children who develop severe RSV disease are otherwise healthy children. Genetic susceptibility to RSV infection is emerging as a complex trait, in which many different host genetic variants contribute to risk for distinct disease manifestations. Initially, host genetic studies focused on severe RSV disease using the candidate gene approach to interrogate common single nucleotide polymorphisms (SNPs). Many studies have reported genetic associations between severe RSV bronchiolitis and SNPs in genes within plausible biological pathways, such as in innate host defense genes (SPA, SPD, TLR4, and VDR), cytokine or chemokine response genes (CCR5, IFN, IL6, IL10, TGFB1), and altered Th1/Th2 immune responses (IL4, IL13). Due to the complexity of RSV susceptibility, genome studies done on a larger scale, such as genome-wide association studies have certainly identified more of the host factors that contribute to the development of severe RSV bronchiolitis or excessive pathology. Furthermore, whole-genome approaches can reveal robust associations between genetic markers and RSV disease susceptibility. Recent introduction of 'exome' genotyping or sequencing, which specifically analyzes the majority of coding variants, should be fruitful in sufficiently large, well-powered studies. The advent of new genomic technologies together with improved computational tools offer the promise of interrogating the host genome in search of genetic factors, rare, uncommon, or common that should give new insights into the underlying biology of susceptibility to or protection from severe RSV infection. Careful assessment of novel pathways and further identification of specific genes could identify new approaches for vaccine development and perhaps lead to effective risk modeling.