Comprehensive evaluation of measles-specific humoral immunity after vaccination is important for determining new and/or additional correlates of vaccine immunogenicity and efficacy.

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.

Older adults experience declining influenza vaccine-induced immunity and are at higher risk of influenza and its complications. For this reason, high dose (e.g., Fluzone) and adjuvanted (e.g., Fluad) vaccines are preferentially recommended for people age 65 years and older. However, T cell transcriptional activity shaping the humoral immune responses to Fluzone and Fluad vaccines in older adults is still poorly understood. We designed a study of 234 older adults (≥65 years old) who were randomly allocated to receive Fluzone or Fluad vaccine and provided blood samples at baseline and at Day 28 after immunization. We measured the humoral immune responses (hemagglutination inhibition/HAI antibody titer) to influenza A/H3N2 and performed mRNA-Seq transcriptional profiling in purified CD4+ T cells, in order to identify T cell signatures that might explain differences in humoral immune response by vaccine type. Given the large differences in formulation (higher antigen dose vs adjuvant), our hypothesis was that each vaccine elicited a distinct transcriptomic response after vaccination. Thus, the main focus of our study was to identify the differential gene expression influencing the antibody titer in the two vaccine groups. Our analyses identified three differentially expressed, functionally linked genes/proteins in CD4+ T cells: the calcium/calmodulin dependent serine/threonine kinase IV (CaMKIV); its regulator the TMEM38B/transmembrane protein 38B, involved in maintenance of intracellular Ca2+ release; and the transcriptional coactivator CBP/CREB binding protein, as regulators of transcriptional activity/function in CD4+ T cells that impact differences in immune response by vaccine type. Significantly enriched T cell-specific pathways/biological processes were also identified that point to the importance of genes/proteins involved in Th1/Th2 cell differentiation, IL-17 signaling, calcium signaling, Notch signaling, MAPK signaling, and regulation of TRP cation Ca2+ channels in humoral immunity after influenza vaccination. In summary, we identified the genes/proteins and pathways essential for cell activation and function in CD4+ T cells that are associated with differences in influenza vaccine-induced humoral immunity by vaccine type. These findings provide an additional mechanistic perspective for achieving protective immunity in older adults.

While waning immunity and SARS-CoV-2 variant immune escape continue to result in high infection rates worldwide, associations between longitudinal quantitative, qualitative, and functional humoral immune responses after SARS-CoV-2 infection remain unclear. In this study, we found significant waning of antibody against Spike S1 (R = -0.32, p = 0.035) and N protein (R = -0.39, p = 0.008), while RBD antibody moderately decreased (R = -0.19, p = 0.203). Likewise, neutralizing antibody titer (ND50) waned over time (R = -0.46, p = 0.001). In contrast, antibody avidity increased significantly over time for Spike S1 (R = 0.62, p = 6.0e-06), RBD (R = 0.54, p = 2.0e-04), and N (R = 0.33, p = 0.025) antibodies. Across all humoral responses, ND50 strongly associated with Spike S1 (R = 0.85, p = 2.7e-13) and RBD (R = 0.78, p = 2.9e-10) antibodies. Our findings provide longitudinal insight into humoral immune responses after infection and imply the potential of Spike S1/RBD antibody titer as surrogate correlates of protection.

Failure to achieve a protected state after influenza vaccination is poorly understood but occurs commonly among aged populations experiencing greater immunosenescence. In order to better understand immune response in the elderly, we studied epigenetic and transcriptomic profiles and humoral immune response outcomes in 50-74 year old healthy participants. Associations between DNA methylation and gene expression reveal a system-wide regulation of immune-relevant functions, likely playing a role in regulating a participant's propensity to respond to vaccination. Our findings show that sites of methylation regulation associated with humoral response to vaccination impact known cellular differentiation signaling and antigen presentation pathways. We performed our analysis using per-site and regionally average methylation levels, in addition to continuous or dichotomized outcome measures. The genes and molecular functions implicated by each analysis were compared, highlighting different aspects of the biologic mechanisms of immune response affected by differential methylation. Both cis-acting (within the gene or promoter) and trans-acting (enhancers and transcription factor binding sites) sites show significant associations with measures of humoral immunity. Specifically, we identified a group of CpGs that, when coordinately hypo-methylated, are associated with lower humoral immune response, and methylated with higher response. Additionally, CpGs that individually predict humoral immune responses are enriched for polycomb-group and FOXP2 transcription factor binding sites. The most robust associations implicate differential methylation affecting gene expression levels of genes with known roles in immunity (e.g. HLA-B and HLA-DQB2) and immunosenescence. We believe our data and analysis strategy highlight new and interesting epigenetic trends affecting humoral response to vaccination against influenza; one of the most common and impactful viral pathogens.

