key: cord-0689370-xo6a58si
authors: Hatzakis, Angelos; Karabinis, Andreas; Roussos, Sotirios; Pantazis, Nikos; Degiannis, Dimitrios; Chaidaroglou, Antigoni; Petsios, Konstantinos; Pavlopoulou, Ioanna; Tsiodras, Sotirios; Paraskevis, Dimitrios; Sypsa, Vana; Psichogiou, Mina
title: Modelling SARS-CoV-2 Binding Antibody Waning 8 Months after BNT162b2 Vaccination
date: 2022-02-13
journal: Vaccines (Basel)
DOI: 10.3390/vaccines10020285
sha: 2570b17ac433e5bbab69b7037e88784de42a24e1
doc_id: 689370
cord_uid: xo6a58si

Several lines of evidence suggest that binding SARS-CoV-2 antibodies such as anti-SARS-CoV-2 RBD IgG (anti-RBD) and neutralising antibodies (NA) are correlates of protection against SARS-CoV-2, and the correlation of anti-RBD and NA is very high. The effectiveness (VE) of BNT162b2 in preventing SARS-CoV-2 infection wanes over time, and this reduction is mainly associated with waning immunity, suggesting that the kinetics of antibodies reduction might be of interest to predict VE. In a study of 97 health care workers (HCWs) vaccinated with the BNT162b2 vaccine, we assessed the kinetics of anti-RBD 30–250 days after vaccination using 388 individually matched plasma samples. Anti-RBD levels declined by 85%, 92%, and 95% at the 4th, 6th, and 8th month from the peak, respectively. The kinetics were estimated using the trajectories of anti-RBD by various models. The restricted cubic splines model had a better fit to the observed data. The trajectories of anti-RBD declines were statistically significantly lower for risk factors of severe COVID-19 and the absence of vaccination side effects. Moreover, previous SARS-CoV-2 infection was associated with divergent trajectories consistent with a slower anti-RBD decline over time. These results suggest that anti-RBD may serve as a harbinger for vaccine effectiveness (VE), and it should be explored as a predictor of breakthrough infections and VE.

Two years after the beginning of the COVID-19 pandemic and one year into the widespread application of COVID-19, there is extensive data on vaccine efficacy/effectiveness (VE) based on pivotal randomised trials and several ongoing vaccine effectiveness studies around the globe. VE assessment is based on several outcomes such as symptomatic disease, hospitalisation, severe disease, death, asymptomatic disease, or any documented infection [1] .

The randomised trials of mRNA vaccines against COVID-19 illness initially provided data suggesting an efficacy of 95% for BNT162b2 and of 94.1% for mRNA-1273 at a median

Participants were vaccinated with 2 doses of BNT162b2 21 days apart. The study was designed to assess immunogenicity at time intervals of 1-2 weeks after the 2nd dose (28-35 days) and 4, 8, and 12 months after the 1st dose. Immunogenicity 1-2 weeks after the 2nd dose was expected to be highest, based on results from phase I/II studies [25] . In our previous study, we found that anti-RBD peaked 9-11 days after the 2nd dose [19] .

Vaccinated health care workers (HCWs) from 2 teaching hospitals, Laiko General Hospital (Hospital 1) and Onassis Cardiac Surgery Center (Hospital 2), participated in the immunogenicity study after signing informed consent (n = 871) [19] . Participants in this study comprised a group of HCWs from Onassis Cardiac Surgery Center who consented to donate an additional sample at the 6th month. All HCWs in this group (n = 97) provided samples at months 1, 4, 6, and 8 since their 1st vaccine dose.

A brief questionnaire was administered to HCWs concerning information about age, gender, education, position within the hospital, body mass index (BMI), history of risk factors for severe COVID-19 (RFS-CoV), previous COVID-19 (Pr-CoV) , and history of self-reported adverse reactions after vaccination (VSEs). VSEs were grouped according to the major symptoms such as local pain, fever, fatigue, and allergic and other systematic reactions. Combinations of VSEs were counted as one VSE in the final analysis [19] . RFS-CoV included hypertension, diabetes mellitus, cardiovascular disease, obesity, malignant neoplasm, chronic renal disease, chronic liver disease, and immunosuppressive treatment [19] .

