Genetic associations with a fever after measles-containing vaccines

ABSTRACT

Children have elevated fever risk 1 to 2 weeks after the first dose of a measles-containing vaccine (MCV), which is likely affected by genetic, immunologic, and clinical factors. Fever after MCV is associated with febrile seizures, though may also be associated with higher measles antibody titers. This exploratory study investigated genetic and immunologic associations with a fever after MCV. Concurrent with a randomized Phase 3 clinical trial of 12–15-month-olds who received their first measles-mumps-rubella (MMR) vaccine in which parents recorded post-vaccination temperatures daily, we consented a subset to collect additional blood and performed human leukocyte antigens (HLA) typing. Association between fever 5–12 days after MMR (“MMR-associated”) and HLA type was assessed using logistic regression. We compared 42-day post-vaccination geometric mean titers (GMT) to measles between children who did and did not have fever using a t-test. We enrolled 86 children and performed HLA typing on 82; 13 (15.1%) had MMR-associated fever. Logistic regressions identified associations between MMR-associated fever and HLA Class I loci A-29:02 (P = .036), B-57:01 (P = .018), C-06:02 (P = .006), C-14:02 (P = .022), and Class II loci DRB1-15 (P = .045). However, Bonferroni's adjustment for multiple comparisons suggests that these associations could have been due to chance. Ninety-eight percent of children had protective antibody titers to measles; however, GMT was higher among those with fever compared with children without fever (P = .006). Fever after the measles vaccine correlated with genetic factors and higher immune response. This study suggests a possible genetic susceptibility to MMR-associated fever.

KEYWORDS:

• MMR
• MMRV
• vaccine
• genetic factors
• fever
• antibody
• titers

Introduction

Fever is a common nonspecific immune response to infection and immunization and may have a beneficial role in response to infection.^1,2 However, in the context of immunization, fever is considered an adverse event that leads to increased medical visits.³ Parental concerns about vaccine adverse events, even mild ones such as fever, have been associated with parental reluctance to vaccinate.^4,5 Parental vaccine hesitancy negatively impacts public health,⁶ as evidenced by 2018–19 measles outbreak in the United States which saw the largest number of cases since 1992 and was mainly attributable to under vaccination.⁷ Elucidating a positive role for fever after vaccination could reassure parents and may help address vaccine hesitancy.

Several lines of evidence suggest that genetic factors influence susceptibility to fever following the first dose of MCV. MCV predictably causes fever 7–10 days after vaccination in 5%–15% of the vaccinees^8–16 and this fever is strongly associated with febrile seizures 7–10 days after MCV.^8,13,15 A large Danish genome-wide association study demonstrated that susceptibility to measles, mumps, rubella (MMR) vaccine-associated febrile seizures have a genetic basis.¹⁷ A Vaccine Safety Datalink (VSD) study evaluating 946,806 children found that children with siblings who had fever 7–10 days after a MCV was more than 3 times as likely themselves to have fever 7–10 days after a MCV.¹⁸ More recently, we reported that parental history of fever and fever with respiratory infection was also associated with an increased risk of fever in their child 7–10 days after MCV.¹⁹

Fever after MMR may be related to higher immune response,²⁰ and genetics influence the immune response to MMR vaccine. An analysis of 5 MCV clinical trials detected an association between post-vaccination fever and higher immune responses, with the level of post-vaccination measles antibody titers positively associated with fever.²⁰ The heritability of antibody titers in response to measles vaccination varies between 49% and 88.5%.^21–23 Studies have reported associations between differential immune responses to MMR and human leukocyte antigen (HLA) alleles,^24–27 as well as with other cellular receptors that influence measles vaccine responses.^28–36

Together these studies provide support for the hypothesis that genetic and biological factors influence susceptibility to fever following measles vaccination. Analyzing biospecimens previously collected from one-year-olds enrolled in a randomized-controlled MMR vaccine clinical trial,¹⁶ the objectives of this exploratory study were to investigate the association between HLA type and fever after MCV and between fever after MCV and antibody titers.

Materials and methods

Study design and population

This study was conducted by separately enrolling subjects who had already been enrolled in a large-randomized, controlled Phase 3 clinical trial evaluating an investigational MMR vaccine (licensed in >100 countries). The Phase 3 trial enrolled subjects aged 12–15 months to receive either a single dose of an investigational MMR vaccine or of the US licensed MMR vaccine (MMR II, Merck & Co) as described.¹⁶ Subjects had blood samples collected before and 42 days after vaccination. As part of the clinical trial, parents were given thermometers, taught how to measure daily temperatures, and instructed to record the highest temperature each day on a diary card prospectively for 42 days after vaccination. The clinical trial enrolled a total of 1736 subjects globally. The incidence of fever 5–12 days post-vaccination was comparable in the clinical trial, with similar rates of fever and other MMR-specific symptoms reported after both MMR vaccines.¹⁶ Kaiser Permanente Northern California (KPNC), as a participating clinical trial site, enrolled 341 of the 1736 subjects.