To assess gene signatures related to humoral response among healthy older subjects following seasonal influenza vaccination, we studied 94 healthy adults (50-74 years old) who received one documented dose of licensed trivalent influenza vaccine containing the A/California/7/2009 (H1N1)-like virus strain. Influenza-specific antibody (HAI) titer in serum samples and next-generation sequencing on PBMCs were performed using blood samples collected prior to (Day 0) and at two timepoints after (Days 3 and 28) vaccination. We identified a number of uncharacterized genes (ZNF300, NUP1333, KLK1 and others) and confirmed previous studies demonstrating specific genes/genesets that are important mediators of host immune responses and that displayed associations with antibody response to influenza A/H1N1 vaccine. These included interferon-regulatory transcription factors (IRF1/IRF2/IRF6/IRF7/IRF9), chemokine/chemokine receptors (CCR5/CCR9/CCL5), cytokine/cytokine receptors (IFNG/IL10RA/TNFRSF1A), protein kinases (MAP2K4/MAPK3), growth factor receptor (TGFBR1). The identification of gene signatures associated with antibody response represents an early stage in the science for which further research is needed. Such research may assist in the design of better vaccines to facilitate improved defenses against new influenza virus strains, as well as better understanding the genetic drivers of immune responses.

There are insufficient system-wide transcriptomic (or other) data that help explain the observed inter-individual variability in antibody titers after measles vaccination in otherwise healthy individuals.

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.

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.

Host antiviral genes are important regulators of antiviral immunity and plausible genetic determinants of immune response heterogeneity after vaccination. We genotyped and analyzed 307 common candidate tagSNPs from 12 antiviral genes in a cohort of 745 schoolchildren immunized with two doses of measles-mumps-rubella (MMR) vaccine. Associations between SNPs/haplotypes and measles virus-specific immune outcomes were assessed using linear regression methodologies in Caucasians and African-Americans. Genetic variants within the DDX58/RIG-I gene, including a coding polymorphism (rs3205166/Val800Val), were associated as single-SNPs (p≤0.017; although these SNPs did not remain significant after correction for false discovery rate/FDR) and in haplotype-level analysis, with measles-specific antibody variations in Caucasians (haplotype allele p-value=0.021; haplotype global p-value=0.076). Four DDX58 polymorphisms, in high LD, demonstrated also associations (after correction for FDR) with variations in both measles-specific IFN-γ and IL-2 secretion in Caucasians (p≤0.001, q=0.193). Two intronic OAS1 polymorphisms, including the functional OAS1 SNP rs10774671 (p=0.003), demonstrated evidence of association with a significant allele-dose-related increase in neutralizing antibody levels in African-Americans. Genotype and haplotype-level associations demonstrated the role of ADAR genetic variants, including a non-synonymous SNP (rs2229857/Arg384Lys; p=0.01), in regulating measles virus-specific IFN-γ Elispot responses in Caucasians (haplotype global p-value=0.017). After correction for FDR, 15 single-SNP associations (11 SNPs in Caucasians and 4 SNPs in African-Americans) still remained significant at the q-value<0.20. In conclusion, our findings strongly point to genetic variants/genes, involved in antiviral sensing and antiviral control, as critical determinants, differentially modulating the adaptive immune responses to live attenuated measles vaccine in Caucasians and African-Americans.

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.

Despite high levels of MMR-II usage in the US, mumps outbreaks continue to occur. Evidence suggests that mumps vaccine-induced humoral immunity wanes over time. Relatively few studies have examined cell-mediated immunity or reported on sex-based differences. To better understand sex-based differences in the immune response to mumps vaccine, we measured neutralizing antibody titers and mumps-specific cytokine/chemokine responses in a cohort of 748 adolescents and young adults after two doses of MMR vaccine. We observed significantly higher neutralizing antibody titers in females than in males (120.8 IU/mL, 98.7 IU/mL, p = 0.038) but significantly higher secretion levels of MIP-1α, MIP-1β, TNFα, IL-6, IFNγ, and IL-1β in males compared to females. These data demonstrate that sex influences mumps-specific humoral and cell-mediated immune response outcomes, a phenomenon that should be considered during efforts to improve vaccines and prevent future outbreaks.

The lack of standardization of the currently used commercial anti-rubella IgG antibody assays leads to frequent misinterpretation of results for samples with low/equivocal antibody concentration. The use of alternative approaches in rubella serology could add new information leading to a fuller understanding of rubella protective immunity and neutralizing antibody response after vaccination.

Although influenza causes significant morbidity and mortality in the elderly, the factors underlying the reduced vaccine immunogenicity and efficacy in this age group are not completely understood. Age and immunosenescence factors, and their impact on humoral immunity after influenza vaccination, are of growing interest for the development of better vaccines for the elderly.

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(-15), FDR<1.00E(-14)), 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(-08) 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.