Serum samples collected after venipuncture were tested for SARS-CoV-2 IgG binding antibodies to nucleocapsid protein (anti-N) and SARS-CoV-2 receptor-binding domain spike protein IgG. Both assays were described in detail in our previous publication [19] .

The first assay was a qualitative one, with an index (sample/calibrator (s/c)) cutoff of 1.4. Samples with an index ≥1.4 are considered positive and <1.4 negative. Reactive results are considered as an indication of natural infection.

The second assay (Abbott SARS-CoV-2 IgG II Quant) or anti-RBD was used to quantify IgG antibodies against the receptor-binding domain (RBD) of the S1 subunit of the spike protein. The linear range is between 21 and 40,000 AU/mL, with a cutoff value of 50 AU/mL. Both assays are based on chemiluminescent microparticle immunoassay (CLIA) [26] .

The correlation coefficient of Abbott anti-RBD with the World Health Organization (WHO) standard is 0.999, and the transformation of Abbott anti-RBD AU/mL to WHO BAU/mL is possible by using the equation BAU/mL = 0.142 × AU/mL [26].

Continuous variables were described using mean and standard deviation, or median, 25th, and 75th percentiles. Antibody titers were presented as geometric mean titers (GMTs) with the corresponding 95% CI. Categorical variables were described using frequencies and percentages. Chi-squared and Mann-Whitney U tests were used to assess differences in demographics and clinical characteristics among the participants of this study (n = 97) and the remaining ones (n = 774).

Exploratory data analysis for anti-SARS-CoV-2 RBD antibodies was based on available measurements within 26-258 days after the first dose of the vaccine or 5-237 days after the second dose. We used log 10 transformation in anti-SARS-CoV-2 RBD antibodies to normalise the distribution.

Longitudinal changes in antibody levels were analysed through three linear mixed models: (a) exponential model (EM), (b) power-law model (PLM), and (c) mixed model using four-knot restricted cubic splines (restricted cubic splines model) (RCSM), assuming that the peak of antibody titers was at Study Day 30 (i.e., 9 days after the second dose) [19] . Linear mixed models have been well studied within the classical/likelihood [27] and Bayesian [28] . frameworks. Typical applications of EM, PLM, and RCSM mixed models can be found in [14] . and [29] . Half-lives (i.e., time after peak required to reach antibody titers equal to 1 2 of the peak levels) were calculated based on the three previous models [14, [27] [28] [29] ].

The general form of the exponential model takes the following form:

where β 0 and β 1 are the fixed effects, intercept, and decay rate (slope), respectively. b 0i and b 1i are the random effects, intercept, and decay rate for each participant, and e i,j is the model errors for participant i at study day j. Hence, log 10 transformation of titers is a linear function of time (i.e., steady decay rate in log scale over time) [14, 30, 31] . The half-life (t 1/2 ) is given by the following equation:

and the 95% confidence interval of half-life was calculated using the delta method.

The general form of the power-law model takes the following form:

Where β 0 and β 1 are the fixed effects, intercept, and decay rate (slope), respectively. b 0i and b 1i are the random effects, intercept, and decay rate for each participant, and e i,j is the model errors for participant i at study day j. Study day was offset by 21 days to account for the 2nd dose regimen. Hence, the log 10 transformation of titers is a linear function of the log 10 transformation of time (i.e., decay rates decrease over time) [14, 31] . The half-life is given from the following form:

and the 95% confidence interval of half-life was calculated using the delta method.

A restricted cubic spline is a set of piecewise cubic functions (polynomials), where the boundaries of these pieces are called knots. The curves pass through all the knots, and both first-and second-order derivatives (the slope and the rate of slope, respectively) are the same for both functions on either side of a knot. Restricted splines are constrained to be linear beyond boundary knots (i.e., before the first knot and after the last one) [32] .

Restricted cubic splines with four knots were used to estimate the time kinetics curves of antibodies. The location of the knots was placed in the percentiles recommended by Harrell [33] . The general form of the mixed model with a four-knot restricted cubic spline takes the following form [32] :

