Abstract
Background Few head-to-head evaluations of immune responses to difference vaccines have been reported.
Methods Surrogate virus neutralization test (sVNT) antibody levels of adults receiving either 2 doses of BNT162b2 (n=366) or CoronaVac (n=360) vaccines in Hong Kong were determined. An age-matched subgroup (BNT162b2 (n=49) vs CoronaVac (n=49)) were tested for plaque reduction neutralizing (PRNT) and spike binding antibody and T cell reactivity in peripheral blood mononuclear cells (PBMC).
Findings One month after the second dose of vaccine, BNT162b2 elicited significantly higher PRNT50, PRNT90, sVNT, spike receptor binding, spike N terminal domain binding, spike S2 domain binding, spike FcR binding and antibody avidity levels than CoronaVac. The geometric mean PRNT50 titres in those vaccinated with BNT162b2 and CoronaVac vaccines were 251.6 and 69.45 while PRNT90 titres were 98.91 and 16.57, respectively. All of those vaccinated with BNT162b2 and 45 (91.8%) of 49 vaccinated with CoronaVac achieved the 50% protection threshold for PRNT90. Allowing for an expected seven-fold waning of antibody titres over six months for those receiving CoronaVac, only 16.3% would meet the 50% protection threshold versus 79.6% of BNT162b2 vaccinees. Age was negatively correlated with PRNT90 antibody titres. Both vaccines induced SARS-CoV-2 specific CD4+ and CD8+ T cell responses at 1-month post-vaccination but CoronaVac elicited significantly higher structural protein-specific CD4+ and CD8+ T cell responses.
Conclusion Vaccination with BNT162b2 induces stronger humoral responses than CoronaVac. CoronaVac induce higher CD4+ and CD8+ T cell responses to the structural protein than BNT162b2.
Summary At a Glance Through the head-to-head comparison, vaccination with BNT162b2 induces significantly higher levels of SARS-CoV-2 specific binding and neutralizing antibody responses when compared to CoronaVac. CoronaVac induce higher CD4+ and CD8+ T cell responses to the structural protein than BNT162b2.
INTRODUCTION
SARS-CoV-2, the cause of COVID-19, emerged in late 2019 leading to a devastating pandemic (1,2). Over 240 million COVID-19 infections including 4.9 million deaths have been reported to the World Health Organization as of October 25th, 2021 (3). Many COVID-19 vaccines were rapidly developed, evaluated and deployed with over 4 billion doses of COVID-19 vaccines administered worldwide so far (3). These include inactivated whole virus, lipid nanoparticle (LNP)-encapsulated mRNA, adenovirus-vectored and protein sub-unit vaccines. Virus neutralizing antibodies correlate with protection from re-infection following natural infection as well as vaccination (4, 5) and in experimental animal models (6). T cell immune responses are consistently elicited after natural infection and correlate with reduced disease severity in humans and reduced viral loads in non-human primates (7-11).
The safety, immunogenicity and efficacy of these vaccines have been evaluated in separate clinical trials (12-14) but there are few “head-to-head” comparisons of different vaccines. Here, we compare the humoral and cellular responses from vaccinees who received either BNT162b2 or CoronaVac.
MATERIALS AND METHODS
Cohort study design and participants
Healthy adults aged between 18-79 years were recruited in Hong Kong SAR, China at vaccination centers at the Chinese University of Hong Kong Medical Centre, Prince of Wales Hospital and Kowloon Bay Vaccination Station between March10 and August 31, 2021. Those with previous COVID-19 infection were excluded. The participants received two doses of either BNT162b2 [Comirnaty(BioNTech)] vaccine (21 day interval between the two doses) or CoronaVac (Sinovac) vaccine (28 day interval between the two doses), according to the manufacturer’s recommendation. Demographic information was collected from each participant prior to vaccination (Supplementary Table 1). Ten ml of heparinized blood were collected from each donor before vaccination and at one month after receiving the second vaccine dose.
Immunological Assays
Humoral and cellular immune responses were examined to compare the levels of SARS-CoV-2 neutralizing, binding and FcR+ binding antibody in plasma and T cell reactivity in peripheral blood mononuclear cells (PBMC) between the two vaccine groups. Details of specimen storage and processing and the serological methods used in ELISA for spike receptor binding domain (RBD), spike N terminal domain (NTD) binding, S2 and N protein antibodies, ELISA for FcγRIIIa-binding and antibody avidity, surrogate virus neutralization test (sVNT), plaque reduction neutralization test (PRNT), and SARS-CoV-2-specific T cells by Intracellular Cytokine Staining (ICS) are provided in the supplemental file.
Statistical analysis
Methods used for statistical analysis are provided in the supplemental file.
Role of the funding source
The sponsor of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.
RESULTS
Adult volunteers (age range 18-79 years) were recruited between March 10 and August 31, 2021; 366 of whom received two doses of BNT162b2 and 360 received two doses of CoronaVac. The mean (SD) age of the two groups were 45.01 (13.16) and 51.78 (9.92) years respectively (p<0.0001). Demographic information is summarized in Supplementary Table 1. We used sVNT, which has high overall concordance to the PRNT (15,16), to examine the level of SARS-CoV-2 specific neutralizing antibody from the plasma samples collected before and at 1 month after the 2nd dose of vaccination. The plasma samples from all vaccinees before vaccination were negative in sVNT (data not shown). The mean % inhibition in the sVNT test in post-vaccine plasma for BNT1626 and CoronaVac was 93.63% (standard deviation (SD) 7.92) and 52.11% (SD 22.06), respectively (p<0.0001) (Figure 1A).
The percentage of inhibition were detected by surrogate virus neutralization test (sVNT) from the plasma collected from adult individuals who received two doses of BNT162b2 (n=366) or received two doses of CoronaVac (n=360) (A). Various antibody responses were further determined from the plasma from an age-matched subgroup of (A) (49 vs 49). (B) The percentage of inhibition from the plasma of pre-vaccination and 1 month after two doses of vaccination were tested by sVNT. The dash line at 30% indicates the negative threshold of the sVNT. Comparison of the (C) PRNT50 and (D) PRNT90 from the plasma collected at 1 month after two doses of vaccination between the BNT162b2 and CoronaVac groups. The levels of (E) RBD-specific (F) NTD-specific and (G) S2-specific IgG antibodies from the plasma of pre-vaccination and 1 month after two doses of vaccination were tested by ELISAs. **** indicates p<0.0001
From the data above, we estimated that a sample size of 49 for each vaccine group provided statistical power higher than 80% with 95% confidence level to discern differences between the two vaccines. To avoid age-related bias, we randomly 49 selected age-matched samples from BNT162b2 or CoronaVac vaccinees for more detailed analysis. The demographic, underlying co-morbidities and other potential risk-factors were shown in Table 1. The mean % inhibition in the sVNT test in post-vaccine plasma for BNT1626 and CoronaVac was 94.8% (standard deviation (SD) 3.45) and 63.9% (SD 16.72), respectively (p=7.71×10−17) (Figure 1B).
Comparison of characteristics between the two vaccine groups (n=98).
PRNT is the gold-standard method to evaluate virus neutralization. All vaccinees developed detectable PRNT50 antibody titers. BNT162b2 vaccinees had significantly higher geometric mean titers (GMT) of 251.60 (±1SD range from 147 to 432) compared to 69.45 (±1SD range from 28 to 172) for the CoronaVac group (p=1.24×10−9) (Figure 1C). The corresponding PRNT90 antibody responses for BNT162b2 (GMT 98.91, ±1SD range from 44 to 221) was significantly higher than for CoronaVac vaccines (GMT 16.57, ±1SD range from 5 to 54) (p=5.82×10−10) (Figure 1D). It has been suggested that the neutralizing antibody titer associated with protection from re-infection in 50% of individuals is 20% of the mean convalescent antibody levels (5). Using RT-PCR confirmed convalescent symptomatic patients sera (16), using the same PRNT methods, we estimated that this 20% convalescent antibody titer threshold for 50% protection from re-infection for PRNT90 was 1:8.75 (95% CI 1:6.6-1:11.6) (Supplementary Figure 1). Thus, we estimated that all 49 vaccinated with BNT162b2 and 45 (91.8%) of 49 of those vaccinated CoronaVac achieved the 50% protection threshold, one month after the second dose of vaccination. Neutralizing antibody levels reportedly fell by 7.3-fold within 6 months of CoronaVac immunization (17). Thus we estimate that only 8 (16.3%) of 49 CoronVac vaccines would retain protective levels of neutralizing antibody 6 months post vaccination while 79.6% of BNT162b2 vaccinees would do so, assuming comparable waning of antibody in BNT162b2 vaccinees.
Adjusting for age, gender and history of other vaccine uptake, subjects who received CoronaVac had a significantly lower PRNT90 value compared to those who received BNT162b2 (Supplementary Table 2). Stratifying regression analysis by choice of vaccine, older age was independently associated with lower PRNT90 values with either vaccine (Supplementary Table 3).
Antibodies against different domains of the spike can facilitate protection (18, 19, 20). Thus, we further tested the levels of IgG antibodies specifically binding to the RBD, NTD and S2 domains of the spike protein, respectively. While both vaccines elicited significant increases in these antibodies as detected by ELISA, BNT162b2 vaccine elicited significantly higher antibody levels to all three antigens (Figure 1E-G). As expected from a whole virus vaccine which also contains the virus N protein, 59.2% (29/49) and 40.8% (20/49) of post-vaccination plasma from CoronaVac group were positive in ELISA against the full length (aa 44-419) and C-terminal domain of the N protein (Supplementary Figure 2). One plasma sample from a BNT162b2 vaccinee had high ELISA binding antibody to N protein (but not C-terminal domain of N) in post-vaccination plasma, but was negative for other viral proteins (ORF8, data not shown), and therefore may represent cross-reactivity in the N antibody assay, perhaps with other seasonal coronaviruses.
Antibodies which bind to Fc Receptors on effector cells such as natural killer (NK) cells may stimulate antibody dependent cell cytotoxicity (ADCC) upon FcgRIIIa binding of virus-antibody complexes leading to down-stream signaling. Antigen specific FcgRIIIa binding significantly corelates with NK cell function of degranulation, killing and cytokine production (21). Furthermore, increased ADCC function is associated with mRNA (22) and protein (23) vaccination and reduced morbidity during infection (24, 25), and is therefore regarded as one of the indicators of immune protection. We found that recipients from both BNT162b2 or CoronaVac elicit increases of FcγRIIIa binding spike antibodies after two doses of vaccination (Figure 2A). However, the antibody level from the BNT162b2 group was significantly higher than the CoronaVac group. CoronaVac recipients had higher levels of FcγRIIIa binding antibodies to N than the BNT162b2 group after vaccination (p <0.0001) as expected, but CoronaVac N-FcγRIIIa responses were still significantly lower than COVID-19 cases (Figure 2B). FcγRIIIa binding is more likely to result in downstream signaling leading to effector function if antibodies have high avidity (26). S-specific FcγRIIIa as well as S-IgG avidity was significantly lower in CoronaVac recipients than BNT162b2 recipients or COVID-19 convalescent patients (Figure 2C, D). The average avidity of N specific FcγRIIIa from CoronaVac recipients were comparable to the patients with COVID-19 (Figure 2E).
The levels of FcγRIIIa binding antibodies and their avidity were detected from the plasma collected from adult individuals who received two doses of BNT162b2 (n=49) or received two doses of CoronaVac (n=49). Recovered COVID-19 cases (n=34, timepoint 56 days post infection +/- 17 days (mean+/- SD)), and healthy adults negative for SARS-CoV-2 (n=40) were served as positive and negative control respectively. The levels of (A) FcγRIIIa binding S antibodies and (B) FcγRIIIa binding N antibodies were tested from the plasma collected before and at 1 month after 2 doses of vaccination. The avidity indexes of (A) S FcγRIIIa, (B) N FcγRIIIa and (C) S IgG were determined as the proportion of antibodies remaining after 3x washes with 8M urea, compared to the total FcγRIIIa-binding antibodies to each protein. **** indicates p<0.0001
We collected PBMCs from the same subgroup where possible, making a cohort of 25 vaccinees who received two doses of BNT162b2 and 30 with two doses of CoronaVac. Overlapping peptides of the structural proteins including S, N, Envelope and Matrix (SNEM) or S only pools were used to induce specific T cell responses in the PBMCs. The structural and S specific CD4+ and CD8+ T cell responses were quantified by flow cytometry (Supplementary Figure 3) from paired samples at pre and post vaccination of CoronaVac and BNT162b2 while RT-PCR confirmed COVID-19 patients were served as positive controls. The average magnitude of the post vaccination CD4+ and CD8+ T cell responses after stimulation of SNEM peptides was significantly higher in CoronaVac group compared those received BNT162b2 (Figure 3A and B). The magnitude of the post vaccination CD4+ T cell response after stimulation of S only peptides was also significantly higher in CoronaVac group but their CD8+ T cell responses was comparable to the BNT162b2 group (Supplementary Figure 4 A and B). Overall, the proportion of subjects that had detectable post vaccination T cell responses (% of IFNγ+ cell is higher than 0.001), termed ‘responders’ was higher in CoronaVac recipients than BNT162b2 recipients, for either SNEM peptides or S peptides alone (Table 2).
PBMCs from pre- (day 0) and post-vaccination (day 30 post after second dose) of BNT162b2 mRNA (pre-vaccine n=25, post-vaccine n=25) and CoronaVac (pre-vaccine n=30, post-vaccine n=30) and recovered COVID-19 cases (n=10, timepoint 59 days post infection +/- 20 days (mean+/- SD)) were stimulated with pooled structural (SNEM) peptides or a DMSO control. The percentage of (A) IFNγ+CD4+ and (B) IFNγ+CD8+ T cell were measured by flow cytometry. Dotted lines represent limit of detection following DMSO background subtraction (IFNγ of CD4+=0.001, IFNγ of CD8+=0.001. (C) The proportion of IFNγ producing IL2 and TNFα CD4+ and CD8+ T cells post vaccination. (D) The phenotype (by CCR7 and CD45RA) of IFNγ responses for T effector memory (TEM), central memory (TCM), terminal effector memory (TeEM) or naïve (TN) CD4+ and CD8+ T cells post vaccination. Bars represent the mean values, and error bars represent SD. Statistical significance was determined by paired t-test between pre-versus post-vaccination timepoint samples, and Kruskal-Wallis test for multiple comparisons between vaccines, and COVID-19 patients. *p<0.05, **<0.01, ***<0.001.
Post vaccination polyfunctional cytokine production, including TNFα and IL-2 by T cells in the vaccine responders was equivalent between the two vaccine types and convalescent COVID-19 samples after stimulation of either the SNEM or S only peptides (Figure 3C; Supplementary Figure 4C). The memory phenotype of IFNγ+ CD4+ T cells in both vaccine groups showed significantly higher percentage of T effector memory (TEM) (p<0.05) and lower in T naïve (TN) (p<0.01) than the COVID-19 recovered patients (Figure 3D). Interestingly, S-specific CD4+ T cell responses memory phenotypes further showed TEM>TCM (T central memory) in the CoronaVac group compared to BNT126b2 which may impact recall at infection and long-term cellular memory (Supplementary Figure 4D).
DISCUSSION
Using a cohort of RT-PCR confirmed convalescent sera from COVID-19 patients, using the assays we have used in the present study, we estimated the 20% convalescent antibody titer threshold for 50% protection from re-infection for PRNT90 was 1:8.75 (95% CI 1:6.6-1:11.6) (5,16). We conclude that all those vaccinated with BNT162b2 and 91.8% of those vaccinated CoronaVac achieved the 50% protection threshold one month after the second dose of vaccination. One month post-second dose likely represents the peak antibody response, beyond which antibody titers are likely to wane. It was reported that neutralizing antibody titres wane 7.3-fold within 6 months after CoronaVac vaccination (17) but comparable data is not available from BNT126b2. If we adjust for a 7.3-fold waning of antibody titres for both vaccines, we estimate that only 8 (16.3%) of 49 receiving CoronaVac vaccines meet the protective threshold while 39 (79.6%) of 49 of those receiving BNT162b2 do so, 6 months post-vaccination. This difference in immunogenicity may explain reported difference in vaccine efficacy between the two vaccines (27-30). Many Phase 3 studies assessed protection within 1 to 3 months after the second dose of vaccine and may not reflect the impact of waning immunity. Furthermore, some of the virus variants (e.g. B.1.351 Beta variant; P1 Gamma variant) circulating in parts of the world lead to an 8 to 10-fold reduction in neutralizing titers and this is likely to further compromise protection afforded by CoronaVac vaccines with its weaker immunogenicity. Our results suggest that booster doses may be required for older CoronaVac recipients and such booster doses. A third dose of CoronaVac vaccine appears to boost antibody levels (17) and our data suggest that these may be needed for older CoronaVac recipients.
The average magnitude of post vaccination responses was higher in CoronaVac subjects for structural and S-specific T cell responses. Since it is likely that T cell responses are important in limiting severity and fatal outcomes (7-11), both vaccines may be effective in preventing such adverse outcomes of COVID-19. This is consistent with high levels of protection against hospitalization and death in CoronaVac vaccines observed in Chile and Turkey despite lower antibody neutralization titers than BNT162b2 (30, 31). In animal models of re-infection, spike-specific CD8+ T cell responses can compensate for inadequate antibody responses (24) and are also highly cross-reactive to different variants of concern (32). Therefore, SARS-CoV-2 specific T cells may provide an additional contribution to the immune correlate of protection. An advantage of inactivated vaccines is that they also contain additional viral antigens, such as the highly abundant and immunogenic N protein which may also elicit T-cell immunity (33) and contributed to the higher responses in CoronaVac subjects here. Sampling at earlier timepoints (day 7-14) post vaccination may reveal differences in response magnitude given the phenotypic expansion of S-specific TEM CD4+ T cells by the BNT162b2 vaccine which may also have greater recall potential at infection (34). Longitudinal follow-up is thus necessary to confirm the long-term T cell responses between the two vaccines.
Proline mutations used in BNT162b2 vaccine to stabilize the spike protein in its pre-fusion state (35) and the serial passages in Vero cells and inactivation procedures used in the CoronaVac vaccine (36) are differences that may contribute to the differences in neutralizing antibody responses elicited by the two vaccines.
Preliminary results from another study has recently reported a marked difference of immunogenicity between BNT162b2 and CoronaVac in healthcare workers (37). However, the median age of the two groups were 37 and 47 years old respectively in that study while our subgroup was designed to be age matched. Moreover, the ADCC and T cell responses were not examined in their cohort.
There were some limitations in our study. The choice of vaccine was not randomized and there might be a selection bias in those opting for each vaccine. Our study only focused on investigating the immunogenicity at 1 month after two doses of vaccination. The durability of immune responses needs to be monitored; and indeed, this cohort will be followed up to address this question in the coming years. Our estimates to adjust for antibody waning was based on reports on CoronaVac, comparable to data for BioNTech was lacking. Thus our assumption of comparable rates of antibody waning for the two vaccines may not be correct. We did not collect plasma prior to the second dose of vaccine to assess the effect of the first dose of the vaccine, or acute phase responses, where earlier responses may account for the final post vaccine differences as our primary study endpoint was neuralization titers after vaccination. Similar comparisons between vaccines in teenagers and older adults will be needed.
Conclusion
Our data showed that the levels of antibodies elicited by CoronaVac were significantly lower than those of BNT162b2 in PRNT50, PRNT90, sVNT, spike RBD ELISA, spike NTD ELISA and spike S2 ELISA assays, and total and high avidity FcgRIIIa spike. CoronaVac induce higher CD4+ and CD8+ T cell responses to the structural protein than BNT162b2.
Data Availability
All data produced in the present work are contained in the manuscript
FIGURE LEGENDS
PRNT90 antibody titres in 70 sera from 50 symptomatic COVID-19 patients collected between 30 to 60 days after illness onset is shown. Geometric mean antibody titers (GMT) is denoted in the solid lines and 20% of GMT corresponding to 50% protective antibody titre is shown by the dashed lines, respectively.
The levels of SARS-CoV-2 (A) full length (N) or (B) C-terminal domain (CTD) of nucleoprotein antibodies were detected from the plasma collected at 1 month after received two doses of BNT162b2 (n=49) or CoronaVac (n=49). Dotted line represents the negative cut-off value.
SARS-CoV-2 Spike specific T cells by Flow Cytometry. (A) Gating strategy for characterization of CD4+ and CD8+ T cells by cytokine production (IFNγ, TNFα, IL2, IL4) or memory phenotype (CCR7, CD45RA). (B) Representative FACS plots of IFNg+ responses (B) and memory phenotypes (C) post-SARS-CoV-2 infection or vaccination with inactivated CoronaVac or mRNA vaccine BNT162b2.
PBMCs from pre- (day 0) and post-vaccination (day 30 post after second dose) of BNT162b2 mRNA (pre-vaccine n=25, post-vaccine n=25) and CoronaVac (pre-vaccine n=30, post-vaccine n=30) and recovered COVID-19 cases (n=10, timepoint 59 days post infection +/- 20 days (mean+/- SD)) were stimulated with S peptides or a DMSO control. The percentage of (A) IFNγ+CD4+ and (B) IFNγ+CD8+ T cell were measured by flow cytometry. Dotted lines represent limit of detection following DMSO background subtraction (IFNγ of CD4+=0.001, IFNγ of CD8+=0.001. (C) The proportion of IFNγ producing IL2 and TNFα CD4+ and CD8+ T cells post vaccination. (D) The phenotype (by CCR7 and CD45RA) of IFNγ responses for T effector memory (TEM), central memory (TCM), terminal effector memory (TeEM) or naïve (TN) CD4+ and CD8+ T cells post vaccination. Bars represent the mean values, and error bars represent SD. Statistical significance was determined by paired t-test between pre-versus post-vaccination timepoint samples, and Kruskal-Wallis test for multiple comparisons between vaccines, and COVID-19 patients. *p<0.05, **<0.01, ***<0.001.
Comparison of characteristics between the two vaccine groups (n=726)..
Multivariate analysis of factors associated with PRNT90 (n=97)
Factors associated with PRNT90 among subjects given the same vaccine (multivariate)
ACKNOWLEDGMENTS
We acknowledge the technical support from Mr Huibin Lv and Mr Ho Lun Lai. We also thank Dr Fung Hong (Chinese University Medical Center), Dr Ken Tsang (Kowloon Bay community vaccination center) and Dr Beatrice Cheng (Prince of Wales Hospital) for allowing us to recruit subjects for this study. The recombinant RBD and NTD proteins were kindly provided as gifts by Prof Ian A. Wilson and Dr Meng Yuan. This project utilized an Invitrogen Attune flow cytometer funded by the Pasteur Foundation Asia.
Footnotes
Funding source: This research was supported by grants from the Health and Medical Research Fund Commissioned Research on the Novel Coronavirus Disease (COVID-19), Hong Kong SAR (COVID1903003) (CKPM, SSN, GCL, MP and DSH), (COVID-190115 and COVID-190126, SAV), Guangdong Province International Scientific and Technological Cooperation Projects (2020A0505100063) (CKPM), the National Research Foundation of Korea (NRF) grant funded through the Korea government (NRF-2018M3A9H4055203) (CKPM), US National Institutes of Health (contract no. HHSN272201400006C) (MP), National Natural Science Foundation of China (NSFC)/Research Grants Council (RGC) Joint Research Scheme (N_HKU737/18) (CKPM and MP) and RGC theme-based research scheme (T11-712/19-N and T11-705/21-N) (DSH). Research was partly supported by Fast Grant ##2161 (Emergent Ventures) to G.K.A., and NIH grants (P01AI120943, R01AI123926 to G.K.A.; R01AI107056 to D.W.L.).
Conflict of Interest We declare no competing interests.
Human ethics approval declaration The study was approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee (Ref no: 2020.229) and all participants provided written consent.