Interferon-induced antiviral genes are crucial players in innate antiviral defense and potential determinants of immune response heterogeneity. We selected 114 candidate single-nucleotide polymorphisms (SNPs) from 12 antiviral genes using an LD tagSNP selection approach and genotyped them in a cohort of 738 school children immunized with two doses of rubella vaccine. Associations between SNPs/haplotypes and rubella virus-specific immune measures were assessed using linear regression methodologies. We identified 23 significant associations (p < 0.05) between polymorphisms within the 2'-5'-oligoadenylate synthetase (OAS) gene cluster, and rubella virus-specific IL-2, IL-10, IL-6 secretion, and antibody levels. The minor allele variants of three OAS1 SNPs (rs3741981/Ser162Gly, rs1051042/Thr361Arg, rs2660), located in a linkage disequilibrium block of functional importance, were significantly associated with an increase in rubella virus-specific IL-2/T(h)1 response (p

The goal of annual influenza vaccination is to reduce mortality and morbidity associated with this disease through the generation of protective immune responses. The objective of the current study was to examine markers of immunosenescence and identify immunosenescence-related differences in gene expression, gene regulation, cytokine secretion, and immunologic changes in an older study population receiving seasonal influenza A/H1N1 vaccination. Surprisingly, prior studies in this cohort revealed weak correlations between immunosenescence markers and humoral immune response to vaccination. In this report, we further examined the relationship of each immunosenescence marker (age, T cell receptor excision circle frequency, telomerase expression, percentage of CD28- CD4+ T cells, percentage of CD28- CD8+ T cells, and the CD4/CD8 T cell ratio) with additional markers of immune response (serum cytokine and chemokine expression) and measures of gene expression and/or regulation. Many of the immunosenescence markers indeed correlated with distinct sets of individual DNA methylation sites, miRNA expression levels, mRNA expression levels, serum cytokines, and leukocyte subsets. However, when the individual immunosenescence markers were grouped by pathways or functional terms, several shared biological functions were identified: antigen processing and presentation pathways, MAPK, mTOR, TCR, BCR, and calcium signaling pathways, as well as key cellular metabolic, proliferation and survival activities. Furthermore, the percent of CD4+ and/or CD8+ T cells lacking CD28 expression also correlated with miRNAs regulating clusters of genes known to be involved in viral infection. Integrated (DNA methylation, mRNA, miRNA, and protein levels) network biology analysis of immunosenescence-related pathways and genesets identified both known pathways (e.g., chemokine signaling, CTL, and NK cell activity), as well as a gene expression module not previously annotated with a known function. These results may improve our ability to predict immune responses to influenza and aid in new vaccine development, and highlight the need for additional studies to better define and characterize immunosenescence.

Human antibody response to measles vaccine is highly variable in the population. Host genes contribute to inter-individual antibody response variation. The killer cell immunoglobulin-like receptors (KIR) are recognized to interact with HLA molecules and possibly influence humoral immune response to viral antigens. To expand on and improve our previous work with HLA genes, and to explore the genetic contribution of KIR genes to the inter-individual variability in measles vaccine-induced antibody responses, we performed a large population-based study in 2,506 healthy immunized subjects (ages 11 to 41 years) to identify HLA and KIR associations with measles vaccine-induced neutralizing antibodies. After correcting for the large number of statistical tests of allele effects on measles-specific neutralizing antibody titers, no statistically significant associations were found for either HLA or KIR loci. However, suggestive associations worthy of follow-up in other cohorts include B*57:01, DQB1*06:02, and DRB1*15:05 alleles. Specifically, the B*57:01 allele (1,040 mIU/mL; p = 0.0002) was suggestive of an association with lower measles antibody titer. In contrast, the DQB1*06:02 (1,349 mIU/mL; p = 0.0004) and DRB1*15:05 (2,547 mIU/mL; p = 0.0004) alleles were suggestive of an association with higher measles antibodies. Notably, the associations with KIR genotypes were strongly nonsignificant, suggesting that KIR loci in terms of copy number and haplotypes are not likely to play a major role in antibody response to measles vaccination. These findings refine our knowledge of the role of HLA and KIR alleles in measles vaccine-induced immunity.

MicroRNAs are important mediators of post-transcriptional regulation of gene expression through RNA degradation and translational repression, and are emerging biomarkers of immune system activation/response after vaccination.

The durability of protective humoral immunity after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccination and infection is largely dependent on the generation and persistence of antigen-specific isotype-switched memory B cells (MBCs) and long-lived plasma cells that reside in the bone marrow and secrete high-affinity neutralizing antibodies. The reactivity of vaccine-induced MBCs to emerging clinically significant SARS-CoV-2 variants of concern (VoCs) is largely unknown. In a longitudinal cohort study (up to 6 months following coronavirus disease 2019 messenger RNA vaccination), we measured MBCs in concert with other functional antibody measures. We found statistically significant differences between the frequencies of MBCs responding to homologous and VoC (Beta, Gamma, and Delta) receptor-binding domains after vaccination that persisted over time. In concert with a waning antibody response, the reduced MBC response to VoCs could translate to a weaker subsequent recall immune response and increased susceptibility to the emerging SARS-CoV-2 variant strains after vaccination.