where b 0i is the random intercept for each participant, e i,j is the model errors for participant i at study day j, S k (t), k = 1, 2, 3 are restricted cubic splines with four-knot terms, s 1 (t), s 2 (t) are restricted cubic splines with three-knot terms, and b 1i , b 2i are the corresponding random effects. The half-life and the corresponding 95% confidence interval of were estimated using a Monte Carlo procedure.

Log-likelihood is a measure of model fit. Higher values indicate better fit. To compare the fit of the evaluated models, we used the likelihood ratio test and the AIC (Akaike information criterion), which takes the following form:

where k is the number of model parameters. Lower values of AIC indicate a better fit.

The best-fitted model was applied to assess relative differences in means of anti-RBD IgG levels by gender, age, risk factors for severe COVID-19 illness, side effects of vaccination, and history of previous COVID-19 infection in the post-vaccination period (5-237 days after the 2nd dose).

All analyses were performed using Stata version 13.0 [33] . All p-values were two-sided, and p < 0.05 was considered statistically significant.

Ninety-seven vaccinated HCWs, donating 388 blood samples, participated in all rounds of blood sampling. Their demographic and clinical characteristics are shown in Table 1 . The majority were females (55.7%), Greek nationals (93.8%), highly educated holding an MSc, Ph.D., or MD degree (40.2%), treating or caring patients (67%) with the mean (SD) age of 50.1 (9.4) years and BMI of 26.4 (5.5) kg/m 2 , including 21.6% obese, 15.5% reporting RFS-CoV, 54.6% VSEs, and 7.2% with Pr-CoV. Overall, the distribution of demographic and clinical characteristics of 97 HCWs was similar to the distribution of 774 of the remaining HCWs, who participated in our previous study [19] , with the exception of HCWs involved in patient care who were more frequent in the studied population (83% vs. 67%) ( Table 1) .

In Supplementary Table S1, the timing of blood samplings is shown. The median (25th-75th) times were 30 (29) (30) (31) (32) , 125 (124-126), 183 (182-185), and 251 (250-252) days from the first vaccine dose to first, second, third, and fourth measurement, respectively. The corresponding times from the second vaccine dose were 9 (8-10), 103 (102-105), 161 (160-163), and 229 (228-231) days, respectively.

Overall, the GMTs (95% CI) at times 30, 125, 183 and 251 were 13,674 (11, 749 ) AU/mL, 2112 (1716-2601) AU/mL, 1036 (834-1288) AU/mL, and 617 (490-778) AU/mL, respectively. The percentage reduction in GMTs of anti-RBD levels from the 30th day to the 125th, 183rd, 251st days was 85%, 92% and 95%, respectively. The percentage reductions were similar when we used the median of anti-RBD levels (86%, 93%, and 96%, respectively) ( Table 2 , Supplementary Figure S1 ). We further explored this significant waning of anti-RBD by studying the pattern of waning. We used three models of reduction: EM, PLM, and RCSM. Based on the AIC criterion, the best fitting was the RCSM with AIC values −202.4 compared with 155.8 and 58.5 for EM and PLM, respectively (Table 3, Figure 1 ). Abbreviation: df, degrees of freedom; AIC, Akaike information criterion; CI, confidence interval.

Comparison of three models (splines, exponential, and power-law) for the predicted trajectories of anti-SARS-CoV-2 RBD IgG antibody levels over time.

According to RCSM, the day where antibody levels reached 50% reduction from the peak levels of the 30th day was day 62 (95% CI 59.9-64.5), while in EM and PLM, this was the 79th (95% CI 76-81) day and the 40th (95% CI 39-41) day, respectively ( Table 3) .

The GMTs of anti-RBD levels trajectories are presented by gender, age, risk factors for severe COVID-19 illness, side effects of vaccination, and history of COVID-19 infection (Table 4 ). According to RCSM, the day where antibody levels reached 50% reduction from the peak levels of the 30th day was day 62 (95% CI 59.9-64.5), while in EM and PLM, this was the 79th (95% CI 76-81) day and the 40th (95% CI 39-41) day, respectively ( Table 3) .