Concurrent with the per-protocol activities for the above clinical trial, KPNC separately consented the parents/legally authorized representatives of enrolled children to obtain and store additional vials of blood after collecting the above per-protocol blood samples. Due to resource and time constraints, only a portion of the 341 children were enrolled in this separate study. The goal of this separate study was to investigate genetic and biologic factors related to developing fever associated with the MMR vaccine using stored biospecimens. This current report focuses on the analyses of these stored biospecimens.

The KPNC and Stanford University Institutional Review Boards approved this study.

Biospecimen collection and storage

We obtained approximately 6 ml of blood in an EDTA tube, inverted to mix 8–10 times, and immediately spun the tube for 10 min at 1,000–3000 G’s. We removed the plasma and the buffy coat layers placed each into separate cryovials which were immediately frozen at −20°C, where they remained for several months. The plasma and buffy coat cryovials were subsequently transferred and stored frozen at −80°C, where they remained until they were processed and analyzed.

DNA extraction

DNA from buffy coat specimens were extracted using the Qiagen DNA extraction kit (catalog number 56504) following standard protocol. DNA samples underwent quality control in three ways: NanoDrop, Qubit, and fragment analyzer. Nanodrop readings were used to assess the overall sample quality. Qubit reading was used for quantification (Thermofisher catalog number Q32854) following standard protocol. Fragment analyzer traces were used to assess sample quality following standard protocol. All samples met the standard quality for library synthesis.

HLA library preparation and sequencing

We used 400 ng of intact and not fragmented genomic DNA to generate HLA libraries using Illumina TruSight HLAv2 library prep kit (Catalog # 20000215) following standard protocol. The TruSight HLA v2 sequencing panel covers all commonly typed HLA loci. A single assay results in the sequencing of 11 HLA loci, 3 of which were Class I (HLA A, HLA B and HLA C) and 8 were Class II (DRB1, DRB3, DRB4, DRB5, DPA1, DPB1, DQA1, and DQB1). All the sequencing was done on MiSeq using a 2 × 150 pair end kit following standard protocol.

Antibody titers

Measles-specific immunoglobulin G (IgG) antibodies (mIU/mL) were measured using commercial ELISA kits as described¹⁶ and were provided by the sponsor of the Phase 3 trial. Measles titers were considered protective if the post-vaccination antibody concentration was ≥200 mIU/mL for children who were seronegative before vaccination (defined as ≤150 mIU/mL) as described.¹⁶

Statistical analyses

For all analyses, to be consistent with the original clinical trial,¹⁶ we considered a temperature ≥100.4°F(38°C) between days 5 and 12 recorded on the daily diary cards as a fever, henceforth known as “MMR-associated fever.”

We described the distribution of selected demographic factors (Table 1) and the proportion of each HLA type among children who did and did not have MMR-associated fever. We assessed genotypes according to a dominant inheritance pattern. Participants were classified as a carrier of the allele if at least one allele copy of the allele of interest was present. We determined allele frequencies across diverse populations in the US using the US National Marrow Donor Program report.³⁷ This report presents allele frequencies across population subgroups, including African American, Asian, European American, and Hispanic cohorts in the US. Based on these reported allele frequencies, we calculated the expected proportion of individuals in each subpopulation carrying at least one copy of the allele of interest (carriers) by assuming Hardy-Weinberg Equilibrium and taking the square of the allele frequency (q^2) and adding it to 2*alleleFreq*(1-alleleFreq).³⁸ We defined the expected population carrier frequency range for each allele as the range of the proportion of individuals who would be expected to carry at least one of these alleles across diverse US subpopulations, where the low end of the range reflects the calculated carrier frequency based on the subpopulation where the allele was reported to be least common and where the high end of the range reflects the calculated carrier frequency based on the subpopulation where the allele was reported to be most common. Thus, we expect the carrier frequency in our cohort to fall within that range regardless of the distribution of specific genetic ancestries. We assessed for an association between MMR-associated fever and HLA type using unadjusted logistic regression. We only report HLA genes where ≥50% of the children had two typed alleles.

Table 1. Racial and ethnic demographics among children with and without MMR-associated fever

	Children without MMR-Associated Fever^a N = 66 (%)	Children with MMR-Associated Fever^b N = 13 (%)
Race
White – Caucasian/European Heritage	32 (48.5)	6 (46.2)
African Heritage/African American	7 (10.6)	0
Asian	14 (21.2)	4 (30.8)
American Indian or Alaskan Native	3 (4.5)	2 (15.4)
Multiple Races	10 (15.2)	1 (7.7)
Ethnicity
Non-Hispanic	58 (87.8)	10 (76.9)
Hispanic	8 (12.2)	3 (23.1)

^aDoes not include 7 subjects for whom fever data was missing (3 White, 2 Asian and 2 African Heritage).

^b14.5% had fever after a single dose of an investigational MMR and 16.1% after the US-licensed MMR vaccine.

For the antibody titer analyses, we compared day 42 post-vaccination measles-specific geometric mean titers (GMT) between children who did and did not have MMR-associated fever using a t-test after natural logarithm transformation to achieve normality. Box plots of geometric means were created using the Surveymeans procedure in SAS.

Analyses were performed using SAS, version 9.4 (SAS Institute, Inc., Cary, North Carolina), and R v3.1.2 (R Core Team, 2012) glm package.

P values were two-sided and those ≤0.05 were considered as statistically significant. We adjusted the P values for the multiplicity of hypotheses tested using Bonferroni correction. We regarded the HLA results as exploratory given the small sample size.

Results

We separately consented and enrolled 86 children and, of these, were able to obtain buffy coat specimens on 83/86 (97%). One child subsequently withdrew consent just for HLA typing, for a final of 82 children on whom we performed HLA typing. Among the total 86 children, 13 (15.1%) had MMR-associated fever (14.5% after investigational MMR and 16.1% after MMR II). Racial and ethnic demographics were broadly similar between the fever and no fever groups, except that none of the children with fever was of Black race (regardless of ethnicity) versus 10.6% of those without fever (Table 1). Comparisons were limited due to the small numbers of children.

HLA sequence data were assessed for a total of 187 alleles. It was available for both alleles for 94% of the children for HLA A, 95% for HLA B, 87% for HLA C, 100% for HLA DPA1, 77% for DPB1, 99% for HLA DQA1, 94% HLA DQB1, and 98% for HLA DRB1. We excluded HLA DRB3, DRB4, and DRB5 alleles because <50% of the children had two typed alleles.

Logistic regressions identified associations between MMR-associated fever and 4 HLA Class I loci and 1 Class II loci, with 1 additional Class II loci having a borderline association (Table 2). None of the children with MMR-associated fever were homozygous for these five loci. These five loci occurred more commonly among children with MMR-associated fever than would be expected based on allele frequencies across diverse US populations. For example, we expected that 3%–7% of our study population would have HLA A-29; however, we observed that 25% of the children with MMR-associated fever had HLA A-29, while only 5% of those without fever had HLA A-29 (P = .036). In contrast, HLA C-4 was expected to occur in 15%–33% of the study population and we noted that 27.1% of the children without fever and 8.3% of the children with MMR-associated fever had the HLA C-4 loci (P = .189; Supplemental Table). However, adjusting for the multiplicity of hypotheses tested (n = 71) using the Bonferroni correction suggests that these associations could be due to chance. Out of the 187 total alleles, there were 66 alleles that were identified in at least two children, with at least one of the two being in the group with MMR-associated fever. These 66 were not associated with MMR-associated fever (Supplemental Table).

Table 2. Association between human leukocyte antigen allele and MMR-associated fever

Allele Associated with Fever	Allele frequency across diverse US subpopulations^a (%)	Expected proportion with at least one allele (%)	Proportion among children without MMR-associated Fever (%)^b	Proportion among children with MMR-associated Fever (%)	Odds Ratio^c (95% Confidence Interval)	P value	P value after Bonferroni correction
HLA Class I
A-29:02	1–4	2–7	3/60 (5)	3/12 (25)	6.3 (1.1, 36.4)	0.036	1
B-57:01	0–4	0–7	2/63 (3.2)	3/12 (25)	10.2 (1.5, 69.4)	0.018	1
C-06:02	5–9	10–17	2/55 (3.6)	4/12 (33.3)	13.3 (2.1, 84.5)	0.006	0.426
C-14:02	1–5	2–10	0/65 (0)	2/12 (16.7)	NE (1.63, NE)^d	0.022	1
HLA Class II
DRB1-04:01	1–9	2–17	4/65 (6.2)	3/12 (25)	5.1 (1.0, 26.5)	0.054	1
DRB1-15	2–15	4–27	18/65 (27.7)	7/12 (58.3)	3.7 (1.0, 13.0)	0.045	0.426

^aWe calculated the expected carrier frequency from allele frequencies³⁷ assuming Hardy-Weinberg Equilibrium.³⁸

^bDenominator includes children for whom both HLA type and fever data were available.

^cComparing children with MMR-associated fever vs those without MMR-associated fever.

^dThe lower bound CI was based on the conditional maximum likelihood estimate of the odds ratio. P value was calculated as an exact P value.

Measles antibody titer data were available for 84 children, of whom 82 (98%) had protective titers to measles 42 days after vaccination. Among the children for whom we had both titer and fever data (n = 77), post-vaccination GMT was higher among those with MMR-associated fever (4609, 95% confidence interval [CI] 3629, 5853) than among children who did not have a fever (2918, 95% CI 2318, 3673, P = .006; Figure 1). Higher post-vaccination GMT and MMR-associated fever were also noted among the 341 children enrolled in the trial at KPNC (GMT among children with fever: 4022 [95% CI 3209, 5042] vs children with no fever: 3082 [95% CI 2792, 3402]; P = .032).

Figure 1. Geometric mean titers (GMT) against measles among children with and without fever after MMR

The horizontal line within each box denotes the median GMT for each group.

Discussion

This study investigated whether genetic factors influence susceptibility to fever and if fever is associated with increased immunogenicity to MCV. In this exploratory study which used biospecimens previously collected during a randomized-controlled MMR vaccine clinical trial, we identified 5 HLA alleles associated with a fever after the MMR vaccine. However, after adjusting for the multiplicity of testing these findings could be due to chance alone because of the small sample size. While the current exploratory study does not allow us to draw firm conclusions regarding specific associations between fever after MMR and HLA types due to its small size, and a stringent Bonferroni correction, it does suggest based on effect size with confidence intervals excluding the null value that specific genetic factors correlate with a fever after MCV. Importantly, our study provides additional support for the growing body of evidence that genetics influences susceptibility to fever after the MMR vaccine.

Genetic factors such as HLA type have not previously been directly correlated with a fever after MCV. Although our findings must be interpreted with caution due to the small size of this study, the HLA alleles identified here as potentially being correlated with susceptibility to fever after the MMR vaccine can serve as important targets in future studies. To put our findings in the context of the US population, we also calculated the expected population carrier frequencies to reflect the range of allele frequencies (from least to most common) that have been reported in populations across racially diverse cohorts in the US. These ranges help interpretation of our findings by providing context as to how our study population compared with alleles which are generally more or less common across diverse communities in the US. For example, based on prior studies reflecting many subpopulations, we expected just 3%–7% of our participants to have the HLA A-29 allele.³⁷ In other words, in US populations where the HLA A-29 allele is reported to be most common, we expect to find just 7% of the people with that allele. However, we found that allele in 25% of our participants who experienced fever, which is substantially higher than in the population with the highest reported allele frequency.

Several candidate HLA alleles we identified in association with a fever after the MMR vaccine have been associated previously with a viral response. For example, HLA-B 57:01 has previously been reported to be strongly associated with a hypersensitivity syndrome characterized by fever, malaise, and gastrointestinal symptoms in up to 8% of those starting the guanosine analog antiviral abacavir.^39–42 Another potential HLA allele that may be worthy of future investigations includes HLA-C 06:02 which was seen more than expected among children with fever (33% vs 10%–17% expected). Finally, although the association between HLA-DRB1-04:01 and fever was borderline, HLA-DRB1-04:05 has been associated with resistance to enteric fever.⁴³

These data implicating genes involved with fever are consistent with prior studies on heritability of fever following the measles vaccine. Studies showing a role of genetics in fever response include a VSD study reporting that children with siblings who had fever 7–10 days after a MCV was more than 3 times more likely themselves to have fever 7–10 days after a MCV.¹⁸ Furthermore, a recent study that parental history of fever and fever with respiratory infection were associated with fever in their child 7–10 days after MCV.¹⁹

We directly demonstrated an association between fever after MMR and higher measles titers, providing an important link between a stronger immune response to MMR vaccine and fever after the MMR vaccine. Our results are consistent with a meta-analysis which assessed pre-licensure measles, mumps, rubella, varicella vaccine clinical trials, and reported an association between post-vaccination fever and higher immune responses, with the level of post-vaccination measles antibody titers positively associated with fever.²⁰ Studies showing lower antibody responses in infants given antipyretics prior to vaccination also lend support to the fever pathway aiding in the immunogenicity of vaccines.⁴⁴ Importantly, a lack of MMR-associated fever did not indicate inadequate immune response to MMR as nearly all subjects, regardless of whether they had a fever, were seroprotected against measles. Rather, genetic variation affecting fever after MMR may reflect an improved immune response to measles and suggests that beneficial immune effects may be initiated through the fever pathways.²

The association between fever after MMR and higher measles titers strengthens the evidence that genetics impact susceptibility to fever after the MMR vaccine. Differential immune responses to the MMR vaccine has been reported in associations with genetic variation in the immune system,^21–23 including HLA alleles,^24–26 cytokine and cytokine receptors,²⁸ viral pattern recognition receptors which influence measles vaccine responses,^30–36 and measles cellular receptors (CD46, SLAM, and others).²⁹ A recent large population-based study evaluated measles titers among adolescents and adults 11–41 years of age previously vaccinated at differing times with MMR (presumably most were second MMR doses) identified HLA alleles of B 57:01, DRB1-15:05, and DQB1-06:02 as being associated with circulating antibody titers to MCV.²⁷ This study suggested associations between DQB1*06:02 and DRB1*15:05 with higher measles antibody titers and between B*57:01 alleles with lower titers. Important differences between this report and our current study was that we only included toddlers receiving their first dose of MMR vaccine thereby focusing exclusively on the primary antibody response to MMR. It is possible that genetic factors which may be associated with the primary measles vaccine immune response may differ from those associated with the maintenance of circulating antibodies to the measles vaccine.

Specific genes implicated with MMR-associated febrile seizures in a prior GWAS have also been associated with immunity to the measles vaccine. In particular, variants IFI44L and CD46 involved with MMR-associated febrile seizures¹⁷ have also been associated with variability in measles antibody titers following vaccination in another GWAS.⁴⁵ Notably, the CD46 protein is the cell-surface receptor for the measles virus, and variation in its expression has been associated with variation in measles-specific IgG antibody levels as well as with IL-6, TNFα, and IFNα responses following measles vaccine.²⁹ The IFI44L protein belongs to the interferon-stimulated genes and is likely involved in innate antiviral immunity.⁴⁶ Based on prior GWAS, genetic risk factors for vaccine-related febrile seizures seem distinct from febrile seizures generally and appear driven by the immune response.¹⁷

Strengths of this study were that we collected biospecimens during a randomized-controlled clinical trial which ensured that all children were between 12 and 15 months of age, were all receiving their first dose of MMR, all antibody titer data were obtained 42 days post-vaccination, and importantly, all post-vaccination fever data were prospectively recorded. Our study included data on genetic, fever, and antibody titer level and we were able to analyze all three factors in the same population.

Our study had notable limitations, the most important of which was its small size. For this reason, we consider the identified HLA alleles as preliminary and need confirmation in a large study. Although the subjects were enrolled as part of a sub-study of the larger randomized clinical trial and reflected the racial distribution of all subjects enrolled within KPNC, it is likely they did not reflect the racial and ethnic background of all subjects who were enrolled in the larger global clinical trial. For example, if an HLA allele was more common in a racial subgroup that also had more fever, we may have identified that the HLA allele as being associated with fever due to the effects of confounding from population structure because of our sample size did not allow for adjustment for race/ethnicity. In addition, we were limited in our ability to assess several HLA alleles (DRB3, DRB4, and DRB5) because these assays did not perform well.

In conclusion, utilizing MCV-associated fever and febrile seizures as a model we can begin to understand the link between population health (surveillance of a large population for vaccine adverse events) and personalized medicine (the genetic basis for susceptibility to MCV-associated fever and febrile seizures), and the implications that these findings might have on the precision vaccinology field. Vaccine-associated fever, potentially due in part to genetic variation, may serve as a proxy for an immunologically meaningful phenotype and may serve to reassure parents and benefit public health by linking it to improved effectiveness.

Disclosure of potential conflicts of interest

NPK report research support from GlaxoSmithKline for the Phase 3 clinical trial from which the patients in the study were recruited. The measles vaccines which are the focus of this study were provided by GlaxoSmithKline as part of this clinical trial. NPK has received research grants from GlaxoSmithKline, Sanofi Pasteur, Merck, Pfizer, and Protein Science. All other authors have no conflicts of interest relevant to this article to disclose.

All authors attest they meet the ICMJE criteria for authorship.

Acknowledgments

We would like to thank Charlie Chao, MA, CCRP for assistance with sample preparation, shipping, and data management.

Supplementary data

Supplemental data for this article can be accessed on the publisher’s website

Additional information

Funding

This work was supported by the U.S. Department of Health and Human Services, under a cooperative agreement with the National Vaccine Program Office, to NPK. OZ was funded in part by grant 1K01AI139275-01 from the National Institute of Allergy and Infectious Diseases.