Abstract
Background Scarce information exists in relation to the comparison of seroconversion and adverse events following immunization (AEFI) with different SARS-CoV-2 vaccines. Our aim was to correlate the magnitude of the antibody response to vaccination with previous clinical conditions and AEFI.
Methods A multicentric comparative study where SARS-CoV-2 spike 1-2 IgG antibodies IgG titers were measured at baseline, 21-28 days after the first and second dose (when applicable) of the following vaccines: BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, ChAdOx1-S, Ad5-nCoV and Ad26.COV2. Mixed model and Poisson generalized linear models were performed.
Results We recruited 1867 subjects [52 (SD 16.8) years old, 52% men]. All vaccines enhanced anti-S1 and anti-S2 IgG antibodies over time (p<0.01). The highest increase after the first and second dose was observed in mRNA-1273 (p<0.001). There was an effect of previous SARS-CoV-2 infection; and an interaction of age with SARS-CoV-2, Gam-COVID-Vac and ChAdOx1-S (p<0.01). There was a negative correlation of Severe or Systemic AEFI (AEs) of naïve SARS-CoV-2 subjects with age and sex (p<0.001); a positive interaction between the delta of antibodies with Gam-COVID-Vac (p=0.002). Coronavac, Gam-COVID-Vac and ChAdOx1-S had less AEs compared to BNT162b (p<0.01). mRNA-1273 had a higher number of AEFIs. The delta of the antibodies showed an association with AEFIs in previously infected individuals (p<0.001).
Conclusions The magnitude of seroconversion is predicted by age, vaccine type and SARS-CoV-2 exposure. AEs are correlated with age, sex, and vaccine type. The delta of the antibody response is positively correlated with AEs in patients previously exposed to SARS-CoV-2.
Registration number NCT05228912
Introduction
Covid-19 is a pneumonia-like disease caused by a coronavirus, named SARS-CoV-2. It is highly contagious, and the World Health Organization declared it in 2020 as a pandemic (1). SARS-CoV-2 causes COVID-19 that has a wide range of clinical presentations, from asymptomatic disease to severe acute respiratory distress syndrome and death(2). SARS-CoV-2 utilizes its surface spike glycoprotein to enter host cells. Each unit of the spike trimer contains an S1 and S2 subunit, with the N-terminal S1 subunit binding to the cellular angiotensin-converting enzyme 2 (ACE2) receptor through an internal receptor-binding domain (RBD)(3).
To prevent SARS-CoV-2 infection, 113 vaccines are being tested in human clinical trials, and 44 have reached the final step of testing(4). Based on the mechanism of action, vaccines can be clustered in four groups: 1) mRNA vaccines, for example, BNT162b2 mRNA and mRNA-1273. They use genetically engineered modified RNA to produce the spike protein that prompts an immune response safely. 2) Viral vector (adenovirus) vaccines, for example, ChAdOx1-S, Ad26.COV2, Ad5-nCoV, and Gam-COVID-Vac. They use a virus that has been genetically engineered so that it cannot cause disease but produces coronavirus proteins to safely generate an immune response. 3) Inactivated virus vaccines, for example, Coronavac, use a form of the virus that has been inactivated or weakened, so it does not cause disease but still generates an immune response. 4) Protein subunits vaccines, for example, NVX-CoV2373. It is a protein-based vaccine that uses harmless fragments of proteins or protein shells that mimic the SARS-CoV-2 to safely generate an immune response(5).
Currently BNT162b2 mRNA, mRNA-1273, ChAdOx1-S, Ad26.COV2, Ad5-nCoV, Gam-COVID-Vac, Coronavac, have been administered to people across 587 countries(4). However, no studies have compared using the same assay and time frame the effectiveness of seroconversion and incidence of adverse event events in response to vaccines in different countries in the same study.
This study aimed to correlate the magnitude of the antibody response after the first and second dose (if applicable) between BNT162b2 mRNA, mRNA-1273, ChAdOx1-S, Ad26.COV2, Ad5-nCoV, Gam-COVID-Vac, and Coronavac by measuring SARS-CoV-2 spike 1-2 IgG antibodies, using the same standardized assay, with previous clinical conditions and systemic and severe adverse events (AEs).
Materials and Methods
This was a multicentric observational study of volunteers who received an approved complete scheme of BNT162b2 mRNA, mRNA-1273, ChAdOx1-S, Ad26.COV2, Ad5-nCoV, Gam-COVID-Vac, or Coronavac COVID-19 vaccine during 2021 in five hospital centers (Hospital Clinica Nova, Humanitas Clinical and Research Center, Fundacion San Francisco Xavier, Ternium Health Center in Rio, Hospital Municipal San Jose, Hospital Interzonal de Agudos San Felipe) from four different countries: Mexico, Italy, Brazil and Argentina. This study followed STROBE guidelines (6). The study was approved by each local Institutional Review Board and conducted per the Code of Ethics of the World Medical Association (Declaration of Helsinki) for experiments that involve humans.
The inclusion criteria were volunteers of both genders, any age, who consented to participate, planned to conclude the immunization regimen of any vaccine, and agreed to be followed up for the duration of the study. The following vaccines were available: BNT162b2 mRNA, mRNA-1273, ChAdOx1-S, Ad26.COV2, Ad5-nCoV, Gam-COVID-Vac, or Coronavac. The exclusion criteria were to have received any SARS-CoV-2 vaccine prior to study entry.
Each country’s Health System defined the vaccines available, the schedule, and dose assignation. On the vaccination day, the research team invited any subject who planned to receive any vaccine scheme, explained the project and asked to sign the informed consent. Inclusion-exclusion criteria were applied, and a plasma sample was collected. The baseline sample was taken before receiving the first dose of any SARS-CoV-2 vaccine (T0). The second (T1) and third samples (T2) were taken after 21 days (+/- 7 days) of the first and second dose, respectively.
At each visit, participants had to answer a questionnaire. The basal questionnaire aimed at obtaining patients’ medical history and previous SARS-CoV-2 infections. The questionnaires applied after the first and second dose of vaccines aimed at recognizing adverse events following immunization (AEFI)(7) and identifying a potential SARS-CoV-2 infection after receiving any vaccine dose. SARS-CoV-2 infection was also monitored by the epidemiology team through PCR testing, which informed the research team of any new infection.
Primary Outcome, IgG determination
Our primary outcome was to correlate the magnitude of the antibody response to vaccination with previous clinical conditions and AEs. To determine the amount of specific anti-S1 and anti-S2 IgG antibodies against SARS-CoV-2 in plasma samples, the laboratory personnel used a chemiluminescence immunoassay (CLIA) developed by DiaSorin, which had a sensitivity of 97.4% (95% CI, 86.8-99.5) and a specificity of 98.5% (95% CI, 97.5-99.2). The results were reported as follows: <12.0 AU/ml was considered negative, 12.0 to 15.0 AU/ml was indeterminate, and > 15 AU/ml was positive. This kit is comparable with other commercial kits and has been used elsewhere (8–11).
The variables we analyzed were age, sex, personal medical history (for example, the presence of type 2 diabetes, hypertension, obstructive pulmonary disease, any heart condition, obesity, cancer, liver steatosis, any autoimmune disease), and confirmed SARS-CoV-2 (through nasal swab or serologic tests). The following AEs were considered of particular interest: fever (>37.5°C), adenopathy, diffuse rash, edema, facial paralysis, orthostatic hypotension, headache, arthralgia, myalgia, nausea, vomit, and diarrhea. Anti-S1 and anti-S2 IgG antibodies against SARS-CoV-2 were measured at baseline, 21-28 days post-first dose (S1), and 21-28 days post-second dose if applicable.
Statistical Methods
The researchers reviewed the quality control and the anonymization of the database. Normality assumption was evaluated with the Shapiro-Wilk test or Kolmogorov. Descriptive statistics such as mean, standard deviation, median, interquartile range, frequencies, and percentages were computed. Kruskal-Wallis test with Dunn’s multiple comparisons test was performed for IgG comparisons. We performed a mixed model where the dependent variable was the delta of the antibodies. The personal variability and time were constructed as a random effect, whereas each type of vaccine, age, history of SARS-CoV-2, and the interaction of each vaccine with age were fixed effects. BNT162b mRNA was the reference vaccine group. For the analysis of AEs, the study population was divided into two separate cohorts by stratifying according to the history of SARS-CoV-2 infection (1050 Naïve SARS-CoV-2 subjects and 471 SARS-CoV-2 history subjects). Poisson generalized linear models (PGLM) were performed on both cohorts with the counts of AEs events after the second dose as the outcome variables. Included regressors were age, sex, body mass index, 1,000 AU/mL variation of IgG levels after the second dose compared to the baseline IgG levels, type of vaccine and the interactions between vaccines and IgG level variation. BNT162b mRNA was the reference vaccine group Missing completely at random values were analyzed through complete case analysis since missing antibody levels were less than 5% and we considered it incorrect to impute any AEFI. A sample size of 1870 patients was calculated, according to the primary aim, by using a mixed model formula with an alpha of 0.05, power of 90%, the effect size of 0.15, and k=7. The statistical programs used were R v.4.0.3 and Python v. 3.8.3. The analyses were two-tailed. A p-value less than 0.05 was statistically significant.
Results
A total of 1867 patients were recruited from all countries: 1352 from Mexico, 42 from Italy, 260 from Brazil, and 213 from Argentina. The most frequent vaccine used was ChAdOx1-S in 666 subjects, Coronavac in 582, BNT162b2 mRNA in 289, Gam-COVID-Vac in 213, mRNA-1273 in 65, Ad26.COV2 in 31, and Ad5-nCoV in 19. The mean (SD) age was 52 (16.8) years, being statistically different across vaccine groups (p<0.05) as some vaccines were proposed to a particular age group. Fifty-two % of subjects were men, 559 (30%) had obesity, and 501 (26.8%) had hypertension. Table 1 shows the medical history of patients divided by vaccine type.
General characteristics and medical history
Table 1 shows the medical history and comparison between vaccine groups, p<0.05 was considered statistically significant.
History of SARS-CoV-2
Positive history of SARS-CoV-2 infection was considered if the volunteer had a confirmatory swab or serological test. Additionally, it was considered if specific anti-S1 and anti-S2 IgG antibodies against SARS-CoV-2 in plasma samples were > 15 AU/ml at baseline. We had a total of 627 positive cases before vaccination, of which 165 (24.8%) received ChAdOx1-S, 246 (42.2%) Coronavac, 107 (37.1%) BNT162b2 mRNA, 55 (25.8%) Gam-COVID-Vac, 31 (47.6%) mRNA-1273, 18 (58%) Ad26.COV2, and 5 (26%) Ad5-nCoV.
Antibody titers
All vaccines showed significant anti-S1 and anti-S2 IgG antibody changes with significant differences across vaccines and depending on SARS-CoV-2 history. In naïve patients, the highest increase [Median (IQR) AU/ml] after the first dose was observed in mRNA-1273 [175.5 (76.7)], then BNT162b mRNA [78.1(55)], and Ad5-nCoV [43.4(55.9)]. In subjects previously exposed to SARS-CoV-2, the highest increase after the first dose was observed in mRNA-1273 [5920(4705)], then BNT162b mRNA [2500(3080)], and ChAdOx1-S [1025(2489)].
In naïve patients, the highest increase [Median (IQR) AU/ml] after the second dose (T2) was observed in mRNA-1273 [1875 (1190)], then BNT162b mRNA [998 (1241)], and Gam-COVID-Vac [501 (55.9)], whereas in SARS-CoV-2 previously exposed subjects the highest increase after the second dose was observed in mRNA-1273 [4950(3060)], then Gam-COVID-Vac [3620 (5356)], and BNT162b mRNA [2630 (2672.5)]. Figure 1 and Figure S1 show the change in antibody levels through time, classified by vaccine type. Table 2 shows the antibody levels classified by vaccine type. The mixed model that showed a positive effect of SARS-CoV-2 previous infection (B=767.468, p<0.001), a positive interaction of age with previous SARS-CoV-2 infection (B=13.03, p=0.003), and an interaction of Gam-COVID-Vac and ChAdOx1-S with age (B=46.5, p<0.001; B=20.5; p=0.006). With respect to BNT162b mRNA, mRNA-1273 showed a higher change in antibody titers, while all the other vaccines had less change (p<0.01). We computed sex in the model, but there was no significant correlation, so we eliminated it from the final model. These results are shown in Table 3 and Figure 2.
IgG antibody response was measured in serum of naïve and SARS-CoV-2 previously exposed (SARS-CoV-2 Exp) subjects at different time points (T0, T1 and T2) and vaccinated with different vaccine types. Samples ≥ 15 AU/mL were considered positive. Log scale on y axis. The box plots show the interquartile range, the horizontal lines show the median values, and the whiskers indicate the minimum-to-maximum range. Each dot corresponds to an individual subject. P-values were determined using 2-tailed Kruskal-Wallis test with Dunn’s multiple comparisons test. P values refer to baseline (T0) when there are no connecting lines.
Antibody differences between vaccine types
Mixed model of antibody changes
Antibody titers classified by age type of vaccine and Time
Figure shows change in antibody concentration through age, divided SARS-CoV-2 history, and vaccine type. Antibodies are expressed in logarithm and age has a Z distribution.
Adverse Events Following Immunization (AEFI)
At least one AEFI was observed after the first dose in 71% of respondents receiving BNT162b2 mRNA, in 93% with mRNA-1273, 38% with Gam-COVID-Vac, 42% with Coronavac, 51% with ChAdOx1-S, 74% with Ad5-nCoV, and 81% with Ad26.COV2. After the second dose, 65%, 88%, 23%, 33%, and 23% of the respondents experienced at least one AEFI after receiving BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, and ChAdOx1-S, respectively. For each vaccine, the majority of the adverse events occurred during the first 24 hrs after injection, either after the first or the second dose. Patients receiving BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, ChAdOx1-S, Ad5-nCoV, and Ad26.COV2 subjectively qualified the AEFI after the first dose as “very mild” or “mild” in 85%, 80%, 95%, 84%, 67%, 93%, and 57% of cases, respectively. The AEFI after the second dose were qualified as “very mild” or “mild” by 82%, 49%, 98%, 89%, and 76% of patients receiving BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, and ChAdOx1-S respectively. Finally, 49% of patients receiving mRNA-1273 qualified for their adverse events after the second dose as “moderate”. Figure 3 shows the percentages of AEs per vaccine. AEs were evaluated by using PGLM. The model for naïve SARS-CoV-2 subjects showed a negative association of the count of adverse events with age (B=-0.0327, p<0.001) and sex (B=-0.8145, p<0.001). A negative effect with respect to BNT162b mRNA was observed for Gam-COVID-Vac (B=-1.8582, p<0.001), Coronavac (B=-0.5736, p=0.047) and ChAdOx1-S (B=-0.7871, p=0.006). Finally, a positive interaction was found between the delta of antibodies due to vaccination and Gam-COVID-Vac (B=0.5604, p=0.002).
Percentage of systemic or severe adverse events after second dose, stratified by vaccine Graph shows the percentage of systemic or severe adverse events grouped by frequency of appearance
The model for subjects previously exposed to SARS-CoV-2 showed a negative association of the count of adverse events with age (B=-0.0190, p=0.013) and sex (B=-1.0042, p<0.001). The delta of the antibodies showed a positive association with the count of adverse events (B=0.1018, p<0.001). A positive effect with respect to BNT162b mRNA was observed for Gam-COVID-Vac (B=1.2737, p=0.038), whereas a negative effect was observed for Coronavac (B=-1.1304, p<0.001). Negative interactions were also found between the delta of antibodies due to vaccination and the vaccines Gam-COVID-Vac (B=-8.5758, p<0.001) and ChAdOx1-S (B=-1.2386, p=0.013). See Tables 4 and 5 and Figure 4 for details.
Number of systemic or severe adverse events following vaccination in SARS-CoV-2 Naïve subjects.
Number of systemic or severe adverse events following vaccination in SARS-CoV-2 previously exposed subjects.
Women and young individuals developed more AEs. There is no effect on Body Mass Index. In the naïve SARS-CoV-2 cohort, the vaccines Gam-COVID-Vac, Coronavac and ChAdOx1-S are related to fewer AEFIs after the second dose than BNT162b mRNA. Also, for naïve SARS-CoV-2 patients who received Gam-COVID-Vac, higher antibody levels after the second dose were related to a greater number of AEs. In people previously exposed to SARS-CoV-2, those receiving Gam-COVID-Vac showed a greater number of AEs compared to BNT162b mRNA, while subjects receiving Coronavac showed significantly fewer events when compared to BNT162b mRNA.
SARS-CoV-2 infection
Subjects were followed up to 28 days after the second dose administration. Between the first and second dose, one (0.34%) patient with BNT162b2 mRNA, 0 with mRNA-1273, 1 (0.46%) with Gam-COVID-Vac, 3 (0.51%) with Coronavac, 4 (0.6%) with ChAdOx1-S, 2 (10%) with Ad5-nCoV and 1 (3%) with Ad26.COV2 got the SARS-CoV-2 infection. After the second dose, 4 (1.3%) with BNT162b2 mRNA, 1 (1.5%) with mRNA-1273, 0 with Gam-COVID-Vac, 18 (3%) with Coronavac, and 11 (1.6%) with ChAdOx1-S became infected with SARS-CoV-2 (p<0.01).
Discussion
This multicentric study compared the 21-28 day seroconversion, the AEFI after first and second doses, and the associated predictors of AEs of seven of the most common vaccines used worldwide: BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, ChAdOx1-S, Ad5-nCoV, and Ad26.COV2(4).
All vaccines showed immunogenicity and seroconversion after the complete vaccination scheme. The greatest change in antibody levels occurred with mRNA-1273, followed by BNT162b2 mRNA and Gam-COVID-Vac. Previous studies have shown a greater change in antibody titers for mRNA vector vaccines(12–14).
Our study found that the previous history of SARS-CoV-2 was related to an increase in change in antibody levels, as confirmed elsewhere (12). Also, we found higher levels in antibodies in older patients when they had previous SARS-CoV-2 infection and received Gam-COVID Vac or ChAdOx1-S. The previous study from Brazil did not find any association of antibody levels with age. We also found that except for the Gam-COVID Vac, in the other vaccines requiring two doses, the first dose was sufficient to induce maximal antibody levels in subjects with a history of SARS-CoV-2, as previously shown for the BNT162b2 mRNA vaccine(9,15–19)
The study of Logunov et al. that tested the efficacy of Gam-COVID Vac reported a higher increase of anti-SARS-CoV-2 IgG levels in the group of 18-30 years and no particular change in other age groups(20); however they excluded patients with previous SARS-CoV-2 infection, so these results are not comparable.
In the study conducted by Ramasamy et al., similar antibody titers were seen 28 days after the boost vaccination of ChAdOx1 across all groups, regardless of age or vaccine dose 21. It is possible that older patients in our group of the study had a more severe disease during SARS-CoV-2 infection and that they developed higher antibody titers, as it has been reported that patients aged higher than 60 can develop this condition (21).
For most patients, at least one AEFI occurred with all types of vaccines during the first 24 hrs. after injection, and most of them qualified the symptoms as “very mild” or “mild” both after the first and second dose, except for mRNA-1273 for which nearly half of the volunteers reported moderate symptoms after the second dose as reported previously(18,20,22–27).
Our models showed a higher risk for women and young people to experience AEs. A real-world study that reported AEFI of patients receiving BNT162b2 mRNA, mRNA-1273, or Ad26COV2 in one or two doses (28), an independent studies conducted with Gam-COVID also confirmed this results (26,29). The study of Ramasamy et al. on ChAdOx1 showed more adverse events in the group of 18-55 years old(22). The study conducted by Scott A et al. showed that severe adverse events (grade 5) were more frequent in women with Ad5-nCoV (25).
On one hand, for the Naïve SARS-CoV-2 cohort, the vaccines Gam-COVID-Vac, Coronavac, and ChAdOx1-S were related to fewer AEs after the second dose than BNT162b mRNA, and to the best of our knowledge, there are no studies that compare AEs across these types of vaccines. Also, the delta of antibodies was related to AEs after the second dose in Gam-COVID-Vac. As previously published, mRNA-1273 had higher AEFIs when compared to BNT162b2 mRN(30).
On the other hand, for subjects previously exposed to SARS-CoV-2, those receiving Gam-COVID-Vac showed a greater number of AEs than BNT162b mRNA, while subjects receiving Coronavac showed significantly fewer events when compared to BNT162b mRNA. In general, an increase in the antibody levels was related to an increase in the number of AEFIs, but for people receiving Gam-COVID-Vac or ChAdOx1-S, greater changes in antibody levels corresponded to fewer AEs. This last result should be taken with caution due to the limited number of patients falling in the category. To the best of our knowledge, this has not been previously shown. The study conducted by Sung Hee-Lim et al. that compared the antibody response with adverse events in patients vaccinated with BNT162b mRNA or ChAdOx1-S reported no association with antibody change and adverse events; however, patients previously infected with SARS-CoV-2 were excluded (31).
This study shows a low rate of infection in all vaccine groups. However, the highest rate of SARS-CoV-2 infection in this follow-up period was for the one-dose vaccines Ad5-nCoV and Ad26.COV2. Also, Coronavac had a high rate of infection. It is important to consider that the patients were followed for a short period.
As a limitation of our study, we recognize that we had a small sample size in some vaccine groups such as Ad5-nCoV; however, to the best of our knowledge, there is no other study where this type of vaccine has been compared to other types of vaccines in a real-world setting. It will be of interest to follow the IgG titers for a longer period and to evaluate the effect of heterologous combinations of these vaccines.
This comparative study of vaccine types shows positive immunogenicity and seroconversion of BNT162b2 mRNA, mRNA-1273, Gam-COVID-Vac, Coronavac, ChAdOx1-S, Ad5-nCoV, and Ad26.COV2. The highest IgG response was for the mRNA vector vaccines and the lower for the inactivated vaccine. Women and young individuals developed more AEs. In the naïve SARS-CoV-2 cohort, the vaccines Gam-COVID-Vac, Coronavac, and ChAdOx1-S were related to fewer AEs after the second dose of BNT162b mRNA. For patients who received Gam-COVID-Vac, higher antibody levels after the second dose were related to a greater number of AEs. In people previously exposed to SARS-CoV-2, those receiving Gam-COVID-Vac showed a greater number of AEs than BNT162b mRNA, while subjects receiving Coronavac showed significantly fewer events when compared to BNT162b mRNA. In general, an increase in antibody levels was related to AEs.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Data Availability Statement
Data are available upon reasonable request to the authors.
Ethics statement
Ethics committee/local Institutional Review Board from Hospital Clinica Nova, Humanitas Clinical and Research Center, Fundacion San Francisco Xavier, Ternium Health Center in Rio, Hospital Municipal San Jose, Hospital Interzonal de Agudos San Felipe gave ethical approval.
Author Approval
all authors all authors read an approved the final version of the manuscript.
Author Contributions
Conceptualization: MERI, MR, MT, MASS, and EA. Formal analysis: MERI, MR, MT, MASS, AGC, CP, RL, S, MM, and RB. Investigation: MERI, AGC, MR, MT, MASS, CP, RL, RS, RB, MM, DRS, YGHR, AGAV, GFDRP, IABF, RGF, EZ, ARAD, MHG, ABB, CA, GPB, MT, and MR. Resources: MR, MT, and MASS Writing – original draft: MERI, AGC, CP, MM, RS, RL, RB, MASS, EA, DRS, YGHR, AGAV, GFDRP, IABF, RGF, EZ, ARAD, MHG, ABB, CA, GPB, MT, MR. Writing – review and editing: MERI, AGC, CP, MM, RS, RL, RB, MASS, EA, DRS, YGHR, AGAV, GFDRP, IABF, RGF, EZ, ARAD, MHG, ABB, CA, GPB, MT, MR. Project administration: MR, MT, MASS. Supervision: MERI, MT, MR, and MASS. Funding acquisition: MT, MR, and MASS. All authors contributed to the article and approved the submitted version.
Funding
This research was conducted using private funding from Techint Group of Companies.
Conflict of Interest
The authors have declared no competing interest.
Acknowledgments
This project has been supported by the Techint Group of Companies.
We thank Erika Bienek (Techint Group Community Relations Director), Walter Bruno (Humanitas Communications Director), Cesar Bueno (Fundación São Francisco Xavier President), Mario Galli (Techint Group Communications Director), Pablo Sturiale (Apsot Director) for their support of the work and discussions.
Footnotes
This version of the manuscript has an update in affiliation of author Riccardo Sarti that should be: 3. IRCCS Humanitas Clinical and Research Center, Rozzano, MI, Italy. 4. Humanitas University, Department of Biomedical Sciences, Pieve Emanuele, MI, Italy.
Abbreviations
- (AEFI)
- Adverse events following immunization
- (AEs)
- Severe adverse events