Abstract
COVID-19 vaccines are safe and highly effective, but some individuals experience unpleasant reactions to vaccination. As the majority of adults in the US have received a COVID-19 vaccine this year, there is an unprecedented opportunity to study the genetics of reactions to vaccination via surveys of individuals who are already part of genetic research studies. Here, we have queried 17,440 participants in the Helix DNA Discovery Project and Healthy Nevada Project about their reactions to COVID-19 vaccination. Our GWAS identifies an association between severe vaccine side effects and HLA-A*03:01. This association was statistically significant only for those who received the Pfizer-BioNTech vaccine (BNT162b2; p=4.70E-11), but showed a trending association in those who received the Moderna vaccine (mRNA-1273; p=0.005) despite similar sample sizes for study. In Pfizer-BioNTech recipients, HLA-A*03:01 was associated with a two-fold increase in risk of severe vaccine side effects. The effect was consistent across ages, sexes, and whether the person had previously had a COVID-19 infection. The reactions experienced by HLA-A*03:01 carriers were driven by associations with chills, fever, fatigue, and in general feeling unwell.
Introduction
Less than one year after the first publication of a SARS-CoV-2 sequence, COVID-19 vaccines were developed, clinically tested and authorized to be administered in the general population. Within months, hundreds of millions of adults worldwide were vaccinated, and rates of hospitalization among vaccinated individuals dropped precipitously; potential side effects were mild1. During the clinical trials of mRNA vaccines (Pfizer-BioNTech BNT162b2 and Moderna mRNA-1273), both local reactions, such as pain at the injection site, and systemic symptoms such as fatigue, fever, chills, and myalgia were observed in some participants 2,3. Both clinical trials showed that only a small fraction of these reactions could be categorized as severe, and most of the severe reactions followed the second dose. These studies raise the question of what factors explain the interindividual variability of reactions following COVID-19 vaccination.
Younger age and a personal history of SARS-CoV-2 infection prior to vaccination are two factors leading to increased reactogenicity to the vaccine4. However, these factors alone do not explain the large interindividual variability in the degree of severity of the reaction following COVID-19 vaccination. The aim of this study is therefore to identify factors associated with severe reactions following COVID-19 vaccination. Similarly to our hypothesis that genetic factors play a role in the severity of COVID-19 disease5,6, we hypothesized that genetic factors help explain differences in reactions following COVID-19 vaccination. To test our hypothesis, we administered online surveys to 17,440 participants from the Helix DNA Discovery project and the Healthy Nevada Project asking about their vaccination status and their reactions following the vaccine. These surveys included questions about 4 local and 20 systemic symptoms, as well as the overall severity of the reaction and its impact on daily activities in the days following vaccination. These participants were previously sequenced with Helix Exome+ assay, which allowed us to perform rare and common variant genetic associations.
Methods
Cohort and survey
We administered an online survey in June and September of 2021 and received responses from 8,125 Helix DNA Discovery Project participants and 9,315 Healthy Nevada Project participants (Table 1)7,8. These are unselected Helix customers and patients in the Renown Health System who chose to consent to participate in research projects and respond to our survey. The survey takes approximately 10 minutes to complete and can be found in the supplement. The participants in this cohort are aged 18 to 89+, 65% are female, and 85% are of European genetic ancestry (Table 1).
Respondents indicated whether they had been vaccinated and which vaccine type they had received. They rated the severity of their vaccine response as indicated in Table 1. They also answered questions about 24 specific side-effects that can occur after COVID-19 vaccination (see Supplemental survey). There were 3,323 individuals who took the survey in both June and September, updating their COVID-19 vaccination status and infection status. For these individuals, the highest severity score and all reactions were used as a phenotype, regardless of in which survey either were reported.
Genotyping
DNA samples were sequenced and analyzed at Helix using the Exome+® assay as previously described9. Imputation of common variants was performed by pre-phasing samples and then imputing. Pre-phasing was performed using reference databases, which include the 1000 Genomes Phase 3 data. This was followed by genotype imputation for all 1000 Genomes Phase 3 sites with MAF>=1% that have genotype quality (GQ) values less than 20. Imputation results were then filtered for quality (GP>=0.95) so that only high precision imputed variant calls were reported10.
HLA genotypes for A, B, C, DPB1, DQA1, DQB1, and DRB1 were imputed using HIBAG with the default recommendations11. Individual genotypes were imputed using the model for the appropriate genetic ancestry for each individual11. Probabilities higher than 0.5 were used as genotype calls.
Genetic analysis
We used Regenie for the genetic analysis12. Briefly, this method builds a whole genome regression model using common variants to account for the effects of relatedness and population stratification, and it accounts for situations where there is an extreme case-control imbalance, which can lead to test statistic inflation with other analysis methods. The covariates we included were age, sex, age*sex, age*age, sex*age*age, 10 PCs, type of COVID-19 vaccine, whether the individual had COVID-19 prior to getting vaccinated, and bioinformatics pipeline version. We only utilized phenotypes with at least 50 cases, with a minimum minor allele count (MAC) cutoff of 5.
As previously described, a representative set of 184,445 coding and noncoding LD-pruned, high-quality common variants were identified for building PCs and the whole genome regression model9. Each genetic ancestry group was analyzed separately: African, East Asian, European, Hispanic, South Asian.
Results
A minority of individuals have a severe reaction following COVID-19 vaccination
We received survey responses from 8,125 Helix participants and 9,315 HNP participants in June and September 2021. As a control, we first checked whether we observed the same trends as reported in the clinical trials of the Pfizer-BioNTech, Moderna and J&J COVID-19 vaccines. The large majority of participants in these clinical trials did not have a personal history of SARS-CoV-2 infection. To assess the frequency of reactions in our cohort, we therefore restricted our initial analysis to participants who reported no previous infection by SARS-CoV-2. Our results were in agreement with what was previously reported2,3,13. For the Pfizer-BioNTech vaccine, we observed that only 2.7% of individuals (152 of 5,590) had a severe or extreme reaction after the first dose, and 10.2% (569 of 5,583) had a severe or extreme reaction after the second dose (Figure S1). For the Moderna vaccine, we observed that 4.4% of individuals (243 of 5,575) had a severe or extreme reaction after the first dose, and this number grew to 19.5% (1,083 of 5,558) after the second dose (Figure S1). For the J&J vaccine, we observed that 17.7% (80 of 452) had a severe reaction after the single dose (Figure S1). Moreover, we observed a higher number of severe reactions in respondents 18-55 years old compared to those over 55 (Figure S1). These analyses supported the observation that there are large differences in reactions to COVID-19 vaccination between individuals, and validated our self-reported survey as a reliable tool to investigate the genetic basis of these differences.
GWAS identifies HLA region associated with severe reaction to COVID-19 vaccination
To test our hypothesis that genetic variation drives some of the differences seen in vaccine reaction, we performed a genome-wide association study (GWAS) using a phenotype of extreme or severe reaction (1,709 cases) compared to no or mild reaction (6,203 controls). Individuals reporting moderate reaction were excluded from the GWAS (Table 1). Our analysis of 12,602,619 SNPs identified 188 genome-wide significant variants on chromosome 6 (p<5e-8; Figure 1). The lead SNP was rs144943243 / chr6:29820015:AAAAT:A, p=3.51E-11. This variant is in a region of AAAT repeats found upstream of HLA-G, with a MAF of 24% in the European genetic ancestry subset of our cohort.
We next analyzed imputed HLA types against the phenotype to identify whether a specific HLA type was associated with the GWAS signal. We identified a significant association with HLA-A*03:01 (p=5.00E-11), which had a MAF of 15% in the European genetic ancestry subset of our cohort. The distribution of HLA-A*03:01 across the globe is shown in Figure S214. A regression including both this HLA type and rs144943243 identified similar signals for each variant, with rs144943243 retaining a better p-value in a joint analysis. The subsequent analyses assessing the impact of the genetic variant by phenotype lead to similar results for rs144943243 and HLA-A*03:01. In this manuscript, we decided to show the results for HLA-A*03:01, a common allele for the well-studied HLA-A gene.
HLA-A*03:01 association is specific to vaccine type
HLA-A*03:01 had an odds ratio (OR) of 1.5 for individuals to experience a severe / extreme vaccine reaction as opposed to a mild or no reaction. We found that this effect was additive, with individuals with 1 copy of HLA-A*03:01 having a phenotype intermediate to those with 0 or 2 copies (Figure 2). We found the effect of this variant to track similarly across age, sex, and whether the person had a history of SARS-CoV-2 infection prior to vaccination (Figure 2, S3). However, our analyses beyond European (p=5.00E-11, OR=1.6) and Hispanic (p=2.93E-04, OR=2.7) genetic ancestries were underpowered to identify associations (Table S1). We also observed that most of the effect seems to occur at the second dose for two-dose vaccines (Figure S4), and that severe reactions subsided within 2 days for 66% of those who experienced severe reactions (Figure S5).
Importantly, we identified that the association signal came almost entirely from reactions to the Pfizer-BioNTech vaccine (p=4.70E-11, OR=2.07), with very little impact on Moderna vaccine response (p=0.005, OR=1.32) despite similar sample sizes for analysis (Pfizer 552 cases vs. 3,142 controls; Moderna 1,018 cases vs. 2,592 controls). While the main phenotype of severe / extreme reaction was more common in Moderna recipients compared to Pfizer-BioNTech recipients, even restricting to a more stringent case definition did not result in statistical significance for the Moderna subset (p=0.008, OR=1.5; 473 extreme reaction cases vs. 1,425 no reaction controls). The sample sizes for those who received J&J or other vaccinations were too low to clearly assess an association. Detailed counts of reaction severity split by genetic ancestry, HLA-A*03:01 status, and vaccine type are reported in Table S1.
Specific vaccine reaction phenotypes are associated with HLA-A*03:01
To understand the severity scores more deeply, we next analyzed the individual symptoms that participants reported as occurring after receiving the vaccine, split into Pfizer-BioNTech and Moderna subsets. We identified that in Pfizer-BioNTech recipients, HLA-A*03:01 was most strongly associated with increased risk of fever (p=1.76E-14), chills (p=6.15E-13), feeling unwell (p=5.73E-10) and fatigue (p=5.60E-07) after receiving the vaccine (Figure 3). The ORs for these reactions ranged from 1.44-2.09, with 41% of those with 1 HLA-A*03:01 copy who received a Pfizer vaccine having at least two of these four symptoms, 27% of those with 2 copies, and 19% of those with 0 copies. Associations with joint pain, headache, swollen lymph nodes, headache and nausea were less predictive but still statistically significant (p<0.001). In contrast, none of these symptoms were significantly associated with HLA-A*03:01 in Moderna recipients (Figure 3).
Discussion
Here, we identified an HLA type, HLA-A*03:01, with a strong association with reactions to COVID-19 vaccines. We find that, all else being equal, individuals with this HLA type who received the Pfizer-BioNTech vaccine are approximately twice as likely to have a severe or extreme reaction following COVID-19 vaccination. Chills and fever were the two specific side effects that were most enriched in individuals carrying 1 or 2 copies of HLA-A*03:01 compared to individuals with two other HLA-A alleles. This association was present in the Helix DNA Discovery cohort, as well as the Healthy Nevada Project cohort. We find this association to trend across age groups, sex, and whether the person had a personal history of COVID-19 prior to vaccination, all of which are known to be associated with severity of vaccine reaction. We find that the effect of this variant was almost entirely driven by participants who had received the Pfizer-BioNTech vaccine. The signal for the association was very weak in those receiving the Moderna vaccine and did not approach genome-wide significance, despite a sample size roughly equal to that of the Pfizer vaccine recipients. This difference is surprising given the similarities between the two mRNA vaccines, but it is possible that the increased prevalence of vaccine reactions in Moderna recipients may make the signal harder to identify in that subgroup (Figure 2, S1, S4). Our sample size for J&J and other vaccines were too low for adequate analysis power.
Replication of the result in other cohorts who were mostly vaccinated with non-mRNA COVID-19 vaccines will be informative to hypothesize about potential mechanisms leading to these severe reactions. An effect could be linked to a common peptide derived from all vaccines showing an association. The alternative where the effect is only seen robustly for the Pfizer-BioNTech mRNA vaccine could lead to the hypothesis that the reaction is driven by a peptide unique to the Pfizer-BioNTech vaccine. Importantly, as we were preparing this manuscript, 23andMe published on their blog (https://blog.23andme.com/23andme-research/reaction-to-covid-vaccine/, accessed 14 November 2021) showing that HLA-A*03:01 was also the strongest genetic association with COVID-19 vaccine response in their cohort with a p-value p=1.5E-130. No additional details regarding the effect by vaccine type, dose, or symptoms were available yet.
An association of HLA-A*03:01 with reactogenicity to other non-COVID-19 vaccines could also be informative. However, to our knowledge, very few studies have been published reporting genetic associations with reactogenicity after non live-attenuated vaccines. In the past, a few studies have investigated the genetics of response to vaccines by looking at the levels of antibodies after a certain time period in vaccine recipients, often identifying associations with HLA genes15,16. Other studies have elucidated why some patients presented with a rare life-threatening disease following vaccinations with live-attenuated vaccines17. For example, genetic defects in the IL-12 dependent IFN-gamma pathway cause BCGitis after BCG vaccination18. The difficulty of collecting appropriate phenotype information on the reaction to a vaccine with standard medical data or electronic health records could be one reason to explain the small number of these studies. Another reason may be the difficulty to justify sequencing individuals to investigate a transient phenotype. The ability for us to ask new questions and gain new biological insights highlight the importance of continuing to survey and engage participants who are enrolled in ongoing genetic research projects.
This study raises many new questions about the potential role of HLA-A*03:01 in severe reactions following vaccination. Lastly, it is important to emphasize that our study defined a severe reaction to the vaccine as a reaction that would interfere with daily routine shortly after receiving the vaccine. These reactions were symptoms such as chills or fever, which cannot be compared with the severity of COVID-19 disease experienced by many individuals. COVID-19 vaccines have consistently been shown to be safe and very effective to prevent hospitalizations and life-threatening disease following SARS-CoV-2 infections2,3,13.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Acknowledgements
Funding was provided to DRI by the Nevada Governor’s Office of Economic Development. Funding was provided to the Renown Institute for Health Innovation by Renown Health and the Renown Health Foundation. We acknowledge the entire Helix Bioinformatics team for their contributions to the production exome sequencing pipeline. We thank C. Clinton, KT Farley and E. Levin for their contribution running the Helix DNA Discovery Project and the work with the IRB. We thank all of the genomic representatives of the Healthy Nevada Project (HNP). We thank Renown Health and DRI marketing for helping to launch the HNP project. We thank the participants in the Helix DNA Discovery Project, and the Healthy Nevada Project.