The GMTs of anti-RBD levels trajectories are presented by gender, age, risk factors for severe COVID-19 illness, side effects of vaccination, and history of COVID-19 infection (Table 4) .

By using the RCSM mixed-effects model, we estimated differences in anti-RBD in the trajectory from 30 up to 251 days from the first dose. In the univariable and multivariable analysis, age and sex were not associated with difference in the decline trajectory. The presence of RFS-CoV was significantly associated with decreasing anti-RBD levels by 43.6% (95% CI 5.1-66.7%, p = 0.031). A lack of VSEs was associated with decreasing levels of anti-RBD throughout the immunogenicity trajectory by 38.6% (95% CI 11.9-57.2%, p = 0.001). The aforementioned effects did not show any statistically significant variation over time (i.e., interactions with time terms were not significant), with the exception of previous SARS-CoV-2 infection, the effect of which on anti-RBD levels increased with time (interaction with time terms p-value = 0.019). The increase ranged from 20.3% to 205.8%, resulting in lesser waning in those vaccinated with Pr-CoV (Table 4, Figure 2 ). 

There is overwhelming evidence that neutralising (NA) and binding antibodies such as anti-spike SARS-CoV-2 IgG (anti-spike) and anti-RBD SARS-CoV-2 IgG are immune correlates of protection (Cor-P). Cor-P are immunological markers that can be used to reliably predict the level of vaccine efficacy in a clinically relevant protection endpoint or the level of protection against breakthrough infections in vaccinated individuals [24, [34] [35] [36] [37] [38] [39] [40] [41] [42] . Cor-P measured longitudinally in a vaccinated cohort may inform on the waning of vaccine 

There is overwhelming evidence that neutralising (NA) and binding antibodies such as anti-spike SARS-CoV-2 IgG (anti-spike) and anti-RBD SARS-CoV-2 IgG are immune correlates of protection (Cor-P). Cor-P are immunological markers that can be used to reliably predict the level of vaccine efficacy in a clinically relevant protection endpoint or the level of protection against breakthrough infections in vaccinated individuals [24, [34] [35] [36] [37] [38] [39] [40] [41] [42] .

Cor-P measured longitudinally in a vaccinated cohort may inform on the waning of vaccine efficacy, the timing of potential anamnestic dose, or the timing to update vaccine composition due to the emergence of new variants. Moreover, this may allow new vaccines to be authorised for use on the basis of immunogenicity and safety alone when large efficacy trials are not feasible. Cor-P have been established against many viral diseases [43] .

Immune correlates of protection are meaningful for protection against symptomatic COVID-19 but not against asymptomatic infection [37, 42] . Studies have shown increasing levels of VE associated with increasing immune marker levels. However, a threshold for protection was not found [37, 42] . Among the immune markers examined, anti-RBD and anti-spike have a similar correlation with NA in predicting VE [37, 42] . This finding has major implications for monitoring the VE, since NA testing is complex and not widely available. Sparse data from the literature suggest a reduction in VE for BNT162b2 against asymptomatic infection from 73% 1 month after the 2nd dose to 48% and 24% at months 3 and 5, respectively [8] . This pattern of VE reduction against infection follows the anti-RBD reduction pattern observed in this and other studies [1, 12, 13, 17, 44] . Therefore, the monitoring of binding antibodies in a vaccinated cohort may serve as a harbinger of waning VE, and it may highlight the need for a booster dose.

The three models applied to assess waning humoral immunity belong to the same family of mixed linear models. The only difference between the three models is related to the treatment of time as entered in the fixed and random effects (i.e., untransformed in the EM, log transformed in the PLM, and through spline functions in the RCSM). Existence, estimation, and properties of the solutions of mixed linear models are well studied within the classical/likelihood based framework and the Bayesian frameworks [27, 28] . Such models have been applied numerous times in biomedical research when the evolution of a continuous marker over time is of interest and repeated measurements per individual are available. Typical examples include the evolution of CD4 cell count or HIV-RNA viral load levels in HIV-infected individuals, prostate-specific antigen (PSA) levels in prostate cancer patients, BMI evolution in children, respiratory parameters in chronic obstructive pulmonary disease (COPD) patients, etc. A typical example of a mixed linear model with restricted cubic splines can be found in the work of Mallon PW et al. [29] .

Previous studies on antibody waning modelling with Moderna mRNA-1273 and Pfizer BNT162b2 vaccines used EM or PLM, respectively [14, 17, 24] . Our current analysis suggests that these models cannot accurately predict the antibody waning pattern, and an RCSM is more appropriate. A mixed-effects model with a restricted cubic spline has the flexibility to capture the more complex relationship between the variables, as well as a rich structure in random effects using a restricted cubic spline with four knots. Thus, better modelling was achieved both for the mean of anti-RBD levels and their variability. However, a large study including more knots may provide more accurate prediction of waning pattern.

The study of antibody trajectories is more appealing than using individual time points to study the demographic and clinical characteristics that might be associated with antibody levels. Based on RCSM, we examined the post-vaccination anti-RBD trajectories 3-251 days after the first dose according to gender, age, RFSC, VSE, and Pr-CoV. The gender and age trajectory differences did not differ in a statistically significant way. Statistically significant differences in the trajectories of anti-RBD were found for RFS-CoV, VSE, and Pr-CoV. Previous studies documented that after completion of the second dose of BNT162b2, those not reporting RFS-CoV and VSE had higher anti-RBD levels [12, 23, 44] . A finding of our study is the increasing anti-RBD difference over time according to previous COVID-19 infection, which suggests a faster decline in vaccinated individuals with no prior infection compared with Pr-CoV [45] . The waning of antibodies in Pr-CoV is consistent with a large study from Qatar, where the cumulative incidence of breakthrough infections was 5.5 times higher in vaccinated individuals with no prior COVID-19 infection compared with Pr-CoV ones [46] . Several studies suggest that the quality of immune response in vaccinated Pr-CoV improves over time in terms of neutralisation capacity and breadth (hybrid immunity), with the improved neutralisation of new variants including VOC Omicron, especially after the use of booster dose [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] .

The limitations of our study include a restricted age range and a small sample size for a thorough assessment of anti-RBD trajectories. Previous studies have shown a sharp reduction in anti-RBD after the age of 65 years [12, 54] . Additionally, humoral immunity studies may not adequately capture SARS-CoV-2 immunity due to the lack of data on cellular immunity such as B-and T-cell SARS-CoV-2-specific markers.

Overall, we conclude that (1) the profound decline of anti-RBD 4-8 months after the first dose is a harbinger of a reduction in vaccine effectiveness to prevent SARS-CoV-2 infection, (2) the use of a more precise PCSM for an assessment of the waning pattern will facilitate modelling studies to evaluate the determinants of humoral immunity patterns, and (3) anti-RBD should be explored as a predictor for breakthrough infections and as a strategy to simplify the assessment of vaccine efficacy for new vaccines [36, 59] . 

None of the authors have relevant conflict of interest to declare. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Duration of effectiveness of vaccines against SARS-CoV-2 infection and COVID-19 disease: Results of a systematic review and meta-regression

Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine

Efficacy and safety of the mRNA-1273 SARS-CoV-2 vaccine

Efficacy of the mRNA-1273 SARS-CoV-2 vaccine at completion of blinded phase

Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine through 6 months

Comparative Effectiveness of BNT162b2 and mRNA-1273 Vaccines in US Veterans

BNT162b2 and mRNA-1273 COVID-19 vaccine effectiveness against the SARS-CoV-2 Delta variant in Qatar

Waning of BNT162b2 vaccine protection against SARS-CoV-2 infection in Qatar

Effectiveness of mRNA BNT162b2 COVID-19 vaccine up to 6 months in a large integrated health system in the USA: A retrospective cohort study

Waning immunity after the BNT162b2 vaccine in Israel

Effectiveness of Covid-19 vaccines against the B. 1.617. 2 (Delta) variant

Waning immune humoral response to BNT162b2 Covid-19 vaccine over 6 months

Dynamics of antibody response to BNT162b2 vaccine after six months: A longitudinal prospective study

Antibody persistence through 6 months after the second dose of mRNA-1273 vaccine for COVID-19

Durability of mRNA-1273 vaccine-induced antibodies against SARS-CoV-2 variants

SARS-CoV-2 neutralization with BNT162b2 vaccine dose 3

Waning of IgG, total and neutralizing antibodies 6 months post-vaccination with BNT162b2 in healthcare workers

Differential kinetics of immune responses elicited by COVID-19 vaccines

Comparative Immunogenicity of BNT162b2 mRNA Vaccine with Natural SARS-CoV-2 Infection

Antibody responses to BNT162b2 mRNA COVID-19 vaccine and their predictors among healthcare workers in a tertiary referral hospital in Japan

Initial observations on age, gender, BMI and hypertension in antibody responses to SARS-CoV-2 BNT162b2 vaccine

Age-related immune response heterogeneity to SARS-CoV-2 vaccine BNT162b2

Does a lack of vaccine side effects correlate with reduced BNT162b2 mRNA vaccine response among healthcare workers and nursing home residents?

Waning of SARS-CoV-2 antibodies targeting the Spike protein in individuals post second dose of ChAdOx1 and BNT162b2 COVID-19 vaccines and risk of breakthrough infections: Analysis of the Virus Watch community cohort

BNT162b2 vaccine induces neutralizing antibodies and poly-specific T cells in humans

Random-effects models for longitudinal data

Illustration of Bayesian inference in normal data models using Gibbs sampling

Weight gain before and after switch from TDF to TAF in a US cohort study

Heterogeneity and longevity of antibody memory to viruses and vaccines

Modeling the long-term antibody response of a hepatitis E vaccine

StataCorp. Stata Statistical Software: Release 13

Correlates of protection against SARS-CoV-2 in rhesus macaques

Neutralizing antibody levels are highly predictive of immune protection from symptomatic SARS-CoV-2 infection

Evidence for antibody as a protective correlate for COVID-19 vaccines

Correlates of protection against symptomatic and asymptomatic SARS-CoV-2 infection

Neutralizing antibodies correlate with protection from SARS-CoV-2 in humans during a fishery vessel outbreak with a high attack rate

Immune correlates of protection by mRNA-1273 vaccine against SARS-CoV-2 in nonhuman primates

COVID-19 breakthrough infections in vaccinated health care workers

New SARS-CoV-2 variants-Clinical, public health, and vaccine implications

Immune correlates analysis of the mRNA-1273 COVID-19 vaccine efficacy clinical trial

Plotkin's Vaccines

BNT162b2 COVID-19 vaccine and correlates of humoral immune responses and dynamics: A prospective, single-centre, longitudinal cohort study in health-care workers

Durability of Antibody Levels after Vaccination with mRNA SARS-CoV-2 Vaccine in Individuals with or without Prior Infection

Association of prior SARS-CoV-2 infection with risk of breakthrough infection following mRNA vaccination in Qatar

Naturally enhanced neutralizing breadth against SARS-CoV-2 one year after infection

Impact of circulating SARS-CoV-2 variants on mRNA vaccine-induced immunity

High genetic barrier to SARS-CoV-2 polyclonal neutralizing antibody escape

Plasma neutralization properties of the SARS-CoV-2 Omicron variant

Booster of mRNA-1273 Vaccine Reduces SARS-CoV-2 Omicron Escape from Neutralizing Antibodies

mRNA booster immunization elicits potent neutralizing serum activity against the SARS-CoV-2 Omicron variant. medRxiv 2021

SARS-CoV-2 B. 1.1. 529 variant (Omicron) evades neutralization by sera from vaccinated and convalescent individuals

Age-dependent immune response to the Biontech/Pfizer BNT162b2 COVID-19 vaccination

Effectiveness of a third dose of the BNT162b2 mRNA COVID-19 vaccine for preventing severe outcomes in Israel: An observational study

Protection against COVID-19 by BNT162b2 booster across age groups

BNT162b2 Vaccine Booster and Mortality Due to COVID-19

Hybrid immunity improves B cells and antibodies against SARS-CoV-2 variants

Neutralising antibody titres as predictors of protection against SARS-CoV-2 variants and the impact of boosting: A meta-